Young Crab-like pulsars and luminous X-ray sources in starbursts and optically dull galaxies

Recent Chandra observations of nearby galaxies have revealed a number of ultraluminous X-ray sources (ULXs) with super-Eddington luminosities, away from the central regions of non-active galaxies. The nature of these sources is still debated. We argue that a fraction of them could be young, Crab-like pulsars, the X-ray luminosity of which is powered by rotation. We use the pulsar birth parameters estimated from radio pulsar data to compute the steady-state pulsar X-ray luminosity distribution as a function of the star formation rate (SFR) in the galaxy. We find that ~10% of optically dull galaxies are expected to have a source with L_x >~ 10^{39} erg/s, while starbursts galaxies should each have several of these sources. We estimate that the X-ray luminosity of a few percents of galaxies is dominated by a single bright pulsar with L_x >~10^{39} erg/s, roughly independently of its SFR. We discuss observational diagnostics that can help distinguish the young pulsar population in ULXs.Comment: 17 pages, 4 figures, accepted to Ap

Deep Radio Imaging of Globular Clusters and the Cluster Pulsar Population

We have obtained deep multifrequency radio observations of seven globular clusters using the Very Large Array and the Australia Telescope Compact Array. Five of these, NGC 6440, NGC 6539, NGC 6544, NGC 6624 and Terzan 5 had previously been detected in a shallower survey for steep spectrum radio sources in globular clusters (Fruchter and Goss 1990). The sixth, the rich globular cluster, Liller 1, had heretofore been undetected in the radio, and the seventh, 47 Tucanae, was not included in our original survey. High resolution 6 and 20 cm images of three of the clusters, NGC 6440, NGC 6539, NGC 6624 reveal only point sources coincident with pulsars which have been discovered subsequent to our first imaging survey. 21 and 18 cm images reveal several point sources within a few core-radii of the center of 47 Tuc. Two of these are identified pulsars, and a third, which is both variable and has a steep spectrum, is also most likely a pulsar previously identified by a pulsed survey. However, the 6, 20 and 90 cm images of NGC 6544, Liller 1 and Terzan 5 display strong steep-spectrum emission which cannot be associated with known pulsars. The image of the rich cluster Terzan 5 displays numerous point sources within $30''$, or 4 core radii of the cluster center. The density of these objects rises rapidly toward the core, where an elongated region of emission is found. The brightest individual sources, as well as the extended emission, possess the steep spectra expected of pulsars. Furthermore, the flux distribution of the sources agrees well with the standard pulsar luminosity function. The total luminosity and number of objects observed suggest that Terzan 5 contains more pulsars than any other Galactic globular cluster.Comment: 33 pages, 6 Postscript figures; Accepted for publication in the Astrophysical Journal; abstract abridged. PDF version also available at http://nemesis.stsci.edu/~fruchter/fg99/fg99.pd

Proper-Motion Measurements with the VLA. II. Observations of Twenty-eight Pulsars

Using the Very Large Array, we have measured the proper motions of twenty-eight radio pulsars. On average, the pulsars studied are fainter and more distant than those studied in earlier work, reducing the selection biases inherent in surveys restricted to the Solar neighborhood. The typical measurement precision achieved is a few milliarcseconds per year, corresponding to a few tens of kilometers per second for a pulsar a kiloparsec away. While our results compare well with higher-precision measurements done using very-long baseline interferometry, we find that several earlier proper motion surveys appear to have reported overly optimistic measurement uncertainties, most likely because of a failure to fully account for ionospheric effects. We discuss difficulties inherent in estimating pulsar velocities from proper motions given poorly constrained pulsar distances. Our observations favor a distribution with 20% of pulsars in a low velocity component (sigma_1D = 99 km/s) and 80% in a high velocity component (sigma_1D = 294 km/s). Furthermore, our sample is consistent with a scale height of pulsar birthplaces comparable to the scale height of the massive stars that are their presumed progenitors. No evidence is found in our data for a significant population of young pulsars born far from the plane. We find that estimates of pulsar ages based on kinematics agree well with the canonical spin-down age estimate, but agreement is improved if braking indexes are drawn from a Gaussian distribution centered at n=3 with width 0.8.Comment: 20 pages. Accepted for publication in the Astronomical Journa

