Neutron star spin–kick velocity correlation effect on binary neutron star coalescence rates and spin–orbit misalignment of the components

We study the effect of the neutron star spin–kick velocity alignment observed in young radio pulsars on the coalescence rate of binary neutron stars. Two scenarios are considered for neutron star formation: when the kick is always present, and when it is small or absent if a neutron star is formed in a binary system as a result of electron-capture degenerate core collapse. The effect is shown to be especially strong for large kick amplitudes and tight alignments, reducing the expected galactic rate of binary neutron star coalescence compared to calculations with randomly directed kicks. The spin–kick correlation also leads to a much narrower neutron star spin–orbit misalignment.     

Pulsar spin–velocity alignment from single and binary neutron star progenitors

The role of binary progenitors of neutron stars (NSs) in the apparent distribution of space velocities and spin–velocity alignment observed in young pulsars is studied. We performed a Monte Carlo synthesis of pulsar populations originated from single and binary stars with different assumptions about the NS natal kick (kick–spin alignment, kick amplitude and kick reduction in electron-capture supernovae in binary progenitors with initial main-sequence masses from the range  8–11 M_⊙  which experienced mass exchange due to Roche lobe overflow). The calculated spin–velocity alignment in pulsars is compared with data inferred from radio polarization measurements. The observed space velocity of pulsars is found to be mostly affected by the natal kick velocity form and its amplitude; the fraction of binaries is not important here for reasonably large kicks. The natal kick–spin alignment is found to strongly affect the spin–velocity correlation of pulsars. Comparison with the observed pulsar spin–velocity angles favours a sizeable fraction of binary progenitors and kick–spin angles  ∼5°–20°  .     

Population statistics study of radio and gamma-ray pulsars in the Galactic plane

We present results of our pulsar population synthesis of normal and millisecond pulsars in the Galactic plane. Over the past several years, a program has been developed to simulate pulsar birth, evolution and emission using Monte Carlo techniques. We have added to the program the capability to simulate millisecond pulsars, which are old, recycled pulsars with extremely short periods. We model the spatial distribution of the simulated pulsars by assuming that they start with a random kick velocity and then evolve through the Galactic potential. We use a polar cap/slot gap model for γ-ray emission from both millisecond and normal pulsars. From our studies of radio pulsars that have clearly identifiable core and cone components, in which we fit the polarization sweep as well as the pulse profiles in order to constrain the viewing geometry, we develop a model describing the ratio of radio core-to-cone peak fluxes. In this model, short period pulsars are more cone-dominated than in our previous studies. We present the preliminary results of our recent study and the implications for observing these pulsars with GLAST and AGILE.      

Velocities and birthplaces of pulsars

The values of the proper motions of 36 pulsars are analysed. It is shown that the birthplaces of pulsars are located within the OB-associations and/or in the galactic arms. Eleven cases of pulsar formation in associations, and fourteen cases of formation in the galactic arms, have been established. More precise parameters of pulsars — such as: distances, ages, spatial velocity components —, and the spatial velocity distribution of pulsars is constructed.     

Spatial Distribution of Pulsars

The spatial distributions of pulsars with ages of less than 10⁶ years and more than 10⁶ years are shown to be quite different. The spatial distributions of pulsars with ages of less than 106 years and of the remnants of supernova outbursts are essentially the same. They lie near the galactic plane in a thin zone of width 400 pc. The overwhelming majority of pulsars with ages exceeding 106 years lie outside this zone. These facts suggest a genetic relationship between pulsars and supernova remnants. The spatial distribution of pulsars with different emission periods also supports this point.     

Evolution of massive binary stars in the LMC and its implications for radio pulsar population

