Skip to content

Isotropic Shielding and Anisotropy

Diagonalizing the Tensor

The raw shielding tensor reported for each atom (see Output File) is a general 3×3 matrix, not necessarily symmetric in the Cartesian frame it's printed in. Diagonalizing it gives three principal values, σ₁₁, σ₂₂, σ₃₃ (Gaussian reports these as Eigenvalues), corresponding to the shielding along the tensor's three principal axes.

For example, for atom 1 (carbon) in the formic acid dimer output:

Eigenvalues:   -64.5819    38.5983    99.7689

Isotropic Shielding

The isotropic shielding is simply the average of the three principal values:

σ_iso = (σ₁₁ + σ₂₂ + σ₃₃) / 3

For atom 1: (−64.5819 + 38.5983 + 99.7689) / 3 = 24.5951 ppm, matching the Isotropic = value Gaussian prints directly. This is the physically relevant quantity for a molecule tumbling freely in solution, since rapid molecular tumbling averages out the orientation-dependent part of the tensor, leaving only this isotropic component to be observed.

Anisotropy

The anisotropy describes how much the shielding varies with orientation, how far the tensor is from being perfectly isotropic. Gaussian computes it as:

Δσ = σ₃₃ − (σ₁₁ + σ₂₂) / 2

using the convention that eigenvalues are ordered so σ₃₃ is the one furthest from the other two. For atom 1: 99.7689 − (−64.5819 + 38.5983)/2 = 99.7689 − (−12.9918) = 112.7607 ppm, matching the reported value.

A large anisotropy (as seen for the bridging oxygen, atom 5, at 523.07 ppm) means the local electronic environment responds very differently depending on the molecule's orientation relative to the field, relevant for solid-state NMR, where molecules are fixed in orientation rather than tumbling freely.

Ask a Question

Have a question about quantum chemistry, or anything discussed on this site? Send it our way and we will reply by email.

Prefer to fill it out directly? Open the form.