Hermitian matrix

All Hermitian matrices are diagonalizable
In mathematics, a **Hermitian matrix** (or **self-adjoint matrix**) is a complex square matrix that is equal to its own conjugate transpose—that is, the element in the i-th row and j-th column is equal to the complex conjugate of the element in the j-th row and i-th column, for all indices i and j:

A
 Hermitian

⟺

a

i
j

=

a

j
i

¯

{\displaystyle A \text{ Hermitian} \quad \iff \quad a\_{ij} = \overline{a\_{ji}}}

or in matrix form:

A
 Hermitian

⟺

# A

A

T

¯

{\displaystyle A \text{ Hermitian} \quad \iff \quad A = \overline {A^\mathsf{T}}}

.
Hermitian matrices can be understood as the complex extension of real symmetric matrices.

If the conjugate transpose of a matrix

A

{\displaystyle A}

 is denoted by

A

H

{\displaystyle A^\mathsf{H}}

, then the Hermitian property can be written concisely as

A
 Hermitian

⟺

# A

A

H

{\displaystyle A \text{ Hermitian} \quad \iff \quad A = A^\mathsf{H}}

Hermitian matrices are named after Charles Hermite, who demonstrated in 1855 that matrices of this form share a property with real symmetric matrices of always having real eigenvalues. Other, equivalent notations in common use are

A

H

=

A
†

=

A
∗

{\displaystyle A^\mathsf{H} = A^\dagger = A^\ast}

, although note that in quantum mechanics,

A
∗

{\displaystyle A^\ast}

 typically means the complex conjugate only, and not the conjugate transpose.

## Alternative characterizations

Hermitian matrices can be characterized in a number of equivalent ways, some of which are listed below:

### Equality with the adjoint

A square matrix

A

{\displaystyle A}

 is Hermitian if and only if it is equal to its adjoint, that is, it satisfies

⟨
w
,
A
v
# ⟩
⟨
A
w
,
v
⟩
,

{\displaystyle \langle w, Av\rangle = \langle Aw, v\rangle,}

for any pair of vectors

v
,
w

{\displaystyle v, w}

, where

⟨
⋅
,
⋅
⟩

{\displaystyle \langle\cdot , \cdot\rangle}

 denotes the inner product operation.
This is also the way that the more general concept of self-adjoint operator is defined.

### Reality of quadratic forms

A square matrix

A

{\displaystyle A}

 is Hermitian if and only if it is such that

⟨
v
,
A
v
⟩
∈

R

,

v
∈
V
.

{\displaystyle \langle v,Av\rangle\in\mathbb R, \quad v\in V.}

### Spectral properties

A square matrix

A

{\displaystyle A}

 is Hermitian if and only if it is unitarily diagonalizable with real eigenvalues.

## Applications

Hermitian matrices are fundamental to the quantum theory of matrix mechanics created by Werner Heisenberg, Max Born, and Pascual Jordan in 1925.

## Examples

In this section, the conjugate transpose of matrix

A

{\displaystyle A}

 is denoted as

A

H

{\displaystyle A^\mathsf{H}}

, the transpose of matrix

A

{\displaystyle A}

 is denoted as

A

T

{\displaystyle A^\mathsf{T} }

 and conjugate of matrix

A

{\displaystyle A}

 is denoted as

A
¯

{\displaystyle \overline{A}}

.

See the following example:

[

2

2
i

4

2
−
i

3

i

4

−
i

1

]

{\displaystyle \begin{bmatrix}
 2 & 2 + i & 4 \\
 2 - i & 3 & i \\
 4 & - i & 1 \\
\end{bmatrix}}

The diagonal elements must be real, as they must be their own complex conjugate.

Well-known families of Pauli matrices, Gell-Mann matrices and their generalizations are Hermitian. In theoretical physics such Hermitian matrices are often multiplied by imaginary coefficients, which results in *skew-Hermitian* matrices (see below).

Here, we offer another useful Hermitian matrix using an abstract example. If a square matrix

A

{\displaystyle A}

 equals the multiplication of a matrix and its conjugate transpose, that is,

# A
B

B

H

{\displaystyle A = BB^\mathsf{H} }

, then

A

{\displaystyle A}

 is a Hermitian positive semi-definite matrix. Furthermore, if

B

{\displaystyle B}

 is row full-rank, then

A

{\displaystyle A}

 is positive definite.

## Properties

* The entries on the main diagonal (top left to bottom right) of any Hermitian matrix are real.

*Proof:* By definition of the Hermitian matrix

H

i
j

=

H
¯

j
i

{\displaystyle H\_{ij} = \overline{H}\_{ji}}

so for *i* = *j* the above follows.
Only the main diagonal entries are necessarily real; Hermitian matrices can have arbitrary complex-valued entries in their off-diagonal elements, as long as diagonally-opposite entries are complex conjugates.
* A matrix that has only real entries is Hermitian if and only if it is symmetric. A real and symmetric matrix is simply a special case of a Hermitian matrix.

*Proof:*

H

i
j

