Measuring Holographic Entanglement Entropy on a Quantum Simulator

3 downloads 0 Views 4MB Size Report
Apr 30, 2017 - 10Perimeter Institute for Theoretical Physics, Waterloo N2L 2Y5, Ontario, ... has recently become a bridge connecting quantum gravity to.
Measuring Holographic Entanglement Entropy on a Quantum Simulator Keren Li,1, 2, ∗ Muxin Han,3, 4, ∗ Guilu Long,1 Yidun Wan,5, 6, † Dawei Lu,2, 7, ‡ Bei Zeng,2, 7, 8, 9, § and Raymond Laflamme2, 9, 10

arXiv:1705.00365v1 [quant-ph] 30 Apr 2017

1

State Key Laboratory of Low-Dimensional Quantum Physics and Department of Physics, Tsinghua University, Beijing 100084, China 2 Institute for Quantum Computing and Department of Physics and Astronomy, University of Waterloo, Waterloo N2L 3G1, Ontario, Canada 3 Department of Physics, Florida Atlantic University, 777 Glades Road, Boca Raton, FL 33431, USA 4 Institut f¨ur Quantengravitation, Universit¨at Erlangen-N¨urnberg, Staudtstr. 7/B2, 91058 Erlangen, Germany 5 Department of Physics and Center for Field Theory and Particle Physics, Fudan University, Shanghai 200433, China 6 Collaborative Innovation Center of Advanced Microstructures, Nanjing, 210093, China 7 Department of Physics, Southern University of Science and Technology, Shenzhen 518055, China 8 Department of Mathematics and Statistics, University of Guelph, Guelph N1G 2W1, Ontario, Canada 9 Canadian Institute for Advanced Research, Toronto M5G 1Z8, Ontario, Canada 10 Perimeter Institute for Theoretical Physics, Waterloo N2L 2Y5, Ontario, Canada Anti-de Sitter/conformal field theory (AdS / CFT) correspondence is one of the most promising realizations of holographic principle towards quantum gravity. The recent development of a discrete version of AdS/CFT correspondence in terms of tensor networks motivates one to simulate and demonstrate AdS/CFT correspondence on quantum simulators. We achieve this goal indeed, in this work, on a six-qubit nuclear magnetic resonance quantum simulator. We demonstrate experimentally the discrete AdS/CFT correspondence, under realistic noises, by measuring the relevant entanglement entropies on the corresponding tensor network state. The fidelity of our experimentally prepared tensor network state is 85.0% via full state tomography and rises to 93.7% if the signaldecay due to decoherence is taken into account. Our experiment serves as the basic module of simulating more complex tensor network states that exhibit 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. PACS numbers:

Anti-de Sitter/conformal field theory (AdS / CFT) correspondence is one of the most prominent approaches towards a quantum theory of gravity for over two decades [1, 2]. 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 actually equivalent to a conformal field theory on the lower-dimensional boundary of the spacetime. Beyond its applications in quantum gravity and unified theory of all fundamental interactions, AdS/CFT is also widely used in strongly coupled quantum systems, quantum chromodynamics, and condense matter physics [3–5]. Furthermore, it has recently become a bridge connecting quantum gravity to quantum information theory [6, 7], which inspires revolutionary ideas of developing quantum gravity using the methods in quantum information and entanglement. Within this new framework, quantum gravity can be understood via exactly solvable quantum models. For instance, a black hole is related to quantum scrambling in the chaotic systems [8–10] and studied in the Sachdev-Ye-Kitaev model [11–13]; the spacetime geometry emerges from entanglement [14, 15], which can be realized by tensor networks and quantum error correcting codes [16–19]. The holographic duality of AdS/CFT offers a dictionary between the observables of the d-dimensional bulk gravity theory and those of the (d−1)-dimensional boundary field theory, i.e. properties of the bulk gravity and geometry can be reconstructed or emerged from the boundary field theory (known as the emergent gravity program [15]), and vice versa. In particular, the bulk geometry can be reconstructed holographically

from the entanglement in the boundary field theory, characterized by the Ryu-Takayanagi (RT) formula S EE (A) =

Armin , 4G N

(1)

