Linear algebra for quantum computing

Linear algebra is the language of quantum computing. Although you don’t need to know it to implement or write quantum programs, it is widely used to describe qubit states, quantum operations, and to predict what a quantum computer does in response to a sequence of instructions.

Just like being familiar with the basic concepts of quantum physics can help you understand quantum computing, knowing some basic linear algebra can help you understand how quantum algorithms work. At the least, you’ll want to be familiar with vectors and matrix multiplication. If you need to refresh your knowledge of these algebra concepts, here are some tutorials that cover the basics:

Vectors and matrices in quantum computing

A qubit can be in a state of 1 or 0 or a superposition of both. Using linear algebra, the state of a qubit is described as a vector and is represented by a single column matrix $\begin{bmatrix} a \\ b \end{bmatrix}$. It is also known as a quantum state vector and must meet the requirement that $|a|^2 + |b|^2 = 1$.

The elements of the matrix represent the probability of the qubit collapsing one way or the other, with $|a|^2$ being the probability of collapsing to zero, and $|b|^2$ being the probability of collapsing to one. The following matrices all represent valid quantum state vectors:

$$\begin{bmatrix} 1 \\ 0 \end{bmatrix}, \begin{bmatrix} 0 \\ 1 \end{bmatrix}, \begin{bmatrix} \frac{1}{\sqrt{2}} \\ \frac{1}{\sqrt{2}} \end{bmatrix}, \begin{bmatrix} \frac{1}{\sqrt{2}} \\ \frac{-1}{\sqrt{2}} \end{bmatrix}, \text{ and }\begin{bmatrix} \frac{1}{\sqrt{2}} \\ \frac{-i}{\sqrt{2}} \end{bmatrix}.$$ Quantum operations can also be represented by a matrix. When a quantum operation is applied to a qubit, the two matrices that represent them are multiplied and the resulting answer represents the new state of the qubit after the operation.

Here are two common quantum operations represented with matrix multiplication.

The X operation is represented by the Pauli matrix $X$,

$$X = \begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix},$$

and is used to flip the state of a qubit from 0 to 1 (or vice-versa), for example

$$\begin{bmatrix}0 &1\\ 1 &0\end{bmatrix}\begin{bmatrix} 1 \\ 0 \end{bmatrix} = \begin{bmatrix} 0 \\ 1 \end{bmatrix}.$$

The H operation is represented by the Hadamard transformation $H$,

$$H = \dfrac{1}{\sqrt{2}}\begin{bmatrix}1 &1\\ 1 &-1\end{bmatrix},$$

and puts a qubit into a superposition state where it has an even probability of collapsing either way, as shown here

$$\frac{1}{\sqrt{2}}\begin{bmatrix}1 &1\\ 1 &-1\end{bmatrix}\begin{bmatrix} 1 \\ 0 \end{bmatrix}=\frac{1}{\sqrt{2}}\begin{bmatrix} 1 \\ 1 \end{bmatrix}.$$

Notice that $|a|^2 =|b|^2 = \frac{1}{2}$, meaning that the probability of collapsing to zero and one state is the same.

A matrix that represents a quantum operation has one requirement – it must be a unitary matrix. A matrix is unitary if the inverse of the matrix is equal to the conjugate transpose of the matrix.

Representing two-qubit states

In the examples above, the state of one qubit was described using a single column matrix $\begin{bmatrix} a \\ b \end{bmatrix}$, and applying an operation to it was described by multiplying the two matrices. However, quantum computers use more than one qubit, so how do you describe the combined state of two qubits?

Note

The real power of quantum computing comes from leveraging multiple qubits to perform computations. For a deeper dive into this topic, see Operations on multiple qubits.

Remember that each qubit is a vector space, so they can't just be multiplied. Instead, you use a tensor product, which is a related operation that creates a new vector space from individual vector spaces, and is represented by the $\otimes$ symbol. For example, the tensor product of two qubit states $\begin{bmatrix} a \\ b \end{bmatrix}$ and $\begin{bmatrix} c \\ d \end{bmatrix}$ is calculated

$$ \begin{bmatrix} a \\ b \end{bmatrix} \otimes \begin{bmatrix} c \\ d \end{bmatrix} =\begin{bmatrix} a \begin{bmatrix} c \\ d \end{bmatrix} \\ b \begin{bmatrix}c \\ d \end{bmatrix} \end{bmatrix} = \begin{bmatrix} ac \\ ad \\ bc \\ bd \end{bmatrix}. $$

The result is a four-dimensional matrix, with each element representing a probability. For example, $ac$ is the probability of the two qubits collapsing to 0 and 0, $ad$ is the probability of 0 and 1, and so on.

Just as a single qubit state $\begin{bmatrix} a \\ b \end{bmatrix}$ must meet the requirement that $|a|^2 + |b|^2 = 1$ in order to represent a quantum state, a two-qubit state $\begin{bmatrix} ac \\ ad \\ bc \\ bd \end{bmatrix}$ must meet the requirement that $|ac|^2 + |ad|^2 + |bc|^2+ |bd|^2 = 1$.

Summary

Linear algebra is the standard language for describing quantum computing and quantum physics. Even though the standard library included with the Microsoft Quantum Development Kit helps you run advanced quantum algorithms without diving into the underlying math, understanding the basics helps you get started quickly and provide a solid foundation to build on.

Next steps