Pauli matrices


In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices that are traceless, Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma (σ), they are occasionally denoted by tau (τ) when used in connection with isospin symmetries.

Wolfgang Pauli (1900–1958), c. 1924. Pauli received the Nobel Prize in physics in 1945, nominated by Albert Einstein, for the Pauli exclusion principle.

These matrices are named after the physicist Wolfgang Pauli. In quantum mechanics, they occur in the Pauli equation, which takes into account the interaction of the spin of a particle with an external electromagnetic field. They also represent the interaction states of two polarization filters for horizontal/vertical polarization, 45 degree polarization (right/left), and circular polarization (right/left).

Each Pauli matrix is Hermitian, and together with the identity matrix I (sometimes considered as the zeroth Pauli matrix σ0 ), the Pauli matrices form a basis for the real vector space of 2 × 2 Hermitian matrices. This means that any 2 × 2 Hermitian matrix can be written in a unique way as a linear combination of Pauli matrices, with all coefficients being real numbers.

Hermitian operators represent observables in quantum mechanics, so the Pauli matrices span the space of observables of the complex two-dimensional Hilbert space. In the context of Pauli's work, σk represents the observable corresponding to spin along the kth coordinate axis in three-dimensional Euclidean space

The Pauli matrices (after multiplication by i to make them anti-Hermitian) also generate transformations in the sense of Lie algebras: the matrices 1, 2, 3 form a basis for the real Lie algebra , which exponentiates to the special unitary group SU(2).[a] The algebra generated by the three matrices σ1, σ2, σ3 is isomorphic to the Clifford algebra of [1] and the (unital) associative algebra generated by 1, 2, 3 functions identically (is isomorphic) to that of quaternions ().

Algebraic properties

Cayley table; the entry shows the value of the row times the column.

All three of the Pauli matrices can be compacted into a single expression:


where the solution to i2 = −1 is the "imaginary unit", and δjk is the Kronecker delta, which equals +1 if j = k and 0 otherwise. This expression is useful for "selecting" any one of the matrices numerically by substituting values of j = 1, 2, 3, in turn useful when any of the matrices (but no particular one) is to be used in algebraic manipulations.

The matrices are involutory:


where I is the identity matrix.

The determinants and traces of the Pauli matrices are


from which we can deduce that each matrix σj has eigenvalues +1 and −1.

With the inclusion of the identity matrix I (sometimes denoted σ0), the Pauli matrices form an orthogonal basis (in the sense of Hilbert–Schmidt) of the Hilbert space   of 2 × 2 Hermitian matrices over  , and the Hilbert space   of all complex 2 × 2 matrices over  .

Commutation and anti-commutation relations


Commutation relations


The Pauli matrices obey the following commutation relations:


where the structure constant εjkl is the Levi-Civita symbol, and Einstein summation notation is used.

These commutation relations make the Pauli matrices the generators of a representation of the Lie algebra  

Anticommutation relations


They also satisfy the anticommutation relations:


where   is defined as   and δjk is the Kronecker delta. I denotes the 2 × 2 identity matrix.

These anti-commutation relations make the Pauli matrices the generators of a representation of the Clifford algebra for   denoted  

The usual construction of generators   of   using the Clifford algebra recovers the commutation relations above, up to unimportant numerical factors.

A few explicit commutators and anti-commutators are given below as examples:

Commutators Anticommutators

Eigenvectors and eigenvalues


Each of the (Hermitian) Pauli matrices has two eigenvalues: +1 and −1. The corresponding normalized eigenvectors are


Pauli vectors


The Pauli vector is defined by[b]   where  ,  , and   are an equivalent notation for the more familiar  ,  , and  .

The Pauli vector provides a mapping mechanism from a vector basis to a Pauli matrix basis[2] as follows:   using Einstein's summation convention.

More formally, this defines a map from   to the vector space of traceless Hermitian   matrices. This map encodes structures of   as a normed vector space and as a Lie algebra (with the cross-product as its Lie bracket) via functions of matrices, making the map an isomorphism of Lie algebras. This makes the Pauli matrices intertwiners from the point of view of representation theory.

Another way to view the Pauli vector is as a   Hermitian traceless matrix-valued dual vector, that is, an element of   that maps  

Completeness relation


Each component of   can be recovered from the matrix (see completeness relation below)   This constitutes an inverse to the map  , making it manifest that the map is a bijection.



The norm is given by the determinant (up to a minus sign)   Then, considering the conjugation action of an   matrix   on this space of matrices,


we find   and that   is Hermitian and traceless. It then makes sense to define   where   has the same norm as   and therefore interpret   as a rotation of three-dimensional space. In fact, it turns out that the special restriction on   implies that the rotation is orientation preserving. This allows the definition of a map   given by


