Skip to content

Dimensionless Variables and Scaling

Dimensionless variables turn a physical Schrödinger problem into an equation whose essential parameters are visible. They remove irrelevant units, reveal natural length and energy scales, improve numerical conditioning, and make limiting regimes easier to recognize.

The goal is not to hide constants. The goal is to identify which combinations of constants actually control the physics.

A formula such as

eikxe^{ikx}

is meaningful because kxkx is dimensionless. Similarly, the Schrödinger phase

e−iEt/ℏe^{-iEt/\hbar}

is meaningful because Et/ℏEt/\hbar is dimensionless. Dimensionless scaling applies the same discipline to whole differential equations.

Scaling helps in four common ways:

  • it shows which parameters are physically independent;
  • it exposes small or large dimensionless quantities for approximation methods;
  • it makes eigenvalue equations cleaner;
  • it keeps numerical variables and matrix entries near natural order-one scales.

Changing units does not change the physics. Choosing natural dimensionless variables changes how clearly the physics is displayed.

Consider a one-dimensional stationary Schrödinger equation

[−ℏ22md2dx2+V(x)]ψ(x)=Eψ(x).\left[ -\frac{\hbar^2}{2m}\frac{d^2}{dx^2} +V(x) \right]\psi(x) =E\psi(x).

Choose a length scale LL and define

ξ=xL.\xi=\frac{x}{L}.

Then

ddx=1Lddξ,d2dx2=1L2d2dξ2.\frac{d}{dx} = \frac{1}{L}\frac{d}{d\xi}, \qquad \frac{d^2}{dx^2} = \frac{1}{L^2}\frac{d^2}{d\xi^2}.

The kinetic term introduces the natural energy associated with length LL:

EL=ℏ22mL2.E_L = \frac{\hbar^2}{2mL^2}.

If the potential has a characteristic size V0V_0 and shape U(ξ)U(\xi),

V(x)=V0U(ξ),V(x)=V_0U(\xi),

then the dimensionless eigenvalue equation is

[−d2dξ2+λU(ξ)]φ(ξ)=εφ(ξ),\left[ -\frac{d^2}{d\xi^2} +\lambda U(\xi) \right]\varphi(\xi) = \varepsilon\varphi(\xi),

where

λ=V0EL,ε=EEL.\lambda=\frac{V_0}{E_L}, \qquad \varepsilon=\frac{E}{E_L}.

The dimensionless parameter λ\lambda is often the real control knob. It compares potential energy to the kinetic energy forced by localization over the scale LL.

The wavefunction must also be scaled if normalization is to keep its meaning. In one dimension, define

ψ(x)=1Lφ(ξ),ξ=xL.\psi(x)=\frac{1}{\sqrt L}\varphi(\xi), \qquad \xi=\frac{x}{L}.

Then

∫−∞∞∣ψ(x)∣2 dx=∫−∞∞∣φ(ξ)∣2 dξ.\int_{-\infty}^{\infty} \lvert\psi(x)\rvert^2\,dx = \int_{-\infty}^{\infty} \lvert\varphi(\xi)\rvert^2\,d\xi.

In three dimensions,

ψ(r)=1L3/2Φ(ξ),ξ=rL,\psi(\mathbf r) = \frac{1}{L^{3/2}}\Phi(\boldsymbol\xi), \qquad \boldsymbol\xi=\frac{\mathbf r}{L},

so that

∫∣ψ(r)∣2 d3r=∫∣Φ(ξ)∣2 d3ξ.\int \lvert\psi(\mathbf r)\rvert^2\,d^3r = \int \lvert\Phi(\boldsymbol\xi)\rvert^2\,d^3\xi.

For radial central-potential problems, the scaling depends on which radial function is used:

R(r)=L−3/2R(s),u(r)=L−1/2w(s),s=rL.R(r)=L^{-3/2}\mathcal R(s), \qquad u(r)=L^{-1/2}w(s), \qquad s=\frac{r}{L}.

These powers are not optional decoration; they preserve the normalization convention.

If the natural energy scale is E0E_0, the corresponding time scale is

t0=ℏE0.t_0=\frac{\hbar}{E_0}.

With

τ=tt0,\tau=\frac{t}{t_0},

the phase Et/ℏEt/\hbar becomes