The Hertzsprung-Russell diagram of the Large Magellanic Cloud compiled recently by Fitzpatrick & Garmany (1990) shows that there are a number of supergiant stars immediately redward of the main sequence although theoretical models of massive stars with normal hydrogen abundance predict that the region 4.5 ≤ logT _eff ≤ 4.3 should be un-populated (“gap”). Supergiants having surface enrichment of helium acquired for example from a previous phase of accretion from a binary companion, however, evolve in a way so that the evolved models and observed data are consistent — an observation first made by Tuchman & Wheeler (1990). We compare the available optical data on OB supergiants with computed evolutionary tracks of massive stars of metallicity relevant to the LMC with and without helium-enriched envelopes and conclude that a large fraction (≈ 60 per cent) of supergiant stars may occur in binaries. As these less evolved binaries will later evolve into massive X-ray binaries, the observed number and orbital period distribution of the latter can constrain the evolutionary scenarios of the supergiant binaries. The distributions of post main sequence binaries and closely related systems like WR + O stars are bimodal-consisting of close and wide binaries in which the latter type is numerically dominating. When the primary star explodes as a supernova leaving behind a neutron star, the system receives a kick and in some cases can lead to runaway O-stars. We calculate the expected space velocity distribution for these systems. After the second supernova explosion, the binaries in most cases, will be disrupted leading to two runaway neutron stars. In between the two explosions, the first born neutron star’s spin evolution will be affected by accretion of mass from the companion star. We determine the steady-state spin and radio luminosity distributions of single pulsars born from the massive stars under some simple assumptions. Due to their great distance, only the brightest radio pulsars may be detected in a flux-limited survey of the LMC. A small but significant number of observable single radio pulsars arising out of the disrupted massive binaries may appear in the short spin period range. Most pulsars will have a low velocity of ejection and therefore may cluster around the OB associations in the LMC.     

Distribution of free electrons in the galactic plane

A method is proposed for constructing a map of the electron density distribution in the galactic plane. Data on the dispersion measures DM of more than 1500 pulsars and their distances, found by an independent method, are used. The independent distances of the pulsars are estimated using an empirical relation of the form L=γP ^α Ṗ ^β W ^δ between the radio luminosity L of the pulsars and their periods P, the rate of change Ṗ. of their periods, and the half width W of their pulses. A map of the electron density distribution in the galactic plane within a ±400 pc layer is provided. __________ Translated from Astrofizika, Vol. 49, No. 2, pp. 277–287 (May 2006).     

Asymmetric supernova explosions and the origin of binary pulsars

We investigate the effect of asymmetric supernova explosions on the orbital parameters of binary systems with a compact component. We relate such explosions to the origin of binary pulsars. The degree of asymmetry of the explosion is represented by the kick velocity gained by the exploding star due to the asymmetric mass ejection. The required kick velocity to produce the observed parameters of the binary pulsar PSR 1913+16 should be larger than 80 km s^–1 if the mass of the exploding star is larger than 4M . We examine the mean survival probability of the binary system (<f>) for various degrees of asymmetry in the explosion. The rare occurrence of a binary pulsar does not necessarily imply that such a probability is low since not all pulsars have originated in a binary system. Assuming the birth rate of pulsars derived by Taylor and Manchester (1977), we derive that <f> would be as high as 0.25. Such values of <f> can be obtained if the mass of the exploding stars is, in general, not large (10M ).     

Relation of pulsars to the remnants of supernova bursts

Based on a large volume of statistical data it is shown that the spatial distributions of radio pulsars in the galaxy with characteristic ages T ≤ 10 ⁶ years and T > 10⁶ years differ significantly. The overwhelming majority of the pulsars with T ≤ 10 ⁶ years lie within a narrow band of width 400 pc around the galactic plane. A large portion of the pulsars with T > 10⁶ years is concentrated outside this zone. In the case of younger pulsars, a larger fraction of them lies within the confines of the above mentioned zone. It is also shown that pulsars with T ≤ 10 ⁶ years and the remnants of supernova explosions have essentially the same spatial distribution. These facts support the existence of a relationship between pulsars and supernova remnants, as well as the acquisition of high spatial velocities by pulsars during their birth. __________ Translated from Astrofizika, Vol. 49, No. 1, pp. 103–110 (February 2006).     

Magnetic Field in Selected Directions of the Galaxy, Direction of the Spiral Arm of Sagittarius

Faraday rotation data on 40 pulsars are used in a detailed study of the magnetic field and its fluctuations in the direction of the spiral arm of Sagittarius. These results mostly agree with standard models for the galactic magnetic field. A magnetic field on the order of 3.2 G is directed from galactic longitude l ₀=55° (toward the sun). However, an asymmetry has been found in the degrees of rotation relative to a plane lying in the southern hemisphere parallel to the galactic plane and at a distance of 50-60 pc from it. All the pulsars with measures of dispersion greater than 30 pc·cm^-3 and lying to the north of this plane have positive measures of rotation which increase linearly with distance, while the pulsars lying to the south of this plane have unusually absolutely low negative measures of rotation. We propose that the spiral arm of Sagittarius lies entirely to the north of this plane, while the negative measures of rotation of the pulsars below this plane are caused by the magnetic field of the halo of the southern hemisphere of the galaxy. The magnetic field in the arm of Sagittarius is regular to a great extent and its fluctuating component is roughly half the regular component.     

