arxiv: v1 [astro-ph.he] 9 Apr 2019

Transcription

1 MNRAS 000, 5 208) Preprint 0 April 209 Compiled using MNRAS LATEX style file v3.0 Electron-positron pairs in hot plasma of accretion column in right X-ray pulsars arxiv: v [astro-ph.he] 9 Apr 209 Alexander A. Mushutkov,,2,3 Igor S. Ognev, 3 and Dmitrij I. Nagirner 4 Leiden Oservatory, Leiden University, NL-2300RA Leiden, The Netherlands 2 Space Research Institute of the Russian Academy of Sciences, Profsoyuznaya Str. 84/32, Moscow 7997, Russia 3 P. G. Demidov Yaroslavl State University, Sovietskaya 4, Yaroslavl, Russia 4 Soolev Astronomical Institute, Saint Petersurg State University, Saint-Petersurg 98504, Russia 0 April 209 INTRODUCTION The majority of accreting X-ray inaries is represented y the systems hostineutron stars NSs). Highly magnetized NSs form a special class of ojects among X-ray inaries powered y accretion - X-ray pulsars XRPs, see e.g. Walter et al. 205). Typical magnetic field strength at the NS surface in XRPs is 0 2 G, which is confirmed independently y a numer of different methods: detection of cyclotron lines see Stauert et al. 209 for review), transitions into the propeller state Tsygankov et al. 206; Lutovinov et al. 207), detection of spin-up and spin-down effects Sugizaki et al. 207). Detected luminosities of XRPs cover a few orders of magnitude from 0 33 erg s and up to 0 4 erg s, where the rightest pulsars elong to the recently discovered class of pulsating ultraluminous X-ray sources ULXs, Bachetti et al., 204; Israel et al., 207). The theoretical explanation of ULX pulsars is still under deates, ut the most of the theories agree that the extreme al.mushtukov@gmail.com AAM) ABSTRACT The luminosity of X-ray pulsars powered y accretion onto magnetized neutron stars covers a wide range over a few orders of magnitude. The rightest X-ray pulsars recently discovered as pulsating ultraluminous X-ray sources reach accretion luminosity aove 0 40 erg s which exceeds the Eddington value more than y a factor of ten. Most of the energy is released within small regions in the vicinity of magnetic poles of accretineutron star - in accretion columns. Because of the extreme energy release within a small volume accretion columns of right X-ray pulsars are ones of the hottest places in the Universe, where the internal temperature can exceed 00 kev. Under these conditions, the processes of creation and annihilation of electron-positron pairs can e influential ut have een largely neglected in theoretical models of accretion columns. In this letter, we investigate properties of a gas of electron-positron pairs under physical conditions typical for accretion columns. We argue that the process of pairs creation can crucially influence oth the dynamics of the accretion process and internal structure of accretion column limiting its internal temperature, dropping the local Eddington flux and increasing the gas pressure. Key words: accretion, accretion discs X-rays: inaries neutrinos stars: neutron radiative transfer accretion onto magnetized NS results in the formation of accretion columns aove the stellar surface, where the matter is confined y a strong magnetic field and produce luminosity well aove the Eddington limit Basko & Sunyaev 976; Wang & Frank 98; Mushtukov et al. 205a). A strong magnetic field of a NS in XRP modifies oth the geometry of accretion flow see Chapter 6 in Frank et al. 2002) and asic properties of matter Harding & Lai 2006). Because of a strong magnetic field, the accreting material reaches NS surface in small regions the typical area is 0 0 cm 2 ) in the vicinity of NS magnetic poles. At extremely high mass accretion rates the radiation pressure is high enough to stop accretion flow aove NS surface in radiation dominated shock Lyuarskii & Syunyaev 982; Mushtukov et al. 205). Below the shock, the matter slowly settles to the stellar surface converting its gravitational and kinetic energy into emission in X-ray energy and. A strong energy release within a small region results in extreme temperatures, which are typically aout a few kev Wang & Frank 98; Mushtukov et al. 205a) and c 208 The Authors