=

H
¯

j
i

{\displaystyle H\_{ij} = \overline{H}\_{ji}}

 by definition. Thus H*ij* = H*ji* (matrix symmetry) if and only if

H

i
j

=

H
¯

i
j

{\displaystyle H\_{ij} = \overline{H}\_{ij}}

 (H*ij* is real).
* Every Hermitian matrix is a normal matrix. That is to say, AAH = AHA.

*Proof:* A = AH, so AAH = AA = AHA.
* The finite-dimensional spectral theorem says that any Hermitian matrix can be diagonalized by a unitary matrix, and that the resulting diagonal matrix has only real entries. This implies that all eigenvalues of a Hermitian matrix A with dimension n are real, and that A has n linearly independent eigenvectors. Moreover, a Hermitian matrix has orthogonal eigenvectors for distinct eigenvalues. Even if there are degenerate eigenvalues, it is always possible to find an orthogonal basis of ℂ*n* consisting of n eigenvectors of A.

* The sum of any two Hermitian matrices is Hermitian.

*Proof:*

(
A
B

)

i
j

=

A

i
j

B

i
j

=

A
¯

j
i

B
¯

j
i

=

(
A
B
)

¯

j
i

,

{\displaystyle (A + B)\_{ij} = A\_{ij} + B\_{ij} = \overline{A}\_{ji} + \overline{B}\_{ji} = \overline{(A + B)}\_{ji},}

 as claimed.
* The inverse of an invertible Hermitian matrix is Hermitian as well.

*Proof:* If

A

−
1

# A
I

{\displaystyle A^{-1}A=I}

, then

# I

I
H

=
(

A

−
1

A

)
H

=

A
H

(

A

−
1

)
H

=
A
(

A

−
1

)
H

{\displaystyle I= I^H = (A^{-1}A)^H=A^H(A^{-1})^H=A (A^{-1})^H}

, so

A

−
1

=
(

A

−
1

)
H

{\displaystyle A^{-1}=(A^{-1})^H}

 as claimed.
* The product of two Hermitian matrices A and B is Hermitian if and only if *AB* = *BA*.

*Proof:* Note that

(
A
B

)

H

=

(
A
B

)

T

¯

=

B

T

A

T

¯

=

B

T

¯

A

T

¯

=

B

H

A

H

=
B
A
.

{\displaystyle (AB)^\mathsf{H} = \overline{(AB)^\mathsf{T}} = \overline{B^\mathsf{T} A^\mathsf{T}} = \overline{B^\mathsf{T}} \overline{A^\mathsf{T}} = B^\mathsf{H} A^\mathsf{H} = BA.}

 Thus

(
A
B

)

H

=
A
B

{\displaystyle (AB)^\mathsf{H} = AB}

 if and only if

A
# B
B
A

{\displaystyle AB = BA}

.
Thus *A**n* is Hermitian if A is Hermitian and n is an integer.
* For an arbitrary complex valued vector v the product

v

H

A
v

{\displaystyle v^\mathsf{H} A v }

 is real because of

v

H

A
# v

(

v

H

A
v

)

H

{\displaystyle v^\mathsf{H} A v = \left(v^\mathsf{H} A v\right)^\mathsf{H} }

. This is especially important in quantum physics where Hermitian matrices are operators that measure properties of a system e.g. total spin which have to be real.

* The Hermitian complex n-by-n matrices do not form a vector space over the complex numbers, ℂ, since the identity matrix *I**n* is Hermitian, but *i* *I**n* is not. However the complex Hermitian matrices *do* form a vector space over the real numbers ℝ. In the 2*n*2-dimensional vector space of complex *n* × *n* matrices over ℝ, the complex Hermitian matrices form a subspace of dimension *n*2. If *E**jk* denotes the n-by-n matrix with a 1 in the *j*,*k* position and zeros elsewhere, a basis (orthonormal w.r.t. the Frobenius inner product) can be described as follows:

E

j
j

 for
1
≤
j
≤
n

(
n
 matrices
)

{\displaystyle E\_{jj} \text{ for } 1 \leq j \leq n \quad (n \text{ matrices}) }

together with the set of matrices of the form

1

2

(

E

j
k

E

k
j

)

 for
1
≤
j
<
k
≤
n

(

n
2

−
n

2

 matrices

)

{\displaystyle \frac{1}{\sqrt{2}}\left(E\_{jk} + E\_{kj}\right) \text{ for } 1 \leq j < k \leq n \quad \left( \frac{n^2-n} 2 \text{ matrices} \right) }

and the matrices

i

2

(

E

j
k

−

E

k
j

)

 for
1
≤
j
<
k
≤
n

(

n
2

−
n

2

 matrices

)

{\displaystyle \frac{i}{\sqrt{2}}\left(E\_{jk} - E\_{kj}\right) \text{ for } 1 \leq j < k \leq n \quad \left( \frac{n^2-n} 2 \text{ matrices} \right) }

where

i

{\displaystyle i}

 denotes the complex number

−
1