where   This map is the concrete realization of the double cover of   by   and therefore shows that   The components of   can be recovered using the tracing process above:




The cross-product is given by the matrix commutator (up to a factor of  )   In fact, the existence of a norm follows from the fact that   is a Lie algebra (see Killing form).

This cross-product can be used to prove the orientation-preserving property of the map above.

Eigenvalues and eigenvectors


The eigenvalues of   are   This follows immediately from tracelessness and explicitly computing the determinant.

More abstractly, without computing the determinant, which requires explicit properties of the Pauli matrices, this follows from   since this can be factorised into   A standard result in linear algebra (a linear map that satisfies a polynomial equation written in distinct linear factors is diagonal) means this implies   is diagonal with possible eigenvalues   The tracelessness of   means it has exactly one of each eigenvalue.

Its normalized eigenvectors are   These expressions become singular for  . They can be rescued by letting   and taking the limit  , which yields the correct eigenvectors (0,1) and (1,0) of  .

Alternatively, one may use spherical coordinates   to obtain the eigenvectors   and  .

Pauli 4-vector


The Pauli 4-vector, used in spinor theory, is written   with components


This defines a map from   to the vector space of Hermitian matrices,


which also encodes the Minkowski metric (with mostly minus convention) in its determinant:


This 4-vector also has a completeness relation. It is convenient to define a second Pauli 4-vector


and allow raising and lowering using the Minkowski metric tensor. The relation can then be written  

Similarly to the Pauli 3-vector case, we can find a matrix group that acts as isometries on   in this case the matrix group is   and this shows   Similarly to above, this can be explicitly realized for   with components


In fact, the determinant property follows abstractly from trace properties of the   For   matrices, the following identity holds:


That is, the 'cross-terms' can be written as traces. When   are chosen to be different   the cross-terms vanish. It then follows, now showing summation explicitly,   Since the matrices are   this is equal to  

Relation to dot and cross product


Pauli vectors elegantly map these commutation and anticommutation relations to corresponding vector products. Adding the commutator to the anticommutator gives


so that,


Contracting each side of the equation with components of two 3-vectors ap and bq (which commute with the Pauli matrices, i.e., apσq = σqap) for each matrix σq and vector component ap (and likewise with bq) yields


Finally, translating the index notation for the dot product and cross product results in



If i is identified with the pseudoscalar σxσyσz then the right hand side becomes  , which is also the definition for the product of two vectors in geometric algebra.

If we define the spin operator as J = ħ/2σ, then J satisfies the commutation relation: Or equivalently, the Pauli vector satisfies: 

Some trace relations


The following traces can be derived using the commutation and anticommutation relations.


If the matrix σ0 = I is also considered, these relationships become


where Greek indices α, β, γ and μ assume values from {0, x, y, z} and the notation   is used to denote the sum over the cyclic permutation of the included indices.

Exponential of a Pauli vector




one has, for even powers, 2p, p = 0, 1, 2, 3, ...


which can be shown first for the p = 1 case using the anticommutation relations. For convenience, the case p = 0 is taken to be I by convention.

For odd powers, 2q + 1, q = 0, 1, 2, 3, ...


Matrix exponentiating, and using the Taylor series for sine and cosine,


In the last line, the first sum is the cosine, while the second sum is the sine; so, finally,



which is analogous to Euler's formula, extended to quaternions.

Note that


while the determinant of the exponential itself is just 1, which makes it the generic group element of SU(2).

A more abstract version of formula (2) for a general 2 × 2 matrix can be found in the article on matrix exponentials. A general version of (2) for an analytic (at a and −a) function is provided by application of Sylvester's formula,[3]


The group composition law of SU(2)


A straightforward application of formula (2) provides a parameterization of the composition law of the group SU(2).[c] One may directly solve for c in  

which specifies the generic group multiplication, where, manifestly,   the spherical law of cosines. Given c, then,  

Consequently, the composite rotation parameters in this group element (a closed form of the respective BCH expansion in this case) simply amount to[4]


(Of course, when   is parallel to  , so is  , and c = a + b.)

Adjoint action


It is also straightforward to likewise work out the adjoint action on the Pauli vector, namely rotation of any angle   along any axis  :  

Taking the dot product of any unit vector with the above formula generates the expression of any single qubit operator under any rotation. For example, it can be shown that  .

Completeness relation


An alternative notation that is commonly used for the Pauli matrices is to write the vector index k in the superscript, and the matrix indices as subscripts, so that the element in row α and column β of the k-th Pauli matrix is σ kαβ.

In this notation, the completeness relation for the Pauli matrices can be written


