Skip to content

Coulomb Potential

The Coulomb potential is the inverse-radius interaction between two point charges. In nonrelativistic quantum mechanics it is the central potential behind hydrogenic atoms, Coulomb scattering, Rydberg scaling, and many useful approximations in atomic physics.

This page sets up the Coulomb problem before the exact Hydrogen Atom solution. It fixes sign conventions, reduced mass, natural length and energy scales, and the distinction between attractive bound-state physics and long-range continuum scattering.

Use e>0e\gt 0 for the elementary positive charge. The charge of an electron is −e-e. For two point charges q1q_1 and q2q_2, the Coulomb potential energy in SI-compatible units is

V(r)=q1q24πϵ0r.V(r) = \frac{q_1q_2}{4\pi\epsilon_0r}.

The sign is physical:

  • if q1q2>0q_1q_2\gt 0, the interaction is repulsive;
  • if q1q2<0q_1q_2\lt 0, the interaction is attractive.

For an electron bound to a point nucleus of charge +Ze+Ze, with Z=1,2,…Z=1,2,\ldots, the relative-coordinate potential is

V(r)=−Ze24πϵ0r.V(r) = - \frac{Ze^2}{4\pi\epsilon_0r}.

It is useful to define

κZ=Ze24πϵ0,\kappa_Z = \frac{Ze^2}{4\pi\epsilon_0},

so the attractive Coulomb potential is

V(r)=−κZr.V(r)=-\frac{\kappa_Z}{r}.

The sign convention matters. Many errors in atomic calculations are really sign errors: a repulsive Coulomb potential has no hydrogenic bound states.

For two particles with masses m1,m2m_1,m_2 and relative separation

r=∣r1−r2∣,r=\lvert\mathbf r_1-\mathbf r_2\rvert,

the nonrelativistic Hamiltonian separates into center-of-mass and relative parts. The relative Hamiltonian is

H^rel=p^ 22μ+q1q24πϵ0r,\hat H_{\mathrm{rel}} = \frac{\hat{\mathbf p}^{\,2}}{2\mu} + \frac{q_1q_2}{4\pi\epsilon_0r},

where the reduced mass is

μ=m1m2m1+m2.\mu = \frac{m_1m_2}{m_1+m_2}.

For an electron and a nucleus of mass MNM_N,

μ=meMNme+MN.\mu = \frac{m_eM_N}{m_e+M_N}.

The infinite-nuclear-mass approximation replaces μ\mu by mem_e. That approximation is often close, but reduced mass is the correct parameter for precision energies and isotope dependence.

For the attractive charge-ZZ Coulomb problem, the natural length is the reduced-mass Bohr scale

aZ=4πϵ0ℏ2μZe2=ℏ2μκZ.a_Z = \frac{4\pi\epsilon_0\hbar^2}{\mu Ze^2} = \frac{\hbar^2}{\mu\kappa_Z}.

For Z=1Z=1, this is the reduced-mass Bohr radius. If one uses mem_e instead of μ\mu, it becomes the usual infinite-proton-mass Bohr radius.

The scaling is simple:

aZ∝1μZ.a_Z\propto\frac{1}{\mu Z}.

Larger nuclear charge pulls the wavefunction inward. Larger reduced mass also shrinks the length scale. This is why muonic atoms are much smaller than ordinary electronic atoms, and why isotope shifts can appear even before relativistic or radiative corrections.

The corresponding Coulomb energy scale is

EZ=ℏ22μaZ2=μZ2e42(4πϵ0)2ℏ2=μκZ22ℏ2.E_Z = \frac{\hbar^2}{2\mu a_Z^2} = \frac{\mu Z^2e^4}{2(4\pi\epsilon_0)^2\hbar^2} = \frac{\mu\kappa_Z^2}{2\hbar^2}.

For Z=1Z=1 and electron mass in the infinite-nuclear-mass approximation, this is the Rydberg energy:

EZ=Z2Ryup to reduced-mass replacement.E_Z=Z^2\mathrm{Ry} \qquad \text{up to reduced-mass replacement}.

The exact nonrelativistic attractive Coulomb bound-state energies are

En=−EZn2,n=1,2,….E_n = - \frac{E_Z}{n^2}, \qquad n=1,2,\ldots .

The full derivation and wavefunctions belong to the Hydrogen Atom page. Here the important lesson is the scale:

∣En∣∝μZ2.\lvert E_n\rvert\propto \mu Z^2.

The scales aZa_Z and EZE_Z are chosen so that the Schrödinger equation has a clean dimensionless form. Let

ρ=raZ,ϵ=EEZ.\rho=\frac{r}{a_Z}, \qquad \epsilon=\frac{E}{E_Z}.

For the attractive Coulomb Hamiltonian,

[−ℏ22μ∇2−κZr]ψ=Eψ,\left[ -\frac{\hbar^2}{2\mu}\nabla^2 - \frac{\kappa_Z}{r} \right]\psi = E\psi,

the dimensionless equation becomes

[−∇ρ2−2ρ]ψ=ϵψ.\left[ -\nabla_\rho^2 - \frac{2}{\rho} \right]\psi = \epsilon\psi.

In a fixed angular-momentum sector, the reduced radial equation becomes

[−d2dρ2+ℓ(ℓ+1)ρ2−2ρ]uℓ=ϵuℓ.\left[ -\frac{d^2}{d\rho^2} + \frac{\ell(\ell+1)}{\rho^2} - \frac{2}{\rho} \right]u_\ell = \epsilon u_\ell.

