Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Linear Algebra

Since I have already studied the basics of the necessary mathematical tools, for now I want to prioritize physical understanding. Thus, I will deepen my understanding of the mathematics to the level necessary to understand, apply, and critically question the physical framework under study. Purely mathematical proofs and developments can be left for when they become necessary or when I decide to delve deeper into the foundations.

For example, it is not necessary, at this point, to delve deeper into the proof of the spectral theorem for Hermitian matrices, since its proof is predominantly mathematical and does not, by itself, add a new physical interpretation to the problem. On the other hand, it is worth deepening mathematical results that establish relationships directly relevant to the physical theory, such as proving that two operators that have the same complete basis of eigenvectors commute, since this relationship allows us to understand mathematically the compatibility between observables.

However, I will be sharing here a purely mathematical review, an almost verbatim reproduction of this material (Delft University of Technology (2026)). A more complete review can be found aqui.

Basic Matrix Arithmetic

An m×nm\times n-matrix is an array of mm rows and nn columns containing (usually real or complex) numbers called the entries (or coefficients) of the matrix. The entries of a matrix AA are denoted by Ai,jA_{i,j}, aija_{ij} or similar, where the first index refers to the row and the second to the column. That is, the entries are laid out like

A=(a1,1a1,2⋯a1,na2,1a2,2⋯a2,n⋮⋮⋱⋮am,1am,2⋯am,n).A=\begin{pmatrix}a_{1,1} & a_{1,2} & \cdots & a_{1,n}\\ a_{2,1} & a_{2,2} & \cdots & a_{2,n}\\ \vdots & \vdots & \ddots & \vdots\\ a_{m,1} & a_{m,2} & \cdots & a_{m,n} \end{pmatrix}.

You will often see notation like A=(ai,j)A=(a_{i,j}) or A=(ai,j)1≤i≤m1≤j≤nA=(a_{i,j})_{{1\le i\le m\atop 1\le j\le n}} to indicate how the entries are denoted and/or how many rows and columns the matrix has.

In this section we focus on the basic operations that can be performed with matrices.

Addition and scalar multiplication

Like vectors, matrices can be added or subtracted entrywise. They can also be multiplied entrywise by a given scalar. For example,

4(103−1)−2(021−1)=(4012−4)−(042−2)=(4−410−2).4\begin{pmatrix}1 & 0\\ 3 & -1 \end{pmatrix}-2\begin{pmatrix}0 & 2\\ 1 & -1 \end{pmatrix}=\begin{pmatrix}4 & 0\\ 12 & -4 \end{pmatrix}-\begin{pmatrix}0 & 4\\ 2 & -2 \end{pmatrix}=\begin{pmatrix}4 & -4\\ 10 & -2 \end{pmatrix}.

Special matrices

A few (types of) matrices occur very often and have their own names:

- the m×nm\times n zero matrix (00⋯000⋯0⋮⋮⋱⋮00⋯0)\begin{pmatrix}0 & 0 & \cdots & 0\\ 0 & 0 & \cdots & 0\\ \vdots & \vdots & \ddots & \vdots\\ 0 & 0 & \cdots & 0 \end{pmatrix}

- the n×nn\times n identity matrix (10⋯001⋯0⋮⋮⋱⋮00⋯1)\begin{pmatrix}1 & 0 & \cdots & 0\\ 0 & 1 & \cdots & 0\\ \vdots & \vdots & \ddots & \vdots\\ 0 & 0 & \cdots & 1 \end{pmatrix}

- diagonal matrices (a1,10⋯00a2,2⋯0⋮⋮⋱⋮00⋯an,n)\begin{pmatrix}a_{1,1} & 0 & \cdots & 0\\ 0 & a_{2,2} & \cdots & 0\\ \vdots & \vdots & \ddots & \vdots\\ 0 & 0 & \cdots & a_{n,n} \end{pmatrix}

- upper triangular matrices (a1,1a1,2⋯a1,n0a2,2⋯a2,n⋮⋱⋱⋮0⋯0an,n)\begin{pmatrix}a_{1,1} & a_{1,2} & \cdots & a_{1,n}\\ 0 & a_{2,2} & \cdots & a_{2,n}\\ \vdots & \ddots & \ddots & \vdots\\ 0 & \cdots & 0 & a_{n,n} \end{pmatrix}

- lower triangular matrices (a1,10⋯0a2,1a2,2⋱⋮⋮⋮⋱0an,1an,2⋯an,n)\begin{pmatrix}a_{1,1} & 0 & \cdots & 0\\ a_{2,1} & a_{2,2} & \ddots & \vdots\\ \vdots & \vdots & \ddots & 0\\ a_{n,1} & a_{n,2} & \cdots & a_{n,n} \end{pmatrix}

Inner product of vectors

Definition Given two vectors v=(v1,…,vn)\boldsymbol{v}=(v_{1},\ldots,v_{n}) and w=(w1,…,wn)\boldsymbol{w}=(w_{1},\ldots,w_{n}), the inner product of v\boldsymbol{v} and w\boldsymbol{w} is the scalar v⋅w=v1w1+⋯+vnwn.\boldsymbol{v}\cdot\boldsymbol{w}=v_{1}w_{1}+\cdots+v_{n}w_{n}.

Properties of the inner product: For vectors v\boldsymbol{v}, w\boldsymbol{w}, x\boldsymbol{x} (all of the same length) and scalars cc the following hold:

Matrix-vector multiplication

Matrices can be multiplied by vectors.

Definition: Given a matrix AA as in 1 and a vector

v=(v1⋮vn)∈Rn\boldsymbol{v}=\begin{pmatrix}v_{1}\\ \vdots\\ v_{n} \end{pmatrix}\in\mathbb{R}^{n}

the (matrix-vector) product of AA and v\boldsymbol{v} is the vector

Av=(a1,1v1+⋯+a1,nvn⋮am,1v1+⋯+am,nvn)∈Rm.A\boldsymbol{v}=\begin{pmatrix}a_{1,1}v_{1}+\cdots+a_{1,n}v_{n}\\ \vdots\\ a_{m,1}v_{1}+\cdots+a_{m,n}v_{n} \end{pmatrix}\in\mathbb{R}^{m}.

Equivalently, we can write

(Av)i=∑j=1nAi,jvj(1≤i≤m).(A\boldsymbol{v})_{i}=\sum^{n}_{j=1}A_{i,j}v_{j}\quad(1\le i\le m).

Properties of matrix-vector multiplication: For matrices AA and BB, vectors v\boldsymbol{v} and w\boldsymbol{w} and scalars cc the following hold whenever the expressions are defined:

Note: Here are two other useful ways to think about matrix-vector products. First, in terms of inner products: the vector AvA\boldsymbol{v} consists of the inner products of each of the rows of AA with the vector v\boldsymbol{v}. Namely, if

A=(a1⋮am),A=\begin{pmatrix}\boldsymbol{a}_{1}\\ \vdots\\ \textbf{a}_{m} \end{pmatrix},

