Measuring holographic entanglement entropy on a quantum simulator

Quantum simulation promises to have wide applications in many fields where problems are hard to model with classical computers. Various quantum devices of different platforms have been built to tackle the problems in, say, quantum chemistry, condensed matter physics, and high-energy physics. Here, we report an experiment towards the simulation of quantum gravity by simulating the holographic entanglement entropy. On a six-qubit nuclear magnetic resonance quantum simulator, we demonstrate a key result of Anti-de Sitter/conformal field theory (AdS/CFT) correspondence—the Ryu-Takayanagi formula is demonstrated by measuring the relevant entanglement entropies on the perfect tensor state. The fidelity of our experimentally prepared the six-qubit state is 85.0% via full state tomography and reaches 93.7% if the signal-decay due to decoherence is taken into account. Our experiment serves as the basic module of simulating more complex tensor network states that exploring AdS/CFT correspondence. As the initial experimental attempt to study AdS/CFT via quantum information processing, our work opens up new avenues exploring quantum gravity phenomena on quantum simulators.


INTRODUCTION
The study of quantum systems requires an exponential amount of resources on conventional computers due to the exponentially growing dimensionality of Hilbert spaces, which makes it impossible to model even with supercomputers. Quantum simulators, conceived by Feymann in 1982, 1 are special purpose devices designed to imitate the behaviors or properties of other less accessible quantum systems. 2,3 Over the past few years, proofof-principle experiments have been realized in simulating quantum phase transitions, 4,5 topological order, 6,7 molecular energies, 8,9 quantum chaos, 10,11 and so on. However, there exists a significant field-quantum gravity-that has never been explored by experimental quantum simulation. Many important ideas such as holographic principle and Anti-de Sitter/conformal field theory (ADS/CFT) correspondence remained unrevealed in experiment. Recent development of a discrete version of AdS/CFT correspondence in terms of tensor networks (TN) motivates us to studying AdS/CFT correspondence on quantum simulators. In this work, we make first steps toward to the simulation of quantum gravity on a 6-qubit nuclear magnetic resonance (NMR) quantum processor, where rank-6 perfect tensor that forms the building block of complex TN is realized with high accuracy.
We start from a basic introduction to AdS/CFT correspondence. AdS/CFT correspondence is one of the most prominent approaches towards a quantum theory of gravity for over two decades. 12,13 It is the most successful realization of the holographic principle to date, by stating that the quantum gravity theory in the bulk anti-de Sitter spacetime is equivalent to a quantum conformal field theory on the lower-dimensional boundary of the spacetime. The AdS/CFT correspondence has recently become a bridge connecting quantum gravity to quantum information theory, 14,15 which inspires revolutionary ideas of developing quantum gravity using the methods in quantum information and entanglement. A key result in this perspective is the holographic entanglement entropy characterized by the Ryu-Takayanagi (RT) formula, which relates the entanglement entropy of the boundary quantum system to the bulk geometry: S EE ðAÞ is the entanglement entropy of a (d − 1)-dimensional boundary region A, while Ar min is the area of the bulk (d − 2)dimensional minimal surface anchored to A. [16][17][18] G N is the Newton constant. See Fig. 1a for a brief illustration.
Recently, a discrete version of AdS/CFT is realized on a type of lattices called TN, [19][20][21][22] making it possible to be demonstrated on a quantum simulator device in practice. In general, TN states are ways of rewriting a many-body wave function in terms of contractions of tensors, aiming at obtaining the ground states of interacting many-body Hamiltonians in a numerically efficient way. As a key observation related to AdS/CFT, the TN state has an emergent bulk dimension built by the layers of tensors, making it an ideal ground for manifesting AdS/CFT in many-body systems. Indeed, the theoretical studies have found that the TN made of perfect tensors (PT) can demonstrate interesting holographic properties. In particular, the entanglement entropy of perfect tensor TN gives a discrete realization the above RT formula. 20 In this work, we demonstrate the RT formula on a quantum simulator that simulates a PT of rank-6. Using a six-qubit quantum register in the NMR system, we create the rank-6 PT and subsequently measure its holographic entanglement entropy. The experimental results demonstrate the RT formula if the decoherence effect is taken into account. As the rank-6 PT serves as the building block to construct the entire TN, our experiment also opens up a new and practical way of studying AdS/CFT and the holographic principle at large.

