Gauge-Including Atomic Orbitals (GIAO)¶
The Core Idea¶
Rather than using one gauge origin for the whole molecule, GIAO gives every basis function its own, local gauge origin, centered on the nucleus that basis function belongs to. Each atomic orbital χ is modified with a field-dependent phase factor:
where R_N is the position of the nucleus the basis function is centered on. This phase factor exactly cancels the gauge-origin dependence that would otherwise appear at the level of an individual basis function.
Why This Fixes the Problem¶
The gauge-dependence artifact discussed in The Gauge Origin Problem comes from a finite basis being unable to fully represent how orbitals should respond to a magnetic field, given a fixed gauge origin. By baking the field-dependence directly into the definition of each basis function, GIAO removes the mismatch at its source: gauge invariance holds essentially exactly, even for a modest, finite basis set — which is exactly why a comparatively small basis like 6-31G still gives internally consistent shielding values in the worked example on this site.
The Practical Cost¶
Because each basis function now explicitly depends on B, the derivatives needed to compute the shielding tensor are more involved than an ordinary energy derivative — this added complexity is what the Coupled-Perturbed HF/DFT equations actually solve. This is also the reason the keyword is specifically NMR=GIAO, not just NMR: it names the particular gauge-fixing scheme being used, since other schemes (e.g. IGAIM, CSGT) exist and solve the same underlying problem differently.