then

A=(a1⋅v⋮am⋅v).A=\begin{pmatrix}\boldsymbol{a}_{1}\cdot\boldsymbol{v}\\ \vdots\\ \boldsymbol{a}_{m}\cdot\boldsymbol{v} \end{pmatrix}.

Second, in terms of linear combinations: the vector AvA\boldsymbol{v} is a linear combination of the columns of AA where the coefficients are the entries of v\boldsymbol{v}. Namely, if

A=(a1 a2 ⋯ an)andv=(v1⋮vn),A=(\boldsymbol{a}_{1}\ \boldsymbol{a}_{2}\ \cdots\ \boldsymbol{a}_{n})\quad\text{and}\quad\boldsymbol{v}=\begin{pmatrix}v_{1}\\ \vdots\\ v_{n} \end{pmatrix},

then

Av=v1a1+v2a2+⋯+vnam.A\boldsymbol{v}=v_{1}\boldsymbol{a}_{1}+\boldsymbol{v}_{2}\boldsymbol{a}_{2}+\cdots+\boldsymbol{v}_{n}\boldsymbol{a}_{m}.

Matrix multiplication

Two matrices AA and BB can be multiplied if the number of columns of AA equals the number of rows of BB.

Definition: If AA is an m×nm\times n-matrix and BB is an n×pn\times p-matrix, then the (matrix) product of AA and BB is the m×pm\times p-matrix ABAB defined as follows: (AB)i,k=∑j=1nAi,jBj,k(1≤i≤m,1≤k≤p).(AB)_{i,k}=\sum^{n}_{j=1}A_{i,j}B_{j,k}\quad(1\le i\le m,1\le k\le p).

The product ABAB can be viewed as pp matrix-vector multiplications: if B=(b1 b2 ⋯bp)B=(\boldsymbol{b}_{1}\ \boldsymbol{b}_{2}\ \cdots\boldsymbol{b}_{p}), then

AB=(Ab1 Ab2 ⋯Abp).AB=(A\boldsymbol{b}_{1}\ A\boldsymbol{b}_{2}\ \cdots A\boldsymbol{b}_{p}).

Properties of matrix multiplication: For matrices AA, BB and CC and scalars cc the following hold whenever the expressions are defined:

Warning: Matrix multiplication is not commutative! That is, in general we have

AB≠BA.AB\ne BA.

Matrix Operations

Very often, an m×nm\times n-matrix AA with real entries represents a linear map from the space Rn\mathbb{R}^{n} of vectors of length nn to the space Rm\mathbb{R}^{m} of vectors of length mm. Similarly, a matrix with complex entries represents a linear map from Cn\mathbb{C}^{n} to Cm\mathbb{C}^{m}. This map is defined by the matrix-vector product.

Definition: Given a real (resp. complex) m×nm\times n-matrix AA, the linear map associated with AA is the map that sends each vector v∈Rn\boldsymbol{v}\in\mathbb{R}^{n} to the vector Av∈RmA\boldsymbol{v}\in\mathbb{R}^{m} (resp. each vector v∈Cn\boldsymbol{v}\in\mathbb{C}^{n} to the vector Av∈CmA\boldsymbol{v}\in\mathbb{C}^{m}). :::

Matrix multiplication is defined so that the product ABAB represents the linear map gotten by first applying BB and then AA. In other words, we have

(AB)v=A(Bv)(AB)\boldsymbol{v}=A(B\boldsymbol{v})

Note that on the left we have one matrix multiplication and one matrix-vector multiplication, while on the right we just have two matrix-vector multiplications. Commutators play an extremely important role in quantum physics.

Definition: Given two n×nn\times n-matrices AA and BB, the commutator of AA and BB is the n×nn\times n-matrix [A,B]=AB−BA.[A,B]=AB-BA.

Definition: The transpose of an m×nm\times n-matrix AA is the n×mn\times m-matrix A⊤A^{\top} defined by (A⊤)i,j=Aj,i.(A^{\top})_{i,j}=A_{j,i}.

Definition: The adjoint (conjugate transpose or Hermitian transpose) of an m×nm\times n-matrix AA is the n×mn\times m-matrix A†A^{\dagger} obtained by applying the complex conjugate and the transpose operation: (A†)i,j=Aj,i‾.(A^{\dagger})_{i,j}=\overline{A_{j,i}}. The adjoint of AA is also known as the ** or ** of AA.

Definition: The trace of an n×nn\times n-matrix AA is the sum of the diagonal entries of AA: Tr(A)=∑i=1nAi,i.\text{Tr}\left(A\right)=\sum^{n}_{i=1}A_{i,i}.

The determinant of an n×nn\times n-matrix AA can be defined in various ways. We will make use of permutations. Using the Einstein summation convention[1] and the Levi-Civita symbol[2] commonly used in physics, we can also write this as

det⁡A=ϵa1…anA1,a1…An,an\det A=\epsilon_{a_{1}\dots a_{n}}A_{1,a_{1}}\dots A_{n,a_{n}}

Properties of the determinant :

  1. A square matrix AA is invertible if and only if det⁡(A)\det(A) is non-zero.

  2. Given two n×nn\times n-matrices AA and BB, we have det⁡(AB)=(det⁡A)(det⁡B).\det(AB)=(\det A)(\det B).

  3. The determinant of a diagonal matrix (and more generally of an upper or lower triangular matrix) is the product of the diagonal entries.

  4. 4. If BB is obtained from AA by scaling a row by a scalar cc, then det⁡B=cdet⁡A\det B=c\det A.

  5. 5. If BB is obtained from AA by swapping two rows, then det⁡B=−det⁡A\det B=-\det A.

  6. 6. If BB is obtained from AA by adding a multiple of a row to another row, then det⁡B=det⁡A\det B=\det A.

Properties 4-6 also hold for columns instead of rows.

Finally, we define the characteristic polynomial of a square matrix; this will be used in diagonalisation.

Definition: The characteristic polynomial of an n×nn\times n-matrix AA is the polynomial of degree nn in a variable tt defined by χA(t)=det⁡(xI−A).\chi_{A}(t)=\det(x\mathbb{I}-A).[3]

Properties that are typically demonstrated in exercises:

Basis transformations

Until now, we have been working with one specific coordinate system in which vectors are identified with their coordinate vectors with respect to the standard basis vectors:

e1=(100⋮0),e2=(010⋮0),…,en=(000⋮1)\boldsymbol{e}_{1}=\begin{pmatrix}1\\ 0\\ 0\\ \vdots\\ 0 \end{pmatrix},\quad\boldsymbol{e}_{2}=\begin{pmatrix}0\\ 1\\ 0\\ \vdots\\ 0 \end{pmatrix},\quad\ldots,\quad\boldsymbol{e}_{n}=\begin{pmatrix}0\\ 0\\ 0\\ \vdots\\ 1 \end{pmatrix}