Birth and Evolution of Isolated Radio Pulsars

We investigate the birth and evolution of Galactic isolated radio pulsars. We begin by estimating their birth space velocity distribution from proper motion measurements of Brisken et al. (2002, 2003). We find no evidence for multimodality of the distribution and favor one in which the absolute one-dimensional velocity components are exponentially distributed and with a three-dimensional mean velocity of 380^{+40}_{-60} km s^-1. We then proceed with a Monte Carlo-based population synthesis, modelling the birth properties of the pulsars, their time evolution, and their detection in the Parkes and Swinburne Multibeam surveys. We present a population model that appears generally consistent with the observations. Our results suggest that pulsars are born in the spiral arms, with a Galactocentric radial distribution that is well described by the functional form proposed by Yusifov & Kucuk (2004), in which the pulsar surface density peaks at radius ~3 kpc. The birth spin period distribution extends to several hundred milliseconds, with no evidence of multimodality. Models which assume the radio luminosities of pulsars to be independent of the spin periods and period derivatives are inadequate, as they lead to the detection of too many old simulated pulsars in our simulations. Dithered radio luminosities proportional to the square root of the spin-down luminosity accommodate the observations well and provide a natural mechanism for the pulsars to dim uniformly as they approach the death line, avoiding an observed pile-up on the latter. There is no evidence for significant torque decay (due to magnetic field decay or otherwise) over the lifetime of the pulsars as radio sources (~100 Myr). Finally, we estimate the pulsar birthrate and total number of pulsars in the Galaxy.Comment: 27 pages, including 15 figures, accepted by Ap

Are Supershells Powered by Multiple Supernovae? Modeling the Radio Pulsar Population Produced by OB Associations

Traditional searches for radio pulsars have targeted individual small regions such as supernova remnants or globular clusters, or have covered large contiguous regions of the sky. None of these searches has been specifically directed towards giant supershells, some of which are likely to have been produced by multiple supernova (SN) explosions from an OB association. Here we perform a Montecarlo simulation of the pulsar population associated with supershells powered by multiple SNe. We predict that several tens of radio pulsars could be detected with current instruments associated with the largest Galactic supershells (with kinetic energies >~ 10^{53} ergs), and a few pulsars with the smaller ones. We test these predictions for some of the supershells which lie in regions covered by past pulsar surveys. For the smaller supershells, our results are consistent with the few detected pulsars per bubble. For the giant supershell GSH 242-03+37, we find the multiple SN hypothesis inconsistent with current data at the 95% level. We stress the importance of undertaking deep pulsar surveys in correlation with supershells. Failure to detect any pulsar enhancement in the largest of them would put serious constraints on the multiple SN origin for them. Conversely, the discovery of the pulsar population associated with a supershell would allow a different/independent approach to the study of pulsar properties.Comment: accepted to ApJ; 17 pages, 2 figures, 1 tabl

The Halo Beaming Model for Gamma-Ray Bursts

We consider a model for gamma-ray bursts (GRBs) from high-velocity neutron stars in the galactic halo. In this model, bursters are born in the galactic disk with large recoil velocities V_r, and GRBs are beamed to within emission cones of half-angle \phi centered on V_r. We describe scenarios for magnetically -channeled GRBs that have such beaming characteristics. We then make detailed comparisons of this halo beaming model (HBM) to BATSE and PVO data for GRB intensity & angular position distributions. Acceptable fits to observations of over 1000 bursts are obtained for \phi = 15 - 30 degrees and for a BATSE sampling depth ~ 180 kpc. Present data favor a truly isotropic (cosmological) model over the HBM, but not by a statistically compelling margin. Bursters born in nearby external galaxies, such as M31, are almost entirely undetectable in the HBM because of misdirected beaming. We analyze several refinements of the basic HBM: gamma-ray intensities that vary with angle from the beam axis; non-standard-candle GRB luminosity functions; and models including a subset of bursters that do not escape from the galaxy. We also discuss the energy budgets for the bursters, the origins of their recoils, and the physics of burst beaming and alignment. One possible physical model is based on the magnetar model of soft gamma repeaters (SGRs). Empirical bounds on the rate of formation and peculiar velocities of SGRs imply that there exist ~ 10^4 to ~ 10^7 aged SGRs in the galactic halo within a distance of 100 kpc. The HBM gives an acceptable fit to observations only if it satisfies certain conditions (e.g. \phi ~ 20 deg) which are possible, but for which there exist no clear & compelling theoretical justifications. The cosmological burster hypothesis is more generic and thus more attractive in this sense. (Abbreviated Abstract).Comment: ApJ accepted, 9 figures, AASTE

