Skip to content

Free Particle

A free particle is a nonrelativistic particle with no potential energy. In one dimension its Hamiltonian is

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

The free particle is the canonical model for continuous spectra, momentum eigenstates, plane waves, wave packets, group velocity, and the local behavior of particles far from interactions.

On the full line, the model consists of the Hilbert space L2(R)L^2(\mathbb R) and the self-adjoint Hamiltonian

H^=p^22m\hat H=\frac{\hat p^2}{2m}

on an appropriate twice-differentiable domain. Saying only that V(x)=0V(x)=0 does not specify every global problem. The same differential expression on a ring, a half-line, or a finite interval has different boundary conditions and therefore a different spectrum.

A spatially constant potential V0V_0 is dynamically equivalent to the free problem up to an energy offset:

H^=p^22m+V0.\hat H=\frac{\hat p^2}{2m}+V_0.

It leaves the spatial eigenfunctions and velocities unchanged and multiplies every evolving state by the common phase e−iV0t/ℏe^{-iV_0t/\hbar}. Potential differences, not an isolated additive constant, affect nongravitational wave mechanics.

The time-independent Schrödinger equation is

−ℏ22md2ψdx2=Eψ.-\frac{\hbar^2}{2m}\frac{d^2\psi}{dx^2} =E\psi.

For E>0E\gt0, define

k=2mEℏ.k=\frac{\sqrt{2mE}}{\hbar}.

Then the general stationary solution is

ψ(x)=Aeikx+Be−ikx.\psi(x)=Ae^{ikx}+Be^{-ikx}.

The two terms represent right-moving and left-moving momentum eigenstates. The corresponding momenta are

p=ℏk,p=−ℏk.p=\hbar k, \qquad p=-\hbar k.

The energy-momentum relation is

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

The other energy cases explain the lower edge of the spectrum. At E=0E=0,

ψ(x)=A+Bx,\psi(x)=A+Bx,

and no nonzero solution is square-integrable on the full line. For E<0E\lt0, writing E=−ℏ2κ2/(2m)E=-\hbar^2\kappa^2/(2m) gives

ψ(x)=Aeκx+Be−κx.\psi(x)=A e^{\kappa x}+B e^{-\kappa x}.

Decay at +∞+\infty requires A=0A=0, whereas decay at −∞-\infty requires B=0B=0. There is therefore no nonzero negative-energy bound state.

The operator gives the same conclusion without solving the equation. For a state in its domain,

⟨ψ∣H^∣ψ⟩=−ℏ22m∫−∞∞ψ∗ψ′′ dx=ℏ22m∫−∞∞∣ψ′∣2 dx≥0,\begin{aligned} \langle\psi|\hat H|\psi\rangle &=-\frac{\hbar^2}{2m} \int_{-\infty}^{\infty} \psi^*\psi''\,dx\\ &=\frac{\hbar^2}{2m} \int_{-\infty}^{\infty} \lvert\psi'\rvert^2\,dx \ge0, \end{aligned}

where the boundary term vanishes on the domain. Thus the full-line free Hamiltonian is a nonnegative operator. Its spectrum is [0,∞)[0,\infty) and is continuous; it has no square-normalizable energy eigenvectors.

A right-moving plane wave has the time-dependent form

ψk(x,t)=Aei(kx−ωt),ω=ℏk22m.\psi_k(x,t)=A e^{i(kx-\omega t)}, \qquad \omega=\frac{\hbar k^2}{2m}.

It is an eigenfunction of momentum:

p^ψk=−iℏddxψk=ℏk ψk.\hat p\psi_k =-i\hbar\frac{d}{dx}\psi_k =\hbar k\,\psi_k.

It is also an energy eigenfunction:

H^ψk=ℏ2k22mψk.\hat H\psi_k =\frac{\hbar^2k^2}{2m}\psi_k.

Plane waves are useful because they diagonalize both p^\hat p and the free Hamiltonian. Their normalization conventions are treated carefully in Plane Waves and Delta Normalization. They are not ordinary normalizable states on the full real line, because ∣ψk∣2\lvert\psi_k\rvert^2 is constant and