In other words, a vector v=(v1v2⋮vn)\boldsymbol{v}=\begin{pmatrix}v_{1}\\ v_{2}\\ \vdots\\ v_{n} \end{pmatrix} means the same as the corresponding linear combination of the standard basis vectors:

v=v1e1+⋯+vnen\boldsymbol{v}=v_{1}\boldsymbol{e}_{1}+\cdots+v_{n}\boldsymbol{e}_{n}

In many situations, it is useful to express vectors in a different coordinate system. For this we need the general concept of a basis.

Definition: A basis of Rn\mathbb{R}^{n} consists of nn vectors (b1,…,bn)(\boldsymbol{b}_{1},\ldots,\boldsymbol{b}_{n}) (or {bi}i=1n\left\{ \boldsymbol{b}_{i}\right\} ^{n}_{i=1}) such that every vector v∈Rn\boldsymbol{v}\in\mathbb{R}^{n} can be written in exactly one way as a linear combination:v=c1b1+⋯+cnbn\boldsymbol{v}=c_{1}\boldsymbol{b}_{1}+\cdots+c_{n}\boldsymbol{b}_{n}, with scalars c1,…,cnc_{1},\ldots,c_{n}. These scalars are called the coordinates of v\boldsymbol{v} relative to {bi}i=1n\left\{ \boldsymbol{b}_{i}\right\} ^{n}_{i=1}. The vector (c1…cn)T\left(\begin{array}{ccc} c_{1} & \dots & c_{n}\end{array}\right)^{T} is called the coordinate vector of v\boldsymbol{v} relative to {bi}i=1n\left\{ \boldsymbol{b}_{i}\right\} ^{n}_{i=1}.

Alternatively, a basis consists of nn vectors {bi}i=1n\left\{ \boldsymbol{b}_{i}\right\} ^{n}_{i=1} with the following two properties:

A priori, we could have defined a basis as a collection of some number (say mm) of vectors that span the space and are linearly independent. However, it can be shown that in fact any two bases contain the same number of vectors. Since the standard basis of Rn\mathbb{R}^{n} consists of nn vectors, the same therefore holds for any basis.

Suppose we have two bases of Rn\mathbb{R}^{n}, say

B=(b1,…,bn)andC=(c1,…,cn).B=(\boldsymbol{b}_{1},\ldots,\boldsymbol{b}_{n})\quad\text{and}\quad C=(\boldsymbol{c}_{1},\ldots,\boldsymbol{c}_{n}).

Definition: The basis transformation matrix (or change-of-basis matrix) from BB to CC is the matrix P=PB→CP=P_{B\to C} such that the jj-th column of PP contains the coordinates of bj\boldsymbol{b}_{j} relative to the basis (c1,…,cn)(\boldsymbol{c}_{1},\ldots,\boldsymbol{c}_{n}).

The definition of PP means:

bj=∑i=1nPi,jci.\boldsymbol{b}_{j}=\sum^{n}_{i=1}P_{i,j}\boldsymbol{c}_{i}.

Given a vector x\boldsymbol{x}, we write x(B)=(x1(B),…,xn(B))\boldsymbol{x}^{(B)}=(x^{(B)}_{1},\ldots,x^{(B)}_{n}) for the coordinate vector of x\boldsymbol{x} with respect to the basis {bi}i=1n\left\{ \boldsymbol{b}_{i}\right\} ^{n}_{i=1}, and we define x(C)\boldsymbol{x}^{(C)} similarly. Then we calculate

x=∑j=1nxj(B)bj=∑j=1nxj(B)∑i=1nPi,jci=∑i=1n(∑j=1nPi,jxj(B))ci=∑i=1n(Px(B))ici,\begin{aligned}\boldsymbol{x} & =\sum^{n}_{j=1}x^{(B)}_{j}\boldsymbol{b}_{j}\\ & =\sum^{n}_{j=1}x^{(B)}_{j}\sum^{n}_{i=1}P_{i,j}\boldsymbol{c}_{i}\\ & =\sum^{n}_{i=1}\left(\sum^{n}_{j=1}P_{i,j}x^{(B)}_{j}\right)\boldsymbol{c}_{i}\\ & =\sum^{n}_{i=1}\bigl(P\boldsymbol{x}^{(B)}\bigr)_{i}\boldsymbol{c}_{i}, \end{aligned}

which shows that

x(C)=PB→Cx(B).\boldsymbol{x}^{(C)}=P_{B\to C}\boldsymbol{x}^{(B)}.

Similarly, given a linear map AA, we write A(B)A^{(B)} for the matrix of AA with respect to the basis {bi}i=1n\left\{ \boldsymbol{b}_{i}\right\} ^{n}_{i=1}, and likewise for A(C)A^{(C)}. Then A(B)A^{(B)} and A(C)A^{(C)} are related by

A(C)=PA(B)P−1A^{(C)}=PA^{(B)}P^{-1}

Namely, consider what happens when we multiply this matrix to a coordinate vector of the form x(C)\boldsymbol{x}^{\left(C\right)}. First, applying P−1P^{-1}[4] to x(C)\boldsymbol{x}^{\left(C\right)} gives x(B)\boldsymbol{x}^{\left(B\right)}. Next, applying A(B)A^{(B)} gives (Ax)(B)(A\boldsymbol{x})^{(B)}. Finally, applying PP gives (Ax)(C)(A\boldsymbol{x})^{(C)}, which is what we wanted to show.

The trace and the determinant are invariant under basis transformation: if AA and PP are square matrices of the same size with PP invertible, then we have

Tr(PAP−1)=Tr(A)anddet⁡(PAP−1)=det⁡A\text{Tr}(PAP^{-1})=\text{Tr}\left(A\right)\quad\text{and}\quad\det(PAP^{-1})=\det A

Because the characteristic polynomial is also defined as a determinant, it too is invariant under basis transformation.

Properties that are typically demonstrated in exercises:

Diagonalisation

A very common example of a basis transformation is diagonalisation of a matrix. For this we first introduce the concept of eigenvalues and eigenvectors. Consider an n×nn\times n-matrix TT with complex entries. We view TT as a linear map on Cn\mathbb{C}^{n}.

Definition: An eigenvector for AA is a non-zero vector v\boldsymbol{v} such that Av=λvA\boldsymbol{v}=\lambda\boldsymbol{v} for some scalar λ\lambda. This λ\lambda is called the eigenvalue of TT acting on v\boldsymbol{v}.

Writing II for the n×nn\times n identity matrix, we see that if v\boldsymbol{v} is an eigenvector with eigenvalue λ\lambda, then

(A−λI)v=0(A-\lambda I)\boldsymbol{v}=0

This means that A−λIA-\lambda I is not invertible[5]; equivalently, its determinant vanishes. Thus, if λ\lambda is an eigenvalue of AA, we have

det⁡(A−λI)=0\det(A-\lambda I)=0

This implies that the eigenvalues of AA are roots of the characteristic polynomial χA\chi_{A} . We can factor the characteristic polynomial as

