System #0

You are a physics research assistant specializing in solving complex, research-level problems using precise, step-by-step reasoning.

Input Problems will be provided in Markdown format.

Output (Markdown format)

Step-by-Step Derivation - Show every non-trivial step in the solution. Justify steps using relevant physical laws, theorems, or mathematical identities.
Mathematical Typesetting - Use LaTeX for all mathematics: $...$ for inline expressions, $$...$$ for display equations.
Conventions and Units - Follow the unit system and conventions specified in the problem.
Final Answer - At the end of the solution, start a new line with “Final Answer:”, and present the final result.

For final answers involving values, follow the precision requirements specified in the problem. If no precision is specified: - If an exact value is possible, provide it (e.g., \$\sqrt(2)\$, \$\pi/4\$). - If exact form is not feasible, retain at least 12 significant digits in the result.
Formatting Compliance - If the user requests a specific output format (e.g., code, table), provide the final answer accordingly.

User #1

Problem setup:

In quantum error correction, you encode quantum states into logical states made of many qubits in order to improve their resilience to errors. In quantum error detection, you do the same but can only detect the presence of errors and not correct them. In this problem, we will consider a single [[4,2,2]] quantum error detection code, which encodes two logical qubits into four physical qubits, and investigate how robust logical quantum operations in this code are to quantum errors.

Our convention is that the four physical qubits in the [[4,2,2]] code are labelled 0,1,2,3. The two logical qubits are labelled A and B. The stabilizers are $XXXX$ and $ZZZZ$, where $X$ and $Z$ are Pauli matrices. The logical $X$ and $Z$ operators on the two qubits are $X_A = XIXI$, $X_B=XXII$, $Z_A = ZZII$, $Z_B = ZIZI$, up to multiplication by stabilizers.

We will consider different state preparation circuits consisting of controlled not $CNOT_{ij}$ gates, where $CNOT_{ij}$ has control qubit $i$ and target qubit $j$. As a simple model of quantum errors in hardware, we will suppose that each $CNOT_{ij}$ gate in the circuit has a two qubit depolarizing error channel following it that produces one of the 15 non-identity two-qubit Paulis with equal probability $p/15$. The probability $p$ indicates the probability of an error in a single two-qubit gate. We will assess the logical infidelity of certain state preparation protocols as a function of the physical infidelity $p$.

Checkpoint 1

Suppose that we wish to prepare a logical two-qubit GHZ state $(|00\rangle_{AB}+|11\rangle_{AB})/\sqrt{2}$ in the [[4,2,2]] code. To do so, we use the following state preparation circuit:

\[ (CNOT_{03}) (H_0) (CNOT_{21}) (H_2). \]

Note that this equation is written in matrix multiplication order, while the quantum operations in the circuit occur in the reverse order (from right-to-left in the above equation). $H$ is a single-qubit Hadamard gate.

What is the physical state fidelity of the final physical 4-qubit state at the end of the circuit as a function of the two-qubit gate error rate $p$?

Assistant #2

Reasoning (model)

My Thinking on Calculating Fidelity in a Noisy Quantum Circuit

My goal is to figure out how well a prepared four-qubit state survives a circuit designed to create it, particularly given the presence of depolarizing noise on the two-qubit CNOT gates. I’m starting with a specific circuit – a Hadamard on qubit 2, a CNOT from qubit 2 to qubit 1, a Hadamard on qubit 0, and a CNOT from qubit 0 to qubit 3. The circuit aims to prepare a logical GHZ state (a superposition of |00> and |11> across two logical qubits), but I need to represent this in terms of the four physical qubits (0, 1, 2, and 3) used in the circuit.

The problem defines a specific encoding for logical qubits A and B using the four physical qubits, and provides stabilizers ($XXXX$ and $ZZZZ$) as well as logical $X$ and $Z$ operators. After some checking, I’ve confirmed that these logical operators are correctly constructed and satisfy the necessary commutation relations with the stabilizers, which means they are a valid representation of the logical qubits. Based on the code structure, I am using logical operators: $X_A = X_0 I_1 X_2 I_3$ $X_B = X_0 X_1 I_2 I_3$ $Z_A = Z_0 Z_1 I_2 I_3$ $Z_B = Z_0 I_1 Z_2 I_3$ I need to find the physical state that represents the logical $|00\rangle_{AB}$ state, which is stabilized by logical Z operators. I carefully construct the initial state in the physical four-qubit space, and calculate the target state for the given state preparation circuit. This calculation involves finding eigenstates of the stabilisers and logical operators. My analysis leads to the physical state:

