5.9 Remarks and exercises

5.9.1 Entangled or not?

Let a joint state of \mathcal{A} and \mathcal{B} be written in a product basis as |\psi\rangle = \sum_{ij} c_{ij}|a_i\rangle\otimes|b_j\rangle. Assume that \mathcal{H}_{\mathcal{A}} and \mathcal{H}_{\mathcal{B}} are of equal dimension.

  1. Show that, if |\psi\rangle is a product state, then \det (c_{ij}) = 0.

  2. Show that the converse (\det(c_{ij})=0\implies|\psi\rangle=|a\rangle|b\rangle) holds only for qubits. Explain why.

  3. Deduce that the state \frac{1}{2}\big(|00\rangle + |01\rangle + |10\rangle + (-1)^k|11\rangle\big) is entangled for odd values of k and unentangled for even values of k. Express the latter case explicitly as a product state.

There is a lot of interesting physics behind the previous innocuous-looking mathematical statement. For example, think again about the state (|00\rangle+|11\rangle)/\sqrt{2}. What happens if you measure just the first qubit? It is equally likely that you get |0\rangle or |1\rangle, right? But after your measurement the two qubits are either in state |00\rangle or in |11\rangle, i.e. they show the same bit value. Now, why might that be disturbing? Well, imagine the second qubit to be light-years away from the first one. It seems that the measurement of the first qubit affects the second qubit right away, which seems to imply faster-than-light communication! This is what Einstein called “spooky action as a distance” in his 1947 letter to Max Born.

But can you actually use this effect to send a message faster than light? What would happen if you tried?

Hopefully you can see that it would not work, since the result of the measurement is random: you cannot choose the bit value you want to send. We shall return to this, and other related phenomena, later on — it is not at all a lost cause!

5.9.2 SWAP circuit

Show that, for any states |\psi_1\rangle and |\psi_2\rangle of two qubits, the circuit below implements the \texttt{SWAP} operation |\psi_1\rangle|\psi_2\rangle \mapsto |\psi_2\rangle|\psi_1\rangle.

(Swapping).

5.9.3 Controlled-NOT circuit

Show that the circuit below gives another implementation of the controlled-\texttt{NOT} gate.

(Controlled-\texttt{NOT}, again).

5.9.4 Measuring with controlled-NOT

The controlled-\texttt{NOT} gate can act as a measurement gate: if you prepare the target in state |0\rangle then the gate acts as |x\rangle|0\rangle\mapsto|x\rangle|x\rangle, and so the target learns the bit value of the control qubit. If you wish, you can think about a subsequent measurement of the target qubit in the computational basis as an observer learning about the bit value of the control qubit.

Take a look at the circuit below, where M stands for measurement in the standard basis.

(?).

Now assume that the top two qubits are in the state |\psi\rangle = \frac{1}{\sqrt3}\big( |01\rangle - |10\rangle + i|11\rangle \big). The measurement M gives two possible outcomes: 0 and 1. What are the probabilities of each outcome, and what is the post-measurement state in each case?

What is this circuit actually measuring?

5.9.5 Arbitrary controlled-U on two qubits

Recall Section 3.5: any unitary operation U on a single qubit can be expressed as U = B^\dagger XBA^\dagger XA for some unitaries A and B, where X\equiv\sigma_x is the Pauli bit-flip operator.

Suppose that you can implement any single-qubit gate, and that you have a bunch of controlled-\texttt{NOT} gates at your disposal. How would you implement any controlled-U operation on two qubits?

5.9.6 Entangled qubits

