The P-31 NMR chemical shifts of adsorbed trimethylphosphine oxide (TMPO) and the configurations of the corresponding TMPOH+ complexes on Bronsted acid sites with varying acid strengths in modeled zeolites have been predicted theoretically by means of density functional theory (DFT) quantum chemical calculations. The configuration of each TMPOH+ complex was optimized at the PW91/DNP level based on an 8T cluster model, whereas the IT chemical shifts were calculated with the gauge including atomic orbital (GIAO) approach at both the HF/TZVP and MP2/TZVP levels. A linear correlation between the P-31 chemical shift of adsorbed TMPO and the proton affinity of the solid acids was observed, and a threshold for superacidity (86 ppm) was determined. This threshold for superacidity was also confirmed by comparative investigations on other superacid systems, such as carborane acid and heteropolyoxometalate H3PW12O40. In conjunction with the strong correlation between the MP2 and the HF P-31 isotropic shifts, the 8T cluster model was extended to more sophisticated models (up to 72T) that are not readily tractable at the GIAO-MP2 level, and a P-31 chemical shift of 86 ppm was determined for TMPO adsorbed on zeolite H-ZSM-5, which is in good agreement with the NMR experimental data.