$|00\rangle_{L} = \frac{1}{\sqrt{2}}(|0000\rangle + |1111\rangle)$, and $|11\rangle_{L} = \frac{1}{\sqrt{2}}(|0110\rangle + |1001\rangle)$.

and for the logical GHZ state:

$|\psi_{target}\rangle = \frac{1}{2}(|0000\rangle + |1111\rangle + |0110\rangle + |1001\rangle)$

The next challenge is accounting for the noise. The CNOT gates in the circuit are followed by two-qubit depolarizing channels, so I need to compute the final state by applying the depolarizing channels. The Hadamard gates are assumed to be perfect. Depolarizing channels scramble the quantum state by introducing a mixture of Pauli operators. For this, I am using the standard formula for fidelity: $\text{Tr}(\rho_{target} \rho_{final})$. I need to figure out what happens to the state, step-by-step, as it goes through the circuit and the two CNOT depolarizing gates.

The noise adds complexities:

The first CNOT and its error.
The Hadamard gate on the 0 qubit.
The second CNOT and its noise.
The target state I need to find the fidelity to. I’m carefully considering how each gate and error channel affects the state. Specifically, I know that the noisy state will become a sum of all possible Pauli errors.

I’ve worked to correctly derive the state after each channel, and then to properly compute the fidelity. The end results I get after all these calculations is a formula: $F = 1 - \frac{8p}{5} + \frac{16p^2}{25}$. This reduces to $F = (1 - \frac{4p}{5})^2$, where $p$ is the error rate. This provides a clear mathematical expression of fidelity as a function of the noise strength, and it’s something I can use to understand the performance of this circuit under depolarizing noise.

The problem asks for the physical state fidelity of the final 4-qubit state after a given state preparation circuit, subjected to depolarizing errors.

1. Determine the target physical state: The goal is to prepare the logical two-qubit GHZ state $|\psi_{GHZ}\rangle_L = \frac{1}{\sqrt{2}}(|00\rangle_{AB}+|11\rangle_{AB})$. First, we need to find the physical representations of the logical basis states $|00\rangle_{AB}$ and $|11\rangle_{AB}$. These states must be in the $+1$ eigenspace of the stabilizers $S_1 = X_0 X_1 X_2 X_3$ and $S_2 = Z_0 Z_1 Z_2 Z_3$. The logical state $|00\rangle_{AB}$ is a $+1$ eigenstate of $Z_A = Z_0 Z_1 I_2 I_3$ and $Z_B = Z_0 I_1 Z_2 I_3$. For a Z-basis state $|q_0 q_1 q_2 q_3\rangle$ to be a $+1$ eigenstate of $Z_0Z_1$ and $Z_0Z_2$: $Z_0Z_1|q_0 q_1 q_2 q_3\rangle = |q_0 q_1 q_2 q_3\rangle \implies (-1)^{q_0+q_1} = 1 \implies q_0+q_1$ is even (i.e., $q_0=q_1$). $Z_0Z_2|q_0 q_1 q_2 q_3\rangle = |q_0 q_1 q_2 q_3\rangle \implies (-1)^{q_0+q_2} = 1 \implies q_0=q_2$. From $S_2 = Z_0Z_1Z_2Z_3 = I$, and $q_0=q_1=q_2$, we have $Z_0Z_0Z_0Z_3 = I \implies Z_0Z_3 = I \implies q_0=q_3$. So, any Z-basis state satisfying these conditions must be of the form $|cccc\rangle$. Now, apply the $S_1 = X_0X_1X_2X_3$ condition. $S_1|0000\rangle = |1111\rangle$ and $S_1|1111\rangle = |0000\rangle$. To be a $+1$ eigenstate of $S_1$, the state must be a superposition: $|00\rangle_{AB} = \frac{1}{\sqrt{2}}(|0000\rangle + |1111\rangle)$.

The logical state $|11\rangle_{AB}$ is a $-1$ eigenstate of $Z_A$ and $Z_B$. $Z_0Z_1|q_0 q_1 q_2 q_3\rangle = -|q_0 q_1 q_2 q_3\rangle \implies q_0 \ne q_1$. $Z_0Z_2|q_0 q_1 q_2 q_3\rangle = -|q_0 q_1 q_2 q_3\rangle \implies q_0 \ne q_2$. From $S_2 = Z_0Z_1Z_2Z_3 = I$, we use $Z_0Z_1 = -I$ and $Z_0Z_2 = -I$: $(Z_0Z_1)(Z_0Z_2)Z_0Z_3 = I \implies (-I)(-I)Z_0Z_3 = I \implies Z_0Z_3 = I \implies q_0 = q_3$. This implies that $q_0=q_3$, $q_0 \ne q_1$, $q_0 \ne q_2$. Example basis states: $|0110\rangle$ or $|1001\rangle$. Check $S_1$ action: $S_1|0110\rangle = |1001\rangle$. So, $|11\rangle_{AB} = \frac{1}{\sqrt{2}}(|0110\rangle + |1001\rangle)$.

