Long-term Simulations of Multi-Dimensional Core-collapse Supernovae: Implications for Neutron Star Kicks

Core-collapse supernovae (CCSNe) are the final stage of massive stars, marking the birth of neutron stars (NSs). The aspherical mass ejection drives a natal kick of the forming NS. In this work, we study the properties of the NS kick based on our long-term hydrodynamics CCSN simulations. We perform two-dimensional (2D) simulations for ten progenitors from a 10.8 to 20 $\, M_{\odot}$ star covering a wide range of progenitor's compactness parameter, and two three-dimensional (3D) simulations for an 11.2 $\, M_{\odot}$ star. Our 2D models present a variety of the explosion energies between $\sim 1.3 \times 10^{50}$ erg and $\sim 1.2 \times 10^{51}$ erg and the NS kick velocities between $\sim 100$ km s$^{-1}$ and $\sim 1500$ km s$^{-1}$. For the 2D exploding models, we find that the kick velocities tend to become higher with the progenitor's compactness. This is because the high progenitor compactness results in the high neutrino luminosity from the proto-neutron star (PNS), leading to more energetic explosions. Since the high-compactness progenitors produce massive PNS, we point out that the NS masses and the kick velocities can be correlated, which is moderately supported by observation. Comparing 2D and 3D models of the 11.2 $\, M_{\odot}$ star, the diagnostic explosion energy in 3D is, as previously identified, higher than that in 2D, whereas the 3D model results in smaller asymmetry in the ejecta distribution and smaller kick velocity than in 2D. Our results confirm the importance of self-consistent CCSN modeling covering a long-term postbounce evolution in 3D for a quantitative prediction of the NS kicks.