Pulsars and interstellar medium: Multiple regression analysis of related parameters

An empirical relation which relates the 408 MHz galactic continuum background temperature (408GCBT) to dispersion measures, position and radio-luminosity of 325 pulsars is obtained by means of multple stepwise regression analysis. This relation showns that pulsars may be considered as galactic probes for the distribution of 408GCBT and interstellar electron density (IED) in interstellar medium (ISM).Peculiar pulsars (O-C±2.5) point out galactic regions where the observedT ₄₀₈ are higher (or lower) andn _e lower (or higher) than the averaged ones.Normal pulsars (–2.5T ₄₀₈ andn _e are in agreement, on, the average.Standard pulsars (O-C±0.05) show galactic regions where observed and computedT ₄₀₈ andn _e are in good agreement. Recent models of pulsar disk systems, suggested by Michel and Dessler (1981) could justify the conclusions drawn for peculiar pulsars having O-C>2.5.     

The galactic surface density derived from K giants

We compare theW velocity dispersions of Brosche, Schwan & Schwarz (2001) with more recent results. The increase with the distance |z | from the galactic plane is confirmed, although perhaps with reduced amplitude. This could be interpreted either as one homogenous population or as a superposition of (at least) two populations with two constant dispersions and two scale heights. For each of the possibilities we propose a simple model. Combined with two observational variants for the velocity variation, we obtain surface densities up to |z | = 250 pc ranging from 5 to 44 M_⊙ pc^–2. Thus the case for considerable dark matter in the neighbourhood of the galactic plane is not supported     

Interstellar electron density

We impose the requirement that the spatial distribution of pulsars deduced from their dispersion measures using a model of the galactic electron density (n _e ) should be consistent with cylindrical symmetry around the galactic centre (assumed to be 10 kpc from the Sun). Using a carefully selected subsample of the pulsars detected by the II Molonglo Survey (II MS), we test a number of simple models and conclude that (i) the effective mean 〈n_e〉) for the whole galaxy is 0.037_-0.012 ^+0.020 cm^-3, (ii) the scale height of electrons is greater than 300 pc and probably about 1 kpc or more, and (iii) there is little evidence for variation of n_e with galactic radius R_GC for R_Gc ≳ 5 kpc. Further, we make a detailed analysis of the contribution to n_e from H II regions. Combining the results of a number of relatively independent calculations, we propose a model for the galactic electron density of the formn _e (z) = 0.030 + 0.020 exp (- |z|/70) cm^-3 where z(pc) is the height above the galactic plane and the second term describes the contribution from H II regions. We believe the statistical uncertainties in the parameters of this model are quite small.     

Pulsar sensitivity studies of the GLAST large area telescope

In this contribution we present our preliminary investigation on pulsar sensitivity of the Large Area Telescope, the main instrument aboard the GLAST mission. In particular we concentrated our attention to pulsars located at low galactic latitudes. We created a set of simulated pulsars having different fluxes in an array of galactic coordinates separated by a distance greater than the LAT Point Spread Function in order to avoid confusion between adjacent sources. Galactic gamma-ray sky background as used during the second LAT Data Challenge (DC2) is also included. We then run an automatic routine for testing periodicity for all the pulsars considering an opportune timing solution. In this way we can obtain a map of the sensitivity of the periodic searches for different fluxes and for various Galactic latitudes. Some assumptions have been made by simulating the pulsar sources, but this study is a first step toward an estimate for pulsed emission sensitivity of the GLAST LAT. The pulsed flux sensitivity profile we generate could also be input to a population synthesis code of Galactic pulsars in order to obtain more accurate predictions of the number of expected pulsar detections by GLAST. On behalf of the GLAST LAT Collaboration.     

The simulation of the magnetic field and spin period evolution of accreting neutron stars