This form makes clear why every attractive point-Coulomb problem has the same dimensionless shape. The physical differences are restored by multiplying lengths by aZa_Z and energies by EZE_Z.

The attractive Coulomb potential has two spectral regimes:

  • E<0E\lt 0: an infinite ladder of normalizable bound states;
  • E>0E\gt 0: continuum states describing ionization and scattering.

The bound energies accumulate at E=0E=0 from below. Large-nn states are weakly bound and spatially extended. The continuum begins at the ionization threshold and is treated as its own overview in Continuum States of the Coulomb Problem.

The repulsive Coulomb potential has no square-normalizable bound states in the two-body point-charge problem. It still has continuum scattering states, but the long-range tail changes the asymptotic form of scattering wavefunctions. This is why Coulomb Scattering requires special conventions beyond ordinary short-range scattering.

For the attractive Coulomb case, the fixed-ℓ\ell effective radial potential is

Vℓ,eff(r)=−κZr+ℏ2ℓ(ℓ+1)2μr2.V_{\ell,\mathrm{eff}}(r) = - \frac{\kappa_Z}{r} + \frac{\hbar^2\ell(\ell+1)}{2\mu r^2}.

For ℓ=0\ell=0, there is no centrifugal barrier. For ℓ>0\ell\gt 0, the 1/r21/r^2 term dominates close to the origin and suppresses small-radius probability. The competition between attraction and the centrifugal barrier gives a qualitative radial scale before the exact polynomial solution is known.

The Effective Radial Potential page develops the turning-point and centrifugal-barrier picture. The Boundary Conditions for Radial Wavefunctions page explains why the ordinary Coulomb singularity is still compatible with regular radial boundary conditions.

The point-Coulomb Hamiltonian is a controlled starting point, not a complete atomic theory. It omits:

  • spin-orbit coupling and relativistic kinetic-energy corrections;
  • the Darwin term and full Dirac-Coulomb structure;
  • Lamb-shift radiative corrections;
  • hyperfine interactions with nuclear spin;
  • finite nuclear size and nuclear polarizability;
  • many-electron screening and correlation;
  • external electric and magnetic fields.

Those effects are smaller than the gross Coulomb binding scale in ordinary hydrogen, but they are essential for precision spectroscopy and for real many-electron atoms. They should be added as corrections to a clearly stated base model, not silently mixed into the Coulomb potential itself.

  • Forgetting that e>0e\gt 0 while the electron charge is −e-e.
  • Using the electron mass when the reduced mass is required.
  • Calling aZa_Z the Bohr radius without stating whether reduced mass and nuclear charge have been included.
  • Treating the unscreened Coulomb continuum as if it had short-range plane-wave asymptotics.
  • Assuming a repulsive Coulomb potential can have hydrogenic bound states.
  • Confusing the Rydberg energy with the Hartree energy; the Hartree is twice the Rydberg in the infinite-proton-mass hydrogen convention.
  • Folding fine structure or Lamb-shift physics into the basic nonrelativistic Coulomb Hamiltonian.
  • 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. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Pergamon, 1977.
  • H. A. Bethe and E. E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms, Springer, 1957.
  1. Check that the charge-ZZ Bohr scale has dimensions of length.
Solution

The Coulomb coupling

κZ=Ze24πϵ0\kappa_Z=\frac{Ze^2}{4\pi\epsilon_0}

has units of energy times length. The quantity

aZ=ℏ2μκZa_Z=\frac{\hbar^2}{\mu\kappa_Z}

therefore has units

(action)2(mass)(energy)(length).\frac{(\text{action})^2}{(\text{mass})(\text{energy})(\text{length})}.

Since energy has units mass times length squared over time squared, this reduces to length. Equivalently, aZa_Z is the length that makes ℏ2/(2μaZ2)\hbar^2/(2\mu a_Z^2) and κZ/aZ\kappa_Z/a_Z comparable.

  1. Derive the dimensionless attractive Coulomb equation using r=aZρr=a_Z\rho.
Solution

With r=aZρr=a_Z\rho,

∇2=1aZ2∇ρ2.\nabla^2=\frac{1}{a_Z^2}\nabla_\rho^2.

The Schrödinger equation is

[−ℏ22μaZ2∇ρ2−κZaZρ]ψ=Eψ.\left[ -\frac{\hbar^2}{2\mu a_Z^2}\nabla_\rho^2 - \frac{\kappa_Z}{a_Z\rho} \right]\psi = E\psi.

Using

EZ=ℏ22μaZ2,κZaZ=2EZ,E_Z=\frac{\hbar^2}{2\mu a_Z^2}, \qquad \frac{\kappa_Z}{a_Z}=2E_Z,

and dividing by EZE_Z gives

[−∇ρ2−2ρ]ψ=ϵψ,ϵ=EEZ.\left[ -\nabla_\rho^2 - \frac{2}{\rho} \right]\psi = \epsilon\psi, \qquad \epsilon=\frac{E}{E_Z}.
  1. If ZZ is doubled while the reduced mass is held fixed, how do the natural length and energy scales change?
Solution

The charge-ZZ length scale is

aZ∝1Z.a_Z\propto\frac{1}{Z}.

Doubling ZZ halves the length scale. The energy scale is

EZ∝Z2.E_Z\propto Z^2.

Doubling ZZ therefore multiplies the binding-energy scale by 44. The dimensionless spectrum is unchanged, but the conversion back to physical units changes.