Skip to content

Periodic Boundary Conditions

Periodic boundary conditions identify the two ends of a finite interval. Instead of putting impenetrable walls at x=0x=0 and x=Lx=L, one declares those endpoints to be the same physical point:

ψ(x+L)=ψ(x).\psi(x+L)=\psi(x).

For the ordinary kinetic-energy Hamiltonian, the derivative must match as well:

ψ′(x+L)=ψ′(x).\psi'(x+L)=\psi'(x).

This turns the interval into a circle of circumference LL. It is also the standard finite-volume regulator for free particles, scattering calculations, many-body systems, numerical diagonalization, and field-theory limits.

Periodic boundary conditions do not describe two hard walls. A particle crossing x=Lx=L reappears at x=0x=0 with the same wavefunction and derivative. There is no reflection from an endpoint, no node forced at the boundary, and no standing-wave requirement.

The Hilbert space may be represented as square-integrable functions on one period, with inner product

⟨ϕ∣ψ⟩=∫0Lϕ∗(x)ψ(x) dx.\langle\phi\vert\psi\rangle = \int_0^L \phi^*(x)\psi(x)\,dx.

The periodicity condition says which functions on the interval represent allowed states on the circle. For the free Hamiltonian

H^=−ℏ22md2dx2,\hat H =-\frac{\hbar^2}{2m}\frac{d^2}{dx^2},

the matching of both ψ\psi and ψ′\psi' makes the kinetic-energy operator compatible with probability conservation and self-adjointness.

A plane wave has the form

ψk(x)=Aeikx.\psi_k(x)=A e^{ikx}.

Periodic boundary conditions require

Aeik(x+L)=Aeikx.Ae^{ik(x+L)}=Ae^{ikx}.

For a nonzero state,

eikL=1.e^{ikL}=1.

Thus

kL=2πn,n∈Z,kL=2\pi n, \qquad n\in\mathbb Z,

and the allowed wavenumbers are

kn=2πnL.k_n=\frac{2\pi n}{L}.

The corresponding momenta are

pn=ℏkn=2πℏnL,n∈Z.p_n=\hbar k_n =\frac{2\pi\hbar n}{L}, \qquad n\in\mathbb Z.

The integer nn can be positive, negative, or zero. Positive and negative values represent opposite circulating directions around the ring.

Choose normalized modes

un(x)=1Leiknx,kn=2πnL.u_n(x) =\frac{1}{\sqrt L}e^{ik_nx}, \qquad k_n=\frac{2\pi n}{L}.

They are orthonormal:

∫0Lum∗(x)un(x) dx=δmn.\int_0^L u_m^*(x)u_n(x)\,dx =\delta_{mn}.

For the free particle, they diagonalize both momentum and energy:

p^ un=−iℏddxun=ℏknun,\hat p\,u_n =-i\hbar\frac{d}{dx}u_n =\hbar k_n u_n,

and

H^un=ℏ2kn22mun.\hat H u_n =\frac{\hbar^2k_n^2}{2m}u_n.

Thus

En=ℏ22m(2πnL)2.E_n =\frac{\hbar^2}{2m} \left( \frac{2\pi n}{L} \right)^2.

The n=0n=0 state is the constant mode with zero kinetic energy. For n≠0n\ne0, the levels nn and −n-n have the same energy but opposite momentum.

Any sufficiently regular periodic wavefunction can be expanded as

ψ(x)=∑n∈Zcnun(x),\psi(x) =\sum_{n\in\mathbb Z} c_n u_n(x),

with coefficients

cn=∫0Lun∗(x)ψ(x) dx.c_n = \int_0^L u_n^*(x)\psi(x)\,dx.

Normalization becomes

∫0L∣ψ(x)∣2 dx=∑n∈Z∣cn∣2.\int_0^L \lvert\psi(x)\rvert^2\,dx = \sum_{n\in\mathbb Z} \lvert c_n\rvert^2.

The completeness relation is

∑n∈Zun(x)un∗(x′)=δL(x−x′),\sum_{n\in\mathbb Z} u_n(x)u_n^*(x') =\delta_L(x-x'),

where δL\delta_L is the periodic delta distribution:

δL(x−x′)=1L∑n∈Zei2πn(x−x′)/L.\delta_L(x-x') = \frac{1}{L} \sum_{n\in\mathbb Z} e^{i2\pi n(x-x')/L}.