∫−∞∞∣ψk(x,t)∣2 dx\int_{-\infty}^{\infty}\lvert\psi_k(x,t)\rvert^2\,dx

diverges for nonzero amplitude.

The free particle on the full line has a continuous spectrum. The momentum label pp can take any real value, and the energy satisfies

E≥0.E\ge 0.

Because E=p2/(2m)E=p^2/(2m), the same positive energy corresponds to two momenta, pp and −p-p, in one dimension. This degeneracy records the two possible directions of motion. For the formal distinction between discrete and continuous spectra, see Discrete and Continuous Spectra.

In the momentum representation, the spectral resolution is especially simple:

H^=∫−∞∞dp p22m∣p⟩⟨p∣.\hat H = \int_{-\infty}^{\infty}dp\, \frac{p^2}{2m} |p\rangle\langle p|.

If states are labeled by energy instead, a direction label must accompany every E>0E\gt0. With pE=2mEp_E=\sqrt{2mE}, one may use ∣E,+⟩|E,+\rangle and ∣E,−⟩|E,-\rangle for the +pE+p_E and −pE-p_E branches. The Jacobian between momentum and energy normalization is discussed in Normalization Conventions.

Physical localized free-particle states are wave packets, built by superposing plane waves:

ψ(x,t)=12π∫−∞∞dk a(k)ei(kx−ω(k)t).\psi(x,t) =\frac{1}{\sqrt{2\pi}} \int_{-\infty}^{\infty} dk\, a(k)e^{i(kx-\omega(k)t)}.

The function a(k)a(k) controls the momentum distribution. A sharply localized packet requires a broad range of kk values.

With the displayed Fourier convention, a normalized packet obeys

∫−∞∞∣a(k)∣2 dk=1.\int_{-\infty}^{\infty}\lvert a(k)\rvert^2\,dk=1.

Free evolution changes only the phase of each momentum component:

a(k,t)=a(k,0)exp⁡(−iℏk22mt).a(k,t) =a(k,0) \exp\left(-i\frac{\hbar k^2}{2m}t\right).

Consequently, ∣a(k,t)∣2\lvert a(k,t)\rvert^2 and every moment of momentum are constant in time. Position-space interference among components still changes because their phases advance at different rates. This dispersive dephasing moves and generally spreads the packet even though its momentum distribution is fixed.

No normalizable free-particle state is exactly stationary. Energy eigenfunctions are generalized plane waves, while a normalizable packet necessarily contains a range of energies and changes shape or position under evolution.

A common regulator is to place the particle in a box of length LL with periodic boundary conditions:

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

Then

kn=2πnL,n∈Z,k_n=\frac{2\pi n}{L}, \qquad n\in\mathbb Z,

and normalized plane waves are

ψn(x)=1Leiknx.\psi_n(x)=\frac{1}{\sqrt L}e^{ik_nx}.

The box discretizes momenta. After computing physical quantities, one often takes L→∞L\to\infty and replaces sums by integrals:

∑n⟶L2π∫dk.\sum_n \longrightarrow \frac{L}{2\pi}\int dk.

This method is especially useful for density-of-states calculations and numerical approximations.

This “box” is a periodic cell, not an infinite square well. Hard walls select standing sine waves and different allowed wave numbers. Periodic boundaries preserve translation invariance on the circle and retain traveling-wave momentum eigenstates.

The finite-volume basis and sum-to-integral conversion are explained in Periodic Boundary Conditions.

On the full line, momentum eigenstates are often delta normalized:

⟨p∣p′⟩=δ(p−p′).\langle p\vert p'\rangle=\delta(p-p').

In position representation,

⟨x∣p⟩=12πℏeipx/ℏ.\langle x\vert p\rangle =\frac{1}{\sqrt{2\pi\hbar}}e^{ipx/\hbar}.

This convention is natural for Fourier transforms and continuum completeness:

∫−∞∞dp ∣p⟩⟨p∣=I^.\int_{-\infty}^{\infty} dp\, \lvert p\rangle\langle p\rvert=\hat I.