The target logical GHZ state is thus:

\[ |\psi_{target}\rangle = \frac{1}{\sqrt{2}}(|00\rangle_{AB} + |11\rangle_{AB}) = \frac{1}{2}(|0000\rangle + |1111\rangle + |0110\rangle + |1001\rangle) \]

2. Verify the circuit prepares the target state from $|0000\rangle$ (physical): The circuit operations (from right to left) are $H_2$, $CNOT_{21}$, $H_0$, $CNOT_{03}$. Starting with $|\psi_{in}\rangle = |0000\rangle$:

$H_2|0000\rangle = \frac{1}{\sqrt{2}}(|0000\rangle + |0010\rangle)$.
$CNOT_{21}\left(\frac{1}{\sqrt{2}}(|0000\rangle + |0010\rangle)\right) = \frac{1}{\sqrt{2}}(|0000\rangle + |0110\rangle)$. (Qubit 1 is target, Qubit 2 is control)
$H_0\left(\frac{1}{\sqrt{2}}(|0000\rangle + |0110\rangle)\right) = \frac{1}{2}((|0000\rangle + |1000\rangle) + (|0110\rangle + |1110\rangle))$.
$CNOT_{03}\left(\frac{1}{2}(|0000\rangle + |1000\rangle + |0110\rangle + |1110\rangle)\right) = \frac{1}{2}(|0000\rangle + |1001\rangle + |0110\rangle + |1111\rangle)$. (Qubit 3 is target, Qubit 0 is control) This final state is precisely $|\psi_{target}\rangle$.

3. Analyze the error model and calculate fidelity: Each $CNOT_{ij}$ gate is followed by a two-qubit depolarizing error channel, $\mathcal{D}(\rho) = (1-p)\rho + \frac{p}{15} \sum_{P \in \mathcal{P}_2^*} P \rho P$, where $\mathcal{P}_2^*$ is the set of 15 non-identity two-qubit Pauli operators. There are two $CNOT$ gates: $CNOT_{21}$ and $CNOT_{03}$. The Hadamard gates are assumed perfect. Let $U_1=H_2$, $U_2=CNOT_{21}$, $U_3=H_0$, $U_4=CNOT_{03}$. The initial state is $\rho_0 = |0000\rangle\langle0000|$. The overall evolution is $\rho_{final} = \mathcal{D}_{03}(U_4 U_3 \mathcal{D}_{21}(U_2 U_1 \rho_0 U_1^\dagger U_2^\dagger) U_3^\dagger U_4^\dagger)$. Let $\rho_{ideal} = U_4 U_3 U_2 U_1 \rho_0 U_1^\dagger U_2^\dagger U_3^\dagger U_4^\dagger = |\psi_{target}\rangle\langle\psi_{target}|$. The physical state fidelity is $F = \text{Tr}(\rho_{ideal} \rho_{final}) = \langle \psi_{target} | \rho_{final} | \psi_{target} \rangle$.