It acts like a delta function on periodic test functions. The mathematical background is ordinary Fourier Series, with the quantum interpretation supplied by the inner product and Hamiltonian.

The spacing between adjacent allowed wavenumbers is

Δk=2πL.\Delta k=\frac{2\pi}{L}.

For a smooth function F(k)F(k), a large periodic box gives the replacement

∑n∈ZF(kn)⟶L2π∫−∞∞F(k) dk\sum_{n\in\mathbb Z} F(k_n) \longrightarrow \frac{L}{2\pi} \int_{-\infty}^{\infty} F(k)\,dk

as L→∞L\to\infty. In momentum variables, since p=ℏkp=\hbar k,

∑nF(pn)⟶L2πℏ∫−∞∞F(p) dp.\sum_n F(p_n) \longrightarrow \frac{L}{2\pi\hbar} \int_{-\infty}^{\infty} F(p)\,dp.

This is the finite-volume origin of the one-dimensional density of states:

number of stateslength in k=L2π.\frac{\text{number of states}}{\text{length in }k} =\frac{L}{2\pi}.

Equivalently, the density per unit length in kk-space is 1/(2π)1/(2\pi).

Periodic boxes are often introduced only to regularize continuum expressions. One computes with discrete, normalized states at finite LL, converts sums to integrals, and then removes the regulator. The factors of LL cancel in physical quantities when states, densities, and probabilities are normalized consistently.

For the free particle,

E=ℏ2k22m.E=\frac{\hbar^2k^2}{2m}.

Counting both k>0k\gt0 and k<0k\lt0 states for E>0E\gt0 gives

dN=2L2π dk=Lπ dk.dN = 2\frac{L}{2\pi}\,dk = \frac{L}{\pi}\,dk.

Since

dEdk=ℏ2km,\frac{dE}{dk} =\frac{\hbar^2k}{m},

the energy density of states is

ρ(E)=dNdE=Lπℏm2E,E>0.\rho(E) =\frac{dN}{dE} = \frac{L}{\pi\hbar} \sqrt{\frac{m}{2E}}, \qquad E\gt0.

The divergence as E→0+E\to0^+ is a one-dimensional feature. In higher dimensions, the power of EE changes because the number of states grows with the volume of a shell in momentum space.

A periodic plane wave can carry nonzero probability current:

jn=ℏknmL.j_n = \frac{\hbar k_n}{mL}.

This is allowed because there is no wall. The current circulates around the ring.

The boundary-condition reason is visible from integration by parts. For the kinetic operator, the boundary form is proportional to

[ϕ∗(x)ψ′(x)−ϕ′∗(x)ψ(x)]0L.\left[ \phi^*(x)\psi'(x) -\phi'^*(x)\psi(x) \right]_{0}^{L}.

If both ϕ\phi and ψ\psi obey the same periodic boundary conditions, the value at LL equals the value at 00, and this expression vanishes. The condition is therefore compatible with a self-adjoint kinetic Hamiltonian. The distinction between a vanishing boundary form and equality with the adjoint domain is developed in Hermitian vs Self-Adjoint Operators.

The infinite square well on 0<x<L0\lt x\lt L uses hard-wall conditions:

ψ(0)=0,ψ(L)=0.\psi(0)=0, \qquad \psi(L)=0.

Its eigenfunctions are standing waves,

2Lsin⁡nπxL,n=1,2,3,…\sqrt{\frac{2}{L}} \sin\frac{n\pi x}{L}, \qquad n=1,2,3,\ldots

not momentum eigenstates. The periodic box instead uses traveling waves,

1Lei2πnx/L,n∈Z.\frac{1}{\sqrt L} e^{i2\pi nx/L}, \qquad n\in\mathbb Z.

Both systems have discrete spectra because the spatial domain is compact, but the physics and boundary conditions differ. The hard-wall box models confinement by impenetrable walls. The periodic box models a ring or a finite-volume regulator without boundaries.

A useful nearby generalization is

ψ(x+L)=eiθψ(x),\psi(x+L)=e^{i\theta}\psi(x),

with the same phase relation for ψ′\psi'. Then

kn=2πn+θL.k_n=\frac{2\pi n+\theta}{L}.

This is not the standard periodic condition unless θ=0\theta=0. Twisted boundary conditions appear in rings threaded by magnetic flux, Bloch theory, and finite-size diagnostics. The present page uses the untwisted case as the canonical starting point.

