Molecular dynamics calculations, based on the diatomics-in-molecules method, have been used to probe the manifestations of spin-orbit (SO) coupling in the experimental absorption bands of atomic sodium isolated in solid xenon. Inclusion of SO coupling of –320 cm−1 in spectral simulations of the 3p 2P 3s 2S transition leads to unequal band spacings which very closely match the asymmetrical bandshape observed for blue single vacancy (SV) site occupancy. This SO value, extracted in a previous MCD study, reveals the dramatic change in the effective SO coupling constant of the Na atom (from the gas phase value of +17 cm−1 ) in solid Xe when it is close to the 12 xenon atoms in the first surrounding sphere. In contrast, the symmetrical three-fold split band of the red tetra vacancy (TV) site in Na/Xe is not affected nearly as much by SO coupling. This reflects a greatly reduced “external heavy atom” effect when the 24 Xe atoms surrounding the Na atom in TV are located at greater distances. The contrasting behavior of sodium in the SV and TV sites suggests a strong dependence of the SO coupling strength on the Na–Xe distance.