Perfect tensors
The TN that we focus on is shown in Fig. 1b, where each hexagon represents a special six-qubit state |ψ〉. |ψ〉 is called a PT, if and only if that any three-qubit subsystem out of six is maximally entangled with the rest. It is shown that, for a TN made by the PT, its entanglement entropy is holographic and gives the discrete RT formula on the lattice. Actually, the entanglement entropy of such TN equals the minimal number of links cut by the virtual surfaces anchored to the boundary, as illustrated in Fig. 1b.
To prove the above statement, we first introduce the form of the rank-6 PT, which is the building block of the TN. Given the single-qubit Hilbert space H ' C 2 , a rank-6 PT |ψ〉 is a state in H 6 , such that for any bipartition of qubits m + k = 6, the entropy of the reduced density matrix is maximal. Assuming m ≥ k, and labeling the orthonormal basis in H m and H k by |α〉 and |i〉 (b) (a) (c) Fig. 1 a A sketch of the RT formula. The hexagonal tiling indicates that the disk is a 2-dimensional ads space. The red solid arc in the bulk is the minimal surface (a line in this case) anchored to the two ends of a chosen boundary region A. b A discretization of a by a tensor network comprised of rank-6 tensors. Each hexagonal node represents a rank-6 tensor state jψi 2 H 6 , and the collection of all such nodes corresponds to the tensor product of all |ψ〉's. Each link ' represents a maximally entangled state . Connecting one leg of the node to a link corresponds to taking the inner product in H. The dangling legs are physical qubits in the many-body system. The red dashed arc illustrates the virtual surface S anchored to region A, which cuts a minimal number of links. c Rank-6 PT from the TN with the minimal number of cuts equal to three. The six legs represent six qubits. Three qubits are at the boundary and the other three are bulk qubits. This is the model realized in our experiment respectively, a PT jψi ¼ P In other words, the reduced density matrix ρ (k) by tracing out m qubits is an identity matrix, whose entanglement entropy is simply k, the number of remained qubits. In this Letter, we use the superscript (k) to represent the k-qubit subsystem.
With the rank-6 PT (explicit form in appendix A (see Supplemental Information for a detailed description of the theory and experiment)) in hand, the TN state illustrated in Fig. 1b is constructed as follows. Each internal link ' represents a two-qubit maximally entangled state where two qubits associate respectively to the two end points of '. If we denote by |ψ(n)〉 the PT associated to the hexagon node n, the total TN state |Ψ〉 in Fig. 1b is written as a (partial) inner product form The inner product takes place at the end points of each internal link ', between one qubit in j'i and the other in |ψ(n)〉. The qubits in |ψ(n)〉 not participating the inner product are boundary qubits corresponding to the dangling legs, and these boundary ones are actually physical qubits, indicating that |Ψ〉 is a state on the boundary. We then pick a boundary region A which collects a subset of the boundary qubits, as shown in Fig. 1b. The reduced density Þis computed by tracing out all boundary qubits outside A. Initially, this partial trace boils down to computing the reduced density matrix of individual tensors closest to the boundary. By applying Eq. (2) and noticing that j'i is maximally entangled, the trace computation can be effectively pushed from the boundary into the bulk, meaning that the partial trace on the boundary is now equivalent to computing the reduced density matrix of the PT inside the bulk (see Supplemental Information for a detailed description of the theory and experiment). Once again, we can apply Eq. (2) and push the trace further inside. This iteration procedure is repeated until the trace reaches S in Fig. 1b, where Eq. (2) is not anymore valid, as the number of qubits participating the trace (number of links cut by S) is less than three for each tensor. Now we have presented a sketch about how to calculate the entanglement entropy of ρ A via Eq. (3), and direct readers to Appendix B (see Supplemental Information for a detailed description of the theory and experiment) for a concise proof using the graphical computation of TN. Firstly, tr ρ A ð Þ is found to be equal to the number of qubits on S, i.e., the same as the number of links cut by S. Moreover, the product ρ 2 A , involving the inner product of boundary qubits in A, gives that ρ 2 A / ρ A . Note that we have ignored all numerical prefactors but they all cancel when calculating trρ n A ðtrρ A Þ n in the entanglement entropy. As a result, the Von Neumann entropy gives (see Supplemental Information for a detailed description of the theory and experiment) S EE ðAÞ ¼ lim n!1 1 1Àn log 2 trρ n A ðtrρ A Þ n ¼ minimal number of cuts by S: The above result is a discrete version of the RT formula in Eq. (1). The "minimal number of cuts" represents the minimal area Ar min (in the unit of Planck scale) in the RT formula. The bulk surface S with minimal area emerges effectively from the entanglement entropy of the TN state. Equation (5) demonstrates explicitly that the bulk geometry are created holographically by the entangled qubits of the boundary many-body system.
It is worth emphasizing that, all descriptions about constructing the TN originate from the PT in Eq. (2). Therefore, this rank-6 PT plays the fundamental role in holographic entanglement entropy, and is a key of emerging bulk gravity from TN states. If we choose S as shown in Fig. 1c by which the minimal number of cuts is three, a rank-6 PT is generated where the boundary and bulk qubits are both three. Here, we demonstrate the emergent gravity program in AdS/CFT for the first time in a six-qubit NMR quantum simulator, by creating the rank-6 PT in Fig. 1c and measuring the relevant entanglement entropies.
Experiment implementation of a rank-6 perfect tensor The six qubits in the NMR quantum register are denoted by the spin-1/2 13 C nuclear spins, labeled as 1 to 6 as shown in Fig. 2a, in 13 C-labeled Dichloro-cyclobutanone dissolved in d 6 -acetone. All experiments were carried out on a Bruker DRX 700 MHZ spectrometer at room temperature. The internal Hamiltonian of this system is where ν j is the resonance frequency of the jth spin and J jk is the Jcoupling strength between spins j and k. All parameters including the relaxation times for each spin are listed in appendix C (see Supplemental Information for a detailed description of the theory and experiment). To control system dynamics, we have external control pulses with four adjustable parameters: the amplitude, frequency, phase, and duration, based on which arbitrary singlequbit rotations can be realized with simulated fidelities over 99.5% (see Supplemental Information for a detailed description of the theory and experiment). A rank-6 PT can be created from |0〉 ⊗n through the circuit as illustrated in Fig. 2b, which involves only Hadamard gates and 4 (a) (b)  Fig. 2 a Molecular structure of the 13 C-labeled six-qubit quantum processor. The six qubits of the rank-6 PT are mapped to 1 to 6, respectively. b Quantum circuit that evolves the system from |0〉 ⊗6 to the PT, constructed by several Hadamard gates (blocks) and controlled-Z operations (lines connecting two dots) controlled-Z gates. Experimentally, this requires an initialization of the system onto |0〉 ⊗n . However, initializing an NMR processor to | 0〉 ⊗n is based upon the pseudo-pure state technique, which leads to an exponential signal attenuation. Here, we adopt a temporal averaging approach that enables the PT preparation directly from the thermal equilibrium of NMR, while skipping the intermediate pseudo-pure state stage to avoid the above problem, as shown in Appendix C (see Supplemental Information for a detailed description of the theory and experiment)).
After the creation, we conducted k-qubit (1 ≤ k ≤ 5) quantum state tomography in the corresponding subspace of the whole system, respectively. For simplicity, the cutting of the links was chosen to be continuous in experiment, i.e., in a cyclic manner. It means that six state tomographies for any given k are performed, e.g., when k = 2, we reconstructed ρ Combined with the fact that S (k) = S (6−k) for a six-qubit pure state, we have S (k) = min{k, 6 − k} for the theoretical PT, as shown by the orange dashed line in Fig. 3a. In experiment however, inevitable errors lead to imperfection and hence impurity in the truly prepared state, so we cannot just measure k ≤ 3 cases to deduce other k's. Therefore, we measured and compare the experimental S (k) for each 1 ≤ k ≤ 5 (red circles) with their theoretical predictions in Fig. 3a. For each k, the mean and error bar of the experimental S (k) value are calculated from the six cyclic tomographic results. When k ≤ 3, the measured entanglement entropies match extremely well with the theory; when k > 3, there are notable discrepancies between theory and experiment, which should be primarily attributed to decoherence errors, as discussed in the following.
The pulse sequence that creates the PT is around 60 ms; this is not a negligible length compared to the T Ã 2 time (~400 ms) of the molecule, meaning that decoherence will induce substantial errors during experiments(see Method). As T Ã 2 relaxation is the dominating factor, the off-diagonal terms in the PT density matrix are mainly affected. To estimate this imperfection, we performed full state tomography 23 on the prepared state and got ρ e . The real part of ρ e is depicted in the right panel of Fig. 3b, by projecting each element onto a two-dimensional plane. As a comparison, the figure of the theoretical PT ρ pt = |ψ〉〈ψ| is placed in the left panel of Fig. 3b. In fact, the diagonal elements of ρ e are almost the same as that of ρ pt , but the off-diagonal are lower due to the T Ã 2 errors. The state fidelity between ρ e and ρ pt , defined as is about 85.0%. Direct observations of ρ e in terms of NMR spectra are also shown in Fig. 3c, where experimental and simulated spectra highly match if the experimental signal is rescaled by 1.25 times to compensate for the decoherence effect.
Although the reconstructed state ρ e is prone to the decoherence errors, the entanglement entropies for the cases k ≤ 3 in Fig. 3a are still in excellent accordance with the theory. The reason is, when we trace out three or more qubits, the reduced density matrix is predicted to be identity according to Eq. (2), so the measured k ≤ 3 reduced density matrices are almost irrelevant to the imperfection of the off-diagonal elements in ρ e . However, when k > 3, the reduced density matrix is no longer the identity, meaning that the imperfect off-diagonal terms in ρ e start to be responsible for calculating S (k) . As a result, in Fig. 3a we have S (4) = 2.91 ± 0.20 and S (5) = 2.32 ± 0.25 (red circles) respectively, which are quite distant from the theoretical curve. After numerically simulating and compensating for the decoherence errors 24,25 during the PT creation, we found that the two entanglement entropies S (4) and S (5) approach much closer to the theory, which are now 2.27 ± 0.46 and 1.37 ± 0.28 (blue squares), respectively. We also calculated the current fidelity between the rescaled experimental state and ρ pt via Eq. (7), and found it improved to 93.7%, which is 8.7% greater than that of ρ e . Experimental results are represented by the red circles, where S (4) and S (5) do not fit very well. If the signal's decay due to decoherence is taken into account, the experimental results are rescaled to the blue squares, which fit much better. As a upper-bound reference, the maximal entropy of a k-qubit subsystem is also plotted (green dotted line) by assuming a six-qubit identity. b Density matrices of the theoretical rank-6 PT ρ pt (left) and the experimentally reconstructed state ρ e (right) on a two-dimensional plane. The rows and columns are labeled by the sixqubit computational basis from |0〉 ⊗6 to |1〉 ⊗6 , respectively. c Direct observation of ρ e in the NMR spectra (red), with probe qubits C 1 (top) and C 4 (bottom), respectively. The simulated spectra of the PT are also shown in blue. For a better visualization, experimental signals are rescaled by 1.25 times to neutralize the decoherence error