In many-particle physics and field theory, one often computes in a periodic box of length LL and later takes a large-volume limit. For a single particle, this limit makes the momentum spacing Δk=2π/L\Delta k=2\pi/L approach zero. For many particles, one usually keeps intensive quantities, such as particle density, fixed while sending the volume to infinity.

The periodic boundary is a regulator, not a claim that the universe is a literal loop. It reduces edge effects and preserves translation symmetry, which is why it is so common in numerical physics and continuum-limit derivations.

  • Confusing periodic boundary conditions with hard-wall boundary conditions.
  • Forgetting that nn runs over all integers, including negative values and zero.
  • Using kn=nπ/Lk_n=n\pi/L instead of kn=2πn/Lk_n=2\pi n/L for a periodic box.
  • Treating sin⁡(nπx/L)\sin(n\pi x/L) as a momentum eigenfunction in the periodic problem.
  • Dropping the factor L/(2π)L/(2\pi) when replacing sums by integrals.
  • Forgetting that derivative matching is part of the kinetic-energy domain.
  • Interpreting finite-volume levels as physical walls when the box was only a regulator.
  • 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.
  • N. W. Ashcroft and N. D. Mermin, Solid State Physics, Holt, Rinehart and Winston, 1976.
  • M. Le Bellac, Quantum and Statistical Field Theory, Oxford University Press, 1991.
  1. Derive the allowed wavenumbers for a periodic box.
Solution

Use ψ(x)=Aeikx\psi(x)=Ae^{ikx} and impose ψ(x+L)=ψ(x)\psi(x+L)=\psi(x):

Aeik(x+L)=Aeikx.Ae^{ik(x+L)}=Ae^{ikx}.

For A≠0A\ne0,

eikL=1.e^{ikL}=1.

Therefore

kL=2πn,n∈Z,kL=2\pi n, \qquad n\in\mathbb Z,

so

kn=2πnL.k_n=\frac{2\pi n}{L}.
  1. Verify the orthonormality of un(x)=L−1/2ei2πnx/Lu_n(x)=L^{-1/2}e^{i2\pi nx/L}.
Solution

Compute

∫0Lum∗(x)un(x) dx=1L∫0Lei2π(n−m)x/L dx.\int_0^L u_m^*(x)u_n(x)\,dx = \frac{1}{L} \int_0^L e^{i2\pi(n-m)x/L}\,dx.

If n=mn=m, the integral gives LL, so the answer is 11. If n≠mn\ne m,

∫0Lei2π(n−m)x/L dx=Li2π(n−m)ei2π(n−m)x/L∣0L=0,\int_0^L e^{i2\pi(n-m)x/L}\,dx = \left. \frac{L}{i2\pi(n-m)} e^{i2\pi(n-m)x/L} \right\rvert_0^L =0,

because ei2π(n−m)=1e^{i2\pi(n-m)}=1. Hence the result is δmn\delta_{mn}.

  1. Convert a large-LL sum over periodic-box wavenumbers into an integral.
Solution

The spacing is

Δk=2πL.\Delta k=\frac{2\pi}{L}.

For a smooth function F(k)F(k),

∑nF(kn)≈1Δk∫−∞∞F(k) dk=L2π∫−∞∞F(k) dk.\sum_n F(k_n) \approx \frac{1}{\Delta k} \int_{-\infty}^{\infty} F(k)\,dk = \frac{L}{2\pi} \int_{-\infty}^{\infty} F(k)\,dk.

In momentum variables, dp=ℏ dkdp=\hbar\,dk, so

∑nF(pn)≈L2πℏ∫−∞∞F(p) dp.\sum_n F(p_n) \approx \frac{L}{2\pi\hbar} \int_{-\infty}^{\infty} F(p)\,dp.
  1. Why can a periodic-box eigenstate have nonzero current while an infinite-square-well energy eigenstate has zero average current?
Solution

A periodic-box eigenstate can be a traveling wave,

un(x)=1Leiknx,u_n(x)=\frac{1}{\sqrt L}e^{ik_nx},

with current

jn=ℏknmL.j_n=\frac{\hbar k_n}{mL}.

There is no wall for the current to hit; it circulates around the identified interval. An infinite-square-well energy eigenstate is a real standing wave, proportional to sin⁡(nπx/L)\sin(n\pi x/L), formed from equal left-moving and right-moving components. The opposite currents cancel, and the hard walls forbid net flux through the endpoints.