Checkpoint-2

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)

  1. Step-by-Step Derivation - Show every non-trivial step in the solution. Justify steps using relevant physical laws, theorems, or mathematical identities.
  2. Mathematical Typesetting - Use LaTeX for all mathematics: $...$ for inline expressions, $$...$$ for display equations.
  3. Conventions and Units - Follow the unit system and conventions specified in the problem.
  4. 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.

  5. 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\).

Checkpoints

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

1. Step–by–Step Derivation

1.1 Circuit and error locations
• The preparation circuit (written left–to–right in time order) is
 H₂ → CNOT₂1 → H₀ → CNOT₀3.
• The hardware model inserts a two–qubit depolarising channel after each CNOT.
 For the CNOT acting on qubits (i,j) the channel is

\[ \mathcal D_{ij}(\rho)= (1-p)\rho+\frac{p}{15}\sum_{P\in\mathcal{P}_{ij}\setminus\{II\}} P\rho P , \]

where \(\mathcal{P}_{ij}\) is the 16-element two–qubit Pauli group on qubits \(i,j\).

Thus two independent errors appear: • \(E_1\) (on qubits 2,1) after CNOT₂1,
\(E_2\) (on qubits 0,3) after CNOT₀3.
Single–qubit Hadamards are assumed perfect.


1.2 Effective Pauli acting on the final state
Because the gates following \(E_1\) act only on qubits 0 and 3, \(E_1\) commutes with all later gates.
Hence the net error on the final state is simply the tensor product

\[ P = E_2\,E_1 = (E_2^{(0,3)})\otimes(E_1^{(1,2)}), \qquad E_1\in\mathcal{P}_{21},\;E_2\in\mathcal{P}_{03}. \]

The two error channels are independent, so
\(\Pr\!\bigl(P=E_2E_1\bigr)=\Pr(E_1)\Pr(E_2).\)


1.3 Ideal target state \(|\psi\rangle\)
After the noiseless circuit the four physical qubits are in the logical Bell (GHZ) state.
It is the stabiliser state generated by

\[ \begin{aligned} S_1 &= XXXX, \\ S_2 &= ZZZZ, \\ S_3 &= X_A X_B = IXXI, \\ S_4 &= Z_A Z_B = IZZI . \end{aligned} \]

The full stabiliser group \(\mathcal S\) has \(2^4=16\) elements.
For a stabiliser state \(|\psi\rangle\),

\[ \langle\psi|P|\psi\rangle = \begin{cases} \pm1 & \text{if }P\in\mathcal S,\\[2pt] 0 & \text{otherwise}. \end{cases} \]

Therefore

\[ F(p)=\langle\psi|\rho|\psi\rangle =\sum_{E_1,E_2}\Pr(E_1)\Pr(E_2)\, \bigl|\langle\psi|E_2E_1|\psi\rangle\bigr|^2 =\Pr\!\bigl(P\in\mathcal S\bigr). \]

1.4 Which Pauli products lie in the stabiliser?
Every stabiliser element has the form

\[ P_{\sigma,\tau} = \sigma\otimes\tau\otimes\tau\otimes\sigma, \qquad \sigma,\tau\in\{I,X,Y,Z\}. \]

Hence \(P=E_2E_1\in\mathcal S\) iff

• the two-qubit Pauli on qubits (0,3) is \(\sigma\sigma\) for some \(\sigma\), and
• the two-qubit Pauli on qubits (1,2) is \(\tau\tau\) for some \(\tau\).


1.5 Probability that a two–qubit depolarising channel outputs \(\sigma\sigma\)

There are 16 two-qubit Paulis.
Exactly four of them have identical single-qubit factors:
\(II,\;XX,\;YY,\;ZZ.\)

\[ \Pr(\text{``same on both qubits''}) = (1-p) \;+\; 3\bigl(p/15\bigr) = 1-p+\tfrac{p}{5} = 1-\tfrac{4p}{5}. \]

1.6 Total fidelity

The two error channels are independent, so

\[ F(p)=\Bigl(1-\tfrac{4p}{5}\Bigr)^2 =1-\frac{8}{5}p+\frac{16}{25}p^{2}. \]

2. Final Answer

Final Answer: The physical four-qubit fidelity after the noisy preparation circuit is

\[ \boxed{F(p)=\left(1-\frac{4p}{5}\right)^{2} =1-\frac{8}{5}p+\frac{16}{25}p^{2}.} \]
User #3

The correct answer is: $ F_{\rm{physical}}=(1-\frac{12}{15}p)^2 $

Now let’s solve the next Checkpoint.

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). \]

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 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

1. Step–by–Step Derivation

We use the same error model as in Sub-problem 1:
after each CNOT a two–qubit depolarising channel acts on the same two qubits as the gate.
For the present circuit there are therefore two independent error channels

  • \(\mathcal D_{21}\) acting on qubits (2,1) with Kraus operator \(E_1\),
  • \(\mathcal D_{03}\) acting on qubits (0,3) with Kraus operator \(E_2\).

The total Pauli error on the four physical qubits is

\[ P \;=\;E_2^{(0,3)}\,E_1^{(1,2)} . \]

Throughout we ignore the overall \(\pm1,\pm i\) phases of Pauli operators.


1.1 Parity vector of a two–qubit Pauli

