A scheme to create and verify scalable entanglement in optical lattice

To achieve scalable quantum information processing, great efforts have been devoted to the creation of large-scale entangled states in various physical systems. Ultracold atom in optical lattice is considered as one of the promising platforms due to its feasible initialization and parallel manipulation. In this work, we propose an efficient scheme to generate and characterize global entanglement in the optical lattice. With only two-layer quantum circuits, the generation utilizes two-qubit entangling gates based on the superexchange interaction in double wells. The parallelism of these operations enables the generation to be fast and scalable. To verify the entanglement of this non-stabilizer state, we mainly design three complementary detection protocols which are less resource-consuming compared to the full tomography. In particular, one just needs two homogenous local measurement settings to identify the entanglement property. Our entanglement generation and verification protocols provide the foundation for the further quantum information processing in optical lattice.


INTRODUCTION
Quantum information and quantum computation 1 , which harvest the intrinsic quantum features, like superposition and entanglement 2 , can show advantages against their classical counterparts.
To build a practical quantum information processor, the computation platform with high scalability is preferred and the qubits should be coupled to each other to form a large-scale entanglement. Hence, many researches have been focusing on the generation of scalable entangled states in different physical systems, e.g, ion trap [3][4][5] , photons 6,7 , Rydberg atoms 8 , and superconducting circuits [9][10][11] . Even though there are significant progresses of the qubit number in various systems 12,13 , the generation of large-scale entanglement is still challenging for the Noisy Intermediate Scale Quantum Devices 14 .
Ultracold atoms in optical lattice 15 could be a practicable system to overcome this challenge due to its feasible initialization and parallel manipulation. By adiabatically increasing the lattice depth, the phase of ultracold atom can be tuned from superfluid (SF) to Mott insulator (MI) 16,17 . Under an unit filling rate, numerous atoms can be confined in the lattice and serves as qubits. Based on this initialization, entangled states in optical lattice have been demonstrated experimentally, for instance, the generation of cluster state with controlled collision gate induced by the spindependent lattice 18,19 , which acts as the resource state for the measurement-based quantum computing 20,21 . The development of superlattice further improves the ability to control ultracold atoms, which is formed by overlapping two different optical lattices to generate a series of double wells. The structure of the double well can be modified to induce different kinds of atomic dynamics, such as superexchange coupling 22 and controlled exchange interaction 23 , which can be both used to realize ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p gate and entangle two atoms in a double well 23,24 . Besides, the entangling operations in the superlattice can be performed in a parallel way based on the periodicity of the lattice system, which is suitable for the fast generation of large amount of entangled pairs and even large-scale entangled states. Note that this kind of parallel operation can also be implemented with multitweezer in Rydberg atom experiments 25,26 .
Every coin has two sides. The periodicity property also induces some restrictions on the quantum operation and measurement. The small lattice spacing, required by the large tunnelling between neighbour sites, creates a challenge for the individual control of atoms, say, some local basis rotations. The tight-focused optical tweezer 27 created with the recently developed highresolution imaging [28][29][30][31][32] could be a solution. However, it is still challenging to perform a few of different single-qubit operations on multiple qubits under a realistic time-scale to ensure the system coherence. As a result, homogeneous operations and measurements on all atoms are preferred in optical lattice experiments.
In this work, we propose a scheme to generate scalable entanglement of ultracold atoms which is suitable for the implementation in the optical superlattice system. It mainly contains two entangling steps: first, entangle the atom-pair in each double well by the ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p gate; second, shift the position of double wells to a single site by changing the phase of superlattice, and then entangle the the new atom-pair with ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p again. In this way, all atoms in the superlattice can be connected with neighbours to form a global entangled state. Our theoretical analysis shows that the final state possesses genuine multipartite entanglement (GME). In addition, the final state is also less sensitive to magnetic noises which can cause decoherence, since it only owns amplitude on the computational basis whose total spin is zero.
In actual experiments, the inevitable noise may degrade the entanglement, which should be verified further. Compared with quantum tomography 33,34 , entanglement witness is a more efficient way [35][36][37] to realize this task based on the pre-knowledge of the preparation. Current entanglement witnesses are usually designed for some structured states, such as permutationinvariant states [38][39][40] and stabilizer states [41][42][43] . However, in our protocol the state generated by the non-Clifford gate ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p is not a standard stabilizer state, which makes the entanglement verification challenging. To overcome this challenge, we first construct an entanglement witness based on the preparation fidelity. By adopting the decomposition using stabilizer-like method, we can lower bound the fidelity with a few spincorrelation measurements. Some correlations among them are inhomogenous and thus require the individual atom addressing whose implementation could be challenging with current techniques. To further ease the experiment realization, we show another complementary protocol which only requires homogenous spin-correlations and only two measurement settings. At last, by reversely evolving the state with the conjugate quantum gates, we provide an intuitive method to indirectly bound the preparation fidelity and thus qualitatively verify the entanglement.