Let’s expand $\rho_{final}$: Define $\rho_A = U_2 U_1 \rho_0 U_1^\dagger U_2^\dagger$. This is the ideal state before the first error. Then $\mathcal{D}_{21}(\rho_A) = (1-p)\rho_A + \frac{p}{15} \sum_{P_{21}^*} P_{21} \rho_A P_{21}$. Next, $\rho_B = U_3 \mathcal{D}_{21}(\rho_A) U_3^\dagger = (1-p) U_3 \rho_A U_3^\dagger + \frac{p}{15} \sum_{P_{21}^*} U_3 P_{21} \rho_A P_{21} U_3^\dagger$. Since $U_3=H_0$ acts on qubit 0 and $P_{21}$ acts on qubits 1 and 2, $U_3$ commutes with $P_{21}$. So $U_3 P_{21} U_3^\dagger = P_{21}$. Let $\rho_C = U_3 \rho_A U_3^\dagger$. This is the ideal state after $H_0$ and before the $CNOT_{03}$. So $\rho_B = (1-p)\rho_C + \frac{p}{15} \sum_{P_{21}^*} P_{21} \rho_C P_{21}$. Finally, $\rho_{final} = \mathcal{D}_{03}(U_4 \rho_B U_4^\dagger) = (1-p) U_4 \rho_B U_4^\dagger + \frac{p}{15} \sum_{P_{03}^*} P_{03} U_4 \rho_B U_4^\dagger P_{03}$. Substitute $\rho_B$: $\rho_{final} = (1-p) U_4 \left[ (1-p)\rho_C + \frac{p}{15} \sum_{P_{21}^*} P_{21} \rho_C P_{21} \right] U_4^\dagger + \frac{p}{15} \sum_{P_{03}^*} P_{03} U_4 \left[ (1-p)\rho_C + \frac{p}{15} \sum_{P_{21}^*} P_{21} \rho_C P_{21} \right] U_4^\dagger P_{03}$. Let $\rho_T = U_4 \rho_C U_4^\dagger = |\psi_{target}\rangle\langle\psi_{target}|$. Note that $U_4=CNOT_{03}$ acts on qubits 0 and 3. $P_{21}$ acts on qubits 1 and 2. Thus $U_4$ commutes with $P_{21}$. So $U_4 P_{21} \rho_C P_{21} U_4^\dagger = P_{21} (U_4 \rho_C U_4^\dagger) P_{21} = P_{21} \rho_T P_{21}$. Thus, $\rho_{final} = (1-p) \left[ (1-p)\rho_T + \frac{p}{15} \sum_{P_{21}^*} P_{21} \rho_T P_{21} \right] + \frac{p}{15} \sum_{P_{03}^*} P_{03} \left[ (1-p)\rho_T + \frac{p}{15} \sum_{P_{21}^*} P_{21} \rho_T P_{21} \right] P_{03}$. The fidelity is $F = \text{Tr}(\rho_T \rho_{final})$. Using $\text{Tr}(\rho_T P_1 \rho_T P_2) = \text{Tr}(\rho_T P_1 \rho_T P_1 P_1 P_2) = \langle \psi_{target} | P_1 | \psi_{target} \rangle \langle \psi_{target} | P_2 | \psi_{target} \rangle$ if $P_1$ and $P_2$ are independent, or more simply, $\text{Tr}(\rho_T P_i \rho_T P_i) = |\langle \psi_{target} | P_i | \psi_{target} \rangle|^2$. $F = (1-p)^2 \text{Tr}(\rho_T \rho_T) + (1-p)\frac{p}{15} \sum_{P_{21}^*} \text{Tr}(\rho_T P_{21} \rho_T P_{21}) + \frac{p}{15}(1-p) \sum_{P_{03}^*} \text{Tr}(\rho_T P_{03} \rho_T P_{03}) + \left(\frac{p}{15}\right)^2 \sum_{P_{03}^*} \sum_{P_{21}^*} \text{Tr}(\rho_T P_{03} P_{21} \rho_T P_{21} P_{03})$. Since $\rho_T$ is a pure state, $\text{Tr}(\rho_T \rho_T) = 1$. The other terms are $|\langle \psi_{target} | P | \psi_{target} \rangle|^2$. $F = (1-p)^2 + (1-p)\frac{p}{15} \sum_{P_{21}^*} |\langle \psi_{target} | P_{21} | \psi_{target} \rangle|^2 + \frac{p}{15}(1-p) \sum_{P_{03}^*} |\langle \psi_{target} | P_{03} | \psi_{target} \rangle|^2 + \left(\frac{p}{15}\right)^2 \sum_{P_1 \in \mathcal{P}_{03}^*} \sum_{P_2 \in \mathcal{P}_{21}^*} |\langle \psi_{target} | P_1 P_2 | \psi_{target} \rangle|^2$.

4. Calculate the overlaps $\langle \psi_{target} | P | \psi_{target} \rangle$: $|\psi_{target}\rangle = \frac{1}{2}(|0000\rangle + |1111\rangle + |0110\rangle + |1001\rangle)$. An overlap $\langle \psi_{target} | P | \psi_{target} \rangle$ is non-zero only if $P|\psi_{target}\rangle = \pm |\psi_{target}\rangle$. Such operators are in the stabilizer group of $|\psi_{target}\rangle$. The stabilizer group $\mathcal{G}$ of $|\psi_{target}\rangle$ is generated by: $g_1 = X_0X_1X_2X_3$ (from $S_1$) $g_2 = Z_0Z_1Z_2Z_3$ (from $S_2$) $g_3 = X_A X_B = X_0 I_1 X_2 I_3 X_0 X_1 I_2 I_3 = I_0 X_1 X_2 I_3$ (physical $X_1X_2$) $g_4 = Z_A Z_B = Z_0 Z_1 I_2 I_3 Z_0 I_1 Z_2 I_3 = I_0 Z_1 Z_2 I_3$ (physical $Z_1Z_2$) Any product of these generators, including $I$, will result in $P|\psi_{target}\rangle = \pm |\psi_{target}\rangle$, meaning $|\langle \psi_{target} | P | \psi_{target} \rangle|^2 = 1$. Otherwise, the overlap is 0. The elements of $\mathcal{G}$ are $I, g_1, g_2, g_3, g_4$ and their products (16 elements).