which identifies the entanglement entropy S EE (A) of a (d−1)dimensional boundary region A with the area Armin of the bulk (d − 2)-dimensional minimal surface anchored to A [14, 20, 21]. G N is the Newton constant. See Fig. 1(a) for a brief illustration. Recently, a discrete version of AdS/CFT is realized on a type of lattices called tensor networks (TN) [17, 18, 22, 23], making it possible to be demonstrated on a quantum simulator device in practice. In general, TN states are ways of rewriting a manybody 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 , a 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. In other words, gravity/geometry is anticipated to be created in the emergent bulk dimension of TN. Indeed, this anticipation is fulfilled by the TN made of perfect tensors, which exhibit interesting holographic properties. In particular, perfect tensors (PT) give exactly the holographic entanglement entropy realizing the RT formula [17]. In this work, we demonstrate the RT formula and hence A d S / CFT on a quantum simulator that simulates the minimal PT of rank-6. Using a six-qubit quantum register in the nuclear magnetic resonance (NMR) system, we create the rank-

2 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 A d S / CFT and the holographic principle at large. Perfect tensors – The TN that we focus on is shown in Fig. 1(b), where each hexagon represents a special six-qubit state |ψi. |ψi 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. 1(b). (a)

(b)

-

-

n



 Bulk qubits

(c)

2

2

3 1

1

3

6

4

4 5 6

5 Boundary qubits

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 |ψi ∈ H ⊗6 , and the collection of all such nodes corresponds to the tensor product of all |ψi’s. Each √ link ` represents a maximally entangled state |`i = (|00i + |11i) / 2. 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.

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 ' C2 , a rank-6 PT |ψi 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 P by |αi and |ii respectively, a PT |ψi = α,i ψαi |αi ⊗ |ii satisfies X

ψ† iα ψα j =

α

1 δi j . 2k

(2)

In other words, the reduced density matrix ρ(k) by tracing out m qubits is an identity matrix, whose entanglement entropy   (k) (k) S (k) EE = −tr ρ log2 ρ

(3)

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 [24]) in hand, the TN state illustrated in Fig. 1(b) is constructed as follows. Each internal link ` represents √ a two-qubit maximally entangled state |`i = (|00i + |11i) / 2, where two qubits associate respectively to the two end points of `. If we denote by |ψ(n)i the PT associated to the hexagon node n, the total TN state |Ψi in Fig. 1(b) is written as a (partial) inner product form O O |Ψi = h`| |ψ(n)i. (4) `

n

The inner product takes place at the end points of each internal link `, between one qubit in |`i and the other in |ψ(n)i. The qubits in |ψ(n)i not participating the inner product are boundary qubits corresponding to the dangling legs, and these boundary ones are actually physical qubits, indicating that |Ψi 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. 1(b). The reduced density matrix ρA = trA¯ (|ΨihΨ|) 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 |`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 [24]. 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. 1(b), 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 [24] 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 ρ2A , involving the inner product of boundary qubits in A, gives that ρ2A ∝ ρA . Note that we have ignored all numerical prefactors but they all cancel when trρn calculating (trρAA)n in the entanglement entropy. As a result, the

3 (a)

Von Neumann entropy gives [24] trρnA

1 log2 n→1 1 − n (trρA )n = minimal number of cuts by S.

S EE (A) = lim

(5)

