Validate the quantity a band calculation is meant to predict
Separate model error, numerical convergence, filling, and comparison with experiment.
On this page
The question for today
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.
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.
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.
Work it through
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.
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.
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.
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 assumed quadratic sequence extrapolates to 1.000 eV, with a 6.25 meV remaining fine-grid error within that model.
The living notebook
Explore the model
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.
- Midpoint after 0 refinements
- 1.5
- Bracket error bound
- 1.5
- Computed residual · x² − c
- 0.25
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.
Pause & reason
Check your understanding
Your turn
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.
Keep the idea
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.