DISCUSSION
RT formula, or explicitly the TN built by the rank-6 PT in Fig. 1b, tells us how to deduce the bulk geometry using the entanglement on the boundary. The implicit condition here is that the global TN state is pure. Otherwise, the information on the boundary cannot uniquely (up to local unitaries) determine the bulk geometry, e.g., it cannot specify whether the TN state is the maximally mixed identity or PT since both give the same entanglement entropies on the boundary (meaning k ≤ 3) as shown in Fig. 3a. In experiments, however, under realistic noises, it is difficult to guarantee the purity of the truly created states because experimental procedures inevitably involve errors-in particular the decoherence that render the TN states mixed. In our experiment of a 6-qubit PT-a build block of a complex TN, we have achieved 85% fidelity, which is already state-of-the-art; however, there is yet some non-negligible decoherence due to the T Ã 2 errors. Therefore, our results successfully test the RT formula up to the decoherence.
The simulation of the holographic entanglement entropy can be generalized to TN s with multiple perfect tensors. In the Section E of the supplemental material, we demonstrate a simulation of the holographic entanglement entropy on a TN with seven tensors. The key to performing the simulation is that measuring the Rényi entropies of the TN s can be reduced to measuring the reduced density matrices of ρ e , and their multiplications and traces, while ρ e is simulated experimentally. The result of the simulation demonstrates agreement with the RT formula, up to the experimental noise in ρ e . The simulation can be generalized to other TNs. In conclusion, our work is an endeavor to demonstrate on a quantum simulator the RT formula (the discrete PT version) in the AdS/CFT correspondence. We utilize a temporal average technique to create the rank-6 PT and perform full state tomography to reconstruct the experimental state. This is also the largest full state characterization in an NMR system to date. Although the imperfection of the created state due to decoherence errors makes the holographic entanglement entropy not exactly agree with the theoretical prediction, we simulate and compensate for such type of errors under the realistic experimental environment, and demonstrate the accordance between theory and experiment thereafter. As the first step towards exploring AdS/CFT correspondence using a quantum simulator, our work provides valid experimental demonstrations about studying quantum gravity in the presence of realistic noises.

Decoherence simulation
To numerically simulate the decoherence effect in our six-qubit system, we made the following assumptions: the environment is Markovian; only the T Ã 2 dephasing mechanism is taken into account since T 1 effect is negligible in our circuit; the dephasing noise is independent between all qubits; the dissipator and the total Hamiltonian commute in each pulse slice as the Δt = 10 μs is small. With these assumptions, we simplified and solved the master equation in two steps for each Δt: evolve the system by the propagator calculated by the internal and control pulse Hamiltonian, and subsequently apply the dephasing factors according to the coherent orders for Δt which is an exponential decay of the off-diagonal elements in the density matrix. For each experiment of the 64 runs, we simulated the above process and obtained the signal's decay due to decoherence. From the experimental result, we then compensated for this decay, and a new state in which the decoherence effect was taken into account was thus achieved. The fidelity now between the rescaled experimental state and the theoretical PT is boosted to 93.7%.