- Step-by-Step Derivation
- Noiseless target state produced by the circuit
Apply the gates from right to left to |0000⟩:
- H2 turns qubit 2 into |+⟩: |+⟩2|0⟩1|0⟩0|0⟩3.
- CNOT21 (control 2, target 1) entangles (2,1) into a Bell pair |Φ+⟩21:
|Φ+⟩21 = (|00⟩21 + |11⟩21)/√2.
- H0 turns qubit 0 into |+⟩.
- CNOT03 (control 0, target 3) entangles (0,3) into |Φ+⟩03.
Hence the ideal final state is
\[
|\psi_{\mathrm{id}}\rangle = |\Phi^+\rangle_{03}\otimes|\Phi^+\rangle_{21}.
\]
- Stabilizer of the target state
A Bell state |Φ+⟩ on a pair (a,b) is stabilized (up to phase) by {I⊗I, X_aX_b, Y_aY_b, Z_aZ_b}.
Therefore
\[
\mathcal{S}_{03} = \langle X_0X_3,\;Z_0Z_3\rangle = \{II,\;X_0X_3,\;Y_0Y_3,\;Z_0Z_3\},
\]
\[
\mathcal{S}_{21} = \langle X_2X_1,\;Z_2Z_1\rangle = \{II,\;X_2X_1,\;Y_2Y_1,\;Z_2Z_1\}.
\]
The total stabilizer (modulo overall phases) is the tensor product
\[
\mathcal{S} = \mathcal{S}_{21}\otimes\mathcal{S}_{03}.
\]
- Noise model and where it acts
After each CNOT, a two-qubit depolarizing channel acts on the same two qubits:
\[
\mathcal{E}(\rho)=(1-p)\rho+\frac{p}{15}\sum_{P\in\mathcal{P}_{2}\setminus\{II\}}P\rho P,
\]
where \(\mathcal{P}_2\setminus\{II\}\) are the 15 non-identity 2-qubit Paulis on that pair.
In our circuit:
- After CNOT21: noise \(Q_{21}\) acts on qubits (2,1).
- After CNOT03: noise \(Q_{03}\) acts on qubits (0,3).
These two channels act on disjoint pairs and are independent, and the operations between them (H0 and CNOT03) do not affect qubits (2,1). Hence the final error operator is simply
\[
E_{\mathrm{tot}}=Q_{21}\otimes Q_{03}.
\]
- Fidelity for a stabilizer state under Pauli noise
For a pure stabilizer state \(|\psi_{\mathrm{id}}\rangle\), and a Pauli error \(E\), the overlap satisfies
\[
1,& E\in\mathcal{S}\\
0,& E\notin\mathcal{S}
\]
Therefore the physical state fidelity is the probability that the net Pauli error lies in the stabilizer:
\[
F=\Pr[E_{\mathrm{tot}}\in\mathcal{S}] = \Pr[Q_{21}\in\mathcal{S}_{21}]\;\Pr[Q_{03}\in\mathcal{S}_{03}],
\]
where we used independence and the fact that \(\mathcal{S}=\mathcal{S}_{21}\otimes\mathcal{S}_{03}\).
- Counting “good” Pauli errors on each pair
On any pair, the depolarizing channel outputs:
- \(II\) with probability \((1-p)\),
- each of the 15 non-identity two-qubit Paulis with probability \(p/15\).
Among the 15 non-identity Paulis, exactly 3 are in the Bell stabilizer: \(XX,\ YY,\ ZZ\). Therefore
\[
\Pr[Q_{\text{pair}}\in\{II,XX,YY,ZZ\}] = (1-p) + 3\cdot\frac{p}{15} = 1-\frac{4p}{5}.
\]
The two pairs are independent, so
\[
F(p) = \left(1-\frac{4p}{5}\right)^2.
\]
-
Mathematical Typesetting
- All equations above are typeset in LaTeX.
-
Conventions and Units
- Standard qubit and Pauli operator conventions; no dimensional units apply.
-
Final Answer:
F(p) = (1 - 4p/5)^2 = 1 - (8/5)p + (16/25)p^2.