0

I'm working my way through the ASPLOS paper "CutQC: Using Small Quantum Computers for Large Quantum Circuit Evaluations". At the heart of the paper is, of course, the technique of circuit cutting. However, I find this early derivation perplexing.

We know that we can Pauli-decompose any $2 \times 2$ matrix as: $$ A = \frac{1}{2}\bigl(Tr(AI) I + Tr(AX)X + Tr(AY)Y + Tr(AZ)Z\bigr). $$

The paper suggests to further decompose the Pauli matrices into their eigenbasis and write as: $$ A = \frac{A_1 + A_2 + A_3 + A_4}{2}, $$ with $$ A_1 = \bigl(Tr(AI) + Tr(AZ)\bigr)|0\rangle\langle 0| \\ A_2 = \bigl(Tr(AI) - Tr(AZ)\bigr)|1\rangle\langle 1| \\ A_3 = Tr(AX)\bigl( 2 |+\rangle\langle +| - |0\rangle\langle 0| - |1\rangle\langle 1|\bigr) \\ A_4 = Tr(AY)\bigl( 2 |+i\rangle\langle +i| - |0\rangle\langle 0| - |1\rangle\langle 1|\bigr) $$

How is this being derived?

glS
  • 24,708
  • 5
  • 34
  • 108
rhundt
  • 988
  • 4
  • 12
  • 1
    If you think about matrices as elements of a vector space, it is exactly the same as writing out your vector in some basis. – Rammus Dec 28 '22 at 20:56
  • I took the liberty of changing the question title because it is about basic linear algebra rather than about CutQC. – MonteNero Dec 29 '22 at 00:49
  • Greatly appreciated. Your answer is awesome! Thanks. – rhundt Dec 29 '22 at 01:21

1 Answers1

4

The Pauli matrices, together with the identity, form an orthonormal basis for a vector space of $2 \times 2$ matrices. As a result, we can express a matrix $A$ as a linear combination of these basis vectors (Pauli + Identity matrices). The coefficients in the linear combination can be expressed with the Hilbert-Schmidt inner product, which is defined as $Tr(U^{\dagger}W)$ for some $2 \times 2$ matrices $U$ and $W$. This is why we can rewrite $A$ as follows: $$ \tag{1} A = \frac{1}{2}\bigl(Tr(AI) I + Tr(AX)X + Tr(AY)Y + Tr(AZ)Z\bigr). $$ In the above equation, daggers (conjugate transposes) disappear because Pauli matrices are Hermitian.

The second part of your question relies on a similar basic fact. Any Hermitian matrix can be rewritten as a linear combination of rank-one projectors given by the eigenvectors of a matrix. This is known as spectral decomposition. Hence, we can express Pauli-Z as $$Z = |0\rangle \langle 0| - |1\rangle \langle 1|.$$ Note that $|0\rangle$ and $|1\rangle$ are eigenvectors of $Z$ with eigenvalues +1 and -1, respectively. Similarly, we can express Pauli-X and Y as \begin{align} X &= |+\rangle \langle +| - |-\rangle \langle -|,\\ Y &= |+i\rangle \langle +i| - |-i\rangle \langle -i|. \end{align} Finally, we use the resolution of the identity: $$\tag{2} I = |0\rangle \langle 0| + |1\rangle \langle 1|.$$

Now, substitute these decompositions into Eq (1) and rearrange all terms to get $A_i$. Getting $A_1$ and $A_2$ should be really trivial at this point, just simple algebra and grouping some terms.

Getting $A_3$ and $A_4$ is a bit trickier but still can be readily seen by inspection.

I will write out $A_3$ as an example. First, I hope it is obvious that $A_3 = Tr(AX)X$. This can be readily seen by looking at the prefactor $Tr(AX)$ in Eq (1) and then observing that $A_3$ has the same prefactor. Hence, all we need to do is to express $X$ in terms of projectors $|+\rangle\langle+|, |0\rangle\langle0|$ and $|1\rangle\langle1|$, this will give us $A_3$. Note that $$2|+\rangle\langle+| = \begin{pmatrix} 1 && 1\\ 1 && 1 \end{pmatrix}. $$ If we add $-I$, which can be decomposed as in Eq (2) we get \begin{align*} &2|+\rangle\langle+| - I \\ &= 2|+\rangle\langle+| - |0\rangle\langle0| - |1\rangle\langle1|\\ &= \begin{pmatrix} 1 && 1\\ 1 && 1 \end{pmatrix} - \begin{pmatrix} 1 && 0\\ 0 && 1 \end{pmatrix} = \begin{pmatrix} 0 && 1\\ 1 && 0 \end{pmatrix} = X. \end{align*} Hence, we have $$A_3 = Tr(AX)X = Tr(AX) \left(2|+\rangle\langle+| - |0\rangle\langle0| - |1\rangle\langle1|\right).$$

You can use a similar derivation to get $A_4$.

MonteNero
  • 2,396
  • 4
  • 22