Skip to content

Hermite Functions

Hermite functions are the normalized shape functions of the one-dimensional quantum harmonic oscillator. They are Hermite polynomials multiplied by a Gaussian and normalized in L2(R)L^2(\mathbb R).

This page is a focused reference for the eigenfunctions. The differential-equation derivation belongs to Differential-Equation Solution, and the polynomial identities belong to Hermite Polynomials.

In the dimensionless coordinate ξ\xi, define

φn(ξ)=1π1/42nn!Hn(ξ)e−ξ2/2,n=0,1,2,…\varphi_n(\xi) = \frac{1}{ \pi^{1/4}\sqrt{2^n n!} } H_n(\xi)e^{-\xi^2/2}, \qquad n=0,1,2,\ldots

Here HnH_n is the physicists’ Hermite polynomial. The first few are

H0(ξ)=1,H1(ξ)=2ξ,H2(ξ)=4ξ2−2.H_0(\xi)=1, \qquad H_1(\xi)=2\xi, \qquad H_2(\xi)=4\xi^2-2.

Thus the first Hermite functions are a Gaussian, a Gaussian times a linear function, and a Gaussian times a quadratic function.

First four normalized Hermite functions with alternating parity and increasing node count

The first four normalized Hermite functions φn(ξ)\varphi_n(\xi). The parity alternates with nn, and the number of nodes is exactly nn.

For the physical oscillator with mass mm and angular frequency ω\omega, the oscillator length is

ℓ=ℏmω.\ell=\sqrt{\frac{\hbar}{m\omega}}.

The dimensionless coordinate is

ξ=xℓ.\xi=\frac{x}{\ell}.

The normalized position-space energy eigenfunctions are scaled Hermite functions:

ψn(x)=1ℓ φn(xℓ).\psi_n(x) = \frac{1}{\sqrt{\ell}}\, \varphi_n\left(\frac{x}{\ell}\right).

Equivalently,

ψn(x)=12nn!(1πℓ2)1/4Hn(xℓ)e−x2/(2ℓ2).\psi_n(x) = \frac{1}{\sqrt{2^n n!}} \left(\frac{1}{\pi\ell^2}\right)^{1/4} H_n\left(\frac{x}{\ell}\right) e^{-x^2/(2\ell^2)}.

The factor ℓ−1/2\ell^{-1/2} is required because dx=ℓ dξdx=\ell\,d\xi. Without it, the physical wavefunction would have the wrong units and the wrong normalization.

The dimensionless Hermite functions obey

∫−∞∞φm(ξ)∗φn(ξ) dξ=δmn.\int_{-\infty}^{\infty} \varphi_m(\xi)^*\varphi_n(\xi)\,d\xi = \delta_{mn}.

Therefore the physical oscillator wavefunctions obey

∫−∞∞ψm(x)∗ψn(x) dx=δmn.\int_{-\infty}^{\infty} \psi_m(x)^*\psi_n(x)\,dx = \delta_{mn}.

The normalization follows from the weighted Hermite-polynomial identity

∫−∞∞e−ξ2Hm(ξ)Hn(ξ) dξ=π 2nn! δmn.\int_{-\infty}^{\infty} e^{-\xi^2}H_m(\xi)H_n(\xi)\,d\xi = \sqrt\pi\,2^n n!\,\delta_{mn}.

The Gaussian weight is essential. The polynomials HnH_n are not orthogonal with respect to the ordinary unweighted integral on the real line.

Hermite polynomials satisfy

Hn(−ξ)=(−1)nHn(ξ).H_n(-\xi)=(-1)^nH_n(\xi).

Because the Gaussian envelope is even,

φn(−ξ)=(−1)nφn(ξ),\varphi_n(-\xi)=(-1)^n\varphi_n(\xi),

and hence

ψn(−x)=(−1)nψn(x).\psi_n(-x)=(-1)^n\psi_n(x).

Even nn gives even functions, and odd nn gives odd functions. The nnth Hermite function has exactly nn real nodes. This agrees with the general one-dimensional bound-state rule that higher energy states have more nodes.

The Hermite function label nn is also the oscillator energy label:

En=ℏω(n+12).E_n=\hbar\omega\left(n+\frac12\right).

The wavefunction ψn\psi_n is the coordinate representation of the number state ∣n⟩\lvert n\rangle:

ψn(x)=⟨x∣n⟩.\psi_n(x)=\langle x\vert n\rangle.

Number-state notation and ladder actions are collected in Number States. This page focuses on the functions themselves.

The dimensionless Hermite functions satisfy useful ladder identities:

ξφn(ξ)=n+12 φn+1(ξ)+n2 φn−1(ξ),\xi\varphi_n(\xi) = \sqrt{\frac{n+1}{2}}\,\varphi_{n+1}(\xi) + \sqrt{\frac{n}{2}}\,\varphi_{n-1}(\xi),

and

dφndξ=n2 φn−1(ξ)−n+12 φn+1(ξ).\frac{d\varphi_n}{d\xi} = \sqrt{\frac{n}{2}}\,\varphi_{n-1}(\xi) - \sqrt{\frac{n+1}{2}}\,\varphi_{n+1}(\xi).

After restoring x=ℓξx=\ell\xi, the position operator has matrix elements

⟨m∣x^∣n⟩=ℓ[n+12 δm,n+1+n2 δm,n−1].\langle m\vert \hat x\vert n\rangle = \ell \left[ \sqrt{\frac{n+1}{2}}\,\delta_{m,n+1} + \sqrt{\frac{n}{2}}\,\delta_{m,n-1} \right].