{\displaystyle \sqrt{-1}}

, called the *imaginary unit*.
* If n orthonormal eigenvectors

u
1

,
…
,

u
n

{\displaystyle u\_1, \dots, u\_n}

 of a Hermitian matrix are chosen and written as the columns of the matrix U, then one eigendecomposition of A is

# A
U
Λ

U

H

{\displaystyle A = U \Lambda U^\mathsf{H}}

 where

U

U

H

=
# I

U

H

U

{\displaystyle U U^\mathsf{H} = I = U^\mathsf{H} U}

 and therefore

# A

∑
j

λ
j

u
j

u
j

H

,

{\displaystyle A = \sum\_j \lambda\_j u\_j u\_j ^\mathsf{H},}

where

λ
j

{\displaystyle \lambda\_j}

 are the eigenvalues on the diagonal of the diagonal matrix

Λ

{\displaystyle \; \Lambda }

.
* The determinant of a Hermitian matrix is real:

*Proof:*

det
(
A
# )
det

(

A

T

)

⇒

det

(

A

H

)

=

det
(
A
)

¯

{\displaystyle \det(A) = \det\left(A^\mathsf{T}\right)\quad \Rightarrow \quad \det\left(A^\mathsf{H}\right) = \overline{\det(A)} }

Therefore if

# A

A

H

⇒

det
(
A
# )

det
(
A
)

¯

{\displaystyle A = A^\mathsf{H}\quad \Rightarrow \quad \det(A) = \overline{\det(A)} }

.
(Alternatively, the determinant is the product of the matrix's eigenvalues, and as mentioned before, the eigenvalues of a Hermitian matrix are real.)
## Decomposition into Hermitian and skew-Hermitian

Additional facts related to Hermitian matrices include:

* The sum of a square matrix and its conjugate transpose

(

A

A

H

)

{\displaystyle \left(A + A^\mathsf{H}\right)}

 is Hermitian.

* The difference of a square matrix and its conjugate transpose

(

A
−

A

H

)

{\displaystyle \left(A - A^\mathsf{H}\right)}

 is skew-Hermitian (also called antihermitian). This implies that the commutator of two Hermitian matrices is skew-Hermitian.

* An arbitrary square matrix C can be written as the sum of a Hermitian matrix A and a skew-Hermitian matrix B. This is known as the Toeplitz decomposition of C.

# C
A
B

with

# A

1
2

(

C

C

H

)

and

# B

1
2

(

C
−

C

H

)

{\displaystyle C = A + B \quad\mbox{with}\quad A = \frac{1}{2}\left(C + C^\mathsf{H}\right) \quad\mbox{and}\quad B = \frac{1}{2}\left(C - C^\mathsf{H}\right)}

## Rayleigh quotient

Main articles: Rayleigh quotient
In mathematics, for a given complex Hermitian matrix *M* and nonzero vector *x*, the Rayleigh quotient

R
(
M
,
x
)

{\displaystyle R(M, x)}

, is defined as:

R
(
M
,
x
)
:=

x

H

M
x

x

H

x

{\displaystyle R(M, x) := \frac{x^\mathsf{H} M x}{x^\mathsf{H} x}}

.
For real matrices and vectors, the condition of being Hermitian reduces to that of being symmetric, and the conjugate transpose

x

H

{\displaystyle x^\mathsf{H}}

 to the usual transpose

x

T

{\displaystyle x^\mathsf{T}}

. Note that

R
(
M
,
c
x
# )
R
(
M
,
x
)

{\displaystyle R(M, c x) = R(M, x)}

 for any non-zero real scalar

c

{\displaystyle c}

. Also, recall that a Hermitian (or real symmetric) matrix has real eigenvalues.

It can be shown that, for a given matrix, the Rayleigh quotient reaches its minimum value

λ
min

{\displaystyle \lambda\_\min}

 (the smallest eigenvalue of M) when

x

{\displaystyle x}

 is

v
min

{\displaystyle v\_\min}

 (the corresponding eigenvector). Similarly,

R
(
M
,
x
)
≤

λ
max

{\displaystyle R(M, x) \leq \lambda\_\max}

 and

R
(
M
,

v
max

# )

λ
max

{\displaystyle R(M, v\_\max) = \lambda\_\max}

.

The Rayleigh quotient is used in the min-max theorem to get exact values of all eigenvalues. It is also used in eigenvalue algorithms to obtain an eigenvalue approximation from an eigenvector approximation. Specifically, this is the basis for Rayleigh quotient iteration.

The range of the Rayleigh quotient (for matrix that is not necessarily Hermitian) is called a numerical range (or spectrum in functional analysis). When the matrix is Hermitian, the numerical range is equal to the spectral norm. Still in functional analysis,

λ
max

{\displaystyle \lambda\_\max}

 is known as the spectral radius. In the context of C\*-algebras or algebraic quantum mechanics, the function that to *M* associates the Rayleigh quotient *R*(*M*, *x*) for a fixed *x* and *M* varying through the algebra would be referred to as "vector state" of the algebra.