Sum for $P_{21}^*$ (acting on $q_1,q_2$): The elements of $\mathcal{G}$ that act only on $q_1,q_2$ are $I, X_1X_2, Z_1Z_2, Y_1Y_2$. The non-identity ones are $X_1X_2, Z_1Z_2, Y_1Y_2$. There are 3 such Paulis in $\mathcal{P}_{21}^*$. For these 3 Paulis, $|\langle \psi_{target} | P | \psi_{target} \rangle|^2 = 1$. For the other $15-3=12$ Paulis, it’s 0. So, $\sum_{P_{21}^*} |\langle \psi_{target} | P_{21} | \psi_{target} \rangle|^2 = 3$.
Sum for $P_{03}^*$ (acting on $q_0,q_3$): The elements of $\mathcal{G}$ that act only on $q_0,q_3$ are $I, g_1 g_3 = (X_0X_1X_2X_3)(X_1X_2) = X_0X_3$, $g_2 g_4 = (Z_0Z_1Z_2Z_3)(Z_1Z_2) = Z_0Z_3$, and $(g_1g_3)(g_2g_4) = X_0X_3Z_0Z_3 = Y_0Y_3$. The non-identity ones are $X_0X_3, Z_0Z_3, Y_0Y_3$. There are 3 such Paulis in $\mathcal{P}_{03}^*$. For these 3 Paulis, $|\langle \psi_{target} | P | \psi_{target} \rangle|^2 = 1$. For the other $15-3=12$ Paulis, it’s 0. So, $\sum_{P_{03}^*} |\langle \psi_{target} | P_{03} | \psi_{target} \rangle|^2 = 3$.
Sum for $P_1 P_2$ where $P_1 \in \mathcal{P}_{03}^*$ and $P_2 \in \mathcal{P}_{21}^*$: The product $P_1 P_2$ is a four-qubit Pauli. We need to identify how many of these are in $\mathcal{G}$. The group of Paulis acting on $(q_0,q_3)$ that stabilize $|\psi_{target}\rangle$ is $G_{03}=\{I, X_0X_3, Z_0Z_3, Y_0Y_3\}$. The group of Paulis acting on $(q_1,q_2)$ that stabilize $|\psi_{target}\rangle$ is $G_{12}=\{I, X_1X_2, Z_1Z_2, Y_1Y_2\}$. The product set $G_{03} \times G_{12}$ consists of $4 \times 4 = 16$ elements. These 16 elements are precisely the stabilizer group $\mathcal{G}$. We are summing over $P_1 \in \mathcal{P}_{03}^*$ and $P_2 \in \mathcal{P}_{21}^*$. This means $P_1 \in G_{03}\setminus\{I\}$ and $P_2 \in G_{12}\setminus\{I\}$. There are $3 \times 3 = 9$ such products. For these 9 products, $|\langle \psi_{target} | P_1 P_2 | \psi_{target} \rangle|^2 = 1$. For the other $15 \times 15 - 9 = 216$ products, it’s 0. So, $\sum_{P_1 \in \mathcal{P}_{03}^*} \sum_{P_2 \in \mathcal{P}_{21}^*} |\langle \psi_{target} | P_1 P_2 | \psi_{target} \rangle|^2 = 9$.

5. Substitute values into the fidelity formula:

\[ F = (1-p)^2 + (1-p)\frac{p}{15}(3) + (1-p)\frac{p}{15}(3) + \left(\frac{p}{15}\right)^2(9) \]

\[ F = (1-p)^2 + \frac{6p(1-p)}{15} + \frac{9p^2}{225} \]

\[ F = (1-p)^2 + \frac{2p(1-p)}{5} + \frac{p^2}{25} \]

Expand and combine terms:

\[ F = (1 - 2p + p^2) + \left(\frac{2p}{5} - \frac{2p^2}{5}\right) + \frac{p^2}{25} \]

\[ F = 1 - 2p + \frac{2p}{5} + p^2 - \frac{2p^2}{5} + \frac{p^2}{25} \]

\[ F = 1 - \frac{10p}{5} + \frac{2p}{5} + p^2 \left(1 - \frac{2}{5} + \frac{1}{25}\right) \]

