Photoelectron spectra of sub-valence 2p resonance ionization in sodium atoms have been studied employing the multi-configuration Dirac–Fock (MCDF) method. Special interest has been focused on the satellite lines arising from the 3s→3p conjugate shake-up transitions. It has been found that the agreement between our theoretical results and the measured spectra is good if the rearrangement of the residual bound-state electron density is taken into account. By means of large-scale MCDF calculations, the competing spin–orbit and Coulomb interactions have been evaluated satisfactorily to reproduce the experimental shake-up satellite spectra.