2 2 A. A. Mushtukov et al. can e as high as aout MeV in the case of ULX pulsars Mushtukov et al. 208). High temperatures in regions of energy release of XRPs result in specific processes which might shape properties of XRPs at extreme mass accretion rates and thus have to e taken into account in theoretical models. In particular, high temperatures result in the creation of electron-positron pairs and possily strong emission of neutrino due to pair annihilation and/or cyclotron neutrino emission Kaminker et al. 992, 994; Mushtukov et al. 208). In this papers, we investigate properties of electronpositron gas under conditions of a strong magnetic field and high temperatures typical for internal regions of optically thick accretion columns. Under the assumption of thermodynamic equilirium, we calculate the numer density of the pairs Section 2.) and the gas pressure along magnetic field lines Section 2.2). We discuss the influence of the pair creation process on internal temperature and dynamics of the accretion process. 2 ELECTRON-POSITRON PAIRS IN ACCRETION COLUMN The accretion columns in right XRPs are confined y a strong magnetic field and supported y the internal radiation and gas pressure Basko & Sunyaev 976). The conditions of matter inside the columns are determined y the local mass density and temperature. The main mechanism of opacity in accretion columns is Compton scattering of photons y electrons/positrons influenced y a strong magnetic field. Because the accretion column is optically thick due to the scattering, the photons created in the accretion flow undergo a numer of scatterings efore they leave the column. A fraction of radiation is truly asored while the photons diffuse towards the edges of the accretion channel. Let consider a case of fully ionized hydrogen plasma. A numer density of protons n p is determined y the mass accretion rate Ṁ, cross section of the accretion channel S and local velocity v = βc of accretion flow n p Ṁ 2Svm p 0 23 Ṁ 20 S 0β cm 3, ) where Ṁ20 = Ṁ/020 gs and S 0 = S/0 0 cm 2. Because the typical geometrical thickness of accretion column is of order of d 0 4 cm, the optical thickness of the column across magnetic field lines due to the scattering can e estimated as τ n pσ d Ṁ 20S 20 β σ /σ T), where σ is ascattering cross section across B-fieldlines and σ T is the Thomson cross section. The numer of scattering experienced y photon in accretion column is N sc τ Ṁ 2 20 S 2 20 β2 σ σ T ) 2. 2) The asorption cross section can e roughly estimated as σ a 23.7σ TZ 4 E/keV) 3, where Z is the atomic numer, E is a photon energy Bethe & Salpeter 957). The accreting material is slowing down from a free-fall velocity v ff to aout v ff /0 in the radiation dominated shock at the top of accretion column Thus, the velocity elow the shock region can eestimated from aove as v < 0.05c and atypical numer of asorption events corresponding to a given numer of scatterings can e estimated as σ a Ṁ N a N sc 30 20Z 2 4 σ S20 2 β/0.05)2 σ σ T ) ) 3 E. 3) The estimated numer of asorption events will e even larger if one will account for the actual chemical composition of accreting material, which is affected y nuclear reactions in accretion column In particular, for the case of caron dominated material Z 6 and an estimated numer of asorption events 3) increases y a factor of 0 3. Thus, the photons emitted in the accretion column will e likely asored on their way to the edges of the accretion channel in the case of sufficiently large mass accretion rates. Under these conditions, the gas of electron/positron pair can e considered to e in thermodynamic equilirium Bisnovatyi-Kogan et al. 97). It worth to note that the outer layers of the accretion column might e far from thermodynamic equilirium and more detailed analyses accounting for detailed alance etween a numer of processes is required see e.g. Guilert & Stepney 985). Another essential condition for estalishing of the equilirium concentrations and distriutions of electrons and positrons is sufficiently long dynamical time scale in accretion column: the dynamical time scale has to e longer than the time scale required for equilirium estalishing. Because the annihilation cross section is of order of σ ann 0 24 cm 2 see e.g. Daugherty & Bussard 980) and typical velocity of electrons and positrons are of order of speed of light at temperatures > 50 kev, the time scale of equilirium estalishing can e estimated as t eq σ annn pc S 0β Ṁ 20 s, while the dynamical time scale in the column can e estimated as t dyn R βc β s. Thus, the time scale of equilirium estalishing is much smaller than the dynamical time scale t eq t dyn ) and pairs in the central parts of accretion column can e considered to e in the equilirium. 2. Numer density of electrons and positrons At low temperature regime, the numer density of electrons is similar to the numer density of protons n n p, while at high temperature regime, the creation of electronpositron pairs ecomes essential and the numer density of leptons oth electrons and positrons) can significantly exceed the numer density of protons. The numer density of electrons n and positrons n + in the equilirium are given y Canuto & Chiu 968; Kaminker et al. 992) n = mc ) 3 4π dp zf 4) dp zf cm 3, where = B/B crit is dimensionless magnetic field strength, B crit m 2 ec 2 /eħ G, the distriution functions MNRAS 000, 5 208)

