We consider the planar Ising model on bounded domains from a conformal invariance point of view. We are interested in the scaling limit of the model at critical temperature. Physics theories, notably Conformal Field Theory, predict the existence of two conformal local fields describing the model in the continuum: the spin and the energy density. We have recently proved the conjectured formulae for the energy field on bounded domains, with an improved precision, using discrete complex analysis techniques, thanks to the introduction of holomorphic spinors, that exhibit a fermionic structure. We relate the correlation functions of the energy to special values of the spinors, and prove convergence of the latter to continuous holomorphic spinors, giving scaling formulae for the correlation functions. Partly based on joint work with Stas Smirnov.