\[ F = 1 - \frac{8p}{5} + p^2 \left(\frac{25 - 10 + 1}{25}\right) \]

\[ F = 1 - \frac{8p}{5} + p^2 \left(\frac{16}{25}\right) \]

This expression can be factored as:

\[ F = \left(1 - \frac{4p}{5}\right)^2 \]

The final answer is $\boxed{\left(1 - \frac{4p}{5}\right)^2}$.

User #3

Checkpoint 2

Suppose that we wish to prepare a logical two-qubit GHZ state $(|00\rangle_{AB}+|11\rangle_{AB})/\sqrt{2}$ in the [[4,2,2]] code. To do so, we use the following state preparation circuit:

\[ (CNOT_{03}) (H_0) (CNOT_{21}) (H_2). \]

What is the logical state fidelity of the final 2-qubit logical state at the end of the circuit as a function of the two-qubit gate error rate $p$, assuming the state is post-selected on all detectable errors in the code?

Assistant #4

Reasoning (model)

Analyzing Logical State Fidelity with Error Detection

Okay, so I’m tackling the fidelity of a prepared GHZ state, but with a crucial condition: we’re only interested in the state after post-selection. That means we’re throwing out any states where errors were detected. My thinking here is to figure out how errors impact the final state, specifically, the logical state fidelity.

First, I need to remind myself of a couple of things. The target logical state I’m aiming for is the GHZ state represented as $|\psi_L\rangle = \frac{1}{\sqrt{2}}(|00\rangle_{AB}+|11\rangle_{AB})$. The physical state I need to prepare is $|\psi_{target}\rangle = \frac{1}{2}(|0000\rangle + |1111\rangle + |0110\rangle + |1001\rangle)$. The actual final state I get, after all the noise, is called $\rho_{final}$.

The code I’m using is a [[4,2,2]] code, which just means I have four qubits to protect two logical qubits, and it can detect two errors. I have two stabilizers: $S_1 = X_0 X_1 X_2 X_3$ and $S_2 = Z_0 Z_1 Z_2 Z_3$. States that are good codewords are in the +1 eigenspace of both stabilizers, meaning that when I act on the state with the stabilizer, the state is unchanged. Now, any error that happens can change the state: $E|\psi_{target}\rangle$. And my job now is to determine what the physical fidelity is.

What I do to detect errors is to measure the syndromes, which will tell me if the state is in the code space, and post-selection is just discarding any state I can detect an error in. The projector onto my code space, the states I’m interested in, looks like $P_{code} = \frac{1}{4}(I+S_1)(I+S_2)$. I can write the post-selected state as $ \rho_{ps} = \frac{P \rho_{final} P}{\text{Tr}(P \rho_{final})}$. The logical fidelity with post-selection is then $F_L = \text{Tr}(|\psi_{target}\rangle\langle\psi_{target}| \rho_{ps})$.

I’ve already figured out $\rho_{final}$ for this problem in subproblem 1. What I had was a lot of operators that I had to include. From here, I will plug in the definition of $P$ into the equation and try to extract information on the state fidelity.

To make the problem easier, I can break it down into four categories: * No error - The main GHZ state (ideal). * $P_{21}$ errors - The error occured during the second CNOT and acts on qubits 1 and 2. * $P_{03}$ errors - The error occurred during the first CNOT and acts on qubits 0 and 3. * $P_{21}$ and $P_{03}$ errors, both occurring on the same state.

The core of the problem comes down to figuring out which of these errors are “undetectable”, meaning they commute with my stabilizer, and what they do to the target state. These “undetectable” errors map to the state, and the states where they map to will go into my fidelity calculation. The other errors that don’t commute will be detected and post-selected out.

I need to classify the operators of the form $P \rho_T P$, where $\rho_T$ is $|\psi_{target}\rangle$, and then multiply by the constant $c_E$ to account for their statistical probability.

The main argument is the probability of the state. These should give me the logical state fidelity of 1, because, by definition, any error that isn’t detected should leave me with a perfect GHZ. After checking my math I had the same answer that I had previously calculated and confirmed that the logical fidelity for this problem is 1.

To determine the logical state fidelity with post-selection, we first need to identify the final state after the noisy circuit, then project it onto the code space (representing no detectable error), normalize it, and finally compute its fidelity with the ideal logical GHZ state.

1. Recall the final physical state from Subproblem 1: The ideal target physical state is $|\psi_{target}\rangle = \frac{1}{2}(|0000\rangle + |1111\rangle + |0110\rangle + |1001\rangle)$. The final density matrix after the noisy circuit, $\rho_{final}$, can be written as:

