  We present in this contribution the predictions on the multiwavelength
spectral energy distribution of our evolutionary population synthesis models
including single and binary stellar systems. The high energy computations
include the emission associated to X-ray binaries and supernovae remnants, as
well as the mechanical energy released into the interstellar medium, which can
be partially reprocessed into thermal X-rays. With these components we compute
the spectral energy distribution of starburst galaxies from X-ray to radio
ranges, and analyze finally the effects of the high energy emission on the H
and He ionizing continuum.