Two entangled qubits in the state \frac{1}{\sqrt{2}}(|00\rangle+|11\rangle) are generated by some source S. One qubit is sent to Alice, and one to Bob, who then both perform measurements in the computational basis.

  1. What is the probability that Alice and Bob will register identical results? Can any correlations they observe be used for instantaneous communication?

  2. Prior to the measurements in the computational basis, Alice and Bob apply unitary operations R_\alpha and R_\beta (respectively), for some fixed values \alpha,\beta\in\mathbb{R}, to their respective qubits:

    where the gate R_\theta is defined by its action on the basis states: \begin{aligned} |0\rangle &\longmapsto \cos\theta|0\rangle + \sin\theta|1\rangle \\|1\rangle &\longmapsto -\sin\theta|0\rangle + \cos\theta|1\rangle. \end{aligned} Show that the state of the two qubits prior to the measurements is \begin{aligned} &\frac{1}{\sqrt{2}}\cos(\alpha-\beta)\big( |00\rangle + |11\rangle \big) \\- &\frac{1}{\sqrt{2}}\sin(\alpha-\beta)\big( |01\rangle - |10\rangle \big). \end{aligned}

  3. What is the probability that Alice and Bob’s outcomes are identical?

  4. What is the geometric interpretation of the operator R_\theta?

5.9.7 Quantum dense coding

5.9.8 Playing with conditional unitaries

The swap gate S on two qubits is defined first on product vectors by S\colon|a\rangle|b\rangle\mapsto|b\rangle|a\rangle and then extended to sums of product vectors by linearity.

  1. Show that the four Bell states \frac{1}{\sqrt{2}}(|00\rangle\pm|11\rangle) and \frac{1}{\sqrt{2}}(|01\rangle\pm|10\rangle) are eigenvectors of S that form an orthonormal basis in the Hilbert space associated to two qubits. Which Bell states span the symmetric subspace (i.e. the space spanned by all eigenvectors with eigenvalue 1), and which the antisymmetric one (i.e. that spanned by eigenvectors with eigenvalue -1)? Can S have any eigenvalues apart from \pm1?

  2. Show that P_\pm = \frac{1}{2}(\mathbf{1}\pm S) are two orthogonal projectors which form the decomposition of the identity and project onto the symmetric and antisymmetric subspaces. Decompose the state vector |a\rangle|b\rangle of two qubits into symmetric and antisymmetric components.

  3. Consider the quantum circuit below, composed of two Hadamard gates, one controlled-S operation (also known as the controlled-swap, or Fredkin gate), and the measurement M in the computational basis. Suppose that the state vectors |a\rangle and |b\rangle are normalised but not orthogonal to one another. Step through the execution of this network, writing down the quantum states of the three qubits after each of the four computational steps. What are the probabilities of observing 0 or 1 when the measurement M is finally performed?

(Symmetric and antisymmetric projection).

  1. Explain why this quantum network implements projections on the symmetric and antisymmetric subspaces of the two qubits.

  2. Two qubits are transmitted through a quantum channel which applies the same randomly chosen unitary operation U to each of them, i.e. U\otimes U. Show that the symmetric and antisymmetric subspaces are invariant under this operation.

  3. Polarised photons are transmitted through an optical fibre. Due to the variation of the refractive index along the fibre, the polarisation of each photon is rotated by the same unknown angle. This makes communication based on polarisation encoding unreliable. However, if you are able to prepare any polarisation state of the two photons then you can still use the channel to communicate without any errors — how?

5.9.9 Tensor products in components

