advanced27 min studyLesson 510 in Physics

Validate the quantity a band calculation is meant to predict

Separate model error, numerical convergence, filling, and comparison with experiment.

On this page

What you’ll understand

  • Distinguish numerical precision from the accuracy of a physical approximation.
  • Use a controlled convergence sequence with a stated asymptotic assumption.
  • Check direct and indirect gaps with full-zone extrema and consistent occupations.

A detailed band plot can conceal several different approximations at once: a restricted orbital basis, an exchange-correlation functional, neglected spin-orbit coupling, finite wavevector sampling, and a structural model. Validation becomes tractable when the target observable is named first. The settings needed for a total-energy difference may not resolve a tiny avoided crossing, effective mass, or Fermi pocket. A successful calculation therefore reports what has converged and what physical assumptions remain to be tested.

Separate physical choices from numerical controls

Start with the Hamiltonian or energy functional, nuclear geometry, electron count, spin configuration, and symmetry constraints. Changing any of these can change the physical question. Numerical controls then approximate the chosen question: basis size or cutoff, wavevector integration, self-consistency tolerance, and interpolation resolution. Compare convergence sequences while holding the physical choices fixed. If a structure relaxes differently at each setting, its changes can mask the numerical trend in the electronic observable.

Convergence of one observable does not imply convergence of every derivative or small energy difference. A coarse path can reproduce broad bandwidths but miss a tiny pocket between path points. A fitted tight-binding model can match energies while misrepresenting curvature or eigenvector geometry. Independent implementations of the same functional test reproducibility and numerical precision; agreement between them does not by itself establish that the functional predicts the experimental excitation spectrum accurately.

Q(h)=Q∞+Ahp+o(hp),Q∞≈rpQ(h/r)−Q(h)rp−1Q(h)=Q_\infty+Ah^p+o(h^p),\qquad Q_\infty\approx\frac{r^pQ(h/r)-Q(h)}{r^p-1}
Q is a specified observable with consistent units, h a dimensionless numerical resolution parameter, r>1 a refinement factor, and p an assumed asymptotic order. The extrapolation is justified only after the observed sequence supports that regime.↔ Scroll wide equations horizontally.

A band gap requires global extrema

The fundamental independent-band gap is the lowest conduction energy anywhere in the zone minus the highest valence energy anywhere in the zone. A direct gap compares the two at the same wavevector. The smallest direct gap can exceed the fundamental gap when the extrema occur at different momenta. A high-symmetry path is a useful diagnostic but cannot guarantee that it contains either extremum; full-zone searches or a justified symmetry argument are needed.

A calculated Kohn–Sham gap, a quasiparticle gap, and an optical onset are different observables. Geometry, temperature, electron-phonon coupling, excitonic effects, and many-body self-energy can affect comparisons. Report which quantity is actually computed before judging disagreement. For metals, verify total occupations, spin channels, smearing dependence, and Fermi-surface convergence rather than using a smooth density-of-states curve as a substitute for an electron-count audit.

Egind=min⁡kEc(k)−max⁡kEv(k),Egdir=min⁡k[Ec(k)−Ev(k)]E_g^{\rm ind}=\min_{\mathbf k}E_c(\mathbf k)-\max_{\mathbf k}E_v(\mathbf k),\qquad E_g^{\rm dir}=\min_{\mathbf k}[E_c(\mathbf k)-E_v(\mathbf k)]
Ec and Ev are conduction and valence energies in the same reference and units. These definitions presume identifiable occupied and unoccupied band manifolds and a global Brillouin-zone search; a negative indirect value signals overlap in this band picture.↔ Scroll wide equations horizontally.

State what the evidence actually establishes

A useful validation record includes the target observable and tolerance, convergence changes, total filling, relevant symmetry or degeneracy checks, and the physical model used for comparison. Distinguish a theorem such as Bloch organization from a controlled weak-potential approximation or an empirical parameter fit. Disagreement may expose a numerical failure, a bad physical approximation, or comparison with a different observable; another digit of printed eigenvalue precision does not identify which explanation applies.

The root-finding lab demonstrates a genuine bound for a much simpler problem. Reset and set target c to 2 and upper bound U to 2, then Step once; predict bracket error bound 0.5. Reset and set c to 2 and U to 4, then Step once; predict bound 1. Its midpoint bounds belong only to bisection of x²−c on a known bracket. The instrument does not certify density-functional errors, band interpolation, or an assumed convergence exponent.

