Evolution of eccentric high-mass X-ray binaries. The case of GX 301-2

The formation of neutron stars is associated with powerful astrophysical transients such as supernovae. In many cases, asymmetries in the supernova explosions are thought to be responsible for the large observed velocities of neutron stars. We aim to study the complete evolutionary history of one particular eccentric high-mass X-ray binary containing a neutron star, GX 301-2, and characterize the natal kick at the time of neutron star formation. We used the publicly available stellar-evolution code MESA to evolve binaries from their initial stages until the core-collapse scenario. We incorporated a natal kick distribution based on observations to continue the evolution during the X-ray binary phase and search for candidates matching current observations of GX 301-2. We find that the range of initial masses is constrained to be less than around $30$ M$_\odot$ depending on the initial mass ratio, as higher initial masses will most likely end up producing a black hole. In the completely conservative mass-transfer scenario under study, only is an interaction between the stars when the donor is still burning Hydrogen in its core, the so-called Case A of mass transfer, able to produce progenitors for GX 301-2. The natal kick study favours kicks of variable strength, which in turn increases the tilt angle between the orbital angular momentum and the spin of the neutron star. We conclude that only a narrow initial progenitor parameter space is able to produce a binary such as GX 301-2. Additionally, the strength of the natal kick can span a wide range of values, but it can be constrained when considering new data concerning the systemic velocity of the binary. Finally, we derive the fraction of the expected number of binaries such as GX 301-2 in the Galaxy to be $\sim 6 \times 10^{-5}$, implying a really low chance of finding a binary similar to GX 301-2.