Properties of proto neutron star PSR J0737-3039A

  • Using five sets of nucleon coupling constants (DD-MEI, GL85, GL97, GM1, and NL2), we find that the radius of the PNS PSR J0737-3039A is R=15.693–18.846 km, the central baryon density is ρc=0.247–0.359 fm3, the central energy density is εc=4.30×1014– 6.49×1014 gcm3, and the central pressure is pc=3.79×1034 – 5.85×1034 dynecm2. With DD-MEI, GL85, GL97, and GM1, baryon octets appear in the PNS PSR J0737-3039A. With NL2, only the baryons n, p, Λ, Σ0, Σ, Ξ0, and Ξ are present. Corresponding to the same baryon density, the relative densities of the same baryon in the PNS PSR J0737-3039A calculated using different nucleon coupling constants differ greatly. The central relative baryon densities of the PNS PSR J0737-3039A calculated using different nucleon coupling constants also differ greatly.
      PCAS:
    • 21.65.Mn(Equations of state of nuclear matter)
  • The 2.8-second pulsar J0737-3039B is the companion of the 23-millisecond pulsar J0737-3039A. The two pulsars form a highly relativistic binary neutron star (NS) system [1], and the mass of the NS J0737-3039A should be related to the mass and metallicity of its companion [2]. The NS J0737-3039A has a mass M=1.34 M, which is typical for NSs [1, 3, 4].

    Morrison et al. constructed a numerical model of the binary pulsar J0737-3039A using its measured mass, spin, and moment of inertia to determine the optimal model for different equations of state (EoS) so as to determine the radius of the pulsar. The NS they used, PSR J0737-3039A, has a gravitational mass M=1.337 M [5]. In 2006, Kramer et al. more accurately determined the mass of the NS PSR J0737-3039A to be M = 1.3381±0.0007 M [6].

    In 2008, Worley et al. placed special emphasis on the PSR J0737-3039A component of the extremely relativistic double NS system PSR J0737-3039. They found its moment of inertia to be between 1.30×1045 gcm2 and 1.63×1045 gcm2 [7]. Using well-known universal relations between NS observables, Landry et al. translated the reported 90%-credible bounds on tidal deformability into a direct constraint on the moment of inertia of PSR J0737-3039A [8].

    In 2020, Avdeev et al. developed a parameterized post-Keplerian formalism for hybrid metric-Palatini f(R)-gravity and obtained the mass ratio MA/MB=1.0714. On this basis, they predicted the mass of the NS PSR J0737-3039A as 1.3374 M<M<1.3440 M in the framework of hybrid f(R)-gravity [9].

    An NS has a high density, strong magnetic field, and high rotation speed [1016]. Several massive NSs, PSR J1614-2230 [17, 18], PSR J0348+0432 [19], and PSR J0740+6620 [20], have been observed in recent years. Although the NS PSR J0737-3039A has a significantly smaller mass than these, its mass is typical, and it is a member of a binary NS system. Therefore, theoretical research on this NS is important.

    NSs are formed by the neutrino radiation cooling of proto neutron stars (PNSs) [21]. Studying PNSs is vital for understanding the evolution of stars. Therefore, we are very interested in the PNS PSR J0737-3039A.

    In this study, we investigate the properties of the PNS PSR J0737-3039A with relativistic mean field (RMF) theory [22] considering a baryon octet.

    The mesons f0(975) (denoted as σ) and ϕ(1020)(denoted as ϕ) represent the interactions between nucleons and hyperons [23]. The Lagrangian density of hadronic matter containing the mesons σ and ϕ is expressed as follows [24]:

    L=B¯ΨB(iγμμmB+gσBσ+gσBσgωBγ0ωgϕBγ0ϕgρBγ0τ3ρ)ΨB12m2σσ213g2σ314g3σ4+12m2ωω2+12m2ρρ212m2σσ2+12m2ϕϕ2+λ=e,μ¯Ψλ(iγμμmλ)Ψλ.

    (1)

    Considering neutrino binding, the baryon partition function of PNS matter is [25, 26]

    lnZB=VTL+B2JB+12π20k2dk{ln[1+e(εB(k)μB)/T]}.

    (2)

    From the above, the total baryon number density is obtained as follows:

    ρ=B2JB+12π2bB0k2nB(k)dk.

    (3)

    The energy density and pressure are respectively

    ε=12m2σσ2+12m2σσ2+13g2σ3+14g3σ4+12m2ωω20+12m2ϕϕ2+12m2ρρ203+B2JB+12π20k2nB(k)dkk2+m2B,

    (4)

    p=12m2σσ212m2σσ213g2σ314g3σ4+12m2ωω20+12m2ϕϕ2+12m2ρρ203+13B2JB+12π20k4k2+m2BnB(k)dk.

    (5)

    Here, nB(k) is the Fermi-Dirac partition function of baryons,

    nB(k)=11+exp[(εB(k)μB)/T],

    (6)

    and mB is the effective mass of baryons,

    mB=mBgσBσgσBσ.

    (7)

    Ignoring the interactions of leptons at finite temperature, we express their partition function as

    lnZL=VTiμ4i24π2[1+2(πTμi)2+715(πTμi)4]+Vλ1π20k2dk{ln[1+e(ελ(k)μλ)/T]},

    (8)

    where the first line represents the contribution of massless neutrinos, and the second line represents the contribution of electrons and μs.

    The lepton number density is

    ρl=1π20k2nl(k)dk,

    (9)

    ρν=π2T2μν+μ3ν6π2

    (10)

    and the contributions of leptons to the energy density and pressure are respectively

    ε=l1π20k2nl(k)dkk2+m2l+ν(7π2T4120+T2μ2ν4+μ4ν8π2),

    (11)

    p=13l1π20k4k2+m2lnl(k)dk+ν1360(7π2T4+30T2μ2ν+15μ4νπ2).

    (12)

    The chemical potentials of baryons are

    μi=μnqi(μeμνe).

    (13)

    To describe neutrino trapping, the content of leptons in a PNS is defined as

    Ylν=Yl+Yνl=ρl+ρνlρ.

    (14)

    We can calculate the mass and radius of a PNS using the Tolman-Oppenheimer-Volkoff (TOV) equations [27, 28]

    dpdr=(p+ε)(M+4πr3p)r(r2M),

    (15)

    M=4πR0εr2dr.

    (16)

    Here, we calculate the parameters of the PNS PSR J0737-3039A with five sets of nucleon coupling constants (DD-MEI [29], GL85 [30], GL97 [24], GM1 [31], and NL2 [32]). For the outer crust with 4.73×1015 <ρ<8.907×103 fm3, we use the EoS of a BPS [24], and for the region of ρ>8.907×103 fm3, we use the EoS for PNS matter corresponding to the nuclear properties above.

    The ratios of the hyperon coupling constant to nucleon coupling constant are defined as xσh=gσh/gσ, xωh=gωh/gω, and xρh=gρh/gρ, with h denoting the hyperons Λ,Σ, and Ξ.

    xρh and xωh are chosen using SU(6) symmetry [33, 34]. Then, xσh is obtained using [24]

    U(N)h=mn(mnmn1)xσh+(gωmω)2ρ0xωh.

    (17)

    In this study, the hyperon-potentials are chosen as U(N)Λ=30 MeV [34, 35, 36], U(N)Σ=30 MeV [34, 35, 36, 37], and U(N)Ξ=14 MeV [38].

    The coupling constants of the mesons σ and ϕ can be chosen as [23]

    gϕΞ=2gϕΛ=2gϕΣ=22gω/3,

    (18)

     gσΛ/gσ=gσΣ/gσ=0.69,

    (19)

     gσΞ/gσ=1.25.

    (20)

    Temperatures in the interiors of PNSs can reach up to 30 MeV. After that, they can be cooled by neutrino radiation [21]. In this study, the temperature of the PNS PSR J0737-3039A is taken as T=20 MeV.

    The mass M of the PNS is calculated as a function of the radius R. In Fig. 1, the coarse curves indicate the masses of the PNSs PSR J0740+6620, PSR J0348+0432, PSR J1614-2230, and PSR J0737-3039A, and the five fine curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2. We find that only the nucleon coupling constant DD-MEI can give the masses of the PNSs PSR J0740+6620, PSR J0348+0432, PSR J1614-2230, and PSR J0737-3039A. NL2 cannot give the mass of the PNS PSR J0740+6620; however, it can give the mass of the PNSs PSR J0348+0432, PSR J1614-2230, and PSR J0737-3039A. GL85, GL97, and GM1 cannot give the masses of the PNSs PSR J0740+6620, PSR J0348+0432, and PSR J1614-2230; however, they can give the masses of the PNS PSR J0737-3039A. All five sets of nucleon coupling constants give the masses of the PNS PSR J0737-3039A.

    As shown in Fig. 1 and Table 1, the radii of the PNS PSR J0737-3039A given by the five sets of nucleon coupling constants are different. NL2 gives the smallest radius, R=15.693 km, and GL85 gives the largest radius, R=18.846 km. In other words, the radius of the PNS PSR J0737–3039A calculated using these five sets of nucleon coupling constants is R=15.693–18.846 km. The NS PSR J0030+0451 has a mass similar to that of the PNS PSR J0737-3039A; hence, it is logical to compare their radii. The mass and radius of the NS PSR J0030+0451 given by Riley et al. are in the ranges M=1.18–1.49 M and R=11.52–13.85 km [39], respectively, and those given by Miller et al. are M=1.30–1.59 M and R=11.96–14.26 km, respectively [40]. The PNS radius calculated in this study is larger than that calculated by Riley et al. and Miller et al., which is reasonable because the radius of a PNS is larger than that of an NS of the same mass [41].

    Figure 1

    Figure 1.  (color online) Mass M of the PNS calculated as a function of the radius R. The coarse curves indicate the masses of the PNSs PSR J0740+6620, PSR J0348+0432, PSR J1614-2230, and PSR J0737-3039A. The five fine curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

    Table 1

    Table 1.  Results calculated in this study. The mass of the PNS PSR J0737-3039A is M=1.338 M. R, ρc, εc, and pc are the radius, central baryon density, central energy density, and central pressure of the PNS PSR J0737-3039A, respectively.
    Parameter R/km ρc/fm3 εc×1015/(gcm)3 pc×1035/(dynecm)2
    DD-MEI 18.006 0.247 0.430 0.379
    GL85 18.846 0.290 0.519 0.427
    GL97 17.353 0.359 0.649 0.585
    GM1 16.798 0.310 0.545 0.506
    NL2 15.693 0.293 0.508 0.506
    DownLoad: CSV
    Show Table

    The central baryon densities ρc of the PNS PSR J0737-3039A calculated with different nucleon coupling constants are different. ρc of the PNS PSR J0737-3039A calculated using GL97 is the largest, ρc=0.359 fm3, whereas ρc calculated using DD-MEI is the smallest, ρc=0.247 fm3. That is, the central baryon density calculated using the five nucleon coupling constants is ρc=0.247–0.359 fm3.

    The effective mass of nucleons, Λ, Σ, and Ξ as a function of the baryon density ρ are shown in Fig. 2. We denote the center effective mass of nucleons, Λ, Σ, and Ξ as mNc, mΛc, mΣc, and mΞc, respectively. As shown in Fig. 2, the effective masses of nucleons, Λ, Σ, and Ξ all decrease with increasing baryon density ρ. At the center of the PNS, the effective mass of the baryon reaches a minimum. Corresponding to the same baryon density ρ, the effective mass calculated using DD-MEI is the smallest, whereas that calculated using GL97 is the largest. As shown in Table 2, the center effective masses of nucleons, Λ, Σ, and Ξ calculated using the above five groups of nucleon coupling constants are mNc=384.521–606.324 MeV, mΛc=778.280–911.405 MeV, mΣc=939.187–1080.822 MeV, and mΞc=1152.637–1218.964 MeV, respectively.

    Figure 2

    Figure 2.  (color online) Effective mass of baryons in the PNS PSR J0737-3039A as a function of baryon density. The five curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

    Table 2

    Table 2.  Effective mass of baryons and the chemical potentials of n and e in the PNS PSR J0737-3039A calculated in this study. mNc, mΛc, mΣc, and mΞc are the center effective masses of nucleons N, Λ, Σ, and Ξ, respectively. μnc and μec are central chemical potentials of n and e, respectively.
    Parameter mNc/MeV mΛc/MeV mΣc/MeV mΞc/MeV μnc/fm1 μec/fm1
    DD-MEI 384.521 778.280 939.187 1152.637 5.340 0.721
    GL85 606.324 911.405 1080.822 1218.964 5.421 0.943
    GL97 573.011 888.351 1071.635 1206.096 5.516 0.989
    GM1 478.656 831.015 1006.070 1179.389 5.418 0.763
    NL2 394.604 780.860 1015.662 1200.968 5.367 0.829
    DownLoad: CSV
    Show Table

    The chemical potentials of neutrons and electrons are shown in Fig. 3. The chemical potentials of neutrons μn and electrons μe increase with increasing baryon density ρ. Different nucleon coupling constants result in different chemical potentials relative to the same baryon density ρ. The central neutron chemical potential of the five nucleon coupling constants is calculated as μnc=5.340–5.516 fm1, and the central electron chemical potential is calculated as μec=0.721–0.989 fm1.

    Figure 3

    Figure 3.  (color online) Chemical potentials of neutrons and electrons in the PNS PSR J0737-3039A as a function of baryon density. The five curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

    The energy density ε and pressure p in the PNS PSR J0737-3039A as a function of the baryon density ρ are shown in Fig. 4. Both the energy density ε and pressure p increase with increasing baryon density ρ. Corresponding to the same baryon density ρ, the energy densities calculated using different nucleon coupling constants have little difference; however, the calculated pressures differ greatly. The central energy density and central pressure calculated for the five sets of nucleon coupling constants are different. GL97 results in the maximum center energy density and center pressure.

    Figure 4

    Figure 4.  (color online) Energy density and pressure in the PNS PSR J0737-3039A as a function of baryon density. The five curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

    As shown in Table 1, the minimum central energy density given with DD-MEI is εc=4.30×1014 gcm3. GL97 gives the largest central energy density, εc=6.49×1014 gcm3. The central energy density calculated from the five nucleon coupling constants is εc=4.30×1014–6.49×1014 gcm3. Similarly, the central pressure calculated using the five nucleon coupling constants is pc=3.79×1034–5.85×1034 dynecm2.

    Figure 5 shows the relative particle density ρi/ρ in the PNS PSR J0737-3039A as a function of the baryon density ρ.

    Figure 5

    Figure 5.  (color online) Relative particle density in the PNS PSR J0737-3039A as a function of baryon density. The five curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

    As shown in Fig. 5 and Table 3, baryons n, p, Λ, Σ+, Σ0, Σ, Ξ0 and Ξ all appear in the PNS PSR J0737-3039A when using the nucleon coupling constants DD-MEI, GL85, GL97, and GM1. However, with NL2, only baryons n, p, Λ, Σ0, Σ, Ξ0, and Ξ are present, while baryon Σ+ is absent.

    Table 3

    Table 3.  Central relative particle densities in the PNS PSR J0737-3039A. ρnc/ρ, ρpc/ρ, ρΛc/ρ, ρΣ+c/ρ, ρΣ0c/ρ, ρΣc/ρ, ρΞ0c/ρ, and ρΞc/ρ are the central relative particle densities of n, p, Λ, Σ+, Σ0, Σ, Ξ0, and Ξ, respectively, in the PNS PSR J0737-3039A.
    Parameter ρnc/ρ% ρpc/ρ% ρΛc/ρ% ρΣ+c/ρ% ρΣ0c/ρ% ρΣc/ρ% ρΞ0c/ρ% ρΞc/ρ%
    DD-MEI 90.7 8.8 0.5 1E-6 2.0E-4 0.04 5.1E-5 0.02
    GL85 81.0 17.5 1.5 2.9E-3 4.0E-4 5.7E-5 1.3E-3 0.05
    GL97 80.5 16.5 2.8 1.3E-3 4.1E-4 1.3E-4 2.4E-3 0.18
    GM1 90.7 8.4 0.8 2E-6 1.8E-4 0.02 1.07E-4 0.05
    NL2 87.5 11.5 1.0 0 1.3E-5 5.2E-4 1.5E-5 5.8E-3
    DownLoad: CSV
    Show Table

    As shown in Fig. 5, corresponding to the same baryon density, the relative particle densities of the same baryon in the PNS PSR J0737-3039A calculated using different nucleon coupling constants differ greatly. As shown in Table 3, the central relative baryon densities of the PNS PSR J0737-3039A calculated using different nucleon coupling constants also differ greatly. For example, for hyperon Λ, the minimum central relative baryon density given using DD-MEI is ρΛc/ρ=0.5%, whereas GL97 gives the largest central baryon density, ρΛc/ρ=2.8%. In other words, the central baryon density of the Λ-hyperon calculated using the five nucleon coupling constants is ρΛc/ρ=0.5%–2.8%. Similarly, the central baryon density of n calculated using five groups of nucleon coupling constants is ρnc/ρ=80.5%–90.7%, the central baryon density of p is ρpc/ρ=8.4%–17.5%, the central baryon density of the Σ0 hyperon is ρΣ0c/ρ=1.3×103%–4.1×104%, the central baryon density of the Σ hyperon is ρΣc/ρ=5.7×105%–0.04%, the central baryon density of the Ξ0 hyperon is ρΞ0c/ρ=1.5×105%–2.4×103%, and the central baryon density of the Ξ hyperon is ρΞc/ρ=5.8×103%–0.18%. The hyperon Σ+ in the PNS PSR J0737-3039A calculated using NL2 does not appear, whereas the central baryon density of Σ+ calculated using DD-MEI, GL85, GL97, and GM1 is ρΣ+c/ρ=1.0×106%–2.9×103%.

    In this study, the influence of nucleon coupling constants on the properties of the PNS PSR J0737-3039A is investigated using RMF theory. Here, five sets of nucleon coupling constants (DD-MEI, GL85, GL97, GM1, and NL2) are used, and the temperature of the PNS PSR J0737-3039A is chosen as T=20 MeV.

    The radius R, central baryon density ρc, central energy density εc, and central pressure pc of the PNS PSR J0737-3039A calculated using different nucleon coupling constants are different. Using these five sets of nucleon coupling constants, we find that the radius of the PNS PSR J0737-3039A is R=15.693–18.846 km, the central baryon density is ρc=0.247–0.359 fm3, the central energy density is εc=4.30×1014–6.49×1014 gcm3, and the central pressure is pc=3.79×1034–5.85×1034 dynecm2. Although the PNS radius should be larger than that of an NS of the same mass, it should not be significantly larger. As shown in Table 1, the radius of the PNS PSR J0737-3039A calculated using the nucleon coupling constant NL2 is the smallest, R=15.693 km, and is closest to the radius calculated by Riley et al. and Miller et al. Therefore, NL2 seems to be the most reasonable coupling constant.

    For four sets of nucleon coupling constants (DD-MEI, GL85, GL97 and GM1), baryons n, p, Λ, Σ+, Σ0, Σ, Ξ0, and Ξ all appear in the PNS PSR J0737-3039A. However, using NL2, only baryons n, p, Λ, Σ0, Σ, Ξ0, and Ξ are present, while baryon Σ+ does not appear.

    Corresponding to the same baryon density, the relative densities of the same baryon in the PNS PSR J0737-3039A calculated using different nucleon coupling constants differ greatly. The central relative baryon densities of the PNS PSR J0737-3039A calculated using different nucleon coupling constants also differ greatly. The central baryon density of n calculated using the five groups of nucleon coupling constants is ρnc/ρ=80.5%–90.7%, the central baryon density of p is ρpc/ρ=8.4%–17.5%, the central baryon density of Λ is ρΛc/ρ=0.5%–2.8%, the central baryon density of the Σ0 hyperon is ρΣ0c/ρ=1.3×103%–4.1×104%, the central baryon density of the Σ hyperon is ρΣc/ρ=5.7×105%–0.04%, the central baryon density of the Ξ0 hyperon is ρΞ0c/ρ=1.5×105%–2.4×103%, and the central baryon density of the Ξ hyperons is ρΞc/ρ=5.8×103%–0.18%. The hyperon Σ+ does not appear with NL2, whereas the central baryon density of hyperon Σ+ calculated using DD-MEI, GL85, GL97, and GM1 is ρΣ+c/ρ=1.0×106%–2.9×103%.

    [1] A. G. Lyne, M. Burgay, and M. Kramer et al., Science 303, 1153 (2004) doi: 10.1126/science.1094645
    [2] Z. L. Deng, Z. F. Gao, X. D. Li et al., ApJ 892, 4 (2020) doi: 10.3847/1538-4357/ab76c4
    [3] S. E. Thorsett and D. Chakrabarty, ApJ 512, 288 (1999) doi: 10.1086/306742
    [4] M. Bailes, S. M. Ord, and H. S. Knight et al., ApJ 595, L49 (2003) doi: 10.1086/378939
    [5] I. A. Morrison, T. W. Baumgarte, and S. L. Shapiro et al., ApJ 617, L135-L138 (2004) doi: 10.1086/427235
    [6] M. Kramer, I. H. Stairs, and R. N. Manchester et al., Science 314, 97-102 (2006) doi: 10.1126/science.1132305
    [7] A. Worley, P. G. Krastev, and B. A. Li, ApJ 685, 390-399 (2008) doi: 10.1086/589823
    [8] P. Landry and B. Kumar, ApJL 868, L22 (2018) doi: 10.3847/2041-8213/aaee76
    [9] N. A. Avdeev, P. I. Dyadina, and S. P. Labazova, J Exp. Theor. Phys. 131, 537-547 (2020) doi: 10.1134/S1063776120100039
    [10] D. H. Wen, W. Chen, Y. G. Lu et al., Mod. Phys. Lett. A 22, 631 (2007) doi: 10.1142/S0217732307023225
    [11] H. Y. Jia, B. X. Sun, J. Meng et al., Chin. Phys. Lett. 18, 1571 (2001) doi: 10.1088/0256-307X/18/12/309
    [12] C. Zhu, Z. F. Gao et al., Mod. Phys. Lett. A 31, 1650070 (2016)
    [13] Y. Lim, A. Bhattacharya, J. W. Holt et al., Phys. Rev. C 104, L032802 (2021) doi: 10.1103/PhysRevC.104.L032802
    [14] P. C. Chu, Y. Zhou, Y. Y. Jiang et al., Eur. Phys. J. C 81, 93 (2021) doi: 10.1140/epjc/s10052-020-08800-3
    [15] Z. F. Gao, H. Shan, and H. Wang, Astron. Nachr. 342, 369-376 (2021) doi: 10.1002/asna.202113936
    [16] H. Wang, Z. F. Gao, H. Y. Jia et al., Universe 6, 63 (2020) doi: 10.3390/universe6050063
    [17] P. B. Demorest, T. Pennucci, S. M. Ransom et al., Nature 467, 1081 (2010) doi: 10.1038/nature09466
    [18] E. Fonseca, T. T. Pennucci, J. A. Ellis et al., ApJ 832, 167 (2016) doi: 10.3847/0004-637X/832/2/167
    [19] J. Antoniadis, P. C. C. Freire, N. Wex et al., Science 340, 448 (2013)
    [20] H. T. Cromartie, E. Fonseca, S.M. Ransom et al., Nat. Astron. 4, 72 (2020) doi: 10.1038/s41550-019-0880-2
    [21] A. Burrows and J. M. Lattier, Astrophy. 307, 178 (1986) doi: 10.1086/164405
    [22] S. G. Zhou, Phys. Scripta 91, 063008 (2016) doi: 10.1088/0031-8949/91/6/063008
    [23] J. Schaffner, C. B. Dover, and A. Gal, Ann. Phys. 235, 35 (1994) doi: 10.1006/aphy.1994.1090
    [24] N. K. Glendenning, Compact Stars: Nuclear Physics, Particle Physics, and General Relativity, Springer-Verlag, New York, Inc, 1997
    [25] N. K. Glendenning, Phys. Lett. B 185, 275 (1987) doi: 10.1016/0370-2693(87)90999-3
    [26] N. K. Glendenning, Nucl. Phys. A 469, 600 (1987) doi: 10.1016/0375-9474(87)90015-7
    [27] R. C. Tolman, Phys. Rev. 55, 364 (1939) doi: 10.1103/PhysRev.55.364
    [28] J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, 374 (1939) doi: 10.1103/PhysRev.55.374
    [29] S. Typel and H. H. Wolter, Nucl. Phys. A 656, 331 (1999) doi: 10.1016/S0375-9474(99)00310-3
    [30] N. K. Glendenning, ApJ 293, 470 (1985) doi: 10.1086/163253
    [31] N. K. Glendenning and S. A. Moszkowski, Phys. Rev. Lett. 67, 2414 (1991) doi: 10.1103/PhysRevLett.67.2414
    [32] S. J. Lee et al., Phys. Rev. Lett. 57, 2916 (1986) doi: 10.1103/PhysRevLett.57.2916
    [33] J. Schaffner and I. N. Mishustin, Phys. Rev. C 53, 1416 (1996) doi: 10.1103/PhysRevC.53.1416
    [34] J. Schaffner-Bielich and A. Gal, Phys. Rev. C 62, 034311 (2000) doi: 10.1103/PhysRevC.62.034311
    [35] S. Weissenborn, D. Chatterjee, and J. Schaffner-Bielich, Nucl. Phys. A 881, 62 (2012) doi: 10.1016/j.nuclphysa.2012.02.012
    [36] A. Gal, E. V. Hungerford, and D. J. Millener, Rev. Mod. Phys. 88, 035004 (2016) doi: 10.1103/RevModPhys.88.035004
    [37] C. J. Batty, E. Friedman, and A. Gal, Phys. Rep. 287, 385 (1997) doi: 10.1016/S0370-1573(97)00011-2
    [38] T. Harada, Y. Hirabayashi, and A. Umeya, Phys. Lett. B 690, 363 (2010) doi: 10.1016/j.physletb.2010.05.053
    [39] T. E. Riley, A. L. Watts, S. Bogdanov et al., ApJL 887, L21 (2019) doi: 10.3847/2041-8213/ab481c
    [40] M. C. Miller, F. K. Lamb, A. J. Dittmann et al., ApJL 887, L24 (2019) doi: 10.3847/2041-8213/ab50c5
    [41] X. F. Zhao, Astrophys. Space Sci. 362, 95 (2017) doi: 10.1007/s10509-017-3073-4
  • [1] A. G. Lyne, M. Burgay, and M. Kramer et al., Science 303, 1153 (2004) doi: 10.1126/science.1094645
    [2] Z. L. Deng, Z. F. Gao, X. D. Li et al., ApJ 892, 4 (2020) doi: 10.3847/1538-4357/ab76c4
    [3] S. E. Thorsett and D. Chakrabarty, ApJ 512, 288 (1999) doi: 10.1086/306742
    [4] M. Bailes, S. M. Ord, and H. S. Knight et al., ApJ 595, L49 (2003) doi: 10.1086/378939
    [5] I. A. Morrison, T. W. Baumgarte, and S. L. Shapiro et al., ApJ 617, L135-L138 (2004) doi: 10.1086/427235
    [6] M. Kramer, I. H. Stairs, and R. N. Manchester et al., Science 314, 97-102 (2006) doi: 10.1126/science.1132305
    [7] A. Worley, P. G. Krastev, and B. A. Li, ApJ 685, 390-399 (2008) doi: 10.1086/589823
    [8] P. Landry and B. Kumar, ApJL 868, L22 (2018) doi: 10.3847/2041-8213/aaee76
    [9] N. A. Avdeev, P. I. Dyadina, and S. P. Labazova, J Exp. Theor. Phys. 131, 537-547 (2020) doi: 10.1134/S1063776120100039
    [10] D. H. Wen, W. Chen, Y. G. Lu et al., Mod. Phys. Lett. A 22, 631 (2007) doi: 10.1142/S0217732307023225
    [11] H. Y. Jia, B. X. Sun, J. Meng et al., Chin. Phys. Lett. 18, 1571 (2001) doi: 10.1088/0256-307X/18/12/309
    [12] C. Zhu, Z. F. Gao et al., Mod. Phys. Lett. A 31, 1650070 (2016)
    [13] Y. Lim, A. Bhattacharya, J. W. Holt et al., Phys. Rev. C 104, L032802 (2021) doi: 10.1103/PhysRevC.104.L032802
    [14] P. C. Chu, Y. Zhou, Y. Y. Jiang et al., Eur. Phys. J. C 81, 93 (2021) doi: 10.1140/epjc/s10052-020-08800-3
    [15] Z. F. Gao, H. Shan, and H. Wang, Astron. Nachr. 342, 369-376 (2021) doi: 10.1002/asna.202113936
    [16] H. Wang, Z. F. Gao, H. Y. Jia et al., Universe 6, 63 (2020) doi: 10.3390/universe6050063
    [17] P. B. Demorest, T. Pennucci, S. M. Ransom et al., Nature 467, 1081 (2010) doi: 10.1038/nature09466
    [18] E. Fonseca, T. T. Pennucci, J. A. Ellis et al., ApJ 832, 167 (2016) doi: 10.3847/0004-637X/832/2/167
    [19] J. Antoniadis, P. C. C. Freire, N. Wex et al., Science 340, 448 (2013)
    [20] H. T. Cromartie, E. Fonseca, S.M. Ransom et al., Nat. Astron. 4, 72 (2020) doi: 10.1038/s41550-019-0880-2
    [21] A. Burrows and J. M. Lattier, Astrophy. 307, 178 (1986) doi: 10.1086/164405
    [22] S. G. Zhou, Phys. Scripta 91, 063008 (2016) doi: 10.1088/0031-8949/91/6/063008
    [23] J. Schaffner, C. B. Dover, and A. Gal, Ann. Phys. 235, 35 (1994) doi: 10.1006/aphy.1994.1090
    [24] N. K. Glendenning, Compact Stars: Nuclear Physics, Particle Physics, and General Relativity, Springer-Verlag, New York, Inc, 1997
    [25] N. K. Glendenning, Phys. Lett. B 185, 275 (1987) doi: 10.1016/0370-2693(87)90999-3
    [26] N. K. Glendenning, Nucl. Phys. A 469, 600 (1987) doi: 10.1016/0375-9474(87)90015-7
    [27] R. C. Tolman, Phys. Rev. 55, 364 (1939) doi: 10.1103/PhysRev.55.364
    [28] J. R. Oppenheimer and G. M. Volkoff, Phys. Rev. 55, 374 (1939) doi: 10.1103/PhysRev.55.374
    [29] S. Typel and H. H. Wolter, Nucl. Phys. A 656, 331 (1999) doi: 10.1016/S0375-9474(99)00310-3
    [30] N. K. Glendenning, ApJ 293, 470 (1985) doi: 10.1086/163253
    [31] N. K. Glendenning and S. A. Moszkowski, Phys. Rev. Lett. 67, 2414 (1991) doi: 10.1103/PhysRevLett.67.2414
    [32] S. J. Lee et al., Phys. Rev. Lett. 57, 2916 (1986) doi: 10.1103/PhysRevLett.57.2916
    [33] J. Schaffner and I. N. Mishustin, Phys. Rev. C 53, 1416 (1996) doi: 10.1103/PhysRevC.53.1416
    [34] J. Schaffner-Bielich and A. Gal, Phys. Rev. C 62, 034311 (2000) doi: 10.1103/PhysRevC.62.034311
    [35] S. Weissenborn, D. Chatterjee, and J. Schaffner-Bielich, Nucl. Phys. A 881, 62 (2012) doi: 10.1016/j.nuclphysa.2012.02.012
    [36] A. Gal, E. V. Hungerford, and D. J. Millener, Rev. Mod. Phys. 88, 035004 (2016) doi: 10.1103/RevModPhys.88.035004
    [37] C. J. Batty, E. Friedman, and A. Gal, Phys. Rep. 287, 385 (1997) doi: 10.1016/S0370-1573(97)00011-2
    [38] T. Harada, Y. Hirabayashi, and A. Umeya, Phys. Lett. B 690, 363 (2010) doi: 10.1016/j.physletb.2010.05.053
    [39] T. E. Riley, A. L. Watts, S. Bogdanov et al., ApJL 887, L21 (2019) doi: 10.3847/2041-8213/ab481c
    [40] M. C. Miller, F. K. Lamb, A. J. Dittmann et al., ApJL 887, L24 (2019) doi: 10.3847/2041-8213/ab50c5
    [41] X. F. Zhao, Astrophys. Space Sci. 362, 95 (2017) doi: 10.1007/s10509-017-3073-4
  • 加载中

