The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. By performing a series of three-dimensional hydrodynamical simulations, we investigate how CCSN explosions affect their binary companion. We find that the amount of removed stellar mass, the resulting impact velocity, and the chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Also, it is foud that the impact effects of CCSN ejecta on the structure of main-sequence (MS) companions, and thus their long term post-explosion evolution, is in general not be dramatic.
Keywords. stars: supernovae: general, stars: kinematics, binaries: close
The interaction of CCSN ejecta with a companion] The interaction of core-collapse supernova ejecta with a stellar companion Zheng-Wei Liu et al.] Zheng-Wei Liu, T. M. Tauris, F. K. Röpke, T. J. Moriya, M. Kruckow, R. J. Stancliffe, R. G. Izzard
The discovery of many low-mass X-ray binaries and millisecond pulsars in tight orbits, i.e. binary neutron stars with orbital periods of less than a few hours, provides evidence for supernova (SN) explosions in close binaries with low-mass companions. The nature of the SN explosion determines whether any given binary system remains bound or is disrupted (Hills 1983). An additional consequence of the SN explosion is that the companion star is affected by the impact of the shell debris ejected from the exploding star (Wheeler et al. 1975). Besides chemical enrichment, such an impact has kinematic effects and may induce significant mass loss and heating of the companion star. Core-collapse supernovae (CCSNe) arise from massive stars. There is growing observational evidence that the fraction of massive stars in close binary systems is large. Sana et al. (2012) found that more than of massive stars are in close binary systems, which supports the idea that binary progenitors contribute significantly to the observed CCSNe.
After a SN explosion occurs in a binary system, the ejected debris is expected to expand freely and eventually impact the companion star. The companion star may be significantly heated and shocked by the SN impact, causing the envelope of the companion star to be partially removed due to the stripping and ablation mechanism (e.g., Wheeler et al. 1975, Marietta et al. 2000, Liu et al. 2012, Liu et al. 2013, Pan et al. 2012, Hirai et al. 2018). In this work, we perform impact simulations using a three-dimensional (3D) smoothed particle hydrodynamics (SPH) method to systematically study, for the first time, the impact of CCSN ejecta on MS companion stars.
2 Results and Conclusions
We use the BEC stellar evolution code to construct the detailed companion structure at the moment of SN explosion. The impact of the SN blast wave on the companion star is followed by means of 3D SPH simulations using the Stellar GADGET code (Pakmor et al. 2012). Figure 1 illustrates the temporal density evolution of the SN ejecta and companion material of our hydrodynamics simulations for a G/K-dwarf companion model. Figure 2 shows the effects of varying the orbital separation parameter , by a factor of about 6, on the total amount of removed companion mass (), the resulting impact velocity (), and total accumulated ejecta mass (), for the and companion star models.
To discuss the effects of an explosion on the companion star in CCSNe of Type Ib/c, we perform populations synthesis calculations for SNe with the binary_c/nucsyn code (Izzard et al. 2004, Izzard et al. 2009). Under an assumption of a Galactic star-formation rate of and the average stellar mass of the Kroupa initial mass function (), the total CCSN rate of the Galaxy predicted from their population synthesis is – (36% SNe Ib/c, 10% SNe IIb, 54% SNe II), consistent with the estimated Galactic CCSN rate of from recent surveys (Li et al. 2011).
In our populations synthesis calculations, most SNe Ib/c have an orbital separation of , which is about a fraction of in our binary population synthesis calculations. Furthermore, with the distributions of in Fig. 3, we can simply estimate , and by applying our power-law relationships stated in Fig. 2. We caution that there are large uncertainties in population synthesis studies, which may influence the results. In this case, these mainly relate to the input physics of common-envelope evolution and the subsequent Case BB roche-lobe overflow from the naked helium star prior to its explosion.
We have investigated the impact of SN ejecta on the companion stars in CCSNe of Type Ib/c using the SPH code Stellar GADGET. Our main results can be summarized as follows (see also Liu et al. 2015):
The dependence of total removed mass (), impact velocity () and the amount of accreted SN ejecta mass () on the pre-SN binary separation () can be fitted with power-law functions. All three quantities are shown to decrease significantly with increasing a, as expected (see Fig. 2).
If our population synthesis is correct, we predict that in most CCSNe less than of the MS companion mass can be removed by the SN impact (i.e. ). In addition, the companion star typically receives an impact velocity, , of a few , and the amount of SN ejecta captured by the companion star after the explosion, , is most often less than .
Because a typical CCSN binary companion is relatively massive and can be located at a large pre-SN distance, we do not expect, in general, that the effects of the SN explosion on the post-impact stellar evolution will be very dramatic.
In the closest pre-SN systems, the MS companion stars are affected more strongly by the SN ejecta impact, leading to , and , depending on the mass of the companion star. In addition, these stars are significantly bloated as a consequence of internal heating by the passing shock wave.
It is possible that the SN-induced high velocity stars (HVSs), or more ordinary, less fast, runaway stars, may be contaminated sufficiently to be identified by their chemical peculiarity as former companion stars to an exploding star if mixing processes are not efficient on a long timescale.
- Hills (1983) Hills, J. G. 1983, ApJ, 267, 322
- Hirai et al. 2018 Hirai, Ryosuke, Podsiadlowski, Ph., & Yamada, S. 2018, ApJ, 864, 119
- Izzard et al. 2004 Izzard, R. G., Tout, C. A., Karakas, A. I., & Pols, O. R. 2004, MNRAS, 350, 407
- Izzard et al. 2009 Izzard, R. G., Glebbeek, E., Stancliffe, R. J., & Pols, O. R. 2009 2009, A&A, 508, 1359
- Li et al. 2011 Li, Weidong, Chornock, R., Leaman, J., Filippenko, A. V., Poznanski, D., Wang, Xiaofeng, Ganeshalingam, M., & Mannucci, F. 2011, MNRAS, 412, 1473
- Liu et al. 2012 Liu, Zheng-Wei, Pakmor, R., Röpke, F. K., Edelmann, P., Wang, B., Kromer, M., Hillebrandt, W., & Han, Z. W. 2012, A&A, 548, A2
- Liu et al. 2013 Liu, Zheng-Wei, Pakmor, R., Seitenzahl, I. R., Hillebrandt, W., Kromer, M., Röpke, F. K., Edelmann, P., Taubenberger, S., Maeda, K., Wang, B.ï¼ & Han, Z. W. 2013, ApJ, 774, 37
- Liu et al. 2015 Liu, Zheng-Wei, Tauris, T. M., Röpke, F. K., Moriya, T. J., Kruckow, M., Stancliffe, R. J., Izzard, R. G. 2015, A&A, 584, A11
- Marietta et al. 2000 Marietta, E., Burrows, A., & Fryxell, B. 2000, ApJS, 128, 615
- Pakmor et al. 2012 Pakmor, R., Edelmann, P., Röpke, F. K., & Hillebrandt, W. 2012, MNRAS, 424, 2222
- Pan et al. 2012 Pan, Kuo-Chuan, Ricker, P. M., & Taam, R. E. 2012, ApJ, 750, 151
- Sana et al. (2012) Sana, H., de Mink, S. E., de Koter, A., Langer, N., Evans, C. J., Gieles, M., Gosset, E., Izzard, R. G., Le Bouquin, J.-B., & Schneider, F. R. N. 2012, Science, 337, 444
- Wheeler et al. 1975 Wheeler, J. C., Lecar, M., & McKee, C. F. 1975, ApJ, 200, 145