Delta-normalized plane waves are not physical localized particles by themselves. They are basis states used to construct normalizable packets.

The phrase “free particle” identifies the interior Hamiltonian, but global geometry and endpoint domains complete the model:

Configuration spaceTypical conditionSpectrum and modes
Full linesquare-integrable packets; generalized plane-wave basiscontinuous E≥0E\ge0
Ring of circumference LLperiodic endpointsdiscrete traveling waves kn=2πn/Lk_n=2\pi n/L
Finite intervalDirichlet hard wallsdiscrete standing waves
Half-lineself-adjoint condition at the origincontinuous sector; some Robin domains also admit a boundary-bound state

Thus a particle can have V=0V=0 throughout its allowed region and nevertheless possess a discrete spectrum. The quantization then comes from topology or boundary conditions, not from a force in the interior. Boundary Conditions develops this distinction.

Free Solutions Inside Constant-Potential Regions

Section titled “Free Solutions Inside Constant-Potential Regions”

Many scattering problems are assembled from regions where V(x)=VjV(x)=V_j is constant. The stationary equation there is

−ℏ22mψ′′=(E−Vj)ψ.-\frac{\hbar^2}{2m}\psi'' =(E-V_j)\psi.

If E>VjE\gt V_j, the local solutions are traveling or standing waves with

kj=2m(E−Vj)ℏ.k_j=\frac{\sqrt{2m(E-V_j)}}{\hbar}.

If E<VjE\lt V_j, the local wave number is imaginary and the solutions are exponential, with

κj=2m(Vj−E)ℏ.\kappa_j=\frac{\sqrt{2m(V_j-E)}}{\hbar}.

The additive constant is unobservable only when it shifts the potential everywhere. Differences among regions change the available kinetic energy and therefore the wavelength, current, reflection, and tunneling behavior. The Potential Step is the first complete matching example.

For a one-dimensional wavefunction, the probability current is

j=ℏ2mi(ψ∗dψdx−ψdψ∗dx).j =\frac{\hbar}{2mi} \left(\psi^*\frac{d\psi}{dx} -\psi\frac{d\psi^*}{dx}\right).

For ψ=Aeikx\psi=Ae^{ikx},

j=ℏkm∣A∣2.j=\frac{\hbar k}{m}\lvert A\rvert^2.

Thus eikxe^{ikx} with k>0k\gt0 carries probability to the right, while e−ikxe^{-ikx} carries probability to the left. This current interpretation is essential for scattering, where reflection and transmission are current ratios.

For a same-energy superposition

ψ(x,t)=(Aeikx+Be−ikx)e−iEt/ℏ,\psi(x,t) = \left(Ae^{ikx}+Be^{-ikx}\right)e^{-iEt/\hbar},

the cross terms cancel in the current, giving

j=ℏkm(∣A∣2−∣B∣2).j =\frac{\hbar k}{m} \left(\lvert A\rvert^2-\lvert B\rvert^2\right).

Equal right- and left-moving intensities form a standing wave with zero net current, even though its probability density is spatially modulated. Conversely, the constant modulus of one exact plane wave should not be interpreted as a normalizable uniform probability distribution on the infinite line; it is a property of a generalized eigenfunction.

The full-line Hamiltonian is translation invariant:

[H^,p^]=0.[\hat H,\hat p]=0.

Momentum is therefore conserved. The Hamiltonian also commutes with parity, which sends pp to −p-p. This symmetry explains why the two propagation directions have equal energy. One may combine them into parity-even and parity-odd generalized eigenfunctions proportional to cos⁡(kx)\cos(kx) and sin⁡(kx)\sin(kx), respectively.

For a spinless particle with no fields, time reversal acts by complex conjugation in position space. It changes eikxe^{ikx} into e−ikxe^{-ikx} while leaving the energy unchanged. Translation, parity, and time-reversal statements depend on the full domain: a hard wall, for example, breaks continuous translation symmetry even where the interior potential vanishes.

The dispersion relation is

ω(k)=ℏk22m.\omega(k)=\frac{\hbar k^2}{2m}.

The group velocity of a wave packet centered at k0k_0 is