By means of the Monte Carlo method, we simulate the evolutionary distribution of accreting neutron stars (NSs) in the magnetic field versus spin period (B‐P) diagram where the accretion induced magnetic‐field decay model is exploited. The simulated results show that by mass accretion the B‐P distribution of the accreting NS would evolve along the equilibrium period line to a region with low field and short period. The B‐P distributions of the simulated accreting NSs are consistent with those of the observed millisecond pulsars (MSPs) after accretion of ∼ 0.1–0.2 M⊙. We also test the effects of the initial magnetic field and the spin period on the evolved B‐P distribution of the accreting NSs. It is shown that the evolved distributions of the simulated samples are independent of the selection of the initial condition when the NS magnetic field decays to a value less than ∼10¹⁰ G. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

How Fast Could a Proto-pulsar Rotate?

According to two estimated relations between the initial period andthe dynamo-generated magnetic dipole field of pulsars, we calculate the statisticaldistributions of pulsar initial periods. It is found that proto-pulsars are very likelyto have rotation periods between 20 ms and 30 ms, and that most of the pulsarsrotate initially at a period ＜ 60 ms. Our result supports the asymmetric neutrinoemission model for pulsar kick.     

An Estimation of the Initial Period of Pulsars

The initial period of a pulsar is an important factor in our understanding of the formation of neutron stars and of the nature of the equation of state of neutron star matter.Up to now this quantity can only be obtained for a few pulsars for which accurate age and braking index are known.Based on the theory of the offcenter dipole emission,in which pulsars obtain theiry high velocities depending on the initial periods,we calculate the initial period using the proper motion data,Because the orbital velocity of the progenitor and asymmetric kick in the supernova explosion may also contribute to the observed velocity of the pusar,the derived values of initial periods are lower limits.For normal pulsars,the initial periods are in the range of 0.6～2.6ms.For the millisecond pulsars,the initial periods are comparable to their current periods,and the ratio between the initial period and the current period increases with the decrease of the current period.For PSR B1937 21 with the shortest period of 1.56ms,the ratio is 0.77.     

Pulsars: Space kinematic characteristics and their connection with supernova remnants

The galactic distribution of pulsars and shell remnants of supernovae (SN) as investigated on the basis of newly-estimated parameters. Special attention was paid to taking into account all possible selection effects and an attempt was made to reveal a statistically-pure ensemble of objects. On the basis of this ensemble we studied the radial andz-distribution of pulsars and supernova remnants (SNR).It is shown that the radial distribution of both objects is considered to have an annular structure with the maximum surface density at a distance of 4.5–5.3 kpc (if the distance of the Sun from the galactic centre is assumed to be equal to 8.5 kpc). The scle-height of the progenitors of SNRs is not more than 110 pc and only 15% of the SNRs, whose progenitors may also be massive runaway stars, are situated at 300 pc. The mean application of the pulsars is not more than 300 pc which agrees with the hypothesis about the genetic connection with type-II SN outbursts at the kinematic aget _k5×10⁶ yr and thez-component of spatial velocity beingV _z=120 km s^–1.The possible precursors of type-I SNRs by the shape of their radial distribution in late spirals and the various model calculations given in the literature, are also discussed.     

How abundant is the population of binary radio pulsars with black holes?

Using a 'scenario machine' we have carried out a population synthesis of radio pulsars with black hole binaries (BH+Psr) in the context of the most widespread assumptions concerning star mass loss during evolution, the mass ratio distribution of binary stars, the kick velocity and the envelope mass loss during collapse. Our purpose is to show that under any plausible parameters for the evolution scenario, the BH+Psr population should be abundant in the Galaxy. It is shown that in all models (including those evolved by Heger et al. and Woosley, Heger & Weaver), the expected number of black holes paired with radio pulsars is sufficient for such systems to be discovered within the next few years.     

Very high energy γ-rays from galactic sources

The field of Very High Energy (VHE) gamma ray astronomy using the Atmospheric Cerenkov Technique has entered an interesting phase with detection of various galactic and extragalactic sources. Among galactic sources, only the Crab nebula has been established as a standard candle. Most observations on pulsars are in agreement as to the necessity for the GeV spectra to steepen at < 200 GeV. While the Imaging method for increase of sensitivity has been successful with many results, an alternate technique —Wavefront Sampling Technique- is also being used by an increasing number of experiments. The recently commissioned experiment at Pachmarhi (PACT) in India is presented as an example of this technique. Preliminary results from this experiment show detection of VHE γ-ray emission from (a) the Crab nebula at a high significance and (b) Crab and Geminga pulsars at > 1.5 TeV which could be the second component of the Outer Gap pulsar models.