In 1947, Bogoliubov suggested a heuristic theory to compute the excitation spectrum of weakly interacting Bose gases. Such a theory predicts a linear excitation spectrum and provides expressions for the thermodynamic functions which are believed to be correct in the dilute limit. Thus far, there are only a few cases where the predictions of Bogoliubov can be obtained by means of rigorous mathematical analysis. A major challenge is to control the corrections beyond Bogoliubov theory, namely, to test the validity of Bogoliubov’s predictions in regimes where the approximations made by Bogoliubov are not valid. In these notes, we discuss how this challenge can be addressed in the case of a system of N interacting bosons trapped in a box with volume one in the Gross-Pitaevskii limit, where the scattering length of the potential is of the order 1/N and N tends to infinity. This is a recent result obtained in Boccato et al. [Commun. Math. Phys. (to be published); preprint arXiv:1812.03086 and Acta Math. 222, 219–335 (2019); e-print arXiv:1801.01389], which extends a previous result obtained in Boccato et al. [Commun. Math. Phys. 359, 975 (2018)], removing the assumption of a small interaction potential.