## 3 Results  
### Scenario 1  
In our first scenario, we show how to superactivate bound entanglement by obtaining distillable two-qubit entanglement between qubits \(D\) and \(E\), as in the original superactivation work [13] via implementing a simple QZD procedure, which does not require teleportation or Bell measurements. Through numerical solutions, we find that choosing \(J_{1}=|1\rangle\langle 1|\otimes|1\rangle\langle 1|\) with \(J_{0}=I^{\otimes 2}-J_{1}\) and \(\theta=\pi/180\) for parties [PERSON], [PERSON], and [PERSON] in all the iterations suffice to achieve the superactivation. As we will discuss in the next section, we do not claim that this choice of parameters is the best, and choosing different values in each iteration and for each party found by intelligent techniques potentially increases the performance of QZD and decreases the number of iterations significantly, which we consider for future work. In this work, we aim to present a proof of concept that due to the power of QZD, even with an arbitrary choice of parameters, bound entanglement can be superactivated by only single-qubit operations and simple threshold measurements.  
With \(J_{0}\) and \(R(\pi/180)\), we implement QZD by repeating the \(\rho\to\rho^{r}\to\rho^{rmrm}\) step and denote the obtained state as \(\rho^{n}\) after the \(n\)th iteration. In each iteration, each of the first six qubits are measured on a \(z\)-basis, and the results are communicated to [PERSON] and [PERSON], potentially leaving the system \(\rho_{DE}^{n}\) of spatially separated and never interacted qubits \(D\) and \(E\) entangled, which is detected by calculating the logarithmic negativity using Equation (11). Although \(LN(\rho_{DE}^{n})\) becomes non-zero after around 70 iterations, showing the superactivation of bound entanglement, we continue the iteration to obtain a more entangled state \(\rho_{DE}\).  
Tracing out the first six qubits from the initial state \(\rho\) in Equation (4), i.e., before starting the QZD procedure, \(\rho_{DE}^{0}\) is left in the maximally mixed state, and if the first six qubits are measured on a \(z\)-basis and are all found in the \(|0\rangle\) state, for example, \(\rho_{DE}^{0}\) is left in the \(|0\rangle\otimes|0\rangle\) state. If a single iteration of QZD is performed with the same \(z\)-measurement results, \(\rho_{DE}^{  
providing the proof of concept for the superactivation of BE via QZD. If QZD is executed for 76 iterations, \(LN(\rho_{DE}^{76})\approx 0.14\) is achieved with state  
\[\rho_{DE}^{76}=\left(\begin{array}{cccc}0.00369126&0.0264018&0.0264018&0.00317 834\\ 0.0264018&0.495455&0.0547474&0.023323\\ 0.0264018&0.0547474&0.495455&0.023323\\ 0.00317834&0.023323&0.023323&0.00539837\end{array}\right). \tag{15}\]  
We find that after 3345 iterations, \(LN(\rho_{DE}^{3345})\approx 0.8\) can be achieved where  
\[\rho_{DE}^{3345}=\left(\begin{array}{cccc}0.470202&-0.0115331&-0.0115331&0. 36229\\ -0.0115331&0.00220691&0.00192791&-0.0257913\\ -0.0115331&0.00192791&0.00220691&-0.0257913\\ 0.36229&-0.0257913&-0.0257913&0.525385\end{array}\right). \tag{16}\]  
For a better visualization, we present the bar plots of elements of \(\rho_{DE}\) after 1, 40, and 76 iterations in Figure 3 and present the logarithmic negativity results of our numerical simulation in Figure 4 for up to 3500 iterations.