∑nkσwkfnkσ=νe,∑kwk=1\sum_{n\mathbf k\sigma}w_{\mathbf k}f_{n\mathbf k\sigma}=\nu_e,\qquad \sum_{\mathbf k}w_{\mathbf k}=1
wk are normalized Brillouin-zone quadrature weights, f occupations per explicitly counted spin channel, and νe electrons per primitive cell. This audit assumes one state per band and spin per cell and avoids an additional hidden degeneracy factor.↔ Scroll wide equations horizontally.

Extrapolate only under a declared law

A dimensionless mesh parameter takes h=1, 1/2, 1/4 and gives a target band-energy difference 1.100, 1.025, 1.00625 eV. Assume Q(h)=Q∞+Ah² over this range. Check the order and estimate the continuum value and fine-grid error.

  1. The first change is 1.100−1.025=0.075 eV and the next is 1.025−1.00625=0.01875 eV. Their ratio is four, consistent with halving h for a quadratic leading error.

  2. Using the two finest values, Q∞≈[4×1.00625−1.025]/3=1.000 eV. The result follows from the supplied asymptotic model; a different functional or structure would define another convergence problem.

  3. The fine-grid value exceeds the extrapolated limit by 0.00625 eV, or 6.25 meV. The last inter-grid change is 18.75 meV, which is not itself the fine-grid error. Neither estimate addresses physical model error relative to an experimental excitation.

The result

The assumed quadratic sequence extrapolates to 1.000 eV, with a 6.25 meV remaining fine-grid error within that model.

Explore the model

Interactive lab

Interactive lab · root-finding

Compare a real bracket bound

Reset and set c=2, U=2, then Step once; predict midpoint 1.5 and error bound 0.5. Reset and set c=2, U=4, then Step once; predict midpoint 1 and bound 1.

Compare a real bracket boundBisection of x² − c on [0,U], capped at 50 iterations for floating-point resolution. The error bound is the maximum distance from the displayed midpoint to either bracket endpoint (half the width in exact arithmetic). Compensated square evaluation preserves the bracket when ordinary multiplication rounds a residual to zero. A truly exact endpoint or midpoint collapses the bracket. This is one monotone quadratic, not a general root solver. Current readings: Midpoint after 0 refinements: 1.5; Bracket error bound: 1.5; Computed residual · x² − c: 0.25. All parameters are adjustable with the labeled controls below.Bisection on x squared minus 2: bracket from 0 to 3 after 0 refinements.00.751.52.253−2.360.192.745.297.84x · bracket stays nonnegativef(x) = x² − cn = 0 · [ 0, 3 ]
Quadratic / midpointCurrent bracketPositive root
Midpoint after 0 refinements
1.5
Bracket error bound
1.5
Computed residual · x² − c
0.25
t = 0 s

Bisection of x² − c on [0,U], capped at 50 iterations for floating-point resolution. The error bound is the maximum distance from the displayed midpoint to either bracket endpoint (half the width in exact arithmetic). Compensated square evaluation preserves the bracket when ordinary multiplication rounds a residual to zero. A truly exact endpoint or midpoint collapses the bracket. This is one monotone quadratic, not a general root solver. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forThese are exact bisection bracket properties for x²−c. They do not supply error bars for electronic bands, physical functionals, or arbitrary numerical extrapolation.

Check your understanding

What does agreement between two well-converged codes using the same approximate functional most directly establish?

Try it yourself

Suppose a model’s only relevant global extrema and direct-gap minima occur at Γ or X. At Γ, Ev=0 and Ec=1.8 eV; at X, Ev=−0.1 and Ec=1.2 eV. Find the fundamental gap and the minimum direct gap.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Maximize the valence edge and minimize the conduction edge independently before imposing equal wavevector for the direct comparison.

Compare with the solution

The highest valence energy is 0 at Γ and the lowest conduction energy is 1.2 eV at X, so the fundamental gap is indirect and equals 1.2 eV. Direct differences are 1.8 eV at Γ and 1.3 eV at X, giving minimum direct gap 1.3 eV under the explicit global-extrema assumption. Without that assumption, the two sampled points would not certify either global minimum.

Take it with you

  • Convergence must be demonstrated for the specific observable of interest.
  • Numerical reproducibility and physical accuracy are separate claims.
  • Gaps, fillings, and surface geometry require appropriate whole-zone checks.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.