3 Pairs in accretion columns of X-ray pulsars 3 µ m e c 2, kev B=0 3 G n p =0 25 cm -3 n p =0 29 cm -3 np =027 cm -3 n - + n +, cm a n p =0 29 cm -3 n p =0 27 cm -3 n p =0 25 cm -3 B=0 3 G µ m e c 2, kev B=0 3 G B=0 4 G B=0 5 G n - + n +, cm B=0 5 G B=0 4 G B=0 3 G Figure. The chemical potential of electrons calculated as a function of temperature for different numer densities of protons at fixed magnetic field strength 0 3 G top panel), and for different magnetic field strength at fixed numer density of protons n p = 0 25 cm 3 ottom panel). of particles f = exp [ Enp z) µ t ] +), 5) dimensionless energy E np z) = + p 2 z + 2n) /2 of electron/positron at nth Landau level, dimensionless temperature t T/m ec 2 ), is the spin degeneracy of Landau level g 0 = and = 2 for n ), and the dimensionless in units of m ec 2 ) chemical potential µ is determined y the condition of electro-neutrality n n + = n p. In the limit of zero magnetic field and relatively low temperatures T m ec 2 ) the numer densities of electrons and positrons can e estimated as Zeldovich & Novikov 97) n + n mec ) 3e /t t 3/2 π e /t t 3/2 cm 3, 6) while at high temperatures T m ec 2 ) the approximations is given y n + n.82 mec ) 3 t t 3 cm 3. 7) π 2 The calculations with the accurate equations 4) and 5) give a result which depends oth on numer density of protons and magnetic field strength see Fig.,2). Influence of a strong magnetic field on the chemical potential and numer density of electron-positron pairs ecomes valuale at B 0 3 G. At lower B-field strength and temperature Figure 2. Numer density of electrons and positrons as a function of temperature calculated for a) different numer densities of protons: 0 25 lack solid line), 0 27 rown dashed lile) and 0 29 cm 3 red dashed-dotted line), and ) different magnetic field strength: 0 3 lack solid line), 0 4 dark-lue dashed line) and 0 5 G lue dashed-dotted line). Black dotted line represents approximation given y 6). T 0keV, the chemical potential and numer density of the pairs can e calculated in approximation of zero field strength. In the case of low numer density of protons n p, the contriution of created electron-positron pairs to the total numer density of leptons is dominant at temperatures of a few tens of kev already see Fig.2a), which is typical temperature for accretion column interior Mushtukov et al. 205a). Note, that the stronger the magnetic field, the larger the numer density of electron-positron pairs see Fig. 2). The fast increase of a numer density of electronpositron pairs with temperature requires that a fraction of energy of accreting material is going into a process of pair creation. Because the energy udget of the accretion process is limited y the mass accretion rate and compactness of a NS, the process of pairs creation limits the increase of internal temperature in the accretion column. The upper limit of a numer density of the pairs can e otained from the assumption that the energy of accretion flow is going entirely in a pair creation. The total energy udget per one aryon due to the accretion process is determined y local free-fall velocity: γ ff γ)m pc 2 where γ ff = β 2 ff) /2 is the Lorentz factor due to a local free-fall velocity and γ is actual Lorentz factor of accreting material at a given height aove stellar surface. Then the numer density of the pairs, which is similar to the numer density of positrons n + in accretion flow, can e limited from aove y γ ff γ)m pn p > 2m en +. 8) MNRAS 000, 5 208)

