Step-by-Step Derivation
- Ideal (error–free) preparation
• The circuit
\[
\bigl(CNOT_{03}\bigr)\, H_{0}\, \bigl(CNOT_{21}\bigr)\, H_{2}
\]
acting on the initial product state $|0000\rangle_{0123}$ produces
\[
\]
(Qubit ordering is $0,1,2,3$.)
The encoded state $|\psi_{\text{ideal}}\rangle$ lies in the
$+1$–eigenspace of the stabilizers $XXXX$ and $ZZZZ$ and represents the logical
two-qubit GHZ state
$(|00\rangle_{AB}+|11\rangle_{AB})/\sqrt2$.
- Error model for a single two–qubit gate
• After every \(CNOT_{ij}\) we apply an independent two-qubit depolarising
channel on the same pair \((i,j)\):
\[
(1-p)\,\rho+\frac{p}{15}\sum_{P\in\mathcal P_{ij}\setminus\{II\}}\!
\]
where $\mathcal P_{ij}$ is the $4\times4$ two-qubit Pauli group on
qubits $i$ and $j$.
• Only the \(CNOT_{21}\) and \(CNOT_{03}\) gates are noisy, so two such
channels act in the whole circuit.
- Fidelity after one noisy two-qubit gate
Let \(\rho=|\psi\rangle\langle\psi|\) be the input to the channel and
\(|\psi\rangle\) any pure state. The output state is
\(\rho'=\mathcal D_{ij}(\rho)\), and its fidelity with \(\rho\) is
\[
=(1-p)+\frac{p}{15}\sum_{P\neq II}
\]
For the particular state \(|\psi_{\text{ideal}}\rangle\) we must evaluate
the overlaps \(\langle\psi_{\text{ideal}}|P|\psi_{\text{ideal}}\rangle\)
for the 15 non-identity Paulis acting on the relevant two qubits.
• Qubits 2 & 1 (after \(CNOT_{21}\))
Among the 15 Paulis on qubits 2 and 1
only the three operators
\[
XX_{21},\;YY_{21},\;ZZ_{21}
\]
leave $|\psi_{\text{ideal}}\rangle$ invariant up to a global phase
(all others map the state to an orthogonal subspace).
Hence
\[
\]
• Qubits 0 & 3 (after \(CNOT_{03}\))
Exactly the same situation holds; the only invariants are
\[
XX_{03},\;YY_{03},\;ZZ_{03},
\]
so the sum of squared overlaps is again 3.
Therefore, for either noisy \(CNOT\) gate,
\[
=(1-p)+\frac{p}{15}\times 3
\]
- Two independent noisy \(CNOT\) gates
The two noise channels act on disjoint pairs of qubits and are
statistically independent.
Writing \(r\equiv 1-\dfrac{4p}{5}\) for the fidelity retention factor of
one gate, the overall process is
\[
\]
The final fidelity with the ideal state is then
\[
=\operatorname{Tr}\!\bigl(|\psi_{\text{ideal}}\rangle\langle\psi_{\text{ideal}}|
\,\rho_{\text{final}}\bigr)
=r^{2}
\]
Expanding,
\[
F_{\text{phys}}(p)=1-\frac{8p}{5}+\frac{16p^{2}}{25}.
\]
Conventions and units: All quantities are dimensionless probabilities.
No numerical approximation was required.
Final Answer:
The physical four-qubit state fidelity after the complete preparation
circuit is
\[
\]