Etℏ=EE0τ.\frac{Et}{\hbar} = \frac{E}{E_0}\tau.

Thus the dimensionless energy ε=E/E0\varepsilon=E/E_0 also controls the rate of phase evolution in dimensionless time.

For a particle in an infinite square well of width LL, the natural kinetic scale is

EL=ℏ22mL2.E_L=\frac{\hbar^2}{2mL^2}.

Using ξ=x/L\xi=x/L on 0<ξ<10<\xi<1, the stationary equation becomes

−d2φdξ2=εφ,φ(0)=φ(1)=0.-\frac{d^2\varphi}{d\xi^2} = \varepsilon\varphi, \qquad \varphi(0)=\varphi(1)=0.

The dimensionless eigenvalues are

εn=n2π2,n=1,2,….\varepsilon_n=n^2\pi^2, \qquad n=1,2,\ldots.

Restoring units gives

En=ℏ2π2n22mL2.E_n = \frac{\hbar^2\pi^2n^2}{2mL^2}.

The scaling makes the main dependence immediate: all energies scale as 1/L21/L^2.

For the harmonic oscillator,

V(x)=12mω2x2.V(x)=\frac12m\omega^2x^2.

The natural length is found by balancing kinetic localization energy with oscillator potential energy:

ℏ22mℓ2∼12mω2ℓ2.\frac{\hbar^2}{2m\ell^2} \sim \frac12m\omega^2\ell^2.

This gives the oscillator length

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

With ξ=x/ℓ\xi=x/\ell and energy in units of ℏω\hbar\omega, the stationary equation becomes

[−12d2dξ2+12ξ2]φ=εφ,ε=Eℏω.\left[ -\frac12\frac{d^2}{d\xi^2} +\frac12\xi^2 \right]\varphi = \varepsilon\varphi, \qquad \varepsilon=\frac{E}{\hbar\omega}.

Some derivations multiply this equation by 22 and use 2E/(ℏω)2E/(\hbar\omega) as the dimensionless eigenvalue. Both conventions are equivalent. The important scales are ℓ\ell and ℏω\hbar\omega.

For the hydrogenic Coulomb problem with reduced mass μ\mu and charge magnitude ee, the Hamiltonian is

H^=−ℏ22μ∇2−e24πϵ0r.\hat H = -\frac{\hbar^2}{2\mu}\nabla^2 - \frac{e^2}{4\pi\epsilon_0r}.

The natural length is the reduced-mass Bohr radius

a0=4πϵ0ℏ2μe2.a_0 = \frac{4\pi\epsilon_0\hbar^2}{\mu e^2}.

The corresponding Hartree-scale energy is

EH=ℏ2μa02=μe4(4πϵ0)2ℏ2.E_{\mathrm H} = \frac{\hbar^2}{\mu a_0^2} = \frac{\mu e^4}{(4\pi\epsilon_0)^2\hbar^2}.

The bound-state energies are then simply

En=−EH2n2.E_n = -\frac{E_{\mathrm H}}{2n^2}.

The Coulomb scaling explains why changing reduced mass or nuclear charge changes atomic sizes and binding energies in systematic ways. The canonical hydrogen solution is developed in Hydrogen Atom.

For a charged particle in a uniform magnetic field, the natural length is the magnetic length

ℓB=ℏ∣q∣B.\ell_B = \sqrt{\frac{\hbar}{\lvert q\rvert B}}.

The corresponding frequency scale is the cyclotron frequency

ωc=∣q∣Bm,\omega_c = \frac{\lvert q\rvert B}{m},

and the energy scale is

ℏωc.\hbar\omega_c.

The magnetic length measures the spatial scale on which magnetic flux and quantum phase compete. It is the natural scale for Landau levels and many quantum Hall calculations.

Dimensionless variables are not only aesthetic. They make numerical calculations more stable and easier to audit.

For example, a discretized Hamiltonian with entries near 10−3410^{-34} in SI units may be mathematically correct but numerically awkward. After scaling by a natural energy, the same matrix may have entries of order 11 to 10310^3, with residuals and tolerances that can be interpreted physically.

A numerical eigenpair residual such as

r=Hv−Evr=H\mathbf v-E\mathbf v

should be compared with a scale:

