3D MHD simulations of runaway pulsars in core collapse supernova remnants

Pulsars are one of the possible final stages in the evolution of massive stars. If a supernova explosion is anisotropic, it can give the pulsar a powerful kick, propelling it to supersonic speeds. The resulting pulsar wind nebula is significantly reshaped by its interaction with the surrounding medium as the pulsar moves through it. First, the pulsar crosses the supernova remnant, followed by the different layers of circumstellar medium formed during different stages of the progenitor star s evolution. We aim to investigate how the evolutionary history of massive stars shapes the bow shock nebulae of runaway kicked pulsars, and how these influences in turn affect the dynamics and non-thermal radio emission of the entire pulsar remnant. We perform three-dimensional magnetohydrodynamic simulations using the PLUTO code to model the pulsar wind nebula generated by a runaway pulsar in the supernova remnant of a red supergiant progenitor, and derive its non-thermal radio emission. The supernova remnant and the pre-supernova circumstellar medium of the progenitor strongly confine and reshape the pulsar wind nebula of the runaway pulsar, bending its two side jets inwards and giving the nebula an arched shape for an observer perpendicular to the jets and the propagation direction, as observed around PSR J1509 5850 and Gemina. We perform the first classical 3D model of a pulsar moving inward through its supernova ejecta and circumstellar medium, inducing a bending of its polar jet that turns into characteristic radio synchrotron signature. The circumstellar medium of young runaway pulsars has a significant influence on the morphology and emission of pulsar wind nebulae, whose comprehension requires a detailed understanding of the evolutionary history of the progenitor star.