In our discussion of tensor products we have so far taken a rather abstract approach. There are, however, situations in which we have to put numbers in, and write tensor products of vectors and matrices explicitly. For example, here is the standard basis of two qubits written explicitly as column vectors:102 \begin{aligned} |00\rangle &\equiv |0\rangle\otimes|0\rangle = \begin{bmatrix}1\\0\end{bmatrix} \otimes \begin{bmatrix}1\\0\end{bmatrix} = \begin{bmatrix}1\\0\\0\\0\end{bmatrix} \\[1em] |01\rangle &\equiv |0\rangle\otimes|1\rangle = \begin{bmatrix}1\\0\end{bmatrix} \otimes \begin{bmatrix}0\\1\end{bmatrix} = \begin{bmatrix}0\\1\\0\\0\end{bmatrix} \\[1em] |10\rangle &\equiv |1\rangle\otimes|0\rangle = \begin{bmatrix}0\\1\end{bmatrix} \otimes \begin{bmatrix}1\\0\end{bmatrix} = \begin{bmatrix}0\\0\\1\\0\end{bmatrix} \\[1em] |11\rangle &\equiv |1\rangle\otimes|1\rangle = \begin{bmatrix}0\\1\end{bmatrix} \otimes \begin{bmatrix}0\\1\end{bmatrix} = \begin{bmatrix}0\\0\\0\\1\end{bmatrix} \end{aligned} Given |a\rangle = \alpha_0|0\rangle + \alpha_1|1\rangle and |b\rangle = \beta_0|0\rangle +\beta_1|1\rangle, we write |a\rangle\otimes|b\rangle as \begin{aligned} |a\rangle\otimes|b\rangle &= \begin{bmatrix}\alpha_0\\\alpha_1\end{bmatrix} \otimes \begin{bmatrix}\beta_0\\\beta_1\end{bmatrix} \\&= \begin{bmatrix}\alpha_0\begin{bmatrix}\beta_0\\\beta_1\end{bmatrix}\\\alpha_1\begin{bmatrix}\beta_0\\\beta_1\end{bmatrix}\end{bmatrix} \\&= \begin{bmatrix}\alpha_0\beta_0\\\alpha_0\beta_1\\\alpha_1\beta_0\\\alpha_1\beta_1\end{bmatrix}. \end{aligned} Note that each element of the first vector multiplies the entire second vector. This is often the easiest way to get the tensor products in practice.

The matrix elements of the tensor product operation A\otimes B

are given by (A\otimes B)_{ik,jl} = A_{ij}B_{kl} where ik\in\{00,01,10,11\} labels the rows, and kl\in\{00,01,10,11\} labels columns, when forming the block matrix: \begin{aligned} A\otimes B &= \begin{bmatrix}A_{00}&A_{01}\\A_{10}&A_{11}\end{bmatrix} \otimes \begin{bmatrix}B_{00}&B_{01}\\B_{10}&B_{11}\end{bmatrix} \\&= \begin{bmatrix}A_{00}B&A_{01}B\\A_{10}B&A_{11}B\end{bmatrix} \\&= \left[ \, \begin{array}{c|c} \begin{matrix}A_{00}B_{00}&A_{00}B_{01}\\A_{00}B_{10}&A_{00}B_{11}\end{matrix} & \begin{matrix}A_{01}B_{00}&A_{01}B_{01}\\A_{01}B_{10}&A_{01}B_{11}\end{matrix} \\\hline \begin{matrix}A_{10}B_{00}&A_{10}B_{01}\\A_{10}B_{10}&A_{10}B_{11}\end{matrix} & \begin{matrix}A_{11}B_{00}&A_{11}B_{01}\\A_{11}B_{10}&A_{11}B_{11}\end{matrix} \end{array} \, \right] \end{aligned}

This product A\otimes B also known as the Kronecker product of matrices, which generalises the outer product of two vectors that we saw in 0.8.

The tensor product induces a natural partition of matrices into blocks. Multiplication of block matrices works pretty much the same as regular matrix multiplication (assuming the dimensions of the sub-matrices are appropriate), except that the entries are now matrices rather than numbers, and so may not commute.

  1. Evaluate the following matrix product of (4\times 4) block matrices: \left[ \, \begin{array}{c|c} \mathbf{1}& X \\\hline Y & Z \end{array} \, \right] \left[ \, \begin{array}{c|c} \mathbf{1}& Y \\\hline X & Z \end{array} \, \right] (where X, Y, and Z are the Pauli matrices).

  2. Using the block matrix form of A\otimes B expressed in terms of A_{ij} and B_{ij} (as described above), explain how the following operations are performed on the block matrix:

    • transposition (A\otimes B)^T;
    • partial transpositions A^T\otimes B and A\otimes B^T;
    • trace \operatorname{tr}(A\otimes B);
    • partial traces (\operatorname{tr}A)\otimes B and A\otimes(\operatorname{tr}B).