χA(x)=(x−λ1)k1⋯(x−λr)kr\chi_{A}(x)=(x-\lambda_{1})^{k_{1}}\cdots(x-\lambda_{r})^{k_{r}}

where λ1,…,λr\lambda_{1},\ldots,\lambda_{r} are the distinct roots of χA(t)\chi_{A}(t) and each kik_{i} is an integer called the multiplicity of λi\lambda_{i} as a root of χA(t)\chi_{A}(t).

The idea of diagonalisation is that it is often easiest to understand a linear map AA if we can find a basis of our space consisting of vectors on which AA acts as multiplication by some scalar.

Definition: An n×nn\times n-matrix AA is diagonalisable if there exist a diagonal[6] matrix DD and an invertible n×nn\times n-matrix PP such that A=PDP−1,A=PDP^{-1},

We can rewrite the above equation as AP=PDAP=PD. Let us write v1,…,vn\boldsymbol{v}_{1},\ldots,\boldsymbol{v}_{n} for the columns of PP and λ1,…,λn\lambda_{1},\ldots,\lambda_{n} for the diagonal entries of DD. Then we compute

AP=(Av1Av2⋯Avn)AP=\begin{pmatrix}A\boldsymbol{v}_{1} & A\boldsymbol{v}_{2} & \cdots & A\boldsymbol{v}_{n}\end{pmatrix}

and

PD=(λ1v1λv2⋯λnvn)PD=\begin{pmatrix}\lambda_{1}\boldsymbol{v}_{1} & \lambda\boldsymbol{v}_{2} & \cdots & \lambda_{n}\boldsymbol{v}_{n}\end{pmatrix}

We therefore see that the equation AP=PDAP=PD comes down to the nn equations Avi=λiviA\boldsymbol{v}_{i}=\lambda_{i}\boldsymbol{v}_{i}. In other words, each vi\boldsymbol{v}_{i} is an eigenvector of AA and λi\lambda_{i} is the corresponding eigenvalue.

The fact that PP is invertible means that the eigenvectors vi\boldsymbol{v}_{i} form a basis of Cn\mathbb{C}^{n}. We obtain the following conclusion.

Theorem: An n×nn\times n-matrix AA is diagonalisable precisely when there exists a basis of Cn\mathbb{C}^{n} consisting of eigenvectors for AA.

This also gives us a way to tell whether AA is diagonalisable, and if so to determine matrices DD and PP as above:

1. Determine the eigenvalues λ1,…,λr\lambda_{1},\ldots,\lambda_{r} of AA by finding the roots of the characteristic polynomial χA\chi_{A}.

2. For each eigenvalue λi\lambda_{i}, say with multiplicity kik_{i}, determine if there exist kik_{i} linearly independent eigenvectors vλi,1,…,vλi,ki\boldsymbol{v}_{\lambda_{i},1},\ldots,\boldsymbol{v}_{\lambda_{i},k_{i}} with eigenvalue λi\lambda_{i}.

3. The matrix AA is diagonalisable if and only if there exist kik_{i} linearly independent eigenvectors with eigenvalue λi\lambda_{i} for all ii. In this case, the matrices

D=(λ1Ik10⋯00λ2Ik2⋯0⋮⋮⋱⋮00⋯λrIkr)D=\begin{pmatrix}\lambda_{1}I_{k_{1}} & 0 & \cdots & 0\\ 0 & \lambda_{2}I_{k_{2}} & \cdots & 0\\ \vdots & \vdots & \ddots & \vdots\\ 0 & 0 & \cdots & \lambda_{r}I_{k_{r}} \end{pmatrix}

and

P=(vλ1,1⋯vλ1,k1vλ2,1⋯vλ2,k2⋯vλr,1⋯vλr,kr)P=\begin{pmatrix}\boldsymbol{v}_{\lambda_{1},1} & \cdots & \boldsymbol{v}_{\lambda_{1},k_{1}} & \boldsymbol{v}_{\lambda_{2},1} & \cdots & \boldsymbol{v}_{\lambda_{2},k_{2}} & \cdots & \boldsymbol{v}_{\lambda_{r},1} & \cdots & \boldsymbol{v}_{\lambda_{r},k_{r}}\end{pmatrix}

give a diagonalisation of AA.

Properties that are typically demonstrated in exercises:

Hilbert Spaces

So far, we have been considering vectors in an nn-dimensional real or complex space. In quantum mechanics, one encounters infinite-dimensional vector spaces as well. In particular, the state of a quantum system is represented mathematically by a (unit) vector in a Hilbert space.

Vector spaces

As is customary in quantum physics and quantum algorithms, we use bra-ket notation for vectors in Hilbert spaces.

Definition: A (real or complex) *vector space* consists of

- a set of vectors, denoted by symbols like ∣α⟩\left|\alpha\right\rangle

- a distinguished vector 0\mathbf{0} called the zero vector

- a way of adding vectors, i.e. given two vectors ∣α⟩\left|\alpha\right\rangle and ∣β⟩\left|\beta\right\rangle we can obtain a third vector ∣α⟩+∣β⟩\left|\alpha\right\rangle +\left|\beta\right\rangle

- a way of multiplying a vector by a scalar: given a scalar cc (a real or complex number, depending on whether we consider a real or complex vector space) and a vector ∣α⟩\left|\alpha\right\rangle we can obtain a vector c∣α⟩c\left|\alpha\right\rangle

Such that a number of properties hold:

- 0∣α⟩=00\left|\alpha\right\rangle =\mathbf{0}

- 1∣α⟩=∣α⟩1\left|\alpha\right\rangle =\left|\alpha\right\rangle

- 0+∣α⟩=∣α⟩\mathbf{0}+\left|\alpha\right\rangle =\left|\alpha\right\rangle

- ∣α⟩+∣β⟩=∣β⟩+∣α⟩\left|\alpha\right\rangle +\left|\beta\right\rangle =\left|\beta\right\rangle +\left|\alpha\right\rangle

- (∣α⟩+∣β⟩)+∣γ⟩=∣α⟩+(∣β⟩+∣γ⟩)(\left|\alpha\right\rangle +\left|\beta\right\rangle )+\left|\gamma\right\rangle =\left|\alpha\right\rangle +(\left|\beta\right\rangle +\left|\gamma\right\rangle )

- c(∣α⟩+∣β⟩)=c∣α⟩+c∣β⟩c(\left|\alpha\right\rangle +\left|\beta\right\rangle )=c\left|\alpha\right\rangle +c\left|\beta\right\rangle

- (c+d)∣α⟩=c∣α⟩+d∣α⟩(c+d)\left|\alpha\right\rangle =c\left|\alpha\right\rangle +d\left|\alpha\right\rangle

- c(d∣α⟩)=(cd)∣α⟩c(d\left|\alpha\right\rangle )=(cd)\left|\alpha\right\rangle

Definition: A basis for a vector space is a collection of vectors {∣ei⟩}i\left\{ \left|e_{i}\right\rangle \right\} _{i} such that every vector can be written in exactly one way as a linear combination of the {∣ei⟩}i\left\{ \left|e_{i}\right\rangle \right\} _{i}. The number of vectors in a basis is called the dimension of the space.