\[ \rho_{final} = (1-p)^2 \rho_T + (1-p)\frac{p}{15} \sum_{P_{21}^*} P_{21} \rho_T P_{21} + (1-p)\frac{p}{15} \sum_{P_{03}^*} P_{03} \rho_T P_{03} + \left(\frac{p}{15}\right)^2 \sum_{P_{03}^*} \sum_{P_{21}^*} P_{03} P_{21} \rho_T P_{21} P_{03} \]

where $\rho_T = |\psi_{target}\rangle\langle\psi_{target}|$. $P_{21}^*$ are the 15 non-identity two-qubit Paulis acting on qubits 1 and 2. $P_{03}^*$ are the 15 non-identity two-qubit Paulis acting on qubits 0 and 3.

2. Project onto the code space (post-selection): The [[4,2,2]] code has stabilizers $S_1 = X_0 X_1 X_2 X_3$ and $S_2 = Z_0 Z_1 Z_2 Z_3$. The projector onto the code space (the $+1$ eigenspace of both stabilizers) is $P_{code} = \frac{1}{4}(I+S_1)(I+S_2) = \frac{1}{4}(I+S_1+S_2+S_1S_2)$. When post-selecting on detectable errors, we are interested in the component of $\rho_{final}$ that remains within the code space. This is given by $P_{code} \rho_{final} P_{code}$.

Let’s analyze the effect of $P_{code}$ on each term $E \rho_T E^\dagger$ in $\rho_{final}$. If a Pauli operator $E$ commutes with all stabilizers (i.e., $[E, S_i]=0$ for all $S_i$), then $E$ maps a code state to another code state. In this case, $P_{code} E \rho_T E^\dagger P_{code} = E \rho_T E^\dagger$. If a Pauli operator $E$ anti-commutes with any stabilizer (i.e., $\{E, S_i\}=0$ for some $S_i$), then $E$ maps a code state out of the code space. In this case, $P_{code} E \rho_T E^\dagger P_{code} = 0$.

Let $\mathcal{N}$ be the set of Pauli operators that commute with both $S_1$ and $S_2$. The projected (unnormalized) state is $\rho_{no\_error} = \sum_{E \in \mathcal{N}} c_E E \rho_T E^\dagger$, where $c_E$ are the coefficients from the expansion of $\rho_{final}$.

From Subproblem 1, we identified the specific Paulis from $P_{21}^*$, $P_{03}^*$, and their products $P_{03}P_{21}$ that commute with $S_1$ and $S_2$. These are: * $P_{21}^*$: The 3 Paulis $X_1X_2, Z_1Z_2, Y_1Y_2$. These form the non-identity part of the stabilizer group $G_{12}$ of the ideal state acting on qubits 1,2. * $P_{03}^*$: The 3 Paulis $X_0X_3, Z_0Z_3, Y_0Y_3$. These form the non-identity part of the stabilizer group $G_{03}$ of the ideal state acting on qubits 0,3. * Products $P_{03}P_{21}$: The $3 \times 3 = 9$ products where $P_{03} \in G_{03}\setminus\{I\}$ and $P_{21} \in G_{12}\setminus\{I\}$.

Let’s check how these undetectable errors affect the logical state, i.e., whether $E \rho_T E^\dagger = \rho_T$ or $E \rho_T E^\dagger = \rho_{L'}$ where $\rho_{L'}$ is a different logical state. The ideal target state is $|\psi_{target}\rangle = \frac{1}{\sqrt{2}}(|00\rangle_{AB}+|11\rangle_{AB})$. * For $E=X_1X_2$: $X_1X_2|\psi_{target}\rangle = \frac{1}{2}(X_1X_2|0000\rangle + X_1X_2|1111\rangle + X_1X_2|0110\rangle + X_1X_2|1001\rangle)$ $= \frac{1}{2}(|0110\rangle + |1001\rangle + |0000\rangle + |1111\rangle) = |\psi_{target}\rangle$. * For $E=Z_1Z_2$: $Z_1Z_2|\psi_{target}\rangle = \frac{1}{2}(Z_1Z_2|0000\rangle + Z_1Z_2|1111\rangle + Z_1Z_2|0110\rangle + Z_1Z_2|1001\rangle)$ $= \frac{1}{2}(|0000\rangle + |1111\rangle + (-1)^{1+1}|0110\rangle + (-1)^{0+0}|1001\rangle) = |\psi_{target}\rangle$. * Since $Y_1Y_2 = -X_1X_2 Z_1Z_2$, $Y_1Y_2|\psi_{target}\rangle = -|\psi_{target}\rangle$, and thus $Y_1Y_2 \rho_T Y_1Y_2 = \rho_T$. * Similarly for $X_0X_3, Z_0Z_3, Y_0Y_3$, they all map $|\psi_{target}\rangle$ to $\pm |\psi_{target}\rangle$, hence $E \rho_T E^\dagger = \rho_T$. * For any product error $E = P_{03} P_{21}$ where $P_{03} \in G_{03}\setminus\{I\}$ and $P_{21} \in G_{12}\setminus\{I\}$, both $P_{03}$ and $P_{21}$ commute with $S_1, S_2$ and act as $\pm I$ on $|\psi_{target}\rangle$. Since $P_{03}$ and $P_{21}$ act on disjoint qubits, they commute. So $E|\psi_{target}\rangle = P_{03}P_{21}|\psi_{target}\rangle = P_{03}(\pm|\psi_{target}\rangle) = \pm (\pm|\psi_{target}\rangle) = \pm |\psi_{target}\rangle$. Thus $E \rho_T E^\dagger = \rho_T$.