5.9.10 Hadamard transforms in components

Consider the Hadamard transform H\otimes H\otimes H on three qubits, which is described by a (2^3\times2^3) matrix. We know that H = \frac{1}{\sqrt{2}}\begin{bmatrix}1&1\\1&-1\end{bmatrix} and so we can calculate that H\otimes H = \frac{1}{2} \left[ \, \begin{array}{c|c} \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}1&1\\1&-1\end{matrix} \\\hline \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}-1&-1\\-1&1\end{matrix} \end{array} \, \right] and thus that H\otimes H\otimes H = \frac{1}{2} \left[ \, \begin{array}{c|c|c|c} \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}1&1\\1&-1\end{matrix} \\\hline \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}-1&-1\\-1&1\end{matrix} & \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}-1&-1\\-1&1\end{matrix} \\\hline \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}-1&-1\\-1&1\end{matrix} & \begin{matrix}-1&-1\\-1&1\end{matrix} \\\hline \begin{matrix}1&1\\1&-1\end{matrix} & \begin{matrix}-1&-1\\-1&1\end{matrix} & \begin{matrix}-1&-1\\-1&1\end{matrix} & \begin{matrix}1&1\\1&-1\end{matrix} \end{array} \, \right]. The rows and columns of H\otimes H\otimes H are labelled by the triples 000,001,\ldots,111. Now, suppose we apply H\otimes H\otimes H to the state |110\rangle:

  1. The output state is a superposition of all binary strings: \sum_{x\in\{0,1\}^3} c_x|x\rangle. Where in the H^{\otimes 3} matrix will you find the coefficients c_x?

Do you want to write down H\otimes H\otimes H\otimes H? Probably not! This is an exponentially growing monster and you may soon run out of space if you actually do try to write it down. Instead, let us try to spot the pattern of the entries \pm1 in these matrices.

Consider again the single-qubit Hadamard gate matrix H=(H_{ab}), where a,b=0,1 are the labels for the rows and the columns. Observe that H_{ab}=(-1)^{ab}/\sqrt{2}. (This may look like a needlessly fancy way of writing the entries of the Hadamard matrix, but it will pay off in a moment).

  1. Using the fact that (A\otimes B)_{ik,jl} = A_{ij}B_{kl}, or any other method, analyse the pattern of the \pm1 in the tensor product of Hadamard matrices. What is the entry H^{\otimes 4}_{0101,1110}?

  2. For any two binary strings a=(a_1,\ldots, a_n) and b =(b_1,\ldots , b_n) of the same length we can define their “scalar” product as a\cdot b = (a_1b_1\oplus \ldots \oplus a_n b_n). Show that, up to the constant (1/\sqrt{2})^n, the entry H^{\otimes n}_{a,b} is (-1)^{a\cdot b} for any n and for any binary strings a and b of length n.

  3. Show that H^{\otimes n} acts as |a\rangle \longmapsto \left(\frac{1}{\sqrt{2}}\right)^n \sum_{b\in\{0,1\}^n} (-1)^{a\cdot b}|b\rangle

  4. A quantum register of 10 qubits holds the binary string 0110101001. The Hadamard Transform is then applied to this register yielding a superposition of all binary strings of length 10. What is the sign in front of the |0101010101\rangle term?

5.9.11 The Schmidt decomposition