The above result is a discrete version of the RT formula in Eq. (1). The “minimal number of cuts” in Eq. (5) represents the minimal area Armin (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. Eq. (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. 1(c) by which the minimal number of cuts is three, a rank-6 PT is generated where the boundary and bulk qubits are both three. In this work, we demonstrate the emergent gravity program in A d S / CFT for the first time in a six-qubit NMR quantum simulator, by creating the rank-6 PT in Fig. 1(c) and measuring the relevant entanglement entropies. Now we turn to the experimental part. 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. 2(a), in 13 C-labeled Dichloro-cyclobutanone dissolved in d6 -acetone. All experiments were carried out on a Bruker DRX 700 MH z spectrometer at room temperature. The internal Hamiltonian of this system is Hint =

6 X j=1

πν j σzj +

6 X π J jk σzj σkz , 2 j 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 [24]; 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. 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 [25] on the prepared state and got ρe . The real part of ρe is depicted in the right panel of Fig. 4(a), by projecting each element onto a two-dimensional plane. As a comparison, the figure of the theoretical PT ρ pt = |ψihψ| is placed in the left panel of Fig. 4(a). 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  q   √ √ F ρ pt , ρe = tr ρ pt ρe ρ pt , (7) is about 85.0%. Direct observations of ρe in terms of NMR

4 5

S iden�ty S PT S exp S exp_nodec

3

�ty Iden

S

(k)

4

2

PT

1 1

2

3

k

4

5

FIG. 3: Entanglement entropy S (k) of the k-qubit subsystem of the rank-6 PT. In theory, S (k) = min{k, 6 − k} as shown by the orange dashed line. 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.

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. Our results further suggest that a larger quantum processor would be able to simulate higher dimensional RT formulae and hence the AdS/CFT correspondence, under either of the two assumptions: the decoherence effect during experiments can be tolerated, or the original TN state is supposed to be pure. The latter has been in fact tested and elaborated in detail in [28].

(a)

1 2/3 1/3 0 -1/3 -2/3 -1

10

(b) NMR signal (a.u.)

spectra are also shown in Fig. 4(b), 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. 3 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. 3 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 [26, 27] 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 . Discussion – RT formula, or explicitly the TN built by the rank-6 PT in Fig. 1(b), 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. 3. 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

1

C1

exp

sim

0 -1

-80 1

-60

-40

-20

0

20

40

C4

60 exp

80 sim

0 -1

-80

-60

-40

-20

0

20

40

60

80

Frequency (Hz)

FIG. 4: (a) Density matrices of the theoretical rank-6 PT ρ pt (left) and the experimentally reconstructed state ρe (right) on a twodimensional plane. The rows and columns are labeled by the sixqubit computational basis from |0i⊗6 to |1i⊗6 , respectively. (b) Direct observation of ρe in the NMR spectra (red), with probe qubits C1 (top) and C4 (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.

Conclusion – Our work is the first endeavor to demonstrate

RT formula in A d S / CFT correspondence using the discrete per-

fect tensor version on a quantum simulator. 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 a spin ensemble 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.

5 Acknowledgements – We thank Jianxin Chen, Markus Grassl, Cheng Guo, Ling-Yan Hung, Zhengfeng Ji, Hengyan Wang, and Nengkun Yu for helpful discussions. This research was supported by CIFAR, NSERC and Industry of Canada. K.L. and G.L. acknowledge National Natural Science Foundation of China under Grants No. 11175094 and No. 91221205. MH acknowledges support from the US National Science Foundation through grant PHY-1602867, and Start-up Grant at Florida Atlantic University, USA. YW thanks the startup grant offered by the Fudan University and the hospitality of IQC and PI during his visit, where this work was partially conducted.

∗ † ‡ §

[1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15]

[16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26]

These authors contributed equally to this work. Electronic address: [email protected] Electronic address: [email protected] Electronic address: [email protected] J. M. Maldacena, Int. J. Theor. Phys. 38, 1113 (1999). O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri, and Y. Oz, Phys. Rept. 323, 183 (2000). U. Gursoy and E. Kiritsis, J. High Energy Phys. 02, 032 (2008). S. Sachdev, Lect. Notes Phys. 828, 273 (2011). H. Liu, J. McGreevy, and D. Vegh, Phys. Rev. D 83, 065029 (2011). J. Maldacena and L. Susskind, Fortsch. Phys. 61, 781 (2013). A. R. Brown, D. A. Roberts, L. Susskind, B. Swingle, and Y. Zhao, Phys. Rev. Lett. 116, 191301 (2016). J. Maldacena, S. H. Shenker, and D. Stanford, J. High Energy Phys. 08, 106 (2016). P. Hosur, X.-L. Qi, D. A. Roberts, and B. Yoshida, J. High Energy Phys. 02, 004 (2016). D. A. Roberts and B. Yoshida, (2016), arXiv:1610.04903 . A. Kitaev, in Talks at KITP, April 7, 2015 and May 27, 2015. S. Sachdev and J. Ye, Phys. Rev. Lett. 70, 3339 (1993). J. Maldacena and D. Stanford, Phys. Rev. D 94, 106002 (2016). S. Ryu and T. Takayanagi, Phys. Rev. Lett. 96, 181602 (2006). M. Van Raamsdonk, in Proceedings, Theoretical Advanced Study Institute in Elementary Particle Physics: New Frontiers in Fields and Strings (TASI 2015): Boulder, CO, USA, June 126, 2015 (2017) pp. 297–351. A. Almheiri, X. Dong, and D. Harlow, J. High Energy Phys. 04, 163 (2015). F. Pastawski, B. Yoshida, D. Harlow, and J. Preskill, J. High Energy Phys. 06, 149 (2015). P. Hayden, S. Nezami, X.-L. Qi, N. Thomas, M. Walter, and Z. Yang, (2016), arXiv:1601.01694 . M. Freedman and M. Headrick, Commun. Math. Phys. 352, 407 (2016). A. Lewkowycz and J. Maldacena, J. High Energy Phys. 08, 090 (2013). X. Dong, A. Lewkowycz, and M. Rangamani, J. High Energy Phys. 11, 028 (2016). B. Swingle, Phys. Rev. D 86, 065007 (2012). M. Han and L.-Y. Hung, Phys. Rev. D 95, 024011 (2017). See Supplemental Information for a detailed description of the theory and experiment . G. M. Leskowitz and L. J. Mueller, Phys. Rev. A 69, 052302 (2004). L. M. Vandersypen, M. Steffen, G. Breyta, C. S. Yannoni, M. H.

Sherwood, and I. L. Chuang, Nature 414, 883 (2001). [27] D. Lu, H. Li, D.-A. Trottier, J. Li, A. Brodutch, A. P. Krismanich, A. Ghavami, G. I. Dmitrienko, G. Long, J. Baugh, and R. Laflamme, Phys. Rev. Lett. 114, 140505 (2015). [28] T. Xin, D. Lu, J. Klassen, N. Yu, Z. Ji, J. Chen, X. Ma, G. Long, B. Zeng, and R. Laflamme, Phys. Rev. Lett. 118, 020401 (2017).

Supplemental Information for “Measuring Holographic Entanglement Entropy on a Quantum Simulator” A.

Perfect tensor

A n-qudit state |ψn  is perfect if for any bipartition, the entropy of the reduced state is maximal. For example, in the case of n = 2, any maximally entangled state is perfect. Moreover, in a two-level (qubit) system, the perfect tensor state only exists in the Hilbert space of two, three, five, and six qubits. Any GHZ state is perfect for n = 3; any state of the five-qubit error correcting code is perfect for n = 5 [1]; for n = 6, the perfect state is unique, up to local unitaries (LU) and permutations [2]. One representation of the six-qubit perfect state is given by the graph state with stabilizers in Table. S1. It can also be expanded in terms of computational basis or Pauli basis only with a difference of LU or permutation as shown in Table. S1. This state is just what we have prepared in experiment. B.

Derivation of the RT formula

Fig. S1 is a graphical presentation of Eq. (2) being the defining property of the perfect tensor (PT). In Fig. S1, the white hexagon node with six legs represents a rank-6 PT  |ψ = ψa1 ···a6 |a1  ⊗ · · · |a6  (1)

=

=

= FIG. S1. Graphical presentation for the properties of PTs. Each figure on the LHS represents a contraction of two rank-6 PTs: an internal line between two PTs implies a common index of the two tensors that is summed over all its possible values.  Each line on |aa| in the the RHS represents a maximally entangled state a corresponding 2-qubit Hilbert space, where each qubit sits at an end of the line.

region 

a1 ,··· ,a6

in H , where |ai  ∈ H is an orthonormal basis. The black hexagon represents the Hermitian conjugate ψ|. Connecting a pair of legs between two hexagon nodes means contracting the a pair of tensor indices between ψa1 ···a6 and its conjugate. Equivalently, it means to have inner product between |ψψ| and the maximally entangled state ai |ai ai | in the two-qubit H ⊗ H∗ space (or identity operator in H):  ai ai |ψψ|a i . Each line on the right is an maximally entangled state a |aa|. Given the tensor network state |Ψ, we pick a boundary region A which collects a subset of the boundary qubits. For convenience, the region A we choose here is the same as the one in Fig. 1b in the main text. The reduced density matrix ρA = trA¯ (|ΨΨ|), tracing out all boundary qubits outside A, can be computed graphically as shown in Fig. S2, by using the rules in Fig. S1. Initially, this partial trace boils down to contracting the indices of tensors closest to the boundary. By applying Fig. S1, the tensor contractions are effectively pushed from the boundary into the bulk (see Fig. S2). Once again, we can apply Fig. S1 and push the contractions further inside. This iteration procedure is repeated until the contractions reach S (the red dashed line in Fig. S2), where Fig. S1 cannot be applied anymore, as the number of qubits participating the trace (number of links cut by S) is less than three for each tensor. Similarly, ρ2A is computed in Fig. S3, which clearly gives 2 ρA ∝ ρA by the contractions of tensors. tr(ρ2A ) and tr(ρA ) ⊗6

=

region  region 

region 

=

region  region 

= ρA

FIG. S2. Graphical computation of the reduced density matrix ρA . In the first step, all boundary qubits outside A are traced out following the rules in Fig. S1. Then, the remaining two hexagon nodes in the middle are further traced out, indicating that the contractions can be pushed inside.

are computed in Fig. S4. Note that in the graphical computation, we have ignored all numerical prefactors. However, trρn the numerical prefactors we ignore all cancel in (trρAA)n when computing the entanglement entropy. As a result, the second Renyi entropy S2 (A) trρ2 − ln2 (trρAA)2 gives S2 (A) = Minimal Number of Cuts by S.

=

(2)

which equals to 3 in the example shown in the graphical computation (dim(H) = 2). The n-th Renyi entropy can be computed similarly, and gives the same result.

2

  

 

   



 

PT

000000  000001  000010

"

     

000011  000100



 

 







 









 



















000110

000111  001000

001001

010010  010011

 010100

011000  011001  011010  011011

011100

010101  010110  010111

101000  101001

100011

011101

100100  100101  100110  100111

101010

101011  101100  101101

101110  101111

110011  110100

110101  110110  110111

111000  111001  111010  111011

 111100  111101  111110

  



001110  001111  010000  010001

 110010

"



 001010  001011  001100  001101

011110  011111  100000  100001  100010

PT

000101

 

110000

110001

111111

IIIIII+IIXYYZ-IIYXXY+IIZZZX+IXIYZY+IXXIXX-IXYZYI+IXZXIZ-IYIXYX-IYXZIY+IYYIZZ IYZYXI +IZIZXZ+IZXXZI-IZYYIX+IZZIYY+XIIZYY+XIXXIX-XIYYZI+XIZIXZ+XXIXXI XXXZZZ-XXYIIY-XXZYYX-XYIYIZ-XYXIYI-XYYZXX+XYZXZY+XZIIZX-XZXYXY +XZYXYZ+XZZZII-YIIYXX-YIXIZY+YIYZIZ-YIZXYI-YXIIYZ-YXXYII-YXYXZX+YXZZXY +YYIZZI-YYXXXZ+YYYYYY+YYZIIX-YZIXIY+YZXZYX +YZYIXI+YZZYZZ+ZIIXZZ +ZIXZXI-ZIYIYX+ZIZYIY+ZXIZIX-ZXXXYY+ZXYYXZ+ZXZIZI-ZYIIXY+ZYXYZX+ZYYXII +ZYZZYZ+ZZIYYI+ZZXIIZ+ZZYZZY-ZZZXXX

TABLE S1. Form of the rank-6 PTS in three representations: the graph state with stabilizers, computational basis, and Pauli basis. The normalization factor is ignored.

region 

ρ2A =

region  region 

region  region 

=

=

dim(H3)

region 

=

region  region 

FIG. S3. Graphical computation of Fig. S2 and the rules in Fig. S1.

C.

ρ2A

based on the result of ρA in

Preparation of the rank-6 perfect tensor state

Our six-qubit processor is per-13 C labeled (1S,4S,5S)7,7-dichloro-6-oxo-2-thiabicyclo[3.2.0]heptane-4-carboxylic acid dissolved in d6-acetone as shown in Fig. S5. It consists of seven 13 C and five 1 H spin-1/2 nuclear spins. The 1 H channel are decoupled throughout the experiments. The experimental preparation of the rank-6 PT is based on the temporal averaging method. As shown in Table. S1, the PT consists of 64 Pauli terms with the same amplitudes when

FIG. S4. Graphical computation of tr(ρ2A ) and tr(ρA ).

decomposed in the Pauli basis. Therefore, we ran 64 independent experiments to get each individual Pauli. For each Pauli term out of 64, we designed a specific circuit mainly according to its Pauli weight. Without loss of generality, let us take the circuits in Fig. S6 which can drive the system from the thermal equilibrium state to IZZIXX and ZZZXXX as examples. Each circuit can in fact be divided into three steps: (a) Eliminate all other 13 C signals but C4 . At room temperature, the thermal equilibrium state of this system is highly mixed, with the form ρeq =

6  1− I +  Zi , 26 i=1

(3)

where  ≈ 10−5 describes the polarization, and I is a 26 × 26

3 C1

C2

C3

C4

C5

C6

C-13 labeled 6-qubit system Dichloro-cyclobutanone

C1

8779

C2

57.58

30020

C3

0.30

0

10333

C4

37.48

1.25

29.02

11928

C5

-1.66

5.54

-3.53

34.57

C6

C3

C5

C4 C1

34381

C6

2.62

-1.25

33.16

21.75

33.16

15745

T1

3.61

7.99

3.72

3.64

7.80

9.89

T2

0.31

0.40

0.25

0.38

0.40

0.40

*

C2

FIG. S5. Molecular structure and Hamiltonian parameters of per-13 C labeled (1S,4S,5S)-7,7-dichloro-6-oxo-2-thiabicyclo[3.2.0]heptane4-carboxylic acid. The diagonal elements are the chemical shifts (in Hz), and the off-diagonal elements are the J-couplings between two spins (in Hz). The relaxation times T1 and T2∗ (in seconds) are also listed at bottom. All parameters are obtained on a Bruker DRX 700 MHz spectrometer at room temperature.

1/2J46

(a) C1

1/2J12

1/2J34

1/4J12

1/4J14

C2

1/4J12

1/4J12

C3

1/4J34

1/4J34

C4

1/4J34

1/4J46

C5

1/4J34

1/4J45

C6

1/4J46

(b)

1/2J46

C1

1/2J12

1/4J12

The final rank-6 PT state was obtained by the summation over all the above 64 Pauli terms.

D.

Six-qubit full state tomography

To verify the goodness of the prepared PT, we performed the six-qubit quantum state tomography. The density matrix was reconstructed through 110 readout pulses [5], where all of them are local. All the readout pulses were realized by 1 ms GRAPE pulses with fidelity over 99.5%. The state fidelity between the experimentally prepared PT and its theoretical form is about 85.0%, as shown in the main text.

E.

1/4J34

C4

(c) Evolve to the target Pauli term. Once we have the diagonal multi-coherence states, local π/2 rotations were applied on the target spins to generate the desired Pauli terms. For instance, in Fig. S6, we applied Ry5,6 (π/2) and Ry4,5,6 (π/2) on IZZIZZ and ZZZZZZ to generate the target Pauli IZZIXX and ZZZXXX, respectively.

1/4J12

1/4J14

C2 C3

rotations plus free evolutions of the internal Hamiltonian to drive the system from IIIZII to the diagonal multi-coherence state such as IZZIZZ and ZZZZZZ in Fig. S6. These single-qubit operations include π and π/2 rotations, and were optimized by the gradient ascent pulse engineering (GRAPE) technique [3, 4]. Each pulse has a 1 ms duration with the fidelity over 99.2%. Moreover, to minimize the discrepancy between the simulated shape and its real performance in experiment, each pulse was rectified via a feedback control technique by a pickup coil attached to the spectrometer.

Decoherence simulation

1/4J46

C5

1/4J45

C6

1/4J46

Ry (90s)

R y (90s)

Selective

Pulse

Gradient Pulse

FIG. S6. Pulse sequences that generate (a) IZZIXX and (b) ZZZXXX from the thermal equilibrium state, respectively. The sequence contains local rotations, free evolutions, and z-gradient pulses to destroy the unwanted coherent terms. The free evolution time 1/(2J) is in the magnitude of 10 ms, determined by the J-coupling strengths of the relevant spins.

identity matrix. We used a collective 1 ms π/2 pulse to rotate all 13 C spins except C4 to the transverse plane, and applied a z-gradient pulse consequently. The remained state now is IIIZII . (b) Generate multi-coherent state. We used single-qubit

To numerically simulate the decoherence effect in our sixqubit system, we made the following assumptions: the environment is Markovian; only the T2∗ dephasing mechanism is taken into account since T1 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%.

4

[1] R. Laflamme, C. Miquel, J. P. Paz, and W. H. Zurek, Phys. Rev. Lett. 77, 198 (1996). [2] E. M. Rains, IEEE Trans. Inf. Theory 45, 266 (1999).

[3] N. Khaneja, T. Reiss, C. Kehlet, T. Schulte-Herbrggen, and S. J. Glaser, J. Magn. Reson. 172, 296 (2005). [4] C. A. Ryan, C. Negrevergne, M. Laforest, E. Knill, and R. Laflamme, Phys. Rev. A 78, 012328 (2008). [5] G. M. Leskowitz and L. J. Mueller, Phys. Rev. A 69, 052302 (2004).

10