Be / X -Ray 双星中的中子星自传演化

1. Be/X-Ray双星中的中子星自传演化 成忠群 南京大学 2013.8.17

2. Contents • 1. Introduction • (1) Observed period gap for BeXBs • (2) Possible interpretation by the authors • 2. What we have did • (1) The theoretical considerations • (2) Our model and detailed calculation process • (3) The results of our calculation • 3. Contrast and discussion • (1) Something we have ignored • (2) Comparing and discussion

3. 1.1 What's BeXBs? HMXBs SGXBs BeXBs Super Giant NS NS Be/star Characteristics of BeXBs: 1. Orbital period 10-100s days 2. High orbital eccentricity 3. Mass transfer driven by stellar wind or equatorial disk 4. Drastically X-ray flux change between outburst and quiet states Be/star Characteristics: 1. Usually 8-18 solar mass 2. Rapid rotation 3. Equatorial disk 4. Infrared excess 5. Emission lines of H and He

4. 1.1 HMXBs on the Corbet diagram. • The bimodal spin period distribution of BeXBs, with two peaks at 10s and 200s and a spin period gap around 40s. • The two types of BeXBs have a percentage of 35% and 65% in the total BeXBs, respectively. Ref. (Christian Knigge, etal 2011)

5. 1.2 The equilibrium spin period of X-ray binaries • For disk accretion neutron star, the spin will reach a equilibrium spin period when the co-rotation radius equal the magnetosphere radius.

6. 1.1 HMXBs on the Corbet diagram. • The bimodal spin period distribution of BeXBs, with two peaks at 10s and 200s and a spin period gap around 40s. • The two types of BeXBs have a percentage of 35% and 65% in the total BeXBs, respectively. Ref. (Christian Knigge, etal 2011)

7. 1.1 Two spin period peaks for BeXBs.

8. 1.2 Is the supernovae responsible for the bimodal? Two types of supernovae Electron-capturing Supernovae (ECS) Core collapse Supernovae (CCS) Lower kick velocity (<50km/s) Less NS mass (<1.3 solar mass) Higher kick velocity (>200km/s) Large NS mass (~1.4 solar mass) Shorter orbital periods and little eccentricity Longer orbital periods and higher eccentricity Lower spin period Higher spin period

9. 1.2 Problems with the explanation. • 1. The equilibrium spin period of wind-fed pulsar is determined by the mass accretion rate and magnetic field. If the two groups of pulsar is formed by the supernova, then the orbital period gap should be more clear than the spin period gap. • 2. The mass accretion rate of pulsar in BeXBs change violently, then there's no reason to neglectthe influence of pulsar spin evolution history on the bimodal distribution.

10. 2.1 The equatorial disk of Be/star • (1). Within the equatorial disk • (Water 1989): • (2). Outside the equatorial disk or it's disappeared:

11. 2.1 Mass transfer rate for BeXBs • Within the equatorial disk: • Outside the equatorial disk or during the disk-loss episodes: • Compare with the X-ray luminosity:

12. 2.1 Specific angular momentum capture rate. • (2). Within the spherical stellar wind, numerical simulation shows (Livio etal 1986, Ruffert 1999): • (1). Within the equatorial disk and the NS orbital plane coplanar with the equatorial disk of Be/star (Wang 1981, Li & van den Heuvel 1996 ):

13. 2.2 Accretion flow configuration. • If the accretion material circularization radius is bigger than the magnetosphere radius, the accretion flow can evolve into accretion disk before contact the NS magnetic field (Ghosh and Lamb 1979). • Otherwise a atmosphere envelope probably formed replace the accretion disk to mediate the energy and angular momentum transfer between the NS magnetic field and accretion flow (Davies & Pringle 1981, Ikhasanov 2007, Shakura etal 2011) • Spherical stellar wind: • Equatorial disk:

14. 2.2 The influence of wind structure change. • The observed Be/star equatorial disk formed and disappeared erratically on a timescale range from months to years, so both the mass and specific angular momentum captured from wind change alternatively from between the two cases. • If the accretion flow can evolve into accretion disk during the equatorial disk existing duration propeller Spin down Disk accretion Spin down propeller Spin up

15. 2.2 The influence of wind structure change. • For the case of no accretion disk formed around NS, the subsonic propeller set in and the co-rotation radius can exceed the magnetosphere radius. The co-rotation radius can be an order of magnitude larger than the magnetosphere radius in spherical accretion (Shakura etal 2011). The spin up torque depend on the coupling between the envelope and the NS magnetic field, the limit spin up torque is when all the captured momentum transfer to the NS. propeller Spin down Subsonic propeller Spin down Subsonic accretion Spin down propeller Spin up not efficient

16. 2.2 The detailed spin evolution paths in our model New born NS: 10-100 ms Propeller: two spin down rates work alternatively Accretion disk No accretion disk

17. 2.2 Parameters input • Binaries population synthesis (BPS) to produce BeXBs with new born NS and particular distribution of orbital period, Be/star mass, other parameters is assumed to be: • NS mass, radius: 1.4 solar mass, 10 km • Magnetic field: Gaussian, no decay, Log B=12.3, deviation: 0.3 • Be/star rotation : around [0.7, 0.8] of critical Keplerian rotation rate • Equatorial disk parameters near Be/star: density: 10^[-12, -11] g/cm^3, • velocity: [2, 20] km/s, \alpha=[0.01, 0.1], n=[2, 4] • Spherical stellar wind: mass lose rate (Nieuwenhuijzen 1990), \beta=0.8 • v_\infinity =[0.4, 1] v_escape. • Time proportion for two wind structure: k = T_disk/T_loss=0.05

18. 2.3 Calculation results

19. 2.3 Calculation results

20. 2.3 Calculation results

21. 2.3 Calculation results

22. 2.3 Calculation results

23. 2.3 Calculation results

24. 3.1 Discussion • High eccentricity and large angle between the orbital plane and equatorial disk plane for BeXBs can be ascribed to the supernovae and kick velocity when the NS was born. A isotropic Maxwellian velocity distribution is imparted to the NS when supernovae happens in the BPS. • For ECS BeXBs : • For CCS BeXBs :

25. 3.1 Discussion

26. 3.1 The contrast of two types of supernovae

27. 3.2 conclusion • We have calculated the spin evolution of NS in a large number of BeXBs with particular orbital and mass distribution get from the BPS. Our result showed that the observed bimodal spin period distribution of BeXBs should be ascribed to the pulsar spin evolution history rather than the supernova channel. • Accretion disk formed BeXBs will stay at a lower equilibrium spin period since the spin up of pulsar is much efficient than the spin down. They are responsible to the observed lower spin period peak of BeXBs. • No accretion disk formed BeXBs will evolve to a higher spin period and then formed the other peak of the bimodal peak

28. Thanks !