As in basis-transformations, any two bases contain the same number of vectors, so the dimension does not depend on the choice of basis.

Exemplo: The simplest examples are the nn-dimensional vector spaces R2\mathbb{R}^{2} and Cn\mathbb{C}^{n} with coordinatewise addition and scalar multiplication. An important special case is the two-dimensional space C2\mathbb{C}^{2}. In bra-ket notation, various symbols are used for the standard basis vectors (10)\binom{1}{0} and (01)\binom{0}{1} in this space, such as:

Hermitian inner products and Hilbert spaces

Definition: A Hermitian inner product assigns to vectors ∣α⟩\left|\alpha\right\rangle and ∣β⟩\left|\beta\right\rangle a complex number ⟨α∣β⟩\left\langle \alpha|\beta\right\rangle such that the following properties hold:

- ⟨α∣β⟩=⟨β∣α⟩‾\left\langle \alpha|\beta\right\rangle =\overline{\left\langle \beta|\alpha\right\rangle };

- ⟨α∣α⟩≥0\left\langle \alpha|\alpha\right\rangle \ge0, and ⟨α∣α⟩=0\left\langle \alpha|\alpha\right\rangle =0 holds precisely when ∣α⟩=0\left|\alpha\right\rangle =\mathbf{0} (note that ⟨α∣α⟩\left\langle \alpha|\alpha\right\rangle is a real number because of the first property);

- ⟨α∣(∣β⟩+∣γ⟩)=⟨α∣β⟩+⟨α∣γ⟩\left\langle \alpha\right|\left(\left|\beta\right\rangle +\left|\gamma\right\rangle \right)=\left\langle \alpha|\beta\right\rangle +\left\langle \alpha|\gamma\right\rangle;

- ⟨α∣(c∣β⟩)=c(⟨α∣β⟩)\left\langle \alpha\right|\left(c\left|\beta\right\rangle \right)=c\left(\left\langle \alpha|\beta\right\rangle \right);

Definition: The norm or length of a vector ∣α⟩\left|\alpha\right\rangle in a Hilbert space is defined as the real number ∥α∥=⟨α∣α⟩\|\alpha\|=\sqrt{\left\langle \alpha|\alpha\right\rangle }.

The properties above imply that for any scalar cc, the inner product of c∣α⟩c\left|\alpha\right\rangle with ∣β⟩\left|\beta\right\rangle equals cˉ⟨α∣β⟩\bar{c}\left\langle \alpha|\beta\right\rangle. In other words, while the inner product is (complex) linear in its second argument, it is only conjugate linear in its first argument. This is the usual convention in physics. It differs from the convention in mathematics, where inner products are usually linear in the first argument and conjugate linear in the second argument.

Definition: A Hilbert space is a complex vector space H\mathcal{H} equipped with a Hermitian inner product ⟨∣⟩\left\langle \enspace|\enspace\right\rangle and satisfying a condition called completeness.

We will not make the notion of completeness precise in these notes. Roughly speaking, it means that we do not have to worry too much about convergence of sequences or series in our Hilbert space, and there is a notion of infinite linear combinations[8].

Example: We take H=Cn\mathcal{H}=\mathbb{C}^{n} and define an inner product by ⟨α∣β⟩=∑i=1nαˉiβi.\left\langle \alpha|\beta\right\rangle =\sum^{n}_{i=1}\bar{\alpha}_{i}\beta_{i}. This satisfies all the properties of a Hermitian inner product.

The dual Hilbert space

With any Hilbert space H\mathcal{H}, we can associate a dual Hilbert space H∗\mathcal{H}^{*}. In the bra-ket formalism used in quantum mechanics, for every “ket vector” ∣α⟩\left|\alpha\right\rangle in H\mathcal{H} we have a dual “bra vector” ⟨α∣\left\langle \alpha\right| in H∗\mathcal{H}^{*}. This correspondence has the following properties:

vector inHcorresponding vector inH∗∣α⟩⟨α∣∣α⟩+∣β⟩⟨α∣+⟨β∣c∣α⟩c‾⟨α∣\begin{array}{cc} \text{vector in}\mathcal{H} & \text{corresponding vector in}\mathcal{H^{*}}\\ \left|\alpha\right\rangle & \left\langle \alpha\right|\\ \left|\alpha\right\rangle +\left|\beta\right\rangle & \left\langle \alpha\right|+\left\langle \beta\right|\\ c\left|\alpha\right\rangle & \overline{c}\left\langle \alpha\right| \end{array}

In mathematical language, H∗\mathcal{H^{*}} is the space of continuous linear maps H→C\mathcal{H}\to\mathbb{C}. The bra vector ⟨α∣\left\langle \alpha\right| then represents the linear map that sends ∣β⟩\left|\beta\right\rangle to ⟨α∣β⟩\left\langle \alpha|\beta\right\rangle.

Example:In the case H=Cn\mathcal{H}=\mathbb{C}^{n} we can identify H\mathcal{H} with Cn\mathbb{C}^{n} as well; for a vector ∣α⟩\left|\alpha\right\rangle in Cn\mathbb{C}^{n}, the dual vector ⟨α∣\left\langle \alpha\right| can then be naturally identified with the corresponding vector in Cn\mathbb{C}^{n}. It can also be convenient to distinguish H\mathcal{H} from H∗\mathcal{H}^{*} by denoting a ket vector ∣α⟩\left|\alpha\right\rangle as a column vector and the corresponding bra vector ⟨α∣\left\langle \alpha\right| by the conjugate transpose of ∣α⟩\left|\alpha\right\rangle, so ⟨α∣\left\langle \alpha\right| is the row vector whose entries are the complex conjugates of those of α\alpha.

In C2\mathbb{C}^{2}, for example, we have

∣α⟩=(1i)⟹⟨α∣=(1−i).\left|\alpha\right\rangle =\begin{pmatrix}1\\ i \end{pmatrix}\Longrightarrow\left\langle \alpha\right|=\begin{pmatrix}1 & -i\end{pmatrix}.

This allows us to view ∣α⟩\left|\alpha\right\rangle as a linear map sending ⟨β∣\left\langle \beta\right| to the scalar obtained by multiplying the row vector ∣α⟩\left|\alpha\right\rangle by the column vector ⟨β∣\left\langle \beta\right|. If you know about continuous linear maps of Hilbert spaces, here is a useful connection between the mathematical definition of H∗\mathcal{H}^{*} and the bra-ket formalism: the Riesz representation theorem states that every continuous linear map H→C\mathcal{H}\to\mathbb{C} is of the form ⟨α∣⟩\left\langle \alpha|\enskip\right\rangle for some α\alpha in H\mathcal{H}. Ou seja, todo mapa linear contínuo f:H→Cf:\mathcal{H}\to\mathbb{C} tem a forma ⟨α∣⟩\left\langle \alpha|\enskip\right\rangle onde ⟨α∣∈H∗\left\langle \alpha\right|\in\mathcal{H}^{*} tem um único vetor correspondente ∣α⟩∈H\left|\alpha\right\rangle \in\mathcal{H}.