4 4 A. A. Mushtukov et al. T max, kev n p =0 29 cm -3 n p =0 27 cm -3 n p =0 25 cm -3 B=0 3 G B=0 4 G β ff Figure 3. The maximal temperature achievale in accretion column as a function of local dimensionless free fall velocity β ff. Different curves are given for various numer density of protons in the accretion flow: n p = 0 25, 0 27 and 0 29 cm 3 down up). Black solid and dark-lue dashed line correspond to local magnetic field strength of 0 3 and 0 4 G respectively. At temperatures aove MeV grey area) the energy losses are strongly affected y neutrino emission. P, dyn cm -2 P, dyn cm a n p =0 29 cm -3 n p =0 27 cm -3 n p =0 25 cm -3 B=0 3 G B=0 5 B=0 4 G G B=0 3 G Thus, a conservative upper limit for the numer density of pairs is n +,max np 2 m p m e γ ff ). 9) The maximal numer density of pairs determines the maximal temperature T max achievale at a given local free-fall velocity β ff and numer density of protons see Fig.3). One can see that the temperatures of a few hundred kev in accretion channel are achievale only in the case of sufficiently high densities see Fig. 3). 2.2 Gas and radiation pressure in accretion column The gas pressure along z-direction due to the electron and positrons can e derived from the distriution function of the particles: P z,ee + = m4 ec 5 3 4π n= dp zβ zp z[f +f +] 0) where the dimensionless velocity along z axis β z = p zγ = p z+p 2 z +2n) /2. Expression 0) can e rewritten as P z,ee ) 0 dp z p 2 z [f p z)+f +p z)] +p 2 z +2n) /2 dyn cm 2. The results of calculations with equation ) are given in Fig. 4. At low temperatures, when the effects of pair creation are negligile and P T, thenthepressure increases rapidly with the increase of the numer density of electron-positron pairs. In the limiting case of high temperatures, P T 4 and the gas pressure ecomes very close to the radiation Figure 4. The pressure of electron-positron gas along magnetic field lines calculated as a function of temperature for a) different numer densities of protons: 0 25 lack solid line), 0 27 rown dashed lile) and 0 29 cm 3 red dashed-dotted line), and ) different magnetic field strength: 0 3 lack solid line), 0 4 darklue dashed line) and 0 5 G lue dashed-dotted line). Black dotted line represents radiation pressure calculated according to 2). pressure, which is given in the equilirium y P rad = 3 at T 4 kev dyn cm 2, 2) where a is the radiative constant. 2.3 The Eddington flux affected y electron-positron pairs The Eddington flux is the photon energy flux which is high enough to compensate the gravitational attraction of thecentral oject. The Eddingtonfluxplaysakeyrole in the theory of super-eddington accretion onto compact ojects and particularly in the theory of the accretion column in right XRPs Basko & Sunyaev 976; Wang & Frank 98; Mushtukov et al. 205a). The Eddington flux is given y F Edd = ρc n +n +)σ eff GM r 2 = µ e m pc σ eff GM r 2, 3) where σ eff is the effective cross section, which is affected y a strong magnetic field and local energy spectrum of photons, and the mean molecular weight is µ e 837n p +2n +. 4) 836 n p +2n + The temperature increase and the corresponding increase of electron/positron numer density result in a drop of the MNRAS 000, 5 208)