Entangling gates in optical superlattices
The behavior of ultracold atoms in optical lattices can be described by the Hubbard model which is characterized by the tunnel coupling between neighbouring sites J, the on-site interaction V, and the effective chemical potential μ. By increasing the depth of the lattice, a phase transition from SF to MI can occur and the atoms start to be localized in each site. To extend the ability to manipulate atoms, a more complicated periodic potential, named superlattice, was proposed and applied in many experiments [44][45][46][47] . It is normally formed by overlapping two distinct optical lattices. The period of the first lattice (denoted as the long lattice) is twice as that of the second (denoted as the short lattice), which induces an array of double wells steadily. The resulting potential shows V t (x, ϕ) = V l cos 2 (πx/a + ϕ) + V s cos 2 (2πx/a) with a being the lattice spacing of the long lattice. The structure of such potential is dependent on both the relative strength of two lattices and the relative phase ϕ. When the phase ϕ is not equal to nπ/2 (n is an integer), all the double wells are biased with a non-zero tilt Δ between the subsites. In addition, the center of each double well can shift a period of short lattice via changing n by an odd number, which is illustrated in Fig. 1 Subsites in each double well are separated by a low barrier, while neighbour double wells are separated by a high barrier. We denote corresponding tunneling strengths as J inter and J inner respectively. As J inter is much smaller than J inner , the hopping event between different double wells can be ignored and the movement of atoms is restricted in each double well. In this case, for two-level bosonic atoms, their dynamics in double wells can be described well with a two-site mode Bose-Hubbard Hamiltonian characterized with V, J inner and Δ as follows.
i¼R;L ðn i;"ni;# þn i;#ni;" Þ þ Δ 2 ðn L;" Àn R;" Àn L;# þn R;# Þ: (1) Here, R(L) denotes right(left) subsite,â y i;σ (â i;σ ) is the creation(annihilation) operator for the boson on site i with the inner state σ, andn i;σ is the number operator. We consider a non-biased double-well system starting with unit filling-one atom per site. Note that the chemical potential μ is fixed here for unit filling and thus has no effect on the dynamics, and it is ignored in Eq. (1). In the limit V ≫ J inner , the large V prevents multiple occupation, so that the system would evolve in a subspace consists of four basis states labeled by inner states of the bosons "; " j i, #; " j i, "; # j i and #; # j i. By using perturbation theory, the above model in this subspace is equivalent to the isotropic Heisenberg Hamiltonian 48 : whereX i ,Ŷ i ,Ẑ i are Pauli operators on subsite i = R, L, and J ex % 2J 2 inner =V is the superexchange coupling between subsites.
In the remaining of this work, we denote spin configuration # j i and " j i as 0 j i and 1 j i. Initialized at 0; 1 j iand driven by this effective Hamiltonian, the state of the system can oscillate between 0; 1 j i and 1; 0 j i with a period of T = 2πℏ/J ex while global phase is also recovered, named the superexchange process. Taking an evolution time of T/8, the product states 0; 1 j i and 1; 0 j i are prepared into two-qubit maximally entangled states 1= ffiffi ffi 2 p ð 0; 1 j iþ i 1; 0 j iÞ and 1= ffiffi ffi 2 p ð 0; 1 j iÀ i 1; 0 j iÞ respectively while 0; 0 j i and 1; 1 j i remain unchanged due to the high energy gap of V. The effective unitary transformation corresponds to a ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p y gate operation, i.e., with respect to the basis 0; 0 j i, 1; 0 j i, 0; 1 j i, 1; 1 j i. Moreover, with an evolution time of 3T/8, one can realize the corresponding ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p gate. In particular, when J inter /J ex ≥ 25, the infidelity of the ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p y operation caused by the tunneling between neighbour double wells would be smaller than 0.1%. For the current gate generation scheme, due to the large ratio V/J inner , the effective J ex is far less than J inner . In addition, such high V/J inner requires a deep potential in each site thus leads to a small J inner . As a result, the period of the superexchange process would be long, e.g, tens of milliseconds for 87 Rb atoms in superlattice 22,24 , which would aggravate the decoherence effect. To ease this problem, an alternative and faster gate generation scheme can be adopted here 49 , which is performed with a small V/J inner . Instead of introducing the large energy ratio, this scheme utilizes the coherent competition of superexchange and atom tunneling to decrease the component with double occupation in state. As V/J inner is set to be a finite value 4= ffiffi ffi 3 p , the undesired component can be eliminated completely at specific time intervals. In particular, with an evolution time of πℏ/V, such elimination leads to a fast ffiffiffiffiffiffiffiffiffiffiffiffi ffi SWAP p y gate realization with high fidelity. Besides, the spin-dependent effect further improves the ability to manipulate the atoms in optical lattices. By adding circular polarization components into one lattice light field of superlattice 45,50 , the tilt Δ can be different for different inner states. Therefore, the energy gap between inner states on right subsite would be different from that on the left which induce two applications here. One is to spin-flip atoms in one subsite of every double well without affecting atoms in the other subsite, enabling the state transfer between the four basis states mentioned above 45 . The other is to generate a relative phase between 0; 1 j i and 1; 0 j i and then transfer the states 1= ffiffi ffi 2 ffiffi ffi 2 p ð 1; 0 j iÀ 0; 1 j iÞ 24 , respectively.