These formulas are often faster than doing integrals directly. They are the wavefunction version of the ladder-operator algebra.

With the unitary dimensionless Fourier convention

(Ff)(k)=12π∫−∞∞e−ikξf(ξ) dξ,(\mathcal F f)(k) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} e^{-ik\xi}f(\xi)\,d\xi,

the Hermite functions are Fourier-transform eigenfunctions:

Fφn=(−i)nφn.\mathcal F\varphi_n = (-i)^n\varphi_n.

This property is one reason Hermite functions are useful in spectral methods and phase-space calculations. The phase depends on the Fourier convention; the eigenfunction structure does not.

The set

{φn:n=0,1,2,…}\{\varphi_n:n=0,1,2,\ldots\}

is a complete orthonormal basis of L2(R)L^2(\mathbb R). A square-integrable function can be expanded as

f(ξ)=∑n=0∞cnφn(ξ),f(\xi) = \sum_{n=0}^{\infty} c_n\varphi_n(\xi),

with coefficients

cn=∫−∞∞φn(ξ)∗f(ξ) dξ.c_n = \int_{-\infty}^{\infty} \varphi_n(\xi)^*f(\xi)\,d\xi.

Completeness is a Hilbert-space statement about convergence in norm. It does not mean every function is well approximated pointwise by a short Hermite expansion.

  • Confusing the Hermite polynomial HnH_n with the normalized Hermite function φn\varphi_n.
  • Dropping the factor ℓ−1/2\ell^{-1/2} when converting to ψn(x)\psi_n(x).
  • Writing the Gaussian envelope as e−x2/ℓ2e^{-x^2/\ell^2} instead of e−x2/(2ℓ2)e^{-x^2/(2\ell^2)}.
  • Forgetting that nn starts at 00.
  • Mixing physicists’ Hermite polynomials with probabilists’ He⁡n\operatorname{He}_n.
  • Treating Fourier-transform phases of Hermite functions as convention-independent.
  • Assuming a short Hermite-function expansion is automatically good for discontinuous or sharply localized functions.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, eds., NIST Handbook of Mathematical Functions, Cambridge University Press, 2010.
  • G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physicists, 7th ed., Academic Press, 2013.
  1. Show that ψn(x)=ℓ−1/2φn(x/ℓ)\psi_n(x)=\ell^{-1/2}\varphi_n(x/\ell) is normalized if φn\varphi_n is normalized.
Solution

Compute

∫−∞∞∣ψn(x)∣2 dx=∫−∞∞1ℓ∣φn(xℓ)∣2dx.\int_{-\infty}^{\infty} \lvert\psi_n(x)\rvert^2\,dx = \int_{-\infty}^{\infty} \frac{1}{\ell} \left\lvert \varphi_n\left(\frac{x}{\ell}\right) \right\rvert^2 dx.

Set ξ=x/ℓ\xi=x/\ell, so dx=ℓ dξdx=\ell\,d\xi. Then

∫−∞∞∣ψn(x)∣2 dx=∫−∞∞∣φn(ξ)∣2 dξ=1.\int_{-\infty}^{\infty} \lvert\psi_n(x)\rvert^2\,dx = \int_{-\infty}^{\infty} \lvert\varphi_n(\xi)\rvert^2\,d\xi =1.
  1. Use the parity of HnH_n to prove the parity of ψn\psi_n.
Solution

The Gaussian factor is even:

e−(−x)2/(2ℓ2)=e−x2/(2ℓ2).e^{-(-x)^2/(2\ell^2)} = e^{-x^2/(2\ell^2)}.

The Hermite polynomial satisfies

Hn(−xℓ)=(−1)nHn(xℓ).H_n\left(-\frac{x}{\ell}\right) = (-1)^nH_n\left(\frac{x}{\ell}\right).

Therefore

ψn(−x)=(−1)nψn(x).\psi_n(-x)=(-1)^n\psi_n(x).
  1. Derive ⟨m∣x^∣n⟩\langle m\vert \hat x\vert n\rangle from the recurrence relation for ξφn\xi\varphi_n.
Solution

Since x^=ℓξ\hat x=\ell\xi,

⟨m∣x^∣n⟩=ℓ∫−∞∞φm(ξ)∗ξφn(ξ) dξ.\langle m\vert \hat x\vert n\rangle = \ell \int_{-\infty}^{\infty} \varphi_m(\xi)^* \xi\varphi_n(\xi)\,d\xi.

Use

ξφn=n+12 φn+1+n2 φn−1.\xi\varphi_n = \sqrt{\frac{n+1}{2}}\,\varphi_{n+1} + \sqrt{\frac{n}{2}}\,\varphi_{n-1}.

Orthonormality gives

⟨m∣x^∣n⟩=ℓ[n+12 δm,n+1+n2 δm,n−1].\langle m\vert \hat x\vert n\rangle = \ell \left[ \sqrt{\frac{n+1}{2}}\,\delta_{m,n+1} + \sqrt{\frac{n}{2}}\,\delta_{m,n-1} \right].
  1. Why are Hermite polynomials alone not acceptable oscillator wavefunctions?
Solution

The polynomial Hn(ξ)H_n(\xi) grows like ξn\xi^n at large ∣ξ∣\lvert\xi\rvert. It is not square-integrable on the real line. The Gaussian factor e−ξ2/2e^{-\xi^2/2} supplies the decay needed for normalizability. The normalized Hermite function is the product of the polynomial, the Gaussian envelope, and the normalization constant.