Orthonormal bases

When working with vector spaces, it is frequently useful to make a choice of basis adapted to the problem at hand. In the context of a Hilbert space, it turns out that a “good” basis is often one that satisfies the following property with respect to the inner product.

Definition: An orthonormal system of vectors in a Hilbert space H\mathcal{H} is a collection of vectors {∣ei⟩}i∈I\left\{ \left|e_{i}\right\rangle \right\} _{i\in I}, where II is some index set, satisfying

⟨ei∣ej⟩={1if i=j,0if i≠j.\left\langle e_{i}|e_{j}\right\rangle =\begin{cases} 1 & \quad\text{if }i=j,\\ 0 & \quad\text{if }i\ne j. \end{cases}

Definition: An orthonormal basis of a Hilbert space H\mathcal{H} is an orthonormal system {∣ei⟩}i∈I\left\{ \left|e_{i}\right\rangle \right\} _{i\in I} in H\mathcal{H} that spans H\mathcal{H}, i.e. every vector in H\mathcal{H} can be written as a (possibly infinite) linear combination of the {∣ei⟩}i∈I\left\{ \left|e_{i}\right\rangle \right\} _{i\in I}.

Example: In the Hilbert space Cn\mathbb{C}^{n} with the standard inner product, the standard basis (e1,…,en)(e_{1},\ldots,e_{n}) is an orthonormal basis.

Linear subspaces

Definition: A linear subspace of a vector space VV is a set of vectors WW contained in VV such that the following properties hold:

- The zero vector of VV lies in WW.

- For all vectors ww and w′w' in WW, the sum w+w′w+w' is also in WW.

- For every vector ww in WW and every scalar λ\lambda, the scalar multiple λw\lambda w is also in WW.

In quantum physics, linear subspaces are important because they encode quantum states that a physical state can collapse to when a measurement is performed. To make this more precise, we need the notions of orthogonal complement of a subspace and of orthogonal projection onto a subspace.

Definition: Orthogonal complement Consider a Hilbert space H\mathcal{H} and a linear subspace WW of H\mathcal{H}. The orthogonal complement of WW in H\mathcal{H}, denoted by W⊥W^{\perp}, is the set of all vectors ∣v⟩\left|v\right\rangle in H\mathcal{H} that are orthogonal to all of WW, i.e.

W⊥={∣v⟩∈H∣⟨v∣w⟩=0for all w∈W}.W^{\perp}=\{\left|v\right\rangle \in\mathcal{H}\mid\left\langle v|w\right\rangle =0\quad\text{for all }w\in W\}.

It is not hard to verify that the orthogonal complement of a linear subspace of H\mathcal{H} is again a linear subspace.

Properties that are typically demonstrated in exercises:

Operators on Hilbert spaces

An operator or linear map on a Hilbert space H\mathcal{H} is a map

T ⁣:H→HT\colon\mathcal{H}\to\mathcal{H}

such that

It is also possible to define linear maps between different Hilbert spaces, but we will not go into this.

In physics, operators are often denoted with a hat, so what we call TT would be written as T^\hat{T}. One reason for this is to distinguish a physical quantity (like momentum or energy) from the corresponding mathematical operator. Since we do not consider physical quantities in this module, we do not use the hat notation.

Operators form a vector space: the sum of two operators TT and UU is defined by the formula

(T+U)∣α⟩)=T∣α⟩+U∣α⟩(T+U)\left|\alpha\right\rangle )=T\left|\alpha\right\rangle +U\left|\alpha\right\rangle

and multiplication of an operator TT by a scalar cc is defined by the formula

(cT)∣α⟩=c(T∣α⟩).(cT)\left|\alpha\right\rangle =c(T\left|\alpha\right\rangle ).

Like matrices (but unlike vectors), operators can be composed: if TT and UU are two operators, then TUTU is the operator defined by

(TU)∣α⟩=T(U∣α⟩).(TU)\left|\alpha\right\rangle =T(U\left|\alpha\right\rangle ).

Matrix entries

In physics, quantities of the form

⟨α∣T∣β⟩\left\langle \alpha\left|T\right|\beta\right\rangle

are often of special interest. These are called matrix entries of TT, especially if ∣α⟩\left|\alpha\right\rangle and ∣β⟩\left|\beta\right\rangle are basis vectors in some fixed orthonormal basis of H\mathcal{H}.

A remark on notation: you will often see notation like ⟨α∣Tβ⟩\left\langle \alpha|T\beta\right\rangle. This is literally meaningless: the operator TT can be applied to vectors, which we denote by ∣β⟩\left|\beta\right\rangle, while the symbol β\beta is just a label. However, in practice it is very convenient and completely unambiguous to write ⟨α∣Tβ⟩\left\langle \alpha|T\beta\right\rangle instead of the strictly speaking more correct ⟨α∣(T∣β⟩)\left\langle \alpha\right|(T\left|\beta\right\rangle ).

Similarly, we will write ⟨Tα∣β⟩\left\langle T\alpha|\beta\right\rangle to mean the inner product of T∣α⟩T\left|\alpha\right\rangle with ∣β⟩\left|\beta\right\rangle. Extending this notation to scalars cc, we can also write

⟨cα∣β⟩=⟨β∣cα⟩‾=cˉ⟨β∣α⟩‾=cˉ⟨α∣β⟩.\left\langle c\alpha|\beta\right\rangle =\overline{\left\langle \beta|c\alpha\right\rangle }=\bar{c}\overline{\left\langle \beta|\alpha\right\rangle }=\bar{c}\left\langle \alpha|\beta\right\rangle .

The adjoint of an operator

Given an operator TT on a Hilbert space H\mathcal{H}, there exists an operator T†T^{\dagger} with the property

⟨T†α∣β⟩=⟨α∣Tβ⟩for all ∣α⟩,∣β⟩∈H.\left\langle T^{\dagger}\alpha|\beta\right\rangle =\left\langle \alpha|T\beta\right\rangle \quad\text{for all }\left|\alpha\right\rangle ,\left|\beta\right\rangle \in\mathcal{H}.

The operator T†T^{\dagger} is called the *adjoint* of TT.

Properties of the adjoint: For all operators TT, UU and scalars cc, we have

1. (cT)†=cˉT†(cT)^{\dagger}=\bar{c}T^{\dagger}

2. (TU)†=U†T†(TU)^{\dagger}=U^{\dagger}T^{\dagger}

3. (T†)†=T(T^{\dagger})^{\dagger}=T

Hermitian operators

Definition: An operator TT on a Hilbert space H\mathcal{H} is called Hermitian if it is equal to its own adjoint, i.e. T=T†.T=T^{\dagger}.