Cited by

1. Li, Y.-H., Ma, W.-Q., Wang, H. Equation of state and anisotropy of pressure of magnetars[J]. Astronomische Nachrichten, 2024. doi: 10.1002/asna.20230180

Figures(5) / Tables(3)

Get Citation
Xian-Feng Zhao, Shui-Rong Zhong and Jian-Li Huo. Properties of proto neutron star PSR J0737-3039A[J]. Chinese Physics C. doi: 10.1088/1674-1137/ac76a6
Xian-Feng Zhao, Shui-Rong Zhong and Jian-Li Huo. Properties of proto neutron star PSR J0737-3039A[J]. Chinese Physics C.  doi: 10.1088/1674-1137/ac76a6 shu
Milestone
Received: 2022-04-30
Article Metric

Article Views(1411)
PDF Downloads(44)
Cited by(1)
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:

Properties of proto neutron star PSR J0737-3039A

    Corresponding author: Xian-Feng Zhao, zhaopioneer.student@sina.com
  • School of Sciences, Southwest Petroleum University, Chengdu 610500, China

Abstract: Using five sets of nucleon coupling constants (DD-MEI, GL85, GL97, GM1, and NL2), we find that the radius of the PNS PSR J0737-3039A is R=15.693–18.846 km, the central baryon density is ρc=0.247–0.359 fm3, the central energy density is εc=4.30×1014– 6.49×1014 gcm3, and the central pressure is pc=3.79×1034 – 5.85×1034 dynecm2. With DD-MEI, GL85, GL97, and GM1, baryon octets appear in the PNS PSR J0737-3039A. With NL2, only the baryons n, p, Λ, Σ0, Σ, Ξ0, and Ξ are present. Corresponding to the same baryon density, the relative densities of the same baryon in the PNS PSR J0737-3039A calculated using different nucleon coupling constants differ greatly. The central relative baryon densities of the PNS PSR J0737-3039A calculated using different nucleon coupling constants also differ greatly.

    HTML

    I.   INTRODUCTION
    • The 2.8-second pulsar J0737-3039B is the companion of the 23-millisecond pulsar J0737-3039A. The two pulsars form a highly relativistic binary neutron star (NS) system [1], and the mass of the NS J0737-3039A should be related to the mass and metallicity of its companion [2]. The NS J0737-3039A has a mass M=1.34 M, which is typical for NSs [1, 3, 4].

      Morrison et al. constructed a numerical model of the binary pulsar J0737-3039A using its measured mass, spin, and moment of inertia to determine the optimal model for different equations of state (EoS) so as to determine the radius of the pulsar. The NS they used, PSR J0737-3039A, has a gravitational mass M=1.337 M [5]. In 2006, Kramer et al. more accurately determined the mass of the NS PSR J0737-3039A to be M = 1.3381±0.0007 M [6].

      In 2008, Worley et al. placed special emphasis on the PSR J0737-3039A component of the extremely relativistic double NS system PSR J0737-3039. They found its moment of inertia to be between 1.30×1045 gcm2 and 1.63×1045 gcm2 [7]. Using well-known universal relations between NS observables, Landry et al. translated the reported 90%-credible bounds on tidal deformability into a direct constraint on the moment of inertia of PSR J0737-3039A [8].

      In 2020, Avdeev et al. developed a parameterized post-Keplerian formalism for hybrid metric-Palatini f(R)-gravity and obtained the mass ratio MA/MB=1.0714. On this basis, they predicted the mass of the NS PSR J0737-3039A as 1.3374 M<M<1.3440 M in the framework of hybrid f(R)-gravity [9].

      An NS has a high density, strong magnetic field, and high rotation speed [1016]. Several massive NSs, PSR J1614-2230 [17, 18], PSR J0348+0432 [19], and PSR J0740+6620 [20], have been observed in recent years. Although the NS PSR J0737-3039A has a significantly smaller mass than these, its mass is typical, and it is a member of a binary NS system. Therefore, theoretical research on this NS is important.

      NSs are formed by the neutrino radiation cooling of proto neutron stars (PNSs) [21]. Studying PNSs is vital for understanding the evolution of stars. Therefore, we are very interested in the PNS PSR J0737-3039A.

      In this study, we investigate the properties of the PNS PSR J0737-3039A with relativistic mean field (RMF) theory [22] considering a baryon octet.

    II.   RMF THEORY OF PNS MATTER
    • The mesons f0(975) (denoted as σ) and ϕ(1020)(denoted as ϕ) represent the interactions between nucleons and hyperons [23]. The Lagrangian density of hadronic matter containing the mesons σ and ϕ is expressed as follows [24]:

      L=B¯ΨB(iγμμmB+gσBσ+gσBσgωBγ0ωgϕBγ0ϕgρBγ0τ3ρ)ΨB12m2σσ213g2σ314g3σ4+12m2ωω2+12m2ρρ212m2σσ2+12m2ϕϕ2+λ=e,μ¯Ψλ(iγμμmλ)Ψλ.

      (1)

      Considering neutrino binding, the baryon partition function of PNS matter is [25, 26]

      lnZB=VTL+B2JB+12π20k2dk{ln[1+e(εB(k)μB)/T]}.

      (2)

      From the above, the total baryon number density is obtained as follows:

      ρ=B2JB+12π2bB0k2nB(k)dk.

      (3)

      The energy density and pressure are respectively

      ε=12m2σσ2+12m2σσ2+13g2σ3+14g3σ4+12m2ωω20+12m2ϕϕ2+12m2ρρ203+B2JB+12π20k2nB(k)dkk2+m2B,

      (4)

      p=12m2σσ212m2σσ213g2σ314g3σ4+12m2ωω20+12m2ϕϕ2+12m2ρρ203+13B2JB+12π20k4k2+m2BnB(k)dk.

      (5)

      Here, nB(k) is the Fermi-Dirac partition function of baryons,

      nB(k)=11+exp[(εB(k)μB)/T],

      (6)

      and mB is the effective mass of baryons,

      mB=mBgσBσgσBσ.

      (7)

      Ignoring the interactions of leptons at finite temperature, we express their partition function as

      lnZL=VTiμ4i24π2[1+2(πTμi)2+715(πTμi)4]+Vλ1π20k2dk{ln[1+e(ελ(k)μλ)/T]},

      (8)

      where the first line represents the contribution of massless neutrinos, and the second line represents the contribution of electrons and μs.

      The lepton number density is

      ρl=1π20k2nl(k)dk,

      (9)

      ρν=π2T2μν+μ3ν6π2

      (10)

      and the contributions of leptons to the energy density and pressure are respectively

      ε=l1π20k2nl(k)dkk2+m2l+ν(7π2T4120+T2μ2ν4+μ4ν8π2),

      (11)

      p=13l1π20k4k2+m2lnl(k)dk+ν1360(7π2T4+30T2μ2ν+15μ4νπ2).

      (12)

      The chemical potentials of baryons are

      μi=μnqi(μeμνe).

      (13)

      To describe neutrino trapping, the content of leptons in a PNS is defined as

      Ylν=Yl+Yνl=ρl+ρνlρ.

      (14)

      We can calculate the mass and radius of a PNS using the Tolman-Oppenheimer-Volkoff (TOV) equations [27, 28]

      dpdr=(p+ε)(M+4πr3p)r(r2M),

      (15)

      M=4πR0εr2dr.

      (16)
    III.   PARAMETERS
    • Here, we calculate the parameters of the PNS PSR J0737-3039A with five sets of nucleon coupling constants (DD-MEI [29], GL85 [30], GL97 [24], GM1 [31], and NL2 [32]). For the outer crust with 4.73×1015 <ρ<8.907×103 fm3, we use the EoS of a BPS [24], and for the region of ρ>8.907×103 fm3, we use the EoS for PNS matter corresponding to the nuclear properties above.

      The ratios of the hyperon coupling constant to nucleon coupling constant are defined as xσh=gσh/gσ, xωh=gωh/gω, and xρh=gρh/gρ, with h denoting the hyperons Λ,Σ, and Ξ.

      xρh and xωh are chosen using SU(6) symmetry [33, 34]. Then, xσh is obtained using [24]

      U(N)h=mn(mnmn1)xσh+(gωmω)2ρ0xωh.

      (17)

      In this study, the hyperon-potentials are chosen as U(N)Λ=30 MeV [34, 35, 36], U(N)Σ=30 MeV [34, 35, 36, 37], and U(N)Ξ=14 MeV [38].

      The coupling constants of the mesons σ and ϕ can be chosen as [23]

      gϕΞ=2gϕΛ=2gϕΣ=22gω/3,

      (18)

       gσΛ/gσ=gσΣ/gσ=0.69,

      (19)

       gσΞ/gσ=1.25.

      (20)

      Temperatures in the interiors of PNSs can reach up to 30 MeV. After that, they can be cooled by neutrino radiation [21]. In this study, the temperature of the PNS PSR J0737-3039A is taken as T=20 MeV.

    IV.   RADIUS OF PNS PSR J0737-3039A
    • The mass M of the PNS is calculated as a function of the radius R. In Fig. 1, the coarse curves indicate the masses of the PNSs PSR J0740+6620, PSR J0348+0432, PSR J1614-2230, and PSR J0737-3039A, and the five fine curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2. We find that only the nucleon coupling constant DD-MEI can give the masses of the PNSs PSR J0740+6620, PSR J0348+0432, PSR J1614-2230, and PSR J0737-3039A. NL2 cannot give the mass of the PNS PSR J0740+6620; however, it can give the mass of the PNSs PSR J0348+0432, PSR J1614-2230, and PSR J0737-3039A. GL85, GL97, and GM1 cannot give the masses of the PNSs PSR J0740+6620, PSR J0348+0432, and PSR J1614-2230; however, they can give the masses of the PNS PSR J0737-3039A. All five sets of nucleon coupling constants give the masses of the PNS PSR J0737-3039A.

      As shown in Fig. 1 and Table 1, the radii of the PNS PSR J0737-3039A given by the five sets of nucleon coupling constants are different. NL2 gives the smallest radius, R=15.693 km, and GL85 gives the largest radius, R=18.846 km. In other words, the radius of the PNS PSR J0737–3039A calculated using these five sets of nucleon coupling constants is R=15.693–18.846 km. The NS PSR J0030+0451 has a mass similar to that of the PNS PSR J0737-3039A; hence, it is logical to compare their radii. The mass and radius of the NS PSR J0030+0451 given by Riley et al. are in the ranges M=1.18–1.49 M and R=11.52–13.85 km [39], respectively, and those given by Miller et al. are M=1.30–1.59 M and R=11.96–14.26 km, respectively [40]. The PNS radius calculated in this study is larger than that calculated by Riley et al. and Miller et al., which is reasonable because the radius of a PNS is larger than that of an NS of the same mass [41].

      Figure 1.  (color online) Mass M of the PNS calculated as a function of the radius R. The coarse curves indicate the masses of the PNSs PSR J0740+6620, PSR J0348+0432, PSR J1614-2230, and PSR J0737-3039A. The five fine curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

      Parameter R/km ρc/fm3 εc×1015/(gcm)3 pc×1035/(dynecm)2
      DD-MEI 18.006 0.247 0.430 0.379
      GL85 18.846 0.290 0.519 0.427
      GL97 17.353 0.359 0.649 0.585
      GM1 16.798 0.310 0.545 0.506
      NL2 15.693 0.293 0.508 0.506

      Table 1.  Results calculated in this study. The mass of the PNS PSR J0737-3039A is M=1.338 M. R, ρc, εc, and pc are the radius, central baryon density, central energy density, and central pressure of the PNS PSR J0737-3039A, respectively.

      The central baryon densities ρc of the PNS PSR J0737-3039A calculated with different nucleon coupling constants are different. ρc of the PNS PSR J0737-3039A calculated using GL97 is the largest, ρc=0.359 fm3, whereas ρc calculated using DD-MEI is the smallest, ρc=0.247 fm3. That is, the central baryon density calculated using the five nucleon coupling constants is ρc=0.247–0.359 fm3.

    V.   EFFECTIVE MASSES OF BARYONS AND THE CHEMICAL POTENTIALS OF NEUTRONS AND ELECTRONS
    • The effective mass of nucleons, Λ, Σ, and Ξ as a function of the baryon density ρ are shown in Fig. 2. We denote the center effective mass of nucleons, Λ, Σ, and Ξ as mNc, mΛc, mΣc, and mΞc, respectively. As shown in Fig. 2, the effective masses of nucleons, Λ, Σ, and Ξ all decrease with increasing baryon density ρ. At the center of the PNS, the effective mass of the baryon reaches a minimum. Corresponding to the same baryon density ρ, the effective mass calculated using DD-MEI is the smallest, whereas that calculated using GL97 is the largest. As shown in Table 2, the center effective masses of nucleons, Λ, Σ, and Ξ calculated using the above five groups of nucleon coupling constants are mNc=384.521–606.324 MeV, mΛc=778.280–911.405 MeV, mΣc=939.187–1080.822 MeV, and mΞc=1152.637–1218.964 MeV, respectively.

      Figure 2.  (color online) Effective mass of baryons in the PNS PSR J0737-3039A as a function of baryon density. The five curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

      Parameter mNc/MeV mΛc/MeV mΣc/MeV mΞc/MeV μnc/fm1 μec/fm1
      DD-MEI 384.521 778.280 939.187 1152.637 5.340 0.721
      GL85 606.324 911.405 1080.822 1218.964 5.421 0.943
      GL97 573.011 888.351 1071.635 1206.096 5.516 0.989
      GM1 478.656 831.015 1006.070 1179.389 5.418 0.763
      NL2 394.604 780.860 1015.662 1200.968 5.367 0.829

      Table 2.  Effective mass of baryons and the chemical potentials of n and e in the PNS PSR J0737-3039A calculated in this study. mNc, mΛc, mΣc, and mΞc are the center effective masses of nucleons N, Λ, Σ, and Ξ, respectively. μnc and μec are central chemical potentials of n and e, respectively.

      The chemical potentials of neutrons and electrons are shown in Fig. 3. The chemical potentials of neutrons μn and electrons μe increase with increasing baryon density ρ. Different nucleon coupling constants result in different chemical potentials relative to the same baryon density ρ. The central neutron chemical potential of the five nucleon coupling constants is calculated as μnc=5.340–5.516 fm1, and the central electron chemical potential is calculated as μec=0.721–0.989 fm1.

      Figure 3.  (color online) Chemical potentials of neutrons and electrons in the PNS PSR J0737-3039A as a function of baryon density. The five curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

    VI.   ENERGY DENSITY AND PRESSURE IN PNS PSR J0737-3039A
    • The energy density ε and pressure p in the PNS PSR J0737-3039A as a function of the baryon density ρ are shown in Fig. 4. Both the energy density ε and pressure p increase with increasing baryon density ρ. Corresponding to the same baryon density ρ, the energy densities calculated using different nucleon coupling constants have little difference; however, the calculated pressures differ greatly. The central energy density and central pressure calculated for the five sets of nucleon coupling constants are different. GL97 results in the maximum center energy density and center pressure.

      Figure 4.  (color online) Energy density and pressure in the PNS PSR J0737-3039A as a function of baryon density. The five curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

      As shown in Table 1, the minimum central energy density given with DD-MEI is εc=4.30×1014 gcm3. GL97 gives the largest central energy density, εc=6.49×1014 gcm3. The central energy density calculated from the five nucleon coupling constants is εc=4.30×1014–6.49×1014 gcm3. Similarly, the central pressure calculated using the five nucleon coupling constants is pc=3.79×1034–5.85×1034 dynecm2.

    VII.   RELATIVE PARTICLE DENSITY IN PNS PSR J0737-3039A
    • Figure 5 shows the relative particle density ρi/ρ in the PNS PSR J0737-3039A as a function of the baryon density ρ.

      Figure 5.  (color online) Relative particle density in the PNS PSR J0737-3039A as a function of baryon density. The five curves represent the results of the five nucleon coupling constants DD-MEI, GL85, GL97, GM1, and NL2.

      As shown in Fig. 5 and Table 3, baryons n, p, Λ, Σ+, Σ0, Σ, Ξ0 and Ξ all appear in the PNS PSR J0737-3039A when using the nucleon coupling constants DD-MEI, GL85, GL97, and GM1. However, with NL2, only baryons n, p, Λ, Σ0, Σ, Ξ0, and Ξ are present, while baryon Σ+ is absent.

      Parameter ρnc/ρ% ρpc/ρ% ρΛc/ρ% ρΣ+c/ρ% ρΣ0c/ρ% ρΣc/ρ% ρΞ0c/ρ% ρΞc/ρ%
      DD-MEI 90.7 8.8 0.5 1E-6 2.0E-4 0.04 5.1E-5 0.02
      GL85 81.0 17.5 1.5 2.9E-3 4.0E-4 5.7E-5 1.3E-3 0.05
      GL97 80.5 16.5 2.8 1.3E-3 4.1E-4 1.3E-4 2.4E-3 0.18
      GM1 90.7 8.4 0.8 2E-6 1.8E-4 0.02 1.07E-4 0.05
      NL2 87.5 11.5 1.0 0 1.3E-5 5.2E-4 1.5E-5 5.8E-3

      Table 3.  Central relative particle densities in the PNS PSR J0737-3039A. ρnc/ρ, ρpc/ρ, ρΛc/ρ, ρΣ+c/ρ, ρΣ0c/ρ, ρΣc/ρ, ρΞ0c/ρ, and ρΞc/ρ are the central relative particle densities of n, p, Λ, Σ+, Σ0, Σ, Ξ0, and Ξ, respectively, in the PNS PSR J0737-3039A.

      As shown in Fig. 5, corresponding to the same baryon density, the relative particle densities of the same baryon in the PNS PSR J0737-3039A calculated using different nucleon coupling constants differ greatly. As shown in Table 3, the central relative baryon densities of the PNS PSR J0737-3039A calculated using different nucleon coupling constants also differ greatly. For example, for hyperon Λ, the minimum central relative baryon density given using DD-MEI is ρΛc/ρ=0.5%, whereas GL97 gives the largest central baryon density, ρΛc/ρ=2.8%. In other words, the central baryon density of the Λ-hyperon calculated using the five nucleon coupling constants is ρΛc/ρ=0.5%–2.8%. Similarly, the central baryon density of n calculated using five groups of nucleon coupling constants is ρnc/ρ=80.5%–90.7%, the central baryon density of p is ρpc/ρ=8.4%–17.5%, the central baryon density of the Σ0 hyperon is ρΣ0c/ρ=1.3×103%–4.1×104%, the central baryon density of the Σ hyperon is ρΣc/ρ=5.7×105%–0.04%, the central baryon density of the Ξ0 hyperon is ρΞ0c/ρ=1.5×105%–2.4×103%, and the central baryon density of the Ξ hyperon is ρΞc/ρ=5.8×103%–0.18%. The hyperon Σ+ in the PNS PSR J0737-3039A calculated using NL2 does not appear, whereas the central baryon density of Σ+ calculated using DD-MEI, GL85, GL97, and GM1 is ρΣ+c/ρ=1.0×106%–2.9×103%.

    VIII.   SUMMARY
    • In this study, the influence of nucleon coupling constants on the properties of the PNS PSR J0737-3039A is investigated using RMF theory. Here, five sets of nucleon coupling constants (DD-MEI, GL85, GL97, GM1, and NL2) are used, and the temperature of the PNS PSR J0737-3039A is chosen as T=20 MeV.

      The radius R, central baryon density ρc, central energy density εc, and central pressure pc of the PNS PSR J0737-3039A calculated using different nucleon coupling constants are different. Using these five sets of nucleon coupling constants, we find that the radius of the PNS PSR J0737-3039A is R=15.693–18.846 km, the central baryon density is ρc=0.247–0.359 fm3, the central energy density is εc=4.30×1014–6.49×1014 gcm3, and the central pressure is pc=3.79×1034–5.85×1034 dynecm2. Although the PNS radius should be larger than that of an NS of the same mass, it should not be significantly larger. As shown in Table 1, the radius of the PNS PSR J0737-3039A calculated using the nucleon coupling constant NL2 is the smallest, R=15.693 km, and is closest to the radius calculated by Riley et al. and Miller et al. Therefore, NL2 seems to be the most reasonable coupling constant.

      For four sets of nucleon coupling constants (DD-MEI, GL85, GL97 and GM1), baryons n, p, Λ, Σ+, Σ0, Σ, Ξ0, and Ξ all appear in the PNS PSR J0737-3039A. However, using NL2, only baryons n, p, Λ, Σ0, Σ, Ξ0, and Ξ are present, while baryon Σ+ does not appear.

      Corresponding to the same baryon density, the relative densities of the same baryon in the PNS PSR J0737-3039A calculated using different nucleon coupling constants differ greatly. The central relative baryon densities of the PNS PSR J0737-3039A calculated using different nucleon coupling constants also differ greatly. The central baryon density of n calculated using the five groups of nucleon coupling constants is ρnc/ρ=80.5%–90.7%, the central baryon density of p is ρpc/ρ=8.4%–17.5%, the central baryon density of Λ is ρΛc/ρ=0.5%–2.8%, the central baryon density of the Σ0 hyperon is ρΣ0c/ρ=1.3×103%–4.1×104%, the central baryon density of the Σ hyperon is ρΣc/ρ=5.7×105%–0.04%, the central baryon density of the Ξ0 hyperon is ρΞ0c/ρ=1.5×105%–2.4×103%, and the central baryon density of the Ξ hyperons is ρΞc/ρ=5.8×103%–0.18%. The hyperon Σ+ does not appear with NL2, whereas the central baryon density of hyperon Σ+ calculated using DD-MEI, GL85, GL97, and GM1 is ρΣ+c/ρ=1.0×106%–2.9×103%.

Reference (41)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return