For a single–qubit Pauli \(P\in\{I,X,Y,Z\}\) define

\[ \alpha(P)= \begin{cases} 1 & \text{if }P\in\{Z,Y\}\\[2pt] 0 & \text{otherwise} \end{cases}, \qquad \beta(P)= \begin{cases} 1 & \text{if }P\in\{X,Y\}\\[2pt] 0 & \text{otherwise}. \end{cases} \]

For a two–qubit Pauli on qubits \((i,j)\) let

\[ s\,\equiv\,\alpha(P_i)+\alpha(P_j)\pmod 2,\qquad t\,\equiv\,\beta(P_i)+\beta(P_j)\pmod 2 . \]

The ordered pair

\[ v(P)=(s,t)\in\{(0,0),(0,1),(1,0),(1,1)\} \]

is called the parity vector of \(P\).


1.2 Detectability condition

The code stabilisers are

\[ S_X=XXXX,\qquad S_Z=ZZZZ . \]

A Pauli commutes with both stabilisers iff its parity vector is \((0,0)\) on every qubit it acts on.
Hence for the total error

\[ P=E_2^{(0,3)}E_1^{(1,2)} \]

to be undetectable (accepted) we must have

\[ v(E_2)=v(E_1). \]

1.3 Statistics of the parity classes

For a single two–qubit depolarising channel

  • probability of getting the identity \(II\) : \(1-p\)
  • probability of any one non-identity Pauli : \(p/15\).

Counting the 16 Pauli operators according to their parity vector gives

\((0,0)\)-class: \(II,\,XX,\,YY,\,ZZ\) (4 operators)
• each of the other three classes: 4 operators.

Therefore

\[ \begin{aligned} P_{(0,0)} &\equiv\Pr\bigl[v(P)=(0,0)\bigr] =(1-p)+3\!\left(\frac{p}{15}\right) = 1-\frac{4p}{5},\\[6pt] P_{\text{other}} &\equiv\Pr\bigl[v(P)=(0,1)\text{ or }(1,0)\text{ or }(1,1)\bigr] =\frac{4p}{15}\quad\text{(for each of the three classes).} \end{aligned} \]

1.4 Acceptance probability

Because \(E_1\) and \(E_2\) are independent,

\[ \Pr(\text{accept}) = P_{(0,0)}^{\,2}+3P_{\text{other}}^{\,2} =\Bigl(1-\frac{4p}{5}\Bigr)^{2}+3\Bigl(\frac{4p}{15}\Bigr)^{2} =1-\frac{8p}{5}+\frac{64p^{2}}{75}. \]

1.5 Logical action of the accepted errors

The physical operators that commute with both stabilisers form the code normaliser.
Within it, the logical Bell state

\[ |\Phi^+\rangle_{AB}=\frac{|00\rangle_{AB}+|11\rangle_{AB}}{\sqrt2} \]

is stabilised by the two logical operators

\[ X_A X_B,\qquad Z_A Z_B . \]

A convenient way to see the logical effect of a parity class is:

\((0,0)\) on qubits \((0,3)\) and \((1,2)\)
\(\;\Longrightarrow\) the four–qubit operator is of the form
\(S\,(X_A X_B)^{u}(Z_A Z_B)^{v}\) with \(u,v\in\{0,1\}\) and \(S\) a code stabiliser.
 These operators stabilise \(|\Phi^+\rangle_{AB}\), so they do not change the logical state.

\((0,1),\,(1,0),\,(1,1)\) on both pairs produce, up to stabilisers, the single-qubit logical Paulis
\(X_A,\;Z_A,\;Y_A\) (or the corresponding operators on \(B\)).
 Each of these maps \(|\Phi^+\rangle_{AB}\) to a state orthogonal to it.

Hence:

  • accepted error with parity \((0,0)\)no logical error (fidelity 1);
  • accepted error with any other parity ⇒ logical error (fidelity 0).

1.6 Logical fidelity conditioned on acceptance

Let

\[ P_{\text{good}} = P_{(0,0)}^{\,2}=\Bigl(1-\frac{4p}{5}\Bigr)^{2}. \]

The conditional logical fidelity is

\[ \begin{aligned} F_{\text{logical}}(p) &=\frac{P_{\text{good}}}{\Pr(\text{accept})}\\[4pt] &=\frac{\bigl(1-\dfrac{4p}{5}\bigr)^{2}} {\bigl(1-\dfrac{4p}{5}\bigr)^{2}+3\bigl(\dfrac{4p}{15}\bigr)^{2}}\\[10pt] &=\frac{\bigl(1-\dfrac{4p}{5}\bigr)^{2}} {1-\dfrac{8p}{5}+\dfrac{64p^{2}}{75}} . \end{aligned} \]

2. Final Answer

Final Answer: With post-selection on all detectable errors, the logical two-qubit state fidelity is

\[ \boxed{% F_{\text{logical}}(p)= \frac{\left(1-\dfrac{4p}{5}\right)^{2}} {\,1-\dfrac{8p}{5}+\dfrac{64p^{2}}{75}} }. \]
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 = ((1 - sp.Rational(4, 5) * p) ** 2) / (1 - sp.Rational(8, 5) * p + sp.Rational(64, 75) * p**2)
    # ---------------------------------------------------------------

    return F_logical