This means that any error that commutes with the stabilizers $S_1, S_2$ (i.e., any undetectable error) happens to map the target state $\rho_T$ to itself. This implies that no logical errors are introduced by the undetectable physical errors in this specific scenario.

Therefore, the unnormalized post-selected state is:

\[ \rho_{no\_error} = \left[ (1-p)^2 \cdot 1 + 3(1-p)\frac{p}{15} \cdot 1 + 3(1-p)\frac{p}{15} \cdot 1 + 9\left(\frac{p}{15}\right)^2 \cdot 1 \right] \rho_T \]

The sum of coefficients in the bracket is precisely the physical fidelity calculated in Subproblem 1:

\[ F_{physical} = (1-p)^2 + \frac{6p(1-p)}{15} + \frac{9p^2}{225} = (1-p)^2 + \frac{2p(1-p)}{5} + \frac{p^2}{25} = 1 - \frac{8p}{5} + \frac{16p^2}{25} = \left(1 - \frac{4p}{5}\right)^2 \]

So, $\rho_{no\_error} = F_{physical} \rho_T$.

3. Normalize the post-selected state and compute logical fidelity: The probability of successfully preparing the state without detecting an error is $P_{success} = \text{Tr}(\rho_{no\_error}) = F_{physical} \text{Tr}(\rho_T) = F_{physical} \cdot 1 = \left(1 - \frac{4p}{5}\right)^2$. The normalized post-selected state is $\rho_{ps} = \frac{\rho_{no\_error}}{P_{success}} = \frac{F_{physical} \rho_T}{F_{physical}} = \rho_T$. The logical state fidelity is the fidelity between this post-selected state $\rho_{ps}$ and the ideal logical GHZ state (which is $\rho_T$):

\[ F_{logical} = \text{Tr}(\rho_T \rho_{ps}) = \text{Tr}(\rho_T \rho_T) = 1 \]

This result means that if no error is detected by the code, the logical state is perfectly prepared. The probability of this occurring is $F_{physical}$.

The final answer is $\boxed{1}$.

User #5

Populate your final answer into the code template provided below. This step is purely for formatting/display purposes. No additional reasoning or derivation should be performed. Do not import any modules or packages beyond what is provided in the template.

import sympy as sp

p = sp.symbols('p')

def answer(p):
    r"""
    Return the expression of logical state fidelity of the final 2-qubit logical state
    at the end of the circuit as a function of the two-qubit gate error rate $p$ in Sympy format.

    Inputs
    ----------
    p: sympy.Symbol, the two-qubit gate error rate, $p$

    Outputs
    ----------
    F_logical: sympy.Expr, the logical state fidelity as a function of $p$
    """

    # ------------------ FILL IN YOUR RESULTS BELOW ------------------
    F_logical = ...  # a SymPy expression of inputs
    # ---------------------------------------------------------------

    return F_logical

Assistant #6

import sympy as sp

p = sp.symbols('p')

def answer(p):
    r"""
    Return the expression of logical state fidelity of the final 2-qubit logical state
    at the end of the circuit as a function of the two-qubit gate error rate $p$ in Sympy format.

    Inputs
    ----------
    p: sympy.Symbol, the two-qubit gate error rate, $p$

    Outputs
    ----------
    F_logical: sympy.Expr, the logical state fidelity as a function of $p$
    """

    # ------------------ FILL IN YOUR RESULTS BELOW ------------------
    F_logical = sp.Integer(1)
    # ---------------------------------------------------------------

    return F_logical