Binding energy produced within the framework of the accretion of millisecond pulsars

  • The role and implication of binding energy through the accretion-induced collapse (AIC) of accreting white dwarfs (WDs) for the production of millisecond pulsars (MSPs) are investigated. The binding energy model is examined due to the dynamic process in closed binary systems, and the possible mass of the companion sufficient to induce their orbital parameters is investigated. The deterministic nature of this interaction has a strong sensitivity to the equation of state of the binary systems (where the compactness of a neutron star is proportional to the amount of binding energy) associated with their initial conditions. This behavior mimics the commonly assumed mass and amount of accreted matter under the instantaneous mass loss ($\Delta M \sim 0.18M_{\odot}$). As a result, this indicates an increase in the MSP's gravitational mass due to angular momentum losses. The outcome of such a system is then a circular binary MSP in which the companion is a low-mass WD, thus distinguishing the binary formation scenarios. In addition, the results of this work could provide constraints on the expected mass and binding energy of a neutron star based on the accretion rate.
  • 加载中
  • [1] R. N. Manchester, G. B. Hobbs, A. Teoh et al., Astron. J. 129, 1993 (2005) doi: 10.1086/428488
    [2] P. C. C. Freire and T. M. Tauris, Mon. Not. R. Astron. Soc. 438, L86 (2014) doi: 10.1093/mnrasl/slt164
    [3] Y. Cai, T. Ali, Y.-h. Zhao et al., Chin. Phys. Lett. 36, 137 (2011) doi: 10.1088/0256-307X/28/8/089702
    [4] M. A. Alpar, A. F. Cheng, M. A. Ruderman et al., Nature 300, 728 (1982) doi: 10.1038/300728a0
    [5] D. Bhattacharya and E. P. J. van den Heuvel, Phys. Rep. 203, 1 (1991) doi: 10.1016/0370-1573(91)90064-S
    [6] A. Taani, C. M. Zhang, M. Al-Wardat et al., Astronomische Nachrichten 333, 53 (2012) doi: 10.1002/asna.201111611
    [7] A. Taani and A. Khasawneh, J. Phys.: Conf. Ser. 869, 012090 (2017) doi: 10.1088/1742-6596/869/1/012090
    [8] Z.B. Dai, P. Szkody, A. Taani et al., A&A 604, 45 (2017)
    [9] T. M. Tauris and E. P. J. van den Heuvel, Compact stellar X-ray sources 39, 623 (2006) doi: 10.48550/arXiv.astro-ph/0303456
    [10] A. Taani, A. Abushattal, and M. K. Mardini, Astronomische Nachrichten 340, 847 (2019) doi: 10.1002/asna.201913713
    [11] A. Taani and J. C. Vallejo, Astronomical Society of Australia ( 34, e024 (2017), arXiv:1704.06412 doi: 10.1017/pasa.2017.17
    [12] A. M. Archibald, I. H. Stairs, and S. M. Ransom, Science 324, 1411 (2009) doi: 10.1126/science.1172740
    [13] K. Nomoto, Astrophys. J. 322, 206 (1987) doi: 10.1086/165716
    [14] J. R. Hurley, C. A. Tout, D. T. Wickramasinghe et al., Mon. Not. R. Astron. Soc. 402, 1437 (2010) doi: 10.1111/j.1365-2966.2009.15988.x
    [15] T. M. Tauris, D. Sanyal, S. C. Yoon et al., Astron. Astrophys. 558, A39 (2013) doi: 10.1051/0004-6361/201321662
    [16] T. M. Tauris, M. Kramer, P. C. C. Freire et al., Astrophys. J. 846, 170 (2017) doi: 10.3847/1538-4357/aa7e89
    [17] J. Schwab and K. A. Rocha, Astrophys. J. 872, 131 (2019) doi: 10.3847/1538-4357/aaffdc
    [18] B. Wang and D. Liu, Research in Astronomy and Astrophysics 20, 135 (2020) doi: 10.1088/1674-4527/20/9/135
    [19] A. Kitsikis, G. Fontaine, and P. Brassard, in 14th European Workshop on White Dwarfs, edited by D. Koester and S. Moehler, vol. 334 (Astronomical Society of the Pacic Conference Series, 2005), p. 65
    [20] B. D. Metzger, A. L. Piro, and E. Quataert, Mon. Not. R. Astron. Soc. 396, 1659 (2009) doi: 10.1111/j.1365-2966.2009.14909.x
    [21] S. Darbha, B. D. Metzger, E. Quataert et al., Mon. Not. R. Astron. Soc. 409, 846 (2010), arXiv:1005.1081 doi: 10.1111/j.1365-2966.2010.17353.x
    [22] M. K. Mardini, V. M. P. Haining Li, and, et al., Astrophys. J. 875, 89 (2019) doi: 10.3847/1538-4357/ab0fa2
    [23] M. K. Mardini, A. T. Vinicius M. Placco, and, et al., Astrophys. J. 882, 27 (2019) doi: 10.3847/1538-4357/ab3047
    [24] M. K. Mardini, N. Ershiadat, M. A. Al-Wardat et al., J. Phys.: Conf. Ser. 1258, 012024 (2019) doi: 10.1088/1742-6596/1258/1/012024
    [25] M. K. Mardini, V. M. Placco, Y. Meiron et al., Astrophys. J. 903, 88 (2020) doi: 10.3847/1538-4357/abbc13
    [26] B. Wang, D. Liu, and H. Chen, Mon. Not. R. Astron. Soc. 510, 6011 (2022) doi: 10.1093/mnras/stac114
    [27] W. G. Newton and B.-A. Li, Phys. Rev. C 80, 065809 (2009) doi: 10.1103/PhysRevC.80.065809
    [28] M. Bagchi, Mon. Not. R. Astron. Soc. 413, L47 (2011) doi: 10.1111/j.1745-3933.2011.01030.x
    [29] R. Iaconi, O. De Marco, J.-C. Passy et al., Mon. Not. R. Astron. Soc. 477, 2349 (2018) doi: 10.1093/mnras/sty794
    [30] A. Taani, S. Karino, L. Song et al., Research in Astronomy and Astrophysics 19, 012 (2019), arXiv:1809.01213 doi: 10.1088/1674-4527/19/1/12
    [31] A. Taani, S. Karino, L. Song et al., J. Phys.: Conf. Ser. 1258, 012029 (2019) doi: 10.1088/1742-6596/1258/1/012029
    [32] A. M. Holgado, arXiv e-prints arXiv:2103.13605 (2021).
    [33] A. M. Holgado, H. O. Silva, P. M. Ricker et al., Astrophys. J. Lett. 910, L22 (2021) doi: 10.3847/2041-8213/abecdd
    [34] E. Alecian and S. M. Morsink, Astrophys. J. 614, 914 (2004) doi: 10.1086/423922
    [35] L. Jiang, C.M. Zhang, A. Taani, et al., Inter. Journal of Modern Physics Conf. Series 23, 95 (2013) doi: 10.1142/S2010194513011124
    [36] J. R. Hurley, C. A. Tout, and O. R. Pols, Mon. Not. R. Astron. Soc. 329, 897 (2002) doi: 10.1046/j.1365-8711.2002.05038.x
    [37] A. Taani, J. C. Vallejo, and M, Abu-Saleem, Journal of High Energy Astrophysics 35, 83 (2022), arXiv:2206.11015 doi: 10.1016/j.jheap.2022.06.002
    [38] L. Rezzolla, K. o. Uryu, and S. Yoshida, Mon. Not. R. Astron. Soc. 327, 888 (2001), arXiv:gr-qc/0107019 doi: 10.1046/j.1365-8711.2001.04796.x
    [39] Y.-C. Wei, A. Taani, Y.-Y. Pan et al., Chin. Phys. Lett. 27, 119801 (2010) doi: 10.1088/0256-307X/27/11/119801
    [40] J. G. Hills, Astron. J. 88, 1857 (1983) doi: 10.1086/113478
    [41] W.-C. Chen, Mon. Not. R. Astron. Soc. 501, 2327 (2020) doi: 10.1093/mnras/staa3701
    [42] A. Taani, A. Abushattal, A. Khasawneh, et al., Jordan J. Phys. 13, 243 (2020), arXiv:2002.03011 doi: 10.47011/13.3.8
    [43] A. G. Riess, A. V. Filippenko, P. Challis et al., Astron. J. 116, 1009 (1998) doi: 10.1086/300499
    [44] B. Paczynski, in Structure and Evolution of Close Binary Systems, vol. 73, edited by P. Eggleton, S. Mitton, and J. Whelan (In IAU Symp., 1976) p. 75.
    [45] J. D. M. Dewi and T. M. Tauris, A&A 360, 1043 (2000), arXiv:astro-ph/0007034
    [46] G. P. Briggs, L. Ferrario, C. A. Tout et al., Mon. Not. R. Astron. Soc. 447, 1713 (2015) doi: 10.1093/mnras/stu2539
    [47] M. K. Mardini, A. Frebel, R. Ezzeddine et al., Mon. Not. R. Astron. Soc. 517, 3993 (2022) doi: 10.1093/mnras/stac2783
    [48] P. Podsiadlowski, N. Langer, A. J. T. Poelarends et al., Astrophys. J. 612, 1044 (2004) doi: 10.1086/421713
    [49] B. Reed and C. J. Horowitz, Phys. Rev. D 102, 103011 (2020) doi: 10.1103/PhysRevD.102.103011
    [50] H. Nagakura, A. Burrows, and D. Vartanyan, Mon. Not. R. Astron. Soc. 506, 1462 (2021) doi: 10.1093/mnras/stab1785
    [51] F. Ozel and P. Freire, Annual Rev. Astron. Astrophys. 54, 401 (2016) doi: 10.1146/annurev-astro-081915-023322
    [52] A. Taani, S. Karino, L. Song et al., Publications of the Astronomical Society of Australia 39, e040 (2022) doi: 10.1017/pasa.2022.32
    [53] J. G. Martinez, K. Stovall, P. C. C. Freire et al., Astrophys. J. 812, 143 (2015) doi: 10.1088/0004-637X/812/2/143
  • 加载中

Figures(6)

Get Citation
Ali Taani. The binding energy produced within the framework of the accretion of millisecond pulsars[J]. Chinese Physics C. doi: 10.1088/1674-1137/acb346
Ali Taani. The binding energy produced within the framework of the accretion of millisecond pulsars[J]. Chinese Physics C.  doi: 10.1088/1674-1137/acb346 shu
Milestone
Received: 2022-11-10
Article Metric

Article Views(1074)
PDF Downloads(30)
Cited by(0)
Policy on re-use
To reuse of subscription content published by CPC, the users need to request permission from CPC, unless the content was published under an Open Access license which automatically permits that type of reuse.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Binding energy produced within the framework of the accretion of millisecond pulsars

    Corresponding author: Ali Taani, ali.taani@bau.edu.jo
  • Physics Department, Faculty of Science, Al Balqa Applied University, 19117 Salt, Jordan

Abstract: The role and implication of binding energy through the accretion-induced collapse (AIC) of accreting white dwarfs (WDs) for the production of millisecond pulsars (MSPs) are investigated. The binding energy model is examined due to the dynamic process in closed binary systems, and the possible mass of the companion sufficient to induce their orbital parameters is investigated. The deterministic nature of this interaction has a strong sensitivity to the equation of state of the binary systems (where the compactness of a neutron star is proportional to the amount of binding energy) associated with their initial conditions. This behavior mimics the commonly assumed mass and amount of accreted matter under the instantaneous mass loss ($\Delta M \sim 0.18M_{\odot}$). As a result, this indicates an increase in the MSP's gravitational mass due to angular momentum losses. The outcome of such a system is then a circular binary MSP in which the companion is a low-mass WD, thus distinguishing the binary formation scenarios. In addition, the results of this work could provide constraints on the expected mass and binding energy of a neutron star based on the accretion rate.

    HTML

    I.   INTRODUCTION
    • Millisecond pulsars (MSPs) are remarkable objects born via type II supernova (SN) explosions and characterized by short spin periods (P$_{\rm spin} \leq$ 20 ms), weak magnetic fields, (B $ \leq 10^{9} $ G) and extremely old age ($ \sim 10^{9} $ yr) based on the recycling process. The ATNF census comprises 445 Galactic MSPs in both the Galactic plane and in globular clusters, see, e.g., [13]. They are frequently found in binaries (about 65%) with white dwarfs (WDs) in circular orbits with companions with masses ranging from $\sim 0.15 M_{\odot} - 0.45M_{\odot}$ , see e.g., [48].

      It has been suggested that MSPs are formed in low-mass X-ray binary systems (LMXBs). This argument called for a recycling process in which a slowly rotating, old Neutron Star (NS) may be spun-up into an MSP via accretion from a binary companion, see, e.g., [4, 5, 911]. This link has been more recently reinforced by the discovery of many pulsars like PSR J1023+0038 [12], whose radio emission is likely to have switched on only recently after the LMXB phase.

      An alternative formation scenario is the accretion induced collapse (AIC) of a ONeMg white dwarf (WD), see, e.g., [2, 7, 1318]. In this scenario, the mass transfer can increase the mass of a WD toward the Chandrasekhar mass. When the star reaches this mass limit, degeneracy pressure can no longer support it because the supporting electron pressure is robbed by inverse β decays [19], consequently losing enough energy in neutrinos. This causes it to collapse and violently release a substantial amount of gravitational energy, which might be observable by gravitational wave observatories such as the LIGO-Virgo Collaboration and GEO [2026].

      The NS equation of state, along with the mass-radius relation, plays a role in modeling the baryonic and gravitational masses as well as binding energy [27]. The majority of efforts have been made to calculate the effect of the binding energy by estimating the amount of evolution driven by mass accretion in order to understand the evolutionary tracks of binary models [2834].

      Reference [32] constrained the baryonic mass and binding energy associated with the equation of state of two double pulsars (J0737-3039B and J1756-2251), finding that they are thought to have formed from an ultra-stripped progenitor. It was also found that these systems are consistent with forming from the AIC.

      This work aims to consider the pre-AIC and post-AIC binding energies of the binary system. Then, taking into account the non-zero velocity of ejected matter, I derive several elementary analytic formulae for evaluating the variation of orbital energy and the possibility of binary survival.

    II.   BINDING ENERGY
    • The impact of kick velocity on the companion depends on the binding energy of the system and the kinetic energy of the kick exerted on the orbit. It is noteworthy to mention here that, after AIC, the mass transfer starts once more due to the evolution of the companion star. Figure 1 shows the evolutionary track of the AIC process as it decreases in mass. However, given a particular initial mass of the WD $ 1.1M_{\odot} $ [35], it takes a time interval shorter than 0.4 Gyr, see, e.g., [36] to trigger the AIC. The mass transfer model used in this work is based on the Roche model and ignores the secondary spin angular momentum. The angular momentum loss due to the gravitational radiation is calculated. This is valid for the WD primary [37]. This further changes the orbital period and the eccentricity [38]. This section covers the binding energy of a dynamical system based on the change in orbital energy. This could occur during the pre-SN orbit and immediately before the explosion, as described by the equation

      Figure 1.  The AIC of binary system evolution is tracked, and the mass accretion-rate is expected to be around ($\Delta M \sim 0.2M_{\odot}$) with a timescale shorter than 0.4 Gyr [see e.g., 36, 8].

      $ E_{{\rm orb},i} = \frac{- GM_{\rm WD} M_{c,i}}{2a_{i}} , $

      (1)

      $ \quad\quad = - \frac{GM_{\rm WD} M_{c,i}}{a_{i}} + \frac{M_{\rm WD}M_{c,i}}{M_{\rm NS}+M_{c,i}} v^{2}_{i} , $

      (2)

      where the $ v_{i} $ is WD's initial velocity relative to its companion (mass = Chandrasekhar mass $M_{\rm ch}$). M$_{\rm NS}$ is the NS mass. G is Newton’s gravitational constant. At the evolution of a binary under instantaneous mass loss, the initial mass of the companion $M_c,\;i$ is assumed to be at least $4M_{\odot}\leq M_{\rm com} \leq6M_{\odot}$, which may be reasonable for AIC binaries where the companion is on the main sequence of helium. $ a_{i} $ is the system's initial semi-major axis, which is approximately (1.5 R$ _{\odot} $). For more information, see, e.g., [39] and references therein.

      The orbital energy after the SNe

      $ E_{{\rm orb},f} = \frac{- GM_{\rm NS} M_{c,f}}{2a_{f}}, $

      (3)

      $ \;\;\quad\quad = - \frac{GM_{\rm NS} M_{c,f}}{a_{f}} + \frac{M_{\rm NS}M_{c,f}}{M_{\rm NS}+M_{c,f}} v^{2}_{f} , $

      (4)

      where the companion $ M_{c,f} $'s final mass is estimated to be at least 1.4M$ _{\odot} $, see, e.g., [9]. The final velocity of WD relative to its companion is given by $ v_{f} $. a$ _{f} $ is the system's final semi-major axis after the AIC process, which is approximately (3R$ _{\odot} $), see i.e, [40, 41]. The geometry of an asymmetric SN in the binary system is illustrated in Fig. 2. This depends on the mass loss and angular momentum variations. The low binding energy of the ejecta and the low explosive energy are, of course, related. The slow rotation is caused by the accretion of gas with small angular momentum during this process.

      Figure 2.  Simple diagram illustrating how the companion star, $ M_c $, receives some of the mass ejected by the exploding WD during the AIC. The top panel corresponds to an instant right before the explosion. The bottom panel corresponds to an instant immediately after it happens. The ejected material has been schematically divided into blobs 1, 2, 3, 4, 5, 6, 7. The sum of the masses of all blobs is asymmetric mass ejection. The motion of each blob is shown with an arrow. Blobs 1, 2 and 7 reach the surface of $ M_c $. Therefore, the masses of blobs 1, 2 and 7 are received and set to be ($\Delta M \sim 0.18M_{\odot}$).

      It is worth noting that different amounts of deposition energy into the common envelope should make a noticeable difference in the two different phases (pre- and post-AIC). This is due to the dissipating mechanism acting on the dynamical orbital evolution. In addition, I also consider the angular momentum that could affect the separation distance for the system based on the amount of accreted mass, see, i.e, [40].

      $ J= \frac{M_{\rm NS} M_{c,f}}{M_{\rm NS}+ M_{c,f}} a_{f} \frac{2\pi}{P} . $

      (5)

      However, this effect could be observed from the gravitational radiation emitted by the binary motion in a circular orbit [34]. In addition, the role of the compactness C = GM/(c$ ^{2} $R), where c is the speed of light and G is the gravitational constant, during accretion is adopted from [32], and points out the need to constrain the equation of state. As a result, a more accreted baryonic mass is transformed into binding energy instead of gravitational mass [33].

      Here, I consider the pre-AIC and post-AIC binding energies of the binary system by assuming the difference in final and initial orbital energies $E_{{\rm orb},f}$ and $E_{{\rm orb},i}$, respectively (see Fig. 3). This could cause the system's orbit to continually circularize and contract due to losing orbital energy.

      Figure 3.  Dynamical evolution of the post-AIC and pre-AIC systems as can be seen from the initial and final binding energies. Note that the typical mass loss taken for this is $ 0.18M_{\odot} $.

      $ \begin{array}{*{20}{l}} \Delta E_{\rm orb}= E_{{\rm orb},f} - E_{{\rm orb},i} .\end{array} $

      (6)

      Following Ref. [42], I employ the efficiency parameter (ζ) during a common envelope evolution [43] to obtain the binding energy, which is dependent on the detailed structure (i.e., radius of the companion at an evolutionary stage) of the dynamical process in binaries [44]. Thus, I have

      $ \zeta = \frac{\Delta E_{\rm orb}}{\Delta E_{\rm binding}}. $

      (7)

      From this, I calculate the $ a_{f} $ of the new binary after the AIC process (see Fig. 4). If the value of $ a_{f} $ is too small, the core overfills its new Roche lobe [45], and the cores then merge. As a result, gravitational waves may be produced during the pre-SN phase [16].

      Figure 4.  Final separation distance as a function of time. An indicator of the presence of a dynamical AIC process in a given system is the strong sensitivity to initial conditions. i.e.. masses and semi-major axes (see system details in the text).

      I also deal with a system that orbits around a common center of mass, see e.g., [36, 46]. Then, $ a_{f} $ grows, but the system remains bounded by $a = \dfrac{r(1+e\; {\rm cos\theta})}{1-{\rm e}^{2}}$. As a result, I can obtain both $\Delta E_{{\rm orb},i}$ and $E_{{\rm orb},f}$.

      The initial binding energy of the system can be given by

      $ E_{{\rm binding},i} \simeq {G\; M_{b}} \left(\frac{1}{2a_{f}} - \frac{1}{2a_{i}} \right), $

      (8)

      where $M_{b} = (M_{\rm WD}+M_{c,i}) = (M_{\rm NS}+M_{c,f}$).

      Hence, the final binding energy can be given as

      $ E_{{\rm binding},f}= E_{{\rm binding},i} + \frac{\Delta E_{\rm orb}}{\zeta}. $

      (9)

      In the picture presented here, I find that the initial binding energy has a significant effect on the final fate of both evolution of the components (see Fig. 5) and the orbital parameters that change during the conservative mass transfer. The majority of the mass is lost in an AIC binary system coming from the WD converting baryonic mass into binding energy during the collapse of the newborn NS, with a mass change of the order of $ 0.18M_{\odot} $ [36]. This can happen due to their larger compactness [32]. If enough orbital energy continues to be lost, this causes a merger of the binary components (like in some CO WDs merging). At this stage of evolution, this phenomenon creates a common envelope phase by reducing the orbital energy by envelope binding energy [45, 47] (see Fig. 6).

      Figure 5.  (color online) Final orbital energy as a function of the final fate of the mass companion. The correlation is roughly linear, as predicted by the H-R diagram.

      Figure 6.  (color online) Final orbital energy as a function of the final binding energy through the AIC process.

      A note should be made concerning the given equation of state, because it will govern the specific binding energy for determining the binary's evolution. In [48] a strong correlation was found between several equations of state and tidal deformability associated with NS binding energies. Furthermore, they discovered that the lower limit of NS mass has a binding energy of 1.52$ \times10^{53} $ ergs. Based on these interesting results. [49] estimated the amount of the gravitational mass of the NS, which is $\sim 1.2 M_{\odot}$. This result is consistent with the NS observations made by [50, 51].

      The result from the mass estimation benefits from the mildly recycled pulsar J0453+1559 [52, 53]. This system consists of a double NS with masses of $ 1.56M_{\odot} $ and $ 1.17M_{\odot} $, an orbital period of 4.07-day, a spin period of 45 ms, and an orbital eccentricity of (e = 0.11). The amount of gravitational binding energy is expected to be around ($\Delta M \sim 0.18M_{\odot}$) as shown in Fig. 6. This can be used to add some additional constraints when given a gravitational mass measurement for double NSs.

    III.   CONCLUSIONS
    • An interesting point should be made about the outcome of the final binding energy within the framework of the AIC process, taking into consideration the non-zero velocity of the ejected matter. However, the equation of state influences how much binding energy a NS has. As a result, the larger the binding energy, the more compact the NS. The binding energy thus dominates the gravitational waves at a given amount of accreted mass associated with the angular momentum loss. This provides more insight into AIC observations in the future. I have also investigated the effect of accretion on the orbital period on an NS's binding energy. With the use of the orbital period along with a donor star mass evolution, I have shown that it is possible to add some constraints on the expected values of both mass and binding energy for NS energy parameters that depend on the accretion rate. However, this depends on the initial binding energy associated with the two systems' initial orbital energies, whether an envelope is present or not, since the final binding energy depends on the initial binding energy. This can help us figure out the accretion-driven evolution modes for binaries. The corresponding increase in orbital energy propels the binary into a longer-period orbit. This causes it to collapse and violently release its gravitational energy, which might be observable by gravitational wave observatories such as the LIGO-Virgo Collaboration.

    ACKNOWLEDGMENTS
    • Part of the content of this manuscript has been presented at the IAU Symposium 366, [doi:10.5281/zenodo.5759007]. Special thanks to Nour AlMusleh for helping the calculations. The author would like to thank the anonymous referees for the careful reading of the manuscript and for all suggestions and comments which allowed us to improve both the quality and the clarity of the paper.

Reference (53)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return