∥r∥E0∥v∥.\frac{\lVert r\rVert}{E_0\lVert\mathbf v\rVert}.

Without E0E_0, a statement such as “the residual is 10−910^{-9}” is incomplete, because the units and physical scale of the residual have not been specified.

Scaling also helps identify stiffness in time-dependent problems. If the largest energy scale in a Hamiltonian is Emax⁡E_{\max}, then the shortest phase-evolution time scale is roughly

ℏEmax⁡.\frac{\hbar}{E_{\max}}.

Time-stepping choices must resolve that scale if those high-energy components matter.

  • Setting ℏ=1\hbar=1 without saying which length, mass, charge, or energy scales are being held fixed.
  • Scaling the coordinate but forgetting the wavefunction normalization factor.
  • Treating a dimensionful quantity inside an exponential, logarithm, or trigonometric function as meaningful.
  • Confusing a unit convention with a physical approximation.
  • Hiding the true small parameter by choosing variables that make it look like 11.
  • Reporting numerical residuals or tolerances without a reference scale.
  • Mixing two scaling conventions in the same derivation.

When scaling a Schrödinger problem:

  1. identify the physical domain and boundary conditions;
  2. choose a natural length LL from geometry, balance of terms, or external fields;
  3. define a kinetic energy scale ℏ2/(2mL2)\hbar^2/(2mL^2) or another natural E0E_0;
  4. rescale the wavefunction so normalization is preserved;
  5. divide the equation by E0E_0;
  6. name the remaining dimensionless parameters;
  7. solve or approximate the dimensionless problem;
  8. restore units only after the dimensionless result is understood.

The best scaling is the one that makes the important physics hard to miss.

  1. Rescale the infinite square well on 0<x<L0<x<L using ξ=x/L\xi=x/L and show that the dimensionless energies are n2π2n^2\pi^2 when EL=ℏ2/(2mL2)E_L=\hbar^2/(2mL^2).
Solution

The stationary equation is

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

With ξ=x/L\xi=x/L,

d2dx2=1L2d2dξ2.\frac{d^2}{dx^2} = \frac{1}{L^2}\frac{d^2}{d\xi^2}.

Dividing by EL=ℏ2/(2mL2)E_L=\hbar^2/(2mL^2) gives

−d2φdξ2=εφ,ε=EEL.-\frac{d^2\varphi}{d\xi^2} = \varepsilon\varphi, \qquad \varepsilon=\frac{E}{E_L}.

The boundary conditions are φ(0)=φ(1)=0\varphi(0)=\varphi(1)=0, so

φn(ξ)=2sin⁡(nπξ),εn=n2π2.\varphi_n(\xi)=\sqrt2\sin(n\pi\xi), \qquad \varepsilon_n=n^2\pi^2.
  1. Balance the kinetic localization scale ℏ2/(2mℓ2)\hbar^2/(2m\ell^2) with the oscillator potential scale mω2ℓ2/2m\omega^2\ell^2/2 and solve for ℓ\ell.
Solution

Set

ℏ22mℓ2=12mω2ℓ2.\frac{\hbar^2}{2m\ell^2} = \frac12m\omega^2\ell^2.

Multiplying by 2mℓ22m\ell^2 gives

ℏ2=m2ω2ℓ4.\hbar^2=m^2\omega^2\ell^4.

Therefore

ℓ2=ℏmω,ℓ=ℏmω.\ell^2=\frac{\hbar}{m\omega}, \qquad \ell=\sqrt{\frac{\hbar}{m\omega}}.
  1. In one dimension, show that ψ(x)=L−1/2φ(x/L)\psi(x)=L^{-1/2}\varphi(x/L) preserves normalization.
Solution

Let ξ=x/L\xi=x/L, so dx=L dξdx=L\,d\xi. Then

∫∣ψ(x)∣2 dx=∫1L∣φ(ξ)∣2L dξ=∫∣φ(ξ)∣2 dξ.\int \lvert\psi(x)\rvert^2\,dx = \int \frac{1}{L} \lvert\varphi(\xi)\rvert^2 L\,d\xi = \int \lvert\varphi(\xi)\rvert^2\,d\xi.

Thus the scaled wavefunction has the same dimensionless norm.

  • 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. J. Higham, Accuracy and Stability of Numerical Algorithms, 2nd ed., SIAM, 2002.