5 Pairs in accretion columns of X-ray pulsars 5 µ e n p =0 25 cm -3 n p =0 27 cm -3 n p =0 29 cm -3 reduces the mean molecular weight of the accreting material see Fig.5) and, thus, drops the Eddington flux, which supports accretion columns. According to our knowledge, these effects were largely neglected in the models of accretion columns developed up to date and have to e taken into account in further numerical models. 0-3 B=0 3 G Figure 5. The mean molecular weight as a function of temperature for different numer densities of protons: 0 25 lack solid line), 0 27 rown dashed lile) and 0 29 cm 3 red dashed-dotted line) at magnetic field strength fixed at 0 3 G. mean molecular weight see Fig. 5). This affects the local Eddington flux and can principally make it smaller y more than 3 orders of magnitude. The accurate calculations of the Eddington flux require accounting for special features of Compton scattering in a strong magnetic field see e.g. Herold 979; Daugherty & Harding 986; Mushtukov et al. 206), which is ehind the scape of this letter. 3 SUMMARY AND DISCUSSION We have investigated properties of electron-positron gas in accretion columns of right X-ray pulsars accounting for specific physical conditions of a strong magnetic field and high temperatures. Because of a large optical thickness of accretion column due to Compton scattering, the gas of electron-positron pairs can e considered to e close to the thermodynamic equilirium in the central parts of accretion column at high mass accretion rates typical for recently discovered ULX pulsars. However, in the outer parts accretion column the gas of pairs can e far from the equilirium and its conditions has to e otained from kinetic equations accounting for a numer of processes see e.g. Guilert & Stepney 985). We have demonstrated that the process of pair creation might have key importance in extreme accretion onto magnetized NSs influencing the dynamical and temperature structure of accretion columns. In the case of high velocity of the accretion flow and relatively low mass density, the numer density of pairs ecomes comparale to the numer densityofprotonsattemperaturesofaoutafewtensofkev already see Fig. 2). Because the energy udget of the accretion process is limited, the process of pair creation puts an upper limit on the internal temperature of the accretion column see Fig.3) and, therefore, on the energy losses due to neutrino production Mushtukov et al. 208). The temperatures of a few hundred kev are achievale only in the case of the sufficiently high mass density of material in the accretion channel and only in the very vicinity of a NS surface. The electron-positron pairs affect the gas pressure, which ecomes comparale to the radiation pressure at temperatures 00keV see Fig.4). The process of pair creation ACKNOWLEDGEMENTS AAM and ISO acknowledge support y the Russian Science Foundation Grant No This work was also supported y the Netherlands Organization for Scientific Research Veni Fellowship AAM). We are grateful to Alexander Kaminker, Vitaly Grigoriev and an anonymous referee for discussion and useful comments. REFERENCES Bachetti M., et al., 204, Nature, 54, 202 Basko M. M., Sunyaev R. A., 976, MNRAS, 75, 395 Bethe H. A., Salpeter E. E., 957, Quantum Mechanics of Oneand Two-Electron Atoms Bisnovatyi-Kogan G. S., Zel dovich Y. B., Syunyaev R. A., 97, Soviet Ast., 5, 7 Canuto V., Chiu H.-Y., 968, Physical Review, 73, 220 Daugherty J. K., Bussard R. W., 980, ApJ, 238, 296 Daugherty J. K., Harding A. K., 986, ApJ, 309, 362 Frank J., King A., Raine D. J., 2002, Accretion Power in Astrophysics: Third Edition Guilert P. W., Stepney S., 985, MNRAS, 22, 523 Harding A. K., Lai D., 2006, Reports on Progress in Physics, 69, 263 Herold H., 979, Phys. Rev. D, 9, 2868 Israel G. L., et al., 207, Science, 355, 87 Kaminker A. D., Gnedin O. Y., Yakovlev D. G., Amsterdamski P., Haensel P., 992, Phys. Rev. D, 46, 433 Kaminker A. D., Gnedin O. Y., Yakovlev D. G., Amsterdamski P., Haensel P., 994, Astronomical and Astrophysical Transactions, 4, 283 Lutovinov A. A., Tsygankov S. S., Krivonos R. A., Molkov S. V., Poutanen J., 207, ApJ, 834, 209 Lyuarskii Y. E., Syunyaev R. A., 982, Soviet Astronomy Letters, 8, 330 Mushtukov A. A., Nagirner D. I., Poutanen J., 206, Phys. Rev. D, 93, Mushtukov A. A., Suleimanov V. F., Tsygankov S. S., Poutanen J., 205a, MNRAS, 454, 2539 Mushtukov A. A., Suleimanov V. F., Tsygankov S. S., Poutanen J., 205, MNRAS, 447, 847 Mushtukov A. A., Tsygankov S. S., Suleimanov V. F., Poutanen J., 208, MNRAS, 476, 2867 Stauert R., et al., 209, A&A, 622, A6 Sugizaki M., Mihara T., Nakajima M., Makishima K., 207, PASJ, 69, 00 Tsygankov S. S., Lutovinov A. A., Doroshenko V., Mushtukov A. A., Suleimanov V., Poutanen J., 206, A&A, 593, A6 Walter R., Lutovinov A. A., Bozzo E., Tsygankov S. S., 205, A&ARv, 23, 2 Wang Y.-M., Frank J., 98, A&A, 93, 255 Zeldovich Y. B., Novikov I. D., 97, Relativistic astrophysics. Vol.: Stars and relativity, University of Chicago Press This paper has een typeset from a TEX/LATEX file prepared y the author. MNRAS 000, 5 208)