An arbitrary vector in the Hilbert space \mathcal{H}_{\mathcal{A}}\otimes \mathcal{H}_{\mathcal{B}} can be expanded in a product basis as |\psi\rangle = \sum_{ij} c_{ij}|a_i\rangle|b_j\rangle. Moreover, for any given joint state |\psi\rangle, we can find orthonormal bases, \{|\tilde{a}_i\rangle\} in \mathcal{H}_{\mathcal{A}} and \{|\tilde{b}_j\rangle\} in \mathcal{H}_{\mathcal{B}}, such that |\psi\rangle can be expressed as |\psi\rangle = \sum_{i} d_{i}|\tilde a_i\rangle|\tilde b_i\rangle, where the coefficients d_i are non-negative numbers. This is known as the Schmidt decomposition of |\psi\rangle.

Any bipartite state can be expressed in this form, but remember that the bases used depend on the state being expanded. Indeed, given two bipartite states |\psi\rangle and |\phi\rangle, we usually cannot perform the Schmidt decomposition using the same orthonormal bases in \mathcal{H}_{\mathcal{A}} and \mathcal{H}_{\mathcal{B}}. The number of terms in the Schmidt decomposition is, at most, the minimum of \dim\mathcal{H}_{\mathcal{A}} and \dim\mathcal{H}_{\mathcal{B}}.

The Schmidt decomposition follows from the singular value decomposition (often abbreviated to SVD): any (n\times m) matrix C can be written as C = UDV where U and V are (respectively) (n\times n) and (m\times m) unitary matrices, and D is an (n\times m) diagonal matrix with real, non-negative elements in descending order d_1\geqslant d_2\geqslant\ldots\geqslant d_{\min\{n,m\}} (and with the rest of the matrix is filled with zeros). The elements d_k are called the singular values of C. We will return to the SVD in more detail later on, in Section 11.10.1.

You can visualize the SVD by thinking of C as representing a linear transformation from m-dimensional to n-dimensional Euclidean space: it maps the unit ball in the m-dimensional space to an ellipsoid in the n-dimensional space; the singular values are the lengths of the semi-axes of that ellipsoid; the matrices U and V carry information about the locations of those axes and the vectors in the first space which map into them. Thus SVD tells us that the transformation C consists of rotating the unit ball (the transformation V), stretching the k-th axis by a factor of d_k (the transformation D), and then rotating the resulting ellipsoid (the transformation U).

Using the index notation C_{ij} = \sum_k U_{ik}d_k V_{kj}, we can thus apply SVD to c_{ij}: \begin{aligned} |\psi\rangle &= \sum_{ij} c_{ij}|a_ib_j\rangle \\&= \sum_{ij} \sum_k U_{ik}d_k V_{kj}|a_ib_j\rangle \\&= \sum_k d_k \left(\sum_i U_{ik}|a_i\rangle\right)\otimes\left(\sum_j V_{kj}|b_j\rangle\right). \end{aligned} The Schmidt decomposition of a separable state of the form |a\rangle\otimes|b\rangle is trivially just this state. The Bell states \Psi^+ and \Phi^+ are already written in their Schmidt form, whereas \Psi^- and \Phi^- can be easily expressed in the Schmidt form. For example, for |\Psi^-\rangle we have d_1 = d_2 = \frac{1}{\sqrt{2}}, and the Schmidt basis is \begin{aligned} |\tilde{a}_1\rangle &= |0\rangle \\|\tilde{a}_2\rangle &= |1\rangle \\|\tilde{b}_1\rangle &= |1\rangle \\|\tilde{b}_2\rangle &= -|0\rangle. \end{aligned} The number of non-zero singular values of c_{ij} is called the rank of c_{ij}, or the rank of the corresponding quantum state, or sometimes, the Schmidt number. You should be able to see that all bipartite states of rank-one are separable.

The Schmidt decomposition is almost unique. The ambiguity arises when we have two or more identical singular values, as, for example, in the case of the Bell states. Then any unitary transformation of the basis vectors corresponding to a degenerate singular value, both in \mathcal{H}_a and in \mathcal{H}_b, generates another set of basis vectors.


  1. We always use the lexicographic order 00<01<10<11.↩︎