The fact that the Pauli matrices, along with the identity matrix I, form an orthogonal basis for the Hilbert space of all 2 × 2 complex Hermitian matrices means that we can express any Hermitian matrix M as   where c is a complex number, and a is a 3-component, complex vector. It is straightforward to show, using the properties listed above, that   where "tr" denotes the trace, and hence that   which can be rewritten in terms of matrix indices as   where summation over the repeated indices is implied γ and δ. Since this is true for any choice of the matrix M, the completeness relation follows as stated above. Q.E.D.

As noted above, it is common to denote the 2 × 2 unit matrix by σ0, so σ0αβ = δαβ. The completeness relation can alternatively be expressed as  

The fact that any Hermitian complex 2 × 2 matrices can be expressed in terms of the identity matrix and the Pauli matrices also leads to the Bloch sphere representation of 2 × 2 mixed states’ density matrix, (positive semidefinite 2 × 2 matrices with unit trace. This can be seen by first expressing an arbitrary Hermitian matrix as a real linear combination of {σ0, σ1, σ2, σ3} as above, and then imposing the positive-semidefinite and trace 1 conditions.

For a pure state, in polar coordinates,   the idempotent density matrix  

acts on the state eigenvector   with eigenvalue +1, hence it acts like a projection operator.

Relation with the permutation operator


Let Pjk be the transposition (also known as a permutation) between two spins σj and σk living in the tensor product space  ,


This operator can also be written more explicitly as Dirac's spin exchange operator,


Its eigenvalues are therefore[d] 1 or −1. It may thus be utilized as an interaction term in a Hamiltonian, splitting the energy eigenvalues of its symmetric versus antisymmetric eigenstates.



The group SU(2) is the Lie group of unitary 2 × 2 matrices with unit determinant; its Lie algebra is the set of all 2 × 2 anti-Hermitian matrices with trace 0. Direct calculation, as above, shows that the Lie algebra   is the three-dimensional real algebra spanned by the set {k}. In compact notation,


As a result, each j can be seen as an infinitesimal generator of SU(2). The elements of SU(2) are exponentials of linear combinations of these three generators, and multiply as indicated above in discussing the Pauli vector. Although this suffices to generate SU(2), it is not a proper representation of su(2), as the Pauli eigenvalues are scaled unconventionally. The conventional normalization is λ = 1/2, so that


As SU(2) is a compact group, its Cartan decomposition is trivial.



The Lie algebra   is isomorphic to the Lie algebra  , which corresponds to the Lie group SO(3), the group of rotations in three-dimensional space. In other words, one can say that the j are a realization (and, in fact, the lowest-dimensional realization) of infinitesimal rotations in three-dimensional space. However, even though   and   are isomorphic as Lie algebras, SU(2) and SO(3) are not isomorphic as Lie groups. SU(2) is actually a double cover of SO(3), meaning that there is a two-to-one group homomorphism from SU(2) to SO(3), see relationship between SO(3) and SU(2).



The real linear span of {I, 1, 2, 3} is isomorphic to the real algebra of quaternions,  , represented by the span of the basis vectors   The isomorphism from   to this set is given by the following map (notice the reversed signs for the Pauli matrices):  

Alternatively, the isomorphism can be achieved by a map using the Pauli matrices in reversed order,[5]


As the set of versors U  forms a group isomorphic to SU(2), U gives yet another way of describing SU(2). The two-to-one homomorphism from SU(2) to SO(3) may be given in terms of the Pauli matrices in this formulation.



Classical mechanics


In classical mechanics, Pauli matrices are useful in the context of the Cayley-Klein parameters.[6] The matrix P corresponding to the position   of a point in space is defined in terms of the above Pauli vector matrix,


Consequently, the transformation matrix Qθ for rotations about the x-axis through an angle θ may be written in terms of Pauli matrices and the unit matrix as[6]


Similar expressions follow for general Pauli vector rotations as detailed above.

Quantum mechanics


In quantum mechanics, each Pauli matrix is related to an angular momentum operator that corresponds to an observable describing the spin of a spin 12 particle, in each of the three spatial directions. As an immediate consequence of the Cartan decomposition mentioned above, j are the generators of a projective representation (spin representation) of the rotation group SO(3) acting on non-relativistic particles with spin 12. The states of the particles are represented as two-component spinors. In the same way, the Pauli matrices are related to the isospin operator.

An interesting property of spin 12 particles is that they must be rotated by an angle of 4π in order to return to their original configuration. This is due to the two-to-one correspondence between SU(2) and SO(3) mentioned above, and the fact that, although one visualizes spin up/down as the north–south pole on the 2-sphere S2, they are actually represented by orthogonal vectors in the two-dimensional complex Hilbert space.

For a spin 12 particle, the spin operator is given by J = ħ/2σ, the fundamental representation of SU(2). By taking Kronecker products of this representation with itself repeatedly, one may construct all higher irreducible representations. That is, the resulting spin operators for higher spin systems in three spatial dimensions, for arbitrarily large j, can be calculated using this spin operator and ladder operators. They can be found in Rotation group SO(3) § A note on Lie algebras. The analog formula to the above generalization of Euler's formula for Pauli matrices, the group element in terms of spin matrices, is tractable, but less simple.[7]

Also useful in the quantum mechanics of multiparticle systems, the general Pauli group Gn is defined to consist of all n-fold tensor products of Pauli matrices.

Relativistic quantum mechanics


In relativistic quantum mechanics, the spinors in four dimensions are 4 × 1 (or 1 × 4) matrices. Hence the Pauli matrices or the Sigma matrices operating on these spinors have to be 4 × 4 matrices. They are defined in terms of 2 × 2 Pauli matrices as


It follows from this definition that the   matrices have the same algebraic properties as the σk matrices.

However, relativistic angular momentum is not a three-vector, but a second order four-tensor. Hence   needs to be replaced by Σμν , the generator of Lorentz transformations on spinors. By the antisymmetry of angular momentum, the Σμν are also antisymmetric. Hence there are only six independent matrices.

The first three are the   The remaining three,   where the Dirac αk matrices are defined as


The relativistic spin matrices Σμν are written in compact form in terms of commutator of gamma matrices as


Quantum information


In quantum information, single-qubit quantum gates are 2 × 2 unitary matrices. The Pauli matrices are some of the most important single-qubit operations. In that context, the Cartan decomposition given above is called the "Z–Y decomposition of a single-qubit gate". Choosing a different Cartan pair gives a similar "X–Y decomposition of a single-qubit gate.

See also



  1. ^ This conforms to the convention in mathematics for the matrix exponential, ⟼ exp(). In the convention in physics, σ ⟼ exp(−), hence in it no pre-multiplication by i is necessary to land in SU(2).
  2. ^ The Pauli vector is a formal device. It may be thought of as an element of  , where the tensor product space is endowed with a mapping   induced by the dot product on  
  3. ^ The relation among a, b, c, n, m, k derived here in the 2 × 2 representation holds for all representations of SU(2), being a group identity. Note that, by virtue of the standard normalization of that group's generators as half the Pauli matrices, the parameters a,b,c correspond to half the rotation angles of the rotation group. That is, the Gibbs formula linked amounts to  .
  4. ^ Explicitly, in the convention of "right-space matrices into elements of left-space matrices", it is  


  1. ^ Gull, S. F.; Lasenby, A. N.; Doran, C. J. L. (January 1993). "Imaginary numbers are not Real – the geometric algebra of spacetime" (PDF). Found. Phys. 23 (9): 1175–1201. Bibcode:1993FoPh...23.1175G. doi:10.1007/BF01883676. S2CID 14670523. Retrieved 5 May 2023 – via
  2. ^ See the spinor map.
  3. ^ Nielsen, Michael A.; Chuang, Isaac L. (2000). Quantum Computation and Quantum Information. Cambridge, UK: Cambridge University Press. ISBN 978-0-521-63235-5. OCLC 43641333.
  4. ^ Gibbs, J.W. (1884). Elements of Vector Analysis. New Haven, CT. p. 67.{{cite book}}: CS1 maint: location missing publisher (link) In fact, however, the formula goes back to Olinde Rodrigues (1840), replete with half-angle: Rodrigues, Olinde (1840). "Des lois géometriques qui regissent les déplacements d' un systéme solide dans l' espace, et de la variation des coordonnées provenant de ces déplacement considérées indépendant des causes qui peuvent les produire" (PDF). J. Math. Pures Appl. 5: 380–440.
  5. ^ Nakahara, Mikio (2003). Geometry, Topology, and Physics (2nd ed.). CRC Press. p. xxii. ISBN 978-0-7503-0606-5 – via Google Books.
  6. ^ a b Goldstein, Herbert (1959). Classical Mechanics. Addison-Wesley. pp. 109–118.
  7. ^ Curtright, T L; Fairlie, D B; Zachos, C K (2014). "A compact formula for rotations as spin matrix polynomials". SIGMA. 10: 084. arXiv:1402.3541. Bibcode:2014SIGMA..10..084C. doi:10.3842/SIGMA.2014.084. S2CID 18776942.


  • The Pauli spin matrices - The Feynman Lectures on Physics
  • Liboff, Richard L. (2002). Introductory Quantum Mechanics. Addison-Wesley. ISBN 0-8053-8714-5.
  • Schiff, Leonard I. (1968). Quantum Mechanics. McGraw-Hill. ISBN 978-0070552876.
  • Leonhardt, Ulf (2010). Essential Quantum Optics. Cambridge University Press. ISBN 978-0-521-14505-3.