Mon. Not. R. Astron. Soc. 388, 393–415 (2008) doi:10.1111/j.1365-2966.2008.13402.x
Populating the Galaxy with pulsars – I. Stellar and binary evolution
Paul D. Kiel, Jarrod R. Hurley, Matthew Bailes and James R. Murray Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Hawthorn, Victoria 3122, Australia
Accepted 2008 April 29. Received 2008 April 2; in original form 2007 December 20
ABSTRACT Downloaded from https://academic.oup.com/mnras/article/388/1/393/1013977 by guest on 28 September 2021 The computation of theoretical pulsar populations has been a major component of pulsar studies since the 1970s. However, the majority of pulsar population synthesis has only regarded isolated pulsar evolution. Those that have examined pulsar evolution within binary systems tend to either treat binary evolution poorly or evolve the pulsar population in an ad hoc manner. Thus, no complete and direct comparison with observations of the pulsar population within the Galactic disc has been possible to date. Described here is the first component of what will be a complete synthetic pulsar population survey code. This component is used to evolve both isolated and binary pulsars. Synthetic observational surveys can then be performed on this population for a variety of radio telescopes. The final tool used for completing this work will be a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Results provided here support the need for further (apparent) pulsar magnetic field decay during accretion, while they conversely suggest the need for a re-evaluation of the assumed typical millisecond pulsar formation process. Results also focus on reproducing the observed P P˙ diagram for Galactic pulsars and how this precludes short time-scales for standard pulsar exponential magnetic field decay. Finally, comparisons of bulk pulsar population characteristics are made to observations displaying the predictive power of this code, while we also show that under standard binary evolutionary assumption binary pulsars may accrete much mass. Key words: binaries: close – stars: evolution – stars: neutron – pulsars: general – Galaxy: stellar content.
produced impart significant kicks to the pulsars, that make their 1 INTRODUCTION survival prospects within binary systems bleak. Since the tentative suggestion that neutron stars (NSs) form from Today, there are well over 100 binary pulsars known, and their violent supernova (SN) events (Baade & Zwicky 1934a,b) and the spin periods and inferred magnetic field strengths offer the oppor- discovery of pulsars by Hewish et al. (1968) the number of observed tunity to attempt models of binary evolution and pulsar spin-up pulsars has risen dramatically. Surveys for radio pulsars have dis- that explain their distribution in the pulsar magnetic field–spin pe- covered over 1500 objects including a rich harvest of binary and riod (Bs–P) diagram. To do this properly, one should take models millisecond pulsars (Manchester et al. 2005; Burgay et al. 2006). of an initial population of zero-age main-sequence (ZAMS) single Precision timing of pulsars in binary systems has not only allowed stars and binaries, trace the binary and stellar evolution, including precise tests of theories of relativistic gravity (e.g. van Straten et al. NS spin-up effects, calculate their Galactic trajectories and initial 2001), but also given insights into their masses and the nature of distribution in the Galaxy, and then perform synthetic surveys as- the binary systems they inhabit (e.g. Bell et al. 1993; Verbiest et al. suming a pulsar luminosity and beaming function. This is what we 2008). Out of the first 300 or so pulsars to be discovered, only three wish to achieve. The large number of assumptions that we require were members of binary systems, despite the progenitor population to complete this effort caution against the absolute predictive power having an extremely high binary fraction (>50 per cent; Duquennoy of such a model. For instance, it is easy to demonstrate that trying & Mayor 1991). This paucity of binary pulsars was an important to use such a model to predict something like the merger rate of clue about the origin and evolution of pulsars. Clearly, something NS–black hole (BH), or double pulsars/NSs based solely on a con- about pulsar generation was contributing to the disruption of bi- sideration of the total number and mass distributions of un-evolved nary systems. We now believe that the SNe in which pulsars are ZAMS binaries would be folly. However, the relative numbers of two or more populations can often only depend on relatively few model assumptions (as shown in Hurley, Tout & Pols 2002, here- after HTP02; Belczynski, Kalogera & Bulik 2002; O’Shaughnessy E-mail: [email protected] et al. 2008). With enough observables in time, we might hope to