It is clear from the definition that all real symmetric matrices are Hermitian. The following two results are of fundamental importance in quantum mechanics. Together, they explain why in the mathematical formalism of quantum mechanics, Hermitian operators correspond to physical observables.

Theorem: The eigenvalues of a Hermitian operator TT are real.

Proof: Suppose λ\lambda is an eigenvalue of TT and ∣α⟩\left|\alpha\right\rangle is a corresponding eigenvector. Then we have:

λ⟨α∣α⟩=⟨α∣Tα⟩=⟨Tα∣α⟩=⟨λα∣α⟩=λˉ⟨α∣α⟩\lambda\left\langle \alpha|\alpha\right\rangle =\left\langle \alpha|T\alpha\right\rangle =\left\langle T\alpha|\alpha\right\rangle =\left\langle \lambda\alpha|\alpha\right\rangle =\bar{\lambda}\left\langle \alpha|\alpha\right\rangle

Since ⟨α∣α⟩\left\langle \alpha|\alpha\right\rangle is non-zero, we can divide by it and obtain λ=λˉ\lambda=\bar{\lambda}.

Theorem: If TT is a Hermitian operator on a Hilbert space H\mathcal{H}, then there exists an orthonormal basis of H\mathcal{H} that consists of eigenvectors for TT.

Unitary operators

Definition: An operator TT on a Hilbert space H\mathcal{H} is called unitary if the adjoint T†T^{\dagger} is a two-sided inverse of TT, i.e.

TT†=IandT†T=ITT^{\dagger}=I\quad\text{and}\quad T^{\dagger}T=I

It follows immediately from the definition that the identity operator on a Hilbert space (sending every vector to itself) is unitary, and if TT is a unitary operator, then its inverse T−1T^{-1} (which equals T†T^{\dagger}) is also unitary. Similarly, if TT and UU are unitary operators, then so is the product TUTU. This means that the collection of all unitary operators on a Hilbert space H\mathcal{H} has the mathematical structure of a group.

Orthogonal projection

Consider a Hilbert space H\mathcal{H} and a linear subspace. One can show that every vector ∣α⟩α∈H\left|\alpha\right\rangle \alpha\in\mathcal{H} can be decomposed uniquely as

∣α⟩=∣α⟩W+∣α⟩W⊥with∣α⟩W in Wand∣α⟩W⊥ in W⊥.\left|\alpha\right\rangle =\left|\alpha\right\rangle _{W}+\left|\alpha\right\rangle _{W^{\perp}}\quad\text{with}\quad\left|\alpha\right\rangle _{W}\text{ in }W\quad\text{and}\quad\left|\alpha\right\rangle _{W^{\perp}}\text{ in }W^{\perp}.

Definition: Orthogonal projection Given a Hilbert space H\mathcal{H} and a linear subspace WW of H\mathcal{H}, the orthogonal projection of H\mathcal{H} onto WW is the map that sends every vector ∣α⟩\left|\alpha\right\rangle to ∣α⟩W\left|\alpha\right\rangle _{W}.

Orthogonal projections onto linear subspaces turn out to be important examples of linear maps.

Example: Let WW be a one-dimensional linear subspace of a Hilbert space H\mathcal{H}. Then for any unit vector ∣e⟩\left|e\right\rangle in WW, the operator ∣e⟩⟨e∣\left|e\right\rangle \left\langle e\right| sending every vector ∣ψ⟩\left|\psi\right\rangle to ∣e⟩⟨e∣ψ⟩=(⟨e∣ψ⟩)∣e⟩\left|e\right\rangle \left\langle e|\psi\right\rangle =(\left\langle e|\psi\right\rangle )\left|e\right\rangle is the orthogonal projection onto WW.

Density operators

The topic of this section is usually not treated in BSc-level mathematics courses. We include it because density matrices are a fundamental mathematical concept in quantum physics and are important examples of Hermitian operators. :::

Density operators on Hilbert spaces, or density matrices, play a central role in situations where quantum states are combined with classical uncertainty. They describe mixed states (ensembles of quantum systems) rather than individual (“pure”) quantum states. Density operators are intensively used in the theory of entanglement and in quantum information theory, for example.

Definition: A density operator on a Hilbert space H\mathcal{H} is an operator ρ ⁣:H→H\rho\colon\mathcal{H}\to\mathcal{H} with the following properties:

1. ρ\rho is Hermitian;

2. ρ\rho is positive semi-definite, i.e. for all ∣α⟩∈H\left|\alpha\right\rangle \in\mathcal{H} we have ⟨α∣ρ∣α⟩≥0\left\langle \alpha\left|\rho\right|\alpha\right\rangle \ge0;

3. ρ\rho has trace 1.

A density matrix is the matrix of a density operator relative to some choice of basis of H\mathcal{H}.

When H\mathcal{H} is finite-dimensional, the trace of ρ\rho is simply the trace of the matrix of ρ\rho with respect to any choice of basis for H\mathcal{H}. When H\mathcal{H} is infinite-dimensional, the definition is somewhat more complicated.

Example: Let ∣e⟩\left|e\right\rangle be a unit vector in H\mathcal{H}. The operator ∣e⟩⟨e∣\left|e\right\rangle \left\langle e\right| sending every vector ∣ψ⟩\left|\psi\right\rangle to ∣e⟩⟨e∣ψ⟩=(⟨e∣ψ⟩)∣e⟩\left|e\right\rangle \left\langle e|\psi\right\rangle =(\left\langle e|\psi\right\rangle )\left|e\right\rangle is a density operator. :::

Definition: A density operator ρ ⁣:H→H\rho\colon\mathcal{H}\to\mathcal{H} is pure when ρ\rho is the orthogonal projection onto a 1-dimensional linear subspace of H\mathcal{H}, and mixed otherwise.

On a finite-dimensional Hilbert space, an operator ρ\rho is a density operator precisely when there exist unit vectors e1,…,ene_{1},\ldots,e_{n} and real numbers p1,…,pn≥0p_{1},\ldots,p_{n}\ge0 with p1+⋯+pn=1p_{1}+\cdots+p_{n}=1 such that ρ\rho can be expressed as the convex combination:

ρ=∑i=1npi∣e⟩⟨e∣\rho=\sum^{n}_{i=1}p_{i}\left|e\right\rangle \left\langle e\right|

of the pure density operators ∣e⟩⟨e∣\left|e\right\rangle \left\langle e\right|. Given a density operator ρ\rho, there are in general multiple ways of expressing ρ\rho as such a convex combination. By definition, ρ\rho is pure precisely when it admits a decomposition as above with n=1n=1. Furthermore, it can be shown that ρ\rho is pure precisely when it satisfies the identity ρ2=ρ\rho^{2}=\rho.

Exponential of an operator

The exponential function can be characterised by the Taylor series

exp⁡(x)=∑n=0∞1n!xn.\exp(x)=\sum^{\infty}_{n=0}\frac{1}{n!}x^{n}.

In quantum physics, it is frequently useful to apply this formula not just to (real or complex) numbers, but to operators as well. One reason is that there is an important application to differential equations.