vg=dωdk∣k0=ℏk0m=p0m.v_g=\left.\frac{d\omega}{dk}\right\rvert_{k_0} =\frac{\hbar k_0}{m} =\frac{p_0}{m}.

This matches the classical velocity. The phase velocity is

vph=ωk=ℏk2m=p2m,v_{\text{ph}}=\frac{\omega}{k} =\frac{\hbar k}{2m} =\frac{p}{2m},

which is not the particle velocity. The group velocity carries the packet envelope.

The free-particle Ehrenfest relations are exact, not merely narrow-packet approximations. In the Heisenberg picture,

dp^Hdt=0,dx^Hdt=p^Hm,\frac{d\hat p_H}{dt}=0, \qquad \frac{d\hat x_H}{dt}=\frac{\hat p_H}{m},

so

p^H(t)=p^,x^H(t)=x^+tmp^.\hat p_H(t)=\hat p, \qquad \hat x_H(t)=\hat x+\frac{t}{m}\hat p.

For any state with the required moments,

⟨x⟩t=⟨x⟩0+tm⟨p⟩0.\langle x\rangle_t =\langle x\rangle_0 +\frac{t}{m}\langle p\rangle_0.

The position variance obeys

(Δx)t2=(Δx)02+tm⟨Δx Δp+Δp Δx⟩0+t2m2(Δp)02.\begin{aligned} (\Delta x)^2_t &=(\Delta x)^2_0 +\frac{t}{m} \left\langle \Delta x\,\Delta p +\Delta p\,\Delta x \right\rangle_0\\ &\quad +\frac{t^2}{m^2}(\Delta p)^2_0. \end{aligned}

An initially position-momentum-correlated packet can contract for a while, but the positive quadratic term dominates at large times whenever Δp≠0\Delta p\ne0. A square-normalizable packet of finite position width cannot have Δp=0\Delta p=0, so asymptotic spreading is unavoidable. Wave Packet Spreading develops this result, while Gaussian Wave Packets gives the standard analytic example.

The Fourier integral is one representation of the full initial-value solution. The same unitary evolution can be written with a position-space kernel. Free-Particle Propagator: First Encounter states that kernel and explains its composition law; the general propagator theory remains in the Dynamics and Formulations volume.

  • Treating a plane wave as a normalizable state on the full line.
  • Claiming that a normalization constant can make an exact full-line plane wave square-integrable.
  • Forgetting that positive energy in one dimension corresponds to left-moving and right-moving momentum states.
  • Identifying the sign of energy with the direction of motion; both momentum signs have the same positive energy.
  • Confusing phase velocity with particle velocity.
  • Ignoring the normalization convention when comparing plane-wave amplitudes.
  • Thinking a free particle must be perfectly delocalized; localized free particles are wave packets.
  • Calling a normalizable packet an energy eigenstate.
  • Forgetting that wave packets spread because ω(k)\omega(k) is nonlinear.
  • Assuming every free packet broadens immediately; a correlated packet can initially contract.
  • Concluding that V=0V=0 guarantees a continuous spectrum without specifying the configuration space and boundary conditions.
  • Treating a constant potential in one region as a globally irrelevant energy shift.
  • 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.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, 2nd ed., World Scientific, 2014.
  • J. R. Taylor, Scattering Theory: The Quantum Theory of Nonrelativistic Collisions, Dover, 2006.
  1. Verify that eikxe^{ikx} is an eigenfunction of both p^\hat p and H^=p^2/(2m)\hat H=\hat p^2/(2m), and find the eigenvalues.
Solution

For momentum,

p^eikx=−iℏddxeikx=ℏkeikx.\hat p e^{ikx} =-i\hbar\frac{d}{dx}e^{ikx} =\hbar k e^{ikx}.

For energy,

H^eikx=−ℏ22md2dx2eikx=ℏ2k22meikx.\hat H e^{ikx} =-\frac{\hbar^2}{2m}\frac{d^2}{dx^2}e^{ikx} =\frac{\hbar^2k^2}{2m}e^{ikx}.
  1. A periodic box has length LL. Show that periodicity implies kn=2πn/Lk_n=2\pi n/L, and derive the one-dimensional continuum replacement for ∑nf(kn)\sum_n f(k_n).
