Ultra-luminous X-ray sources (ULXs) are X-ray sources located off-galactic centre and have luminosities exceeding the Eddington limit of a stellar-mass black hole ($L_X>10^{39}\;{\rm erg\,s}^{-1}$). With the discovery of X-ray pulsations in some of these objects (e.g. M82 X-2), we know a part of the ULX population has a neutron star (NS) accretor. We present systematic modelling of intermediate-mass X-ray binaries (IMXBs; donor-star mass range $2.0$ to $8.0\;$M$_{\odot}$) to explain the formation of this sub-population of ULXs. Using MESA, we explore the allowed initial parameter space of binary systems consisting of a neutron star and an intermediate-mass donor star, that could explain the observed properties of ULXs. Our simulations take into account beaming effects and also include stellar rotation, tides, general angular momentum losses, and a detailed and self-consistent calculation of the mass-transfer rate. Exploring the initial parameters that lead to the formation of neutron-star ULXs, we study the conditions that lead to dynamical stability of these systems, which depends strongly on the response of the donor star to mass loss. Using two values for the initial NS mass ($1.3\;$M$_{\odot}$ and $2.0\;$M$_{\odot}$) we find that IMXBs can produce NS-ULXs with typical time-averaged isotropic-equivalent X-ray luminosities between $10^{39}$--$10^{41}\;{\rm erg\,s}^{-1}$ on a timescale of up to $\sim\!1.0\;$Myr for the lower luminosities, assuming geometrical beaming. Finally, we estimate their likelihood of detection, the types of white-dwarf remnants left behind by the donors, as well as the total amount of mass accreted by the neutron stars. We also compare our results to the observed pulsating ULXs and infer their initial parameters. Our results suggest that a large subset of the observed pulsating ULX population can be explained by IMXBs in a super-Eddington mass-transfer phase.

Relativistic supernovae constitute a sub-class of type Ic supernovae (SNe). Their non-thermal, radio emission differs notably from that of regular type Ic supernovae as they have a fast expansion speed (with velocities $\sim$ 0.6-0.8 c) which can not be explained by a "standard", spherical SN explosion but advocates for a quickly evolving, mildly relativistic ejecta associated with the SN. In this paper, we compute the synchrotron radiation emitted by the cocoon of a long gamma-ray burst jet (GRB). We show that the energy and velocity of the expanding cocoon, and the radio non-thermal light curves and spectra are consistent with those observed in relativistic SNe. Thus, the radio emission from this events is not coming from the SN shock front, but from the mildly relativistic cocoon produced by the passage of a GRB jet through the progenitor star. We also show that the cocoon radio emission dominates the GRB emission at early times for GRBs seen off-axis, and the flux can be larger at late times compared with on-axis GRBs if the cocoon energy is at least comparable with respect to the GRB energy.

We model the late evolution and mass loss history of rapidly rotating Wolf-Rayet stars in the mass range $5\,\rm{M}_{\odot}\dots 100\,\rm{M}_{\odot}$. We find that quasi-chemically homogeneously evolving single stars computed with enhanced mixing retain very little or no helium and are compatible with Type\,Ic supernovae. The more efficient removal of core angular momentum and the expected smaller compact object mass in our lower mass models lead to core spins in the range suggested for magnetar driven superluminous supernovae. Our more massive models retain larger specific core angular momenta, expected for long-duration gamma-ray bursts in the collapsar scenario. Due to the absence of a significant He envelope, the rapidly increasing neutrino emission after core helium exhaustion leads to an accelerated contraction of the whole star, inducing a strong spin-up, and centrifugally driven mass loss at rates of up to $10^{-2}\,\rm{M}_{\odot}~\rm{yr^{-1}}$ in the last years to decades before core collapse. Since the angular momentum transport in our lower mass models enhances the envelope spin-up, they show the largest relative amounts of centrifugally enforced mass loss, i.e., up to 25\% of the expected ejecta mass. Our most massive models evolve into the pulsational pair-instability regime. We would thus expect signatures of interaction with a C/O-rich circumstellar medium for Type~Ic superluminous supernovae with ejecta masses below $\sim 10\,\rm{M}_{\odot}$ and for the most massive engine-driven explosions with ejecta masses above $\sim 30\,\rm{M}_{\odot}$. Signs of such interaction should be observable at early epochs of the supernova explosion, and may be related to bumps observed in the light curves of superluminous supernovae, or to the massive circumstellar CO-shell proposed for Type~Ic superluminous supernova Gaia16apd.