Definition: If TT is an operator on a Hilbert space H\mathcal{H}, then the exponential of TT is the operator

exp⁡(T)=∑n=0∞1n!Tn.\exp(T)=\sum^{\infty}_{n=0}\frac{1}{n!}T^{n}.

In the case where H=Cn\mathcal{H}=\mathbb{C}^{n}, an operator TT corresponds to an n×nn\times n matrix AA. The same formula as above then defines the matrix exponential exp⁡(A)\exp(A) of AA.

Warning: Since operators do not commute in general, the exponential does not have the property exp⁡(T+U)=exp⁡(T)exp⁡(U)\exp(T+U)=\exp(T)\exp(U) as one might expect from the case of ordinary exponentials.

Properties that are typically demonstrated in exercises:

Tensor products

In the mathematical formalism of quantum mechanics, the state of a system is a (unit) vector in a Hilbert space. Likewise, a composite of two quantum systems can be described by the tensor product of the corresponding Hilbert spaces.

In this section, H1\mathcal{H}_{1} and H2\mathcal{H}_{2} are two Hilbert spaces.

Definition: The tensor product of H1\mathcal{H}_{1} and H2\mathcal{H}_{2} is a Hilbert space H1⊗H1\mathcal{H}_{1}\otimes\mathcal{H}_{1}. This space contains vectors of the form ∣α⟩∣β⟩=∣α⟩⊗∣β⟩\left|\alpha\right\rangle \left|\beta\right\rangle =\left|\alpha\right\rangle \otimes\left|\beta\right\rangle, where ∣α⟩\left|\alpha\right\rangle is in H1\mathcal{H}_{1} and ∣β⟩\left|\beta\right\rangle is in H2\mathcal{H}_{2}. These vectors are called pure tensors. Importantly, H1⊗H2\mathcal{H}_{1}\otimes\mathcal{H}_{2} also contains linear combinations of such pure tensors (possibly with infinitely many terms); these are called (non-pure) tensors.

The inner product between two pure tensors ∣α⟩∣β⟩\left|\alpha\right\rangle \left|\beta\right\rangle and ∣α′⟩∣α′⟩\left|\alpha'\right\rangle \left|\alpha'\right\rangle is by definition ⟨α∣α′⟩⟨β∣β′⟩\left\langle \alpha|\alpha'\right\rangle \left\langle \beta|\beta'\right\rangle. (Note that the first inner product is taken in H1\mathcal{H}_{1} and the second in H2\mathcal{H}_{2}!)

Note: In quantum mechanics, a non-pure tensor in H1⊗H2\mathcal{H}_{1}\otimes\mathcal{H}_{2} represents an entangled state of a composite quantum system.

We can easily write down a basis for a tensor product H1⊗H2\mathcal{H}_{1}\otimes\mathcal{H}_{2} using bases of H1\mathcal{H}_{1} and H2\mathcal{\mathcal{H}}_{2}. It turns out that the “obvious” definition works well with inner products.

Property: If {∣ei⟩}i\left\{ \left|e_{i}\right\rangle \right\} _{i} is an orthonormal basis of H1\mathcal{H}_{1} and {∣fi⟩}i\left\{ \left|f_{i}\right\rangle \right\} _{i} is an orthonormal basis of H2\mathcal{H}_{2}, then {∣ei⟩∣ej⟩}i,j\left\{ \left|e_{i}\right\rangle \left|e_{j}\right\rangle \right\} _{i,j} is an orthonormal basis of H1⊗H2\mathcal{H}_{1}\otimes\mathcal{H}_{2}.

If AA is an operator (i.e. a linear map) on H1\mathcal{H}_{1} and BB is a linear map on H2\mathcal{H}_{2}, then there is a linear map A⊗BA\otimes B on H1⊗H2\mathcal{H}_{1}\otimes\mathcal{H}_{2}. This is defined on pure tensors by

(A⊗B)(∣α⟩∣β⟩)=(A∣α⟩)(B∣β⟩).(A\otimes B)(\left|\alpha\right\rangle \left|\beta\right\rangle )=(A\left|\alpha\right\rangle )(B\left|\beta\right\rangle ).

However, just like not every tensor is a pure tensor, there are many operators on H1⊗H2\mathcal{H}_{1}\otimes\mathcal{H}_{2} that cannot be written in the form T⊗UT\otimes U (where TT and UU are operators on H1\mathcal{H}_{1} and H2\mathcal{H}_{2}, respectively).

Footnotes
  1. Einstein notation is a notational convention that implies summation over a set of indexed terms in a formula, thus achieving brevity. There are essentially three rules of Einstein summation notation, namely: 1) Repeated indices are implicitly summed over. 2) Each index can appear at most twice in any term. 3) Each term must contain identical non-repeated indices (Wolfram MathWorld).

  2. More generally, in nn dimensions, the Levi-Civita symbol is defined by: +1 if (a1,…,an)\left(a_{1},\dots,a_{n}\right)is an even permutation from (1,2,…,n)\left(1,2,\dots,n\right), -1 if is an odd permutation and 0 otherwhise.

  3. An alternative definition of the characteristic polynomial, which agrees with the one above up to a factor (−1)n(-1)^{n}, is χA(t)=det⁡(A−xI).\chi_{A}(t)=\det(A-x\mathbb{I}).

  4. Note that if x(C)=Px(B)\boldsymbol{x}^{(C)}=P\boldsymbol{x}^{(B)}and given that P−1P=IP^{-1}P=\mathbb{I}, then P−1x(C)=x(B)P^{-1}\boldsymbol{x}^{(C)}=\boldsymbol{x}^{(B)}, that is, P−1=PC→BP^{-1}=P_{C\rightarrow B}.

  5. If B=(A−λI)B=(A-\lambda I) is invertible, it means that B−1B=IB^{-1}B=I ,so v=B−10=0\boldsymbol{v}=B^{-1}0=0 This is contradictory, since by definition v\boldsymbol{v} is nonzero. If it is not invertible, then det(B)=0\text{det}\left(B\right)=0.

  6. A diagonal matrix is one that has elements only on the main diagonal Dii≠0D_{ii}\neq0, for any j≠0j\neq0, we have Dij=0D_{ij}=0.

  7. A matrix is scalar when it is a multiple cc of the identity matrix II , that is, D=cID=cI .

  8. In short, we can think of a complete space as one in which a sequence of vectors cannot converge to a vector outside the space.

  9. The orthogonal complement of WW is the set of vectors that are orthogonal to every vector in WW. For example, if WW has as its basis the vectors u=(1,0,0)\boldsymbol{u}=\left(1,0,0\right) and v=(0,1,0)v=\left(0,1,0\right), then W⊥W^{\perp} is the subspace corresponding to the zz-axis.

References
  1. Delft University of Technology. (n.d.). Mathematics for Quantum Physics. Retrieved August 31, 2026, from https://interactivetextbooks.tudelft.nl/mqp-v2/