Solution

For ψ(x)=eikx\psi(x)=e^{ikx}, periodicity requires

eik(x+L)=eikx.e^{ik(x+L)}=e^{ikx}.

Thus eikL=1e^{ikL}=1, so kL=2πnkL=2\pi n for n∈Zn\in\mathbb Z. Therefore kn=2πn/Lk_n=2\pi n/L. Adjacent values have spacing Δk=2π/L\Delta k=2\pi/L. When ff varies slowly on that spacing,

∑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.
  1. Use positivity of the free Hamiltonian to show that it has no nonzero square-integrable eigenstate with E≤0E\le0 on the full line.
Solution

For an eigenstate in the Hamiltonian domain, integration by parts gives

E⟨ψ∣ψ⟩=⟨ψ∣H^∣ψ⟩=ℏ22m∫−∞∞∣ψ′(x)∣2 dx≥0.E\langle\psi|\psi\rangle =\langle\psi|\hat H|\psi\rangle =\frac{\hbar^2}{2m} \int_{-\infty}^{\infty}\lvert\psi'(x)\rvert^2\,dx \ge0.

Therefore a nonzero eigenstate cannot have E<0E\lt0. If E=0E=0, the integral of ∣ψ′∣2\lvert\psi'\rvert^2 must vanish, so ψ\psi is constant almost everywhere. The only constant function in L2(R)L^2(\mathbb R) is the zero function. Hence there is no nonzero square-integrable eigenstate at zero energy either.

  1. For
ψ(x,t)=(Aeikx+Be−ikx)e−iEt/ℏ,\psi(x,t) = \left(Ae^{ikx}+Be^{-ikx}\right)e^{-iEt/\hbar},

calculate the probability current. Under what condition is this state a standing wave with zero net current?

Solution

Insert the wavefunction into

j=ℏ2mi(ψ∗ψ′−ψψ′∗).j=\frac{\hbar}{2mi} \left(\psi^*\psi'-\psi\psi'^*\right).

Terms proportional to A∗Be−2ikxA^*B e^{-2ikx} and its conjugate cancel. The remaining terms give

j=ℏkm(∣A∣2−∣B∣2).j=\frac{\hbar k}{m} \left(\lvert A\rvert^2-\lvert B\rvert^2\right).

The net current vanishes when ∣A∣=∣B∣\lvert A\rvert=\lvert B\rvert. The relative phase controls where the standing-wave nodes lie but does not change the zero-current condition.

  1. Let
C0=⟨Δx Δp+Δp Δx⟩0.C_0 = \left\langle \Delta x\,\Delta p +\Delta p\,\Delta x \right\rangle_0.

Use free Heisenberg evolution to find the time at which the position variance is smallest when C0<0C_0\lt0. Find that minimum variance and explain why it cannot be negative.

Solution

The exact variance is

(Δx)t2=(Δx)02+C0mt+(Δp)02m2t2.(\Delta x)^2_t =(\Delta x)^2_0 +\frac{C_0}{m}t +\frac{(\Delta p)^2_0}{m^2}t^2.

Differentiating with respect to time gives

tmin⁡=−mC02(Δp)02,t_{\min} =-\frac{mC_0}{2(\Delta p)^2_0},

which is positive when C0<0C_0\lt0. Substitution yields

(Δx)min⁡2=(Δx)02−C024(Δp)02.(\Delta x)^2_{\min} =(\Delta x)^2_0 -\frac{C_0^2}{4(\Delta p)^2_0}.

The Schrödinger–Robertson uncertainty relation implies

(Δx)02(Δp)02−C024≥ℏ24,(\Delta x)^2_0(\Delta p)^2_0 -\frac{C_0^2}{4} \ge\frac{\hbar^2}{4},

so

(Δx)min⁡2≥ℏ24(Δp)02>0.(\Delta x)^2_{\min} \ge \frac{\hbar^2}{4(\Delta p)^2_0} \gt0.