Surface features of P-doped silica: a comparison between IR spectroscopy and theoretical modelling