Entanglement generation protocol
Our protocol is expected to be performed on one-dimensional atomic chains along X direction with a short lattice isolating the atoms. On each site, the filling is initialized into unit through the SF-MI phase transition. Assuming the inner state of atoms are prepared into 0 j i and taking 10-qubit system as example, the protocol can be divided into the following steps which are illustrated in Fig. 2. 1. Turn on the spin-dependent superlattice 50 by ramping up a spin-dependent long lattice along X. With selective spin-flip by coupled field, the Néel state is prepared: 2. Turn off the spin-dependent effect and perform ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p y gate on the atom pairs (1, 2), (3, 4),..., (9, 10) in each double well, 3. Turn on the spin-dependent effect again and introduce a relative phase π/2 on 10 j i. As a result, the state of the system is the tensor product of five Bell states, 4. Change the relative phase ϕ by π/2 to trap the atom pair (2, 3), (4, 5),..., (8,9) in different double wells.

Perform
ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p y gates on the atom pairs (2, 3), (4, 5),..., (8,9) in each double well, the final state is During the entire process atomic motions along other directions are frozen by deep lattice potentials. Finally, onedimensional entangled systems are generated. The gate operations entangle every pair of neighbour atoms so that the target state is GME by itself. As the gate operation is not present, the depth of each site should be deep enough to avoid any undesirable hopping events.
The target state only owns non-zero projection on a few bases whose total spin numbers are zero, and the magnetic field fluctuation can only introduce a global phase. Consequently, our target can be robust to the magnetic field fluctuation, and thus has longer coherence time, compared to for example the GHZ state. In this basis collection, each basis can find another basis whose spin on each qubit is anti-parallel to its own, e.g., 0; 1; 1; 0; ::: j iand 1; 0; 0; 1; ::: j i , and the projections on each pair of bases show same modules but different arguments.
According to the definition of GME 2 , it is inferred that a pure state whose subsystems are all mixed states must be a GME state. Based on this, the purity of all subsystems of the target state in the case of 10-qubit are calculated and found to be less Y. Zhou et al. than 1, which verifies that this 10-qubit target state possesses GME. This conclusion can be extended to the target state with a different number of qubits. We remark that Ref. 51 also studies the generation of entanglement by ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p in the cold atom system. The state generated there is like a GHZ state, which needs much more gate operations, compared to the cluster-like state here. In addition, there is no valid method given there to verify the entanglement.
In the following sections, we exhibit theoretical methods to detect the multipartite entanglement in the cold atom system. The prepared state shown in the Section 2.2 is generated by parallel ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p y gates and thus not a stabilizer state, such as the GHZ state and the cluster state. This fact makes the entanglement detection task challenging and the methods proposed in literature are not suitable in this scenario.
To overcome this challenge, we propose three complementary methods. The first one is based on the fidelity estimation to the target state to detect the strongest form of entanglementgenuine multipartite entanglement. By evaluating the non-Pauli stabilizer after the entangling gate evolution, the method only needs constant number of measurement setting with respect to the system size to lower bound the fidelity. To further release the experiment efforts, the second method adopts homogeneous measurements, with only two measurement settings. Thus it is very efficient to realize and can tell whether the prepared state is separable or not considering any bipartition of the whole system. The third one is more intuitional but can reveal GME with less experiment requirements. In particular, we indirectly estimate the fidelity between the prepared state and the target state by evolving the entangling process forward and backward.
All these methods need the spin measurement on single atom, however, the fluorescence imaging used in high-resolution experiments can not resolve different spins. As the state we consider here only has one atom on each site, it could be done by measuring the atom distribution after removing one spin component in which the occupied site and unoccupied site represent different spins 52 or splitting different spin components with gradient along perpendicular direction 53 .
Before showing these three methods, we give some related definitions about multipartite entanglement.
A pure state is (bi-)separable if it is in a tensor product form Ag is a bipartition of the qubits in the system. A mixed state is separable if it can be written as a mixing of pure separable states, There is another restricted way for the extension to mixed states. A state is P 2 -separable, if it is a mixing of pure separable states with a same partition P 2 , and we denote the state set as Since S P2 b & S b , GME is a stronger claim than full entanglement. We also remark that the recently demonstrated entanglement in the IBM cloud quantum computing 54 is actually the full entanglement defined here. By Def. 1, for a state with full entanglement, it is possible to prepare it by mixing bi-separable states with different bipartitions 55 . On the other hand, GME describes the strongest form of quantum entanglement, that is, all the qubits in the system are indeed entangled with each other. GME is essential in various multipartite quantum information tasks, such as quantum cryptography 56 , quantum nonlocality 57 , quantum networks 58,59 , quantum metrology 60 , and measurementbased quantum computing 61 .

