The dynamical charge-density response of bulk bismuth has been studied within time-dependent density functional perturbation theory, explicitly accounting for spin-orbit coupling. The use of the Liouville-Lanczos approach allows us to calculate electron energy loss spectra for excitation energies as large as 100 eV. Effects of 5d semicore electronic states, spin-orbit coupling, exchange and correlation, local fields, and anisotropy are thoroughly investigated. The account of the 5d states in the calculation turns out to be crucial to correctly describe the loss spectra above 10 eV and, in particular, the position and shape of the bulk-plasmon peak at 14.0 eV at vanishing transferred momentum. Our calculations reveal the presence of interband transitions at 16.3 eV, which had never been discussed before. The origin of the peak at 5.8 eV is revisited as due to mixed interband and collective excitations. Finally, our study supplements the lack of experiments at finite transferred momenta.