G328.4+0.2 : A large and luminous Crab-like supernova remnant

We report on radio continuum and HI observations of the radio source G328.4+0.2 using the Australia Telescope Compact Array. Our results confirm G328.4+0.2 to be a filled-center nebula with no surrounding shell, showing significant linear polarization and an almost flat spectral index. These results lead us to conclude that G328.4+0.2 is a Crab-like, or ``plerionic'', supernova remnant (SNR), presumably powered by an unseen central pulsar. HI absorption towards G328.4+0.2 puts a lower limit on its distance of 17.4 +/- 0.9 kpc, making it the largest (D=25 pc) and most luminous (L_R = 3e35 erg/s) Crab-like SNR in the Galaxy. We infer G328.4+0.2 to be significantly older than the Crab Nebula, but powered by a pulsar which is fast spinning (P<20 ms) and which has a comparatively low magnetic field (B<1e12 G). We propose G328.4+0.2, G74.9+1.2 and N157B as a distinct group of large-diameter, high-luminosity Crab-like SNRs, all powered by fast-spinning low-field pulsars.Comment: 7 pages, 3 embedded EPS figures, uses emulateapj.sty. Accepted to ApJ. Abstract corrected so that distance is now in kpc, not pc

Thin accretion disc with a corona in a central magnetic field

We study the steady-state structure of an accretion disc with a corona surrounding a central, rotating, magnetized star. We assume that the magneto-rotational instability is the dominant mechanism of angular momentum transport inside the disc and is responsible for producing magnetic tubes above the disc. In our model, a fraction of the dissipated energy inside the disc is transported to the corona via these magnetic tubes. This energy exchange from the disc to the corona which depends on the disc physical properties is modified because of the magnetic interaction between the stellar magnetic field and the accretion disc. According to our fully analytical solutions for such a system, the existence of a corona not only increases the surface density but reduces the temperature of the accretion disc. Also, the presence of a corona enhances the ratio of gas pressure to the total pressure. Our solutions show that when the strength of the magnetic field of the central neutron star is large or the star is rotating fast enough, profiles of the physical variables of the disc significantly modify due to the existence of a corona.Comment: Accepted for publication in Astrophysics & Space Scienc

Pulsar Wind Nebulae in Evolved Supernova Remnants

For pulsars similar to the one in the Crab Nebula, most of the energy input to the surrounding wind nebula occurs on a timescale of less than 1000 years; during this time, the nebula expands into freely expanding supernova ejecta. On a timescale 10,000 years, the interaction of the supernova with the surrounding medium drives a reverse shock front toward the center of the remnant, where it crushes the PWN (pulsar wind nebula). One- and two-dimensional, two-fluid simulations of the crushing and re-expansion phases of a PWN show that (1) these phases are subject to Rayleigh-Taylor instabilities that result in the mixing of thermal and nonthermal fluids, and (2) asymmetries in the surrounding interstellar medium give rise to asymmetries in the position of the PWN relative to the pulsar and explosion site. These effects are expected to be observable in the radio emission from evolved PWN because of the long lifetimes of radio emitting electrons. The scenario can explain the chaotic and asymmetric appearance of the Vela X PWN relative to the Vela pulsar without recourse to a directed flow from the vicinity of the pulsar. The displacement of the radio nebulae in G327.1--1.1, MSH15--56 (G326.3--1.8), G0.9+0.1, and W44 relative to the X-ray nebulae may be due to this mechanism. On timescales much greater than the nebular crushing time, the initial PWN may be mixed with thermal gas and become unobservable, so that even the radio emission is dominated by recently injected particles.Comment: 22 pages, 9 figures; submitted to Ap