Entanglement detection based on fidelity estimation
In this section, we show an entanglement detection protocol based on the fidelity value between the prepared state and the target state. Proposition 1. The operator W Ψ can witness genuine multipartite entanglement near Ψ j i, with hW Ψ i ! 0 for any separable state in S b .
According to Proposition 1, if the fidelity of the prepared state ρ pre with the target state Ψ j i, i.e., Trðρ pre Ψ j i Ψ h jÞ, exceeds 5 8 , ρ pre possesses GME. However, it is generally difficult to evaluate the quantity Trðρ pre Ψ j i Ψ h jÞ by the direct projection on Ψ j i, as it is an entangled state.
Alternatively, one needs to decompose the density operator Ψ into the summation of many local measurements in the form N i¼1 O i , which is easier to implement in experiments. Here O i is Hermitian operator of the ith qubit. The number of local measurements characterizes the experiment effort to estimate the fidelity.
Here, in order to reduce the measurement effort, instead of direct decomposing Ψ, we give a lower bound of the fidelity by using the stabilizer-like operator for the non-stabilizer state Ψ.
Proposition 2. For the target state Ψ j i Ψ h j and its N independent stabilizers S 0 , the following inequality holds: where A≥B indicates that (A−B) is positive semidefinite. The stabilizer is determined by the evolution of ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p gates for k = 1,..., N/2, are the stabilizers for the Bell pairs.
Due to ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p y is not a Clifford gate, the corresponding stabilizer S 0 is not in a tensor product of Pauli operators, but the summation of them. One can directly get S 0 by the evolution of the parallel ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p y on S i . Due to the locality of the evolution, S i can only be transformed by the gate with overlapping support. For example, Both are the linear combinations of 16 Pauli tensors. For other stabilizers, similar forms could be obtained in this way. In total, the number of these Pauli tensors is almost 32N which scales linearly with the qubit number. In fact, the measurement effort can be further reduced, as several Pauli tensors can be grouped and measured in one local measurement setting (LMS) simultaneously. For the Pauli tensors from S 0 ¼ UX 2kÀ1 X 2k U y , we introduce the following expression for the example case N = 10, where one can periodically select one of the three Pauli tensors in every brace to construct one LMS, such that these LMSs cover all the possible Pauli operators from S 0 . That is, select same tensor in (2,3) and (6,7) and the same Pauli operator in (4,5) and (8,9). Thus, only 9 LMSs are needed here. Following the same way, we can also find another 9 LMSs for the Pauli operators from S 0 ¼ UZ 2kÀ1 Z 2k U y . As a result, totally only a constant number of 18 LMSs is needed to obtain all the expectation values of the stabilizers.
To evaluate the robustness of our witness, we apply the white noise model ρ pre ¼ pI=2 N þ ð1 À pÞ Ψ j i Ψ h j, and it shows that the noise tolerance is p ¼ 3 4N . As N = 10, it equals 7.5%. In fact, by utilizing the trade-off between the robustness and the measure budget 62,63 , one can enhance the noise tolerance of the witness further. In Ref. 64 , some of us further generalize the witness here by utilizing more local measurement settings, and the witnesses constructed there would be more suitable to realize in other quantum systems, such as superconducting-qubit. In particular, a numerical algorithm is developed in Ref. 64 to search for the witness with the optimal noise tolerance under given measurement settings. For instance, one can reach the noise tolerance p ¼ 3 16ð1À2 ÀN=2 Þ with about 2⋅3 (N/2−1) settings.

Entanglement detection with homogeneous measurements
The entanglement witness shown in Section 2.3 based on the fidelity estimation can detect GME, with a few number of LMSs. However, in each of the LMS, the Pauli operators may be not the same, such as X 1 Y 2 Z 3 ,.... This kind of inhomogeoneus measurement needs additional local basis rotation, which is challenging for the current cold atom system.
In this section, we reduce the measurement efforts further by only considering the homogenous measurement, say O ⊗N . In particular, we only need two LMSs, X ⊗N and Z ⊗N . Note that X ⊗N and Y ⊗N are symmetric for our generation, thus we only need to measure one of them. Detailed illustration on this point is given in Method.
Instead of detecting GME, the protocol here aims to detect the full entanglement property defined in Eq. (8), i.e., not separable with respective to any bipartition. First, let us show the following Lemma which plays an important role of the detection.
Lemma 1. For a k-qubit quantum state ρ with k being even, if it is P 2 -separable, with P 2 ¼ fA; Ag and the number of qubits contained in A and A being odd, the following inequality holds, jhX k ij þ jhY k ij þ jhZ k ij 1; where hOi ¼ TrðρOÞ is the expectation value.
The proof of Lemma is based on the anticommutative relation 65,66 on subsystems A and A respectively. As k = 2, it becomes to the common criterion for the Bell state. As a result, the violation of the bound in Eq. (16) indicates that the underlying even-qubit state is non-separable regarding any odd-odd bipartitions.
Second, we give the following observation which shows the relation of the entanglement property of the state and its reduced density matrix (RDM).
Observation 1 Suppose a quantum state ρ 0 B;B is created from ρ A;A by local operation and classical communications(LOCC) Λ with respective to the bipartiton fA; Ag, i.e., If ρ A;A is separable for the bipartition fA; Ag, ρ 0 B;B is also separable for the bipartition fB; Bg; in other words, if ρ 0 B;B is entangled for fB; Bg, ρ A;A is also entangled for fA; Ag.
The observation holds as LOCC can not create entanglement. Specifically, the partial trace operation is an example of LOCC, where B; B are subsystems of A; A respectively.
By Definition 1, an N-qubit quantum state is fully entangled, if it cannot been written in the following form: for all possible bipartitions fA; Ag of the whole system. In the following, we employ Lemma 1 and Observation 1 time and time again on RDMs of the 1-D prepared state, and finally certify its full entanglement. The RDM of the prepared state ρ pre , for example, on the first two qubits, is denoted by ρ 12 .
Note that the given expectation values are from the perfect target state. The practical value in the experiment may deviate from it, but can still reveal the entanglement property if the deviation is not too large. Detailed noise analysis are shown in Method.
To certify the full entanglement, we prove it by contradiction. Specifically, we first assume that the state can be written as right hand side of Eq. (18) for some fA; Ag, and then show that all the qubit are either in A or in A, which is actually not separable and leads to the contradiction.
Here we take N = 6 for simplicity, and it can be directly generalized to any even qubit number N. The procedure is listed as follows: 1. For the RDM ρ 12 , one has

Entanglement detection with reverse evolution
Compared with the previous two entanglement detection protocols, the third one is more intuitional and aims to reveal GME. As shown in Section 2.3, it is not easy to obtain the fidelity to the target non-stabilizer state. The first method should use several non-homogenous LMSs to lower bound the fidelity.
In this section, we indirectly estimate the fidelity between the prepared state and the target state by evolving the entangling process backwards. Since the fidelity to both the Bell state and product state can be measured in the experiment, we can indirectly estimate the fidelity to Ψ j i and certify GME according to Lemma 1. The reverse evolution follows the state preparation process shown in Section 2.2.
• Adjust the relative phase between 01 j i and 10 j i in each of the Bell pairs in Φ 4 j i.
In the above reverse protocol, the state is evolved backwards finally to the initial product state. Note that here we list the perfect reverse evolution, that is, Φ 4 j i ¼ Φ 2 j i and Φ 5 j i ¼ Φ 1 j i shown in the generation protocol in Section 2.2. The practical states actually deviate from them due to noises, denoted by ρ ϕ i for i = 1, 2,..., 5, and the fidelity decreases due to the imperfection of the underlying gates.
One can estimate the fidelity to product states Φ 1 j i; Φ 5 j i with the Z ⊗10 measurement setting, entangled-pair states Φ 0 j i; Φ 0 j i with XYXY...XY and Z ⊗10 , and Bell-pair states Φ 2 j i; Φ 4 j i with X ⊗10 and Z ⊗10 . In Fig. 3, we show an illustration. By fitting the fidelity values from the experiment, one can indirectly estimate the fidelity of the prepared state ρ pre to the target state Ψ j i, which can be compared with the theoretical bound for GME, say 5 8 , denoted by the dotted horizon line in the figure. Note that ffiffiffiffiffiffiffiffiffiffiffiffi ffi SWAP p ¼ ffiffiffiffiffiffiffiffiffiffiffiffi ffi SWAP p y3 , i.e., it takes three times of evolution time, and one may needs to consider this point in the fitting. We remark that this entanglement detection method is intuitive, and one could make it more rigorous by assuming more specific noise models.

DISCUSSION
In this work, we propose an experimental scheme to generate and characterize large-scale entanglement of cold atoms confined in the optical superlattice. The generation scheme utilizes the entangling gate induced by the superexchange interaction, and is robust to decoherence. To characterize the entanglement, we propose several complementary methods considering the experimental implementation feasibility. Moreover, it is straightforward to generate our scheme to the high dimension scenarios, and it is also interesting to construct some efficient entanglement verification tools for them. In summary, our entanglement generation and verification protocols are well tailored for the cold atom system, and lay the foundation for the further applications, such as measurement-based quantum computing.

Proof of Proposition 1
Proof. To prove that Eq. (10) is a legal GME witness, we need to show that the maximal fidelity between the separable state and our target state is upper bounded by 5/8, that is Trðρ s ΨÞ 5=8: By the convexity of S b , we can reduce the maximization to the pure separable state where we maximize over Φ s from a given bipartition P 2 and also all possible bipartitions. The second line is due to the fact that the maximization result over Φ s is just the square of the largest Schmidt coefficient of Ψ with respective to P 2

67
, where the Schmidt decomposition In fact, we only need to consider the bipartitions where A contains the first n A and A contains the last n A ¼ N À n A qubits on the 1-D chain. That is, there is only one boundary between A and A. Other choices with more boundaries can lead to a smaller λ 1 . Consider the one boundary scenario, there are two types of bipartations depending on that the boundary is between {2k−1, 2k} or {2k, 2k + 1}, as shown in Fig. 4.
Based on the causality of the quantum gate and the fact that quantum gates inside the subsystem A and A do not change the Schmidt number, we can determine λ 1 for both of cases. For the {2k − 1, 2k} case, λ 1 ¼ 1 2 , Fig. 3 Illustration of reverse revolution. At point 3, the state preparation is accomplished; At point 4 and point 5, the system is under reverse evolution. In the perfect case, the state can return back to Bell pairs and Néel state. Here, due to the noise and decoherence, we denote the corresponding mixed state as ρ ϕ i , and the fidelity is Trðρ ϕ i Φ i Þ to the perfect case. Note that one can estimate the fidelity to Φ 1 , Φ 5 with one measurement settings Z ⊗10 , and the one to Φ 2 , Φ 4 with two settings X ⊗10 , Z ⊗10 .
Here for denotation simplicity we take k = 1 without loss of generality, and ϕ Bell j i¼ 1 ffiffi 2 p 01 j iþ 10 j i. By calculation one finds that the reduced density matrix of the first two-qubit is i.e., the mixture of the Bell state and the maximal mixed state. As a result, ρ 12 is a Bell-diagonal state and thus its four eigenvalues on the Bell basis is In summary, we prove that F max ¼ 5 8 .

Proof of Proposition 2
Proof. We first show that S 0 stabilizes the target Ψ j i by definition. Before the final layer of ffiffiffiffiffiffiffiffiffiffiffiffi SWAP p y gates, the state is the product of a few of Bell pairs, denoted by Ψ Bell j i, which is stabilized by S i defined in Eq. (12). That is, Then we prove the inequality in Eq. (11). We remark that Eq. (11) holds for any generalized stabilizer state. The N independent stabilizers S i commute with each other and their common eigenstates together determine an orthogonal basis denoted by Ψ It is clear that both operators on the left and right of Eq. (11) are diagonal in this basis, thus we can prove it by checking for the matrix elements for every Ψ x !. For x ! ¼ 11 Á Á Á 1, one has 1 ! 1 2 N À ðN=2 À 1Þ ¼ 1. For x ! which contains only one zero, such as x ! ¼ 01 Á Á Á 1, one has 0 ! 1 2 ðN À 1 À 1Þ À ðN=2 À 1Þ ¼ 0. For the x ! containing more zeros, one can prove the inequality similarly.

Symmetry of the target state
Remember that in principle we should apply three measurement settings {X ⊗N , Y ⊗N , Z ⊗N } to detect the full entanglement in Section 2.4. Here we show rigorously as follows that one does not need to worry about the direction on the X−Y plane in every measurement. In fact, one only needs to measure 〈X ⊗N 〉 by the symmetry of the target state.
In the experiment, actually there is no reference pulse to decide the direction of the measurement, that is, one chooses a random angle X θ ¼ cosðθÞX þ sinðθÞY for the measurement. Note that there is a corresponding unitary (rotation around Z), such that X θ ¼ U θ XU y θ . As a result, suppose there is a standard X direction, the actual measurement shows, where ρ is the prepared state in the experiment, and ρ sym R 2π θ¼0 U yN θ ρU N θ dθ. In other words, measuring the state ρ in random direction is equivalent to "twirling" the state, and the resulting state is symmetric with respective to the measurement direction.
Specifically, X = X θ=0 and Y ¼ X θ¼ π 2 , and one also has TrðZ N ρ sym Þ ¼ TrðZ N ρÞ; (31) since Z commutes with U θ defined in Eq. (28). It is not hard to see that the entanglement of ρ sym is not stronger than ρ, since ρ sym is obtained from ρ using LOCC, that is, the twirling can not increase the entanglement. Thus if one can detect the entanglement property of ρ sym , this property should also holds for ρ.
On the other hand, the target state is a symmetric state on X−Y plane, i.e., U N θ ΨU yN θ ¼ Ψ; 8θ. In fact, one can write the rotation unitary as Since our target state Ψ j i only has non-zero projection on the computational bases whose total spin is zero, i.e., half number of 0 and 1, the rotation unitary introduce the same phase for these bases, and the state is unchanged.

Proof of Lemma 1
Proof. Without loss of generality, we assume that the subsystem A contains the first k 1 qubits, and |A| = k 1 , jAj ¼ k À k 1 ¼ k 2 are both odd numbers. Since the left hand of Eq. (16) is a convex function of the state, thus we only need to consider the pure separable state in the form Ψ s j i ¼ Ψ k1 j i Ψ k2 j i, and the expectation value can be written apart as, where the second line is due to Cauchy-Schwarz inequality. Since k 1 is an odd number, one can check that X k1 ; Y k1 and Z k1 anticommute with each other, and thus jhX k1 ij 2 þ jhY k1 ij 2 þ jhZ k1 ij 2 1 65,66 . Similarly one has jhX k2 ij 2 þ jhY k2 ij 2 þ jhZ k2 ij 2 1. As a result, we finish the proof.

Proof of Observation 1
Proof. The Observation is right by the definition of LOCC. Here we show the case for the partial trace operation by contradiction. Suppose ρ A; is still separable, which contradicts with the assumption that it is entangled.

Influence from operation errors
The fidelity of the state creation is affected by the error of each operation, caused by the decoherence effect and noises. For the large-scale entangled state, the degradation of fidelity and also the entanglement detection is more apparent compared to few-body states. In the meantime, the read-out error also affects the entanglement detection. As a result, it is meaningful to analyze the influence of various errors in the experimental implementation.
Here we consider the errors in the following operations: the initial preparation of Néel state; the final projective measurement of spins; and the fidelity of the intermediate quantum gates. We show their influences on the entanglement detection of a ten-qubit system, i.e., N = 10, by using the witness method in Section 2.4, which is more practical for experiments. As shown in the procedure of entanglement detection in Section 2.4, we should measure the witness value in Eq. (16) for a few of reduced density matrices of subsystems and also the whole system. Here we take the subsystem as f1; 2g; f1; 3g; f1; 2; 3; 4g; f1; 2; 3; 5g; f1 À 6g; f1 À 5; 7g; f1 À 8g; f1 À 7; 9g (34) which contain n = 2, 4, 6, 8 qubits. In the first step of Section 2.2, starting from the all-zero state 00 Á Á Á 0 j i , one should flip the spins on the odd sites to prepare the Néel state in Eq. (4). There is the probability of a non-flip for the odd site, and also the probability of a wrong flip for the even site. For simplicity, we assume both error probabilities are equal, and denote the probability of the correct flip as P SF for each site. In this case, the system is initially prepared into a mixture of different product states weighted with corresponding probabilities. In Fig. 5a, we plot the witness values of Eq. (16) for different subsystem size n with respective to P SF . Similarly, in the final measurement, the probability to measure the spin of single atom correctly is denoted as P MS . That is, there is probability 1−P MS to recognize 0 j i as 1 j i or vice versa. The witness value with respective to P MS is shown in Fig. 5b. From these two figures, we found that the witness values are the same for the subsystem with the same qubit number, for example, n = 2 case with ρ 1,2 and ρ 1,3 . When there is no error, say P SF = 1 and P MS = 1, all the witness values return 1.5 for all the subsystems, except the whole system with the value 3. One can see that even though there is some error, the witness can detect the entanglement as the value is larger than 1. For larger subsystems, the values decay faster with P SF and P MS .
At last, in Fig. 5c, we study the influence of the fidelity of the entangling step between the above two steps by taking P SF = 0.98 and P MS = 0.985. Here we assume that this operation has a probability P ES to be performed perfectly while it contributes zero to the witness with probability 1−P ES , that is, essentially outputs a maximal mixed state. As shown in Fig. 5c, the full entanglement could be verified when the fidelity of entangling step exceeds 0.85.