Search

We present updated analyses of pulse profiles and their arrival-times from PSR B1534+12, a 37.9-ms radio pulsar in orbit with another neutron star. A high-precision timing model is derived from twenty-two years of timing data, and accounts for all astrophysical processes that systematically affect pulse arrival-times. Five "post-Keplerian" parameters are measured that represent relativistic corrections to the standard Keplerian quantities of the pulsar's binary orbit. These relativistic parameters are then used to test general relativity by comparing the measurements with their predicted values. We conclude that relativity theory is confirmed to within 0.17% of its predictions. Furthermore, we derive the following astrophysical results from our timing analysis: a distance of $d_{\rm GR} = 1.051 \pm 0.005$ kpc to the pulsar-binary system, by relating the "excess" orbital decay to Galactic parameters; evidence for pulse "jitter" in PSR B1534+12 due to short-term magnetospheric activity; and evolution in pulse-dispersion properties. As a secondary study, we also present several analyses on pulse-structure evolution and its connection to relativistic precession of the pulsar's spin axis. The precession-rate measurement yields a value of $\Omega_1^{\rm spin}$ = $0.59^{+0.12}_{-0.08}$ $^{\circ}$/year (68% confidence) that is consistent with expectations, and represents an additional test of relativistic gravity., Comment: Accepted for publication in ApJ on 6 April 2014. v3: added references, a structure-function analysis/discussion, edited MCMC figure and timing-error analysis, corrected minor errors in calculations for distance and MCMC results, and used most up-to-date planetary ephemeris for timing

We report a 5.4\sigma\ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 \mu s, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density, Comment: 5 figures, 2 tables, accepted for publication in ApJ

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs., Comment: To appear in Science Magazine on February 15th; contact authors: Stefan Funk, Takaaki Tanaka, Yasunobu Uchiyama

We report the detection of pulsed gamma-ray emission from the fast millisecond pulsars (MSPs) B1937+21 (also known as J1939+2134) and B1957+20 (J1959+2048) using 18 months of survey data recorded by the \emph{Fermi} Large Area Telescope (LAT) and timing solutions based on radio observations conducted at the Westerbork and Nan\c{c}ay radio telescopes. In addition, we analyzed archival \emph{RXTE} and \emph{XMM-Newton} X-ray data for the two MSPs, confirming the X-ray emission properties of PSR B1937+21 and finding evidence ($\sim 4\sigma$) for pulsed emission from PSR B1957+20 for the first time. In both cases the gamma-ray emission profile is characterized by two peaks separated by half a rotation and are in close alignment with components observed in radio and X-rays. These two pulsars join PSRs J0034-0534 and J2214+3000 to form an emerging class of gamma-ray MSPs with phase-aligned peaks in different energy bands. The modeling of the radio and gamma-ray emission profiles suggests co-located emission regions in the outer magnetosphere., Comment: Accepted for publication in the Astrophysical Journal

Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power pulsar (Edot = 3.1 * 10^36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise typical of very young pulsars. In the very-high-energy domain, H.E.S.S. has reported the detection of the extended source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1.4 GHz, we phase-fold more than two years of gamma-ray data acquired by Fermi-LAT as well as those collected with XMM-Newton, and perform gamma-ray spectral modeling. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 1.5 +/- 0.3stat +/- 0.3syst with an exponential cut-off at 0.8 +/- 0.3stat +/- 0.3syst GeV and an integral photon flux above 100 MeV of (6.5 +/- 1.6stat +/- 2.3syst) * 10^-8 cm^-2 s^-1. The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin (2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of L_gamma = (2.13 +/- 0.25stat +/- 0.83syst) * 10^34 erg/s, consistent with an L_gamma \propto sqrt(Edot) relationship. The Fermi non-detection of the pulsar wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV emission., Comment: Accepted for publication by A&A; 7 pages, 5 figures

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively., Comment: 17 pages, 6 figures, Accepted for publication in ApJ

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced observations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 0.13+-0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+1.3-1.1 kpc. We detail our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. Based on our new model-independent estimates for distance and transverse velocity, we investigate the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distribution. We find that young pulsars are moving away from the Galactic plane, as expected, and that age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable exceptions. Given its present trajectory, the pulsar B2045-16 was plausibly born in the open cluster NGC 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the parallax distances by a factor of two, while in others the estimates agree with or are larger than the parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The VLBA astrometric results for the recycled pulsar J1713+0747 are consistent with two independent estimates from pulse timing, enabling a consistency check between the different reference frames., Comment: 16 pages, 9 figures, 4 tables; results unchanged; revised version accepted by ApJ

We describe a comprehensive pulsar monitoring campaign for the Large Area Telescope (LAT) on the {\em Fermi Gamma-ray Space Telescope} (formerly GLAST). The detection and study of pulsars in gamma rays give insights into the populations of neutron stars and supernova rates in the Galaxy, into particle acceleration mechanisms in neutron star magnetospheres, and into the "engines" driving pulsar wind nebulae. LAT's unprecedented sensitivity between 20 MeV and 300 GeV together with its 2.4 sr field-of-view makes detection of many gamma-ray pulsars likely, justifying the monitoring of over two hundred pulsars with large spin-down powers. To search for gamma-ray pulsations from most of these pulsars requires a set of phase-connected timing solutions spanning a year or more to properly align the sparse photon arrival times. We describe the choice of pulsars and the instruments involved in the campaign. Attention is paid to verifications of the LAT pulsar software, using for example giant radio pulses from the Crab and from PSR B1937+21 recorded at Nan\c{c}ay, and using X-ray data on PSR J0218+4232 from XMM-Newton. We demonstrate accuracy of the pulsar phase calculations at the microsecond level. Data Table 1 is only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ ., Comment: 13 pages, 6 figures. Accepted for publication in Astronomy & Astrophysics

We investigate the age constraints that can be placed on the double pulsar system using models for the spin-down of the first-born 22.7-ms pulsar A and the 2.77-s pulsar B with characteristic ages of 210 and 50 Myr respectively. Standard models assuming dipolar spin-down of both pulsars suggest that the time since the formation of B is ~50 Myr, i.e. close to B's characteristic age. However, adopting models which account for the impact of A's relativistic wind on B's spin-down we find that the formation of B took place either 80 or 180 Myr ago, depending the interaction mechanism. Formation 80 Myr ago, closer to B's characteristic age, would result in the contribution from J0737-3039 to the inferred coalescence rates for double neutron star binaries increasing by 40%. The 180 Myr age is closer to A's characteristic age and would be consistent with the most recent estimates of the coalescence rate. The new age constraints do not significantly impact recent estimates of the kick velocity, tilt angle between pre and post-supernova orbital planes or pre-supernova mass of B's progenitor., Comment: 5 pages, 2 figures, 1 table, accepted for publication by MNRAS

Recent timing observations of the double pulsar J0737-3039A/B have shown that its transverse velocity is extremely low, only 10 km/s, and nearly in the Plane of the Galaxy. With this new information, we rigorously re-examine the history and formation of this system, determining estimates of the pre-supernova companion mass, supernova kick and misalignment angle between the pre- and post-supernova orbital planes. We find that the progenitor to the recently formed `B' pulsar was probably less than 2 MSun, lending credence to suggestions that this object may not have formed in a normal supernova involving the collapse of an iron core. At the same time, the supernova kick was likely non-zero. A comparison to the history of the double-neutron-star binary B1534+12 suggests a range of possible parameters for the progenitors of these systems, which should be taken into account in future binary population syntheses and in predictions of the rate and spatial distribution of short gamma-ray burst events., Comment: To appear in MNRAS Letters. Title typo fix only; no change to paper

When massive stars exhaust their fuel they collapse and often produce the extraordinarily bright explosions known as core-collapse supernovae. On occasion, this stellar collapse also powers an even more brilliant relativistic explosion known as a long-duration gamma-ray burst. One would then expect that long gamma-ray bursts and core-collapse supernovae should be found in similar galactic environments. Here we show that this expectation is wrong. We find that the long gamma-ray bursts are far more concentrated on the very brightest regions of their host galaxies than are the core-collapse supernovae. Furthermore, the host galaxies of the long gamma-ray bursts are significantly fainter and more irregular than the hosts of the core-collapse supernovae. Together these results suggest that long-duration gamma-ray bursts are associated with the most massive stars and may be restricted to galaxies of limited chemical evolution. Our results directly imply that long gamma-ray bursts are relatively rare in galaxies such as our own Milky Way., Comment: 27 pages, 4 figures, submitted to Nature on 22 August 2005, revised 9 February 2006, online publication 10 May 2006. Supplementary material referred to in the text can be found at http://www.stsci.edu/~fruchter/GRB/locations/supplement.pdf . This new version contains minor changes to match the final published version

The highest velocity neutron stars establish stringent constraints on natal kicks, asymmetries in supernova core collapse, and the evolution of close binary systems. Here we present the first results of a long-term pulsar astrometry program using the VLBA. We measure a proper motion and parallax for the pulsar B1508+55, leading to model-independent estimates of its distance (2.37+0.23-0.20 kpc) and transverse velocity (1083+103-90 km/s), the highest velocity directly measured for a neutron star. We trace the pulsar back from its present Galactic latitude of 52.3 degrees to a birth site in the Galactic plane near the Cyg OB associations, and find that it will inevitably escape the Galaxy. Binary disruption alone is insufficient to impart the required birth velocity, and a natal kick is indicated. A composite scenario including a large kick along with binary disruption can plausibly account for the high velocity., Comment: 5 pages, including 2 figures; accepted by ApJL; associated NRAO press release at http://www.nrao.edu/pr/2005/fastpulsar/

Determination of pulsar parallaxes and proper motions addresses fundamental astrophysical questions. We have recently finished a VLBI astrometry project to determine the proper motions and parallaxes of 27 pulsars, thereby doubling the total number of pulsar parallaxes. Here we summarise our astrometric technique and present the discovery of a pulsar moving in excess of 1000 km/s. As an example of the application of high precision pulsar astrometry we also infer the identification of 2 pulsars originating from a disrupted binary in the Cygnus Superbubble., Comment: 3 pages, 2 figures, to appear in the proceedings of the "Stellar End Products" workshop, 13-15 April 2005, Granada, Spain (for publication in MmSAI vol.77)

We study the internal structural properties of 37 long duration gamma-ray burst (GRB) host galaxies imaged with the Hubble Space Telescope. Fitting exponential disk and r^1/4 models to the surface brightness profiles of eight z < 1.2 bright host galaxies, we find that the disk model is slightly preferred for most hosts, although two galaxies are fit best with an r^1/4 profile. We furthermore measure the central concentrations and asymmetries of all 37 host galaxies using the CAS system, and compare with field galaxies in the Hubble Deep Field. Our first main conclusion is that GRB hosts exhibit a surprisingly high diversity of galaxy types. A significant fraction (68%) of host galaxies are situated in a region of the concentration-asymmetry diagram occupied by spirals or peculiar/merging galaxies. Twelve hosts (32%) are situated in the region occupied by elliptical galaxies, having high concentration indices. We also find some evidence for evolution in GRB host galaxy morphology, such that hosts at z > 1 have a relatively high light concentration, indicating that these systems are perhaps progenitors of massive galaxies, or are compact blue star forming galaxies. We find that GRB hosts at z > 1 are different from the general field population at z > 1 in terms of light concentration at >99.5% confidence, yet have sizes similar to the general z > 1 population. This is the opposite of the effect at z < 1 where GRB hosts are smaller than average. We argue that GRB hosts trace the starburst population at high redshift, as similarly concentrated galaxies at z > 1 are undergoing a disproportionate amount of star formation for their luminosities. Our results show that GRBs are perhaps ideal tracers of typical galaxies undergoing star formation at any epoch., Comment: Astrophysical Journal, in press, Oct. 10th, 2005 issue

(Abridged) We report on HST observations of the late-time afterglow and host galaxy of GRB021004 (z=2.33). Although this GRB is one of the best observed so far, there is substantial disagreement between different interpretations of data sets on this burst in the literature. We have observed the field of GRB021004 with the HST at multiple epochs from 3 days until almost 10 months after the burst. With STIS spectroscopy we cover the spectral region 600-1700AA in the restframe. From the limit on the flux recovery bluewards of the Lyman-limit we constrain the HI column density to be above 1x10^18 cm-2 Based on ACS and NICMOS imaging we find that the afterglow evolved achromatically within the errors. The color changes observed by other authors during the first four days must be related to a stochastic phenomenon superimposed on an afterglow component with a constant spectral shape. This achromaticity implies that the cooling break has remained on the blue side of the optical part of the spectrum for at least two weeks after the burst. The optical-to-X-ray slope is consistent with being the same at 1.4 and 52.4 days after the burst. This indicates that the cooling frequency is constant and hence, according to fireball models, that the circumburst medium has a constant density profile. The late-time slope of the lightcurve is in the range alpha_2 = 1.8-1.9, and is inconsistent with a single power-law. This could be due to a late-time flattening caused by the transition to non-relativistic expansion or due to excess emission (a lightcurve `bump') 7 days after burst. The host is a starburst galaxy with no evidence for dust and with strong Ly-alpha emission. The afterglow was located very close (~100 pc) to the center of the host suggesting that the progenitor was associated with a circumnuclear starburst., Comment: 8 figures, to appear in the Astrophysical Journal

We have recently measured the angle between the spin and orbital angular momenta of PSR B1534+12 to be either 25+/-4 deg or 155+/-4 deg. This misalignment was almost certainly caused by an asymmetry in the supernova explosion that formed its companion neutron star. Here we combine the misalignment measurement with measurements of the pulsar and companion masses, the orbital elements, proper motion, and interstellar scintillation. We show that the orbit of the binary in the Galaxy is inconsistent with a velocity kick large enough to produce a nearly antialigned spin axis, so the true misalignment must be ~25 deg. Similar arguments lead to bounds on the mass of the companion star immediately before its supernova: 3+/-1 Msun. The result is a coherent scenario for the formation of the observed binary. After the first supernova explosion, the neutron star that would eventually become the observed pulsar was in a Be/X-ray type binary system with a companion of at least 10--12 Msun. During hydrogen (or possibly helium) shell burning, mass transfer occurred in a common envelope phase, leaving the neutron star in a roughly half-day orbit with a helium star with mass above ~3.3 Msun. A second phase of mass transfer was then initiated by Roche lobe overflow during shell helium burning, further reducing both the helium star mass and orbital period before the second supernova. Scenarios that avoid Roche lobe overflow by the helium star require larger helium star masses and predict space velocities inconsistent with our measurements. The companion neutron star experienced a velocity kick of 230+/-60 km/s at birth, leading to a systemic kick to the binary of 180+/-60 km/s., Comment: 9 pages, submitted to ApJ. Abstract shortened. Version with high-resolution figures available at http://www.astro.ubc.ca/people/stairs/papers/tds04_orig.ps.gz

In relativistic gravity, a spinning pulsar will precess as it orbits a compact companion star. We have measured the effect of such precession on the average shape and polarization of the radiation from PSR B1534+12. We have also detected, with limited precision, special-relativistic aberration of the revolving pulsar beam due to orbital motion. Our observations fix the system geometry, including the misalignment between the spin and orbital angular momenta, and yield a measurement of the precession timescale consistent with the predictions of General Relativity., Comment: 4 pages, accepted to PRL. Version with high-resolution figure 2 available at http://www.astro.ubc.ca/people/stairs/papers/sta04b.ps.gz

We present the discovery and monitoring of the optical transient (OT) associated with GRB 020410. The fading OT was found by Hubble Space Telescope (HST) observations taken 28 and 65 days after burst at a position consistent with the X-ray afterglow. Subsequent re-examination of early ground based observations revealed that a faint OT was present 6 hours after burst, confirming the source association with GRB 020410. A deep non-detection after one week requires that the OT re-brightened between day 7 and day 28, and further late time HST data taken approximately 100 days after burst imply that it is very red.We compare both the flux and color of the excess with supernova models and show that the data are best explained by the presence of a Type Ib/c supernova at a redshift z ~ 0.5, which occured roughly coincident with the day of GRB., Comment: 23 Pages, 9 figures, submitted to ApJ

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 Journal

The Monogem ring is a bright, diffuse, 25-degree-diameter supernova remnant easily visible in soft X-ray images of the sky. Projected within the ring is a young radio pulsar, PSR B0656+14. An association between the remnant and pulsar has been considered, but was seemingly ruled out by the direction and magnitude of the pulsar proper motion and by a distance estimate that placed the pulsar twice as far from Earth as the remnant. Here we show that in fact the pulsar was born very close to the center of the expanding remnant, both in distance and projection. The inferred pulsar and remnant ages are in good agreement. The conclusion that the pulsar and remnant were born in the same supernova explosion is nearly inescapable. The remnant distance and age are in remarkable concordance with the predictions of a model for the primary cosmic ray energy spectrum in which the `knee' feature is produced by a single dominant source., Comment: 4 pages, to appear in the Astrophys. J. Lett. Full size color figure can be found at http://www.thorsett.org/research

We present the result of astrometric observations of the radio pulsar PSR B0656+14, made using the Very Long Baseline Array. The parallax of the pulsar is pi = 3.47 +- 0.36 mas, yielding a distance of 288 +33 -27 pc. This independent distance estimate has been used to constrain existing models of thermal x-ray emission from the neutron star's photosphere. Simple blackbody fits to the x-ray data formally yield a neutron star radius R_inf ~ 7-8.5 km. With more realistic fits to a magnetized hydrogen atmosphere, any radius between ~13 and ~20 km is allowed., Comment: 7 pages including 1 figure. Submitted to ApJL. AAStex

We have continued our long term study of the double-neutron-star binary pulsar PSR B1534+12, using new instrumentation to make very high precision measurements at the Arecibo Observatory. We have significantly improved our solution for the astrometric, spin, and orbital parameters of the system, as well as for the five "post-Keplerian" orbital parameters that can be used to test gravitation theory. The results are in good agreement with the predictions of general relativity. With the assumption that general relativity is the correct theory of gravity in the classical regime, our measurements allow us to determine the masses of the pulsar and its companion neutron star with high accuracy: 1.3332 \pm 0.0010 Msun and 1.3452 \pm 0.0010Msun, respectively. The small but significant mass difference is difficult to understand in most evolutionary models, as the pulsar is thought to have been born first from a more massive progenitor star and then undergone a period of mass accretion before the formation of the second neutron star. PSR B1534+12 has also become a valuable probe of the local interstellar medium. We have now measured the pulsar distance to be 1.02 \pm 0.05 kpc, giving a mean electron density along this line of sight of 0.011/cc. We continue to measure a gradient in the dispersion measure, though the rate of change is now slower than in the first years after the pulsar's discovery., Comment: 23 pages, 6 figures, to be published in ApJ

We have made the first direct interferometric proper motion measurements of the radio pulsar B1757-24, which sits at the tip of the "beak" of the putative "Duck" supernova remnant. The peculiar morphology of this radio complex has been used to argue alternately that the pulsar's space motion was either surprisingly high or surprisingly low. In fact, we show that the pulsar's motion is so small that it and its associated nonthermal nebula G5.27-0.9 (the "head") are almost certainly unrelated to the much larger G5.4-1.2 (the "wings")., Comment: 10 pages, 1 figure, accepted for publication in ApJL

We determined the distances to nine pulsars by parallax measurements using the NRAO Very Long Baseline Array, doubling the number of pulsars with accurate distance measurements. Broadband phase modeling was used to calibrate the varying dispersive effects of the ionosphere and remove the resulting phase errors from the phase-referenced VLBI data. The resulting parallaxes have a typical accuracy of 100 microarcseconds or better, yielding distances measurements as accurate as 2%. We also report new proper motion measurements of these pulsars, accurate to 0.4 mas per year or better., Comment: 36 pages, Latex, 15 figures. Accepted for publication in ApJ, 2002 June 1

The pulsar velocity distribution provides information about the binary history of pulsar progenitors as well as asymmetries of the supernova events in which pulsars are born. Studies of local pulsars present a biased view of this distribution, since they preferentially select low velocity pulsars that have remained near their birthplaces in the Galactic plane. Using the VLA, we have studied the proper motions of a large sample of distant pulsars. These pulsars are generally faint, and the expected proper motions are small. In this paper, we describe the data analysis techniques that we have developed to allow precise astrometric measurements of faint sources with the VLA. These techniques include ``gating'' the VLA correlator to increase the signal-to-noise ratio of the pulsar by gathering data only during the pulse. Wide-field imaging techniques, including multiband imaging to account for bandwidth smearing, were used to allow identification of multiple in-beam background sources for astrometric calibration. We present the analysis of three pulsars, and demonstrate that astrometric accuracy of about ten milliarcseconds can be obtained for individual sources with our technique, allowing measurement of proper motions with errors of only a few milliarcseconds per year over our seven year baseline., Comment: Accepted for publication in Astronomical Journal, 19 pages, 5 figures

We present HST/STIS observations of the optical counterpart (OT) of the gamma-ray burster GRB 000301C obtained on 2000 March 6, five days after the burst. CCD clear aperture imaging reveals a R ~ 21.50+/-0.15 source with no apparent host galaxy. An 8000 s, 1150 < lambda/A < 3300 NUV-MAMA prism spectrum shows a relatively flat continuum (in f_lambda) between 2800 and 3300 A, with a mean flux 8.7 (+0.8,-1.6)+/- 2.6 10^(-18) ergs/s/cm^2/A, and a sharp break centered at 2797+/-25 A. We interpret it as HI Lyman break at z = 2.067+/-0.025 indicating the presence of a cloud with a HI column density log(HI) > 18 on the line-of-sight to the OT. This value is conservatively a lower limit to the GRB redshift. However, the facts that large N(HI) system are usually considered as progenitors of present day galaxies and that other OTs are found associated with star forming galaxies strongly suggest that it is the GRB redshift. In any case, this represents the largest direct redshift determination of a gamma-ray burster to date. Our data are compatible with an OT spectrum represented by a power-law with an intrinsic index \alpha = 1.2((f_nu \propto nu^-alpha) and no extinction in the host galaxy or with alpha = 0.5 and extinction by a SMC-like dust in the OT rest-frame with A_V = 0.15. The large N(HI) and the lack of detected host is similar to the situation for damped Ly-alpha absorbers at z > 2., Comment: Replaced by final version. 10 p., 2 fig. Scheduled to appear in ApJ 555 n2 Jul 10, 2001. Minor changes, both redshift and mean near UV flux are revised with slightly larger values, due to a wrong offset sign in the wavelength calibration

A new technique has been developed to remove the ionosphere's distorting effects from low frequency VLBI data. By fitting dispersive and non-dispersive components to the phases of multi-frequency data, the ionosphere can be effectively removed from the data without the use of {\em a priori} calibration information. This technique, along with the new gating capability of the VLBA correlator, was used to perform accurate astrometry on pulsar B0950+08, resulting in a much improved measurement of this pulsar's proper motion ($\mu_{\alpha} = -1.6 \pm 0.4$ mas/yr, $\mu_{\delta} = 29.5 \pm 0.5$ mas/yr) and parallax ($\pi = 3.6 \pm 0.3$ mas). This puts the pulsar at a distance of $280 \pm 25$ parsecs, about twice as far as previous estimates, but in good agreement with models of the electron density in the local bubble., Comment: 5 pages, Latex with AASTEX. Accepted for publication in ApJ

Digital signal recorders are becoming widely used in several subfields of centimetre-wavelength radio astronomy. We review the benefits and design considerations of such systems and describe the Princeton Mark IV instrument, an implementation designed for coherent-dedispersion pulsar observations. Features of this instrument include corrections for the distortions caused by coarse quantization of the incoming signal, as well algorithms which effectively excise both narrowband and broadband radio-frequency interference. Observations at 430 MHz using the Mark IV system in parallel with a system using a 250 kHz filter bank and incoherent dedispersion demonstrated timing precision improvement by a factor of 3 or better for typical millisecond pulsars., Comment: 20 pages, 8 figures. To appear in MNRAS. Corrected typo in and text near equation 2

We present Arecibo observations of PSR B1534+12 which confirm previous suggestions that the pulse profile is evolving secularly. This effect is similar to that seen in PSR B1913+16, and is almost certainly due to general relativistic precession of the pulsar's spin axis., Comment: 4 pages, 2 figures. Corrected a reference. To appear in proceedings of IAU Colloquium 177: Pulsar Astronomy - 2000 and Beyond

Giant pulses are short, intense outbursts of radio emission with a power-law intensity distribution that have been observed from the Crab Pulsar and PSR B1937+21. We have undertaken a systematic study of giant pulses from PSR B1937+21 using the Arecibo telescope at 430, 1420, and 2380 MHz. At 430 MHz, interstellar scattering broadens giant pulses to durations of $\sim50 \mu$secs, but at higher frequencies the pulses are very short, typically lasting only $\sim1$-$2 \mu$secs. At each frequency, giant pulses are emitted only in narrow ($\lsim10 \mu$s) windows of pulse phase located $\sim 55$-$70 \mu$sec after the main and interpulse peaks. Although some pulse-to-pulse jitter in arrival times is observed, the mean arrival phase appears stable; a timing analysis of the giant pulses yields precision competitive with the best average profile timing studies. We have measured the intensity distribution of the giant pulses, confirming a roughly power-law distribution with approximate index of -1.8, contributing $\gsim0.1%$ to the total flux at each frequency. We also find that the intensity of giant pulses falls off with a slightly steeper power of frequency than the ordinary radio emission., Comment: 21 pages, 10 Postscript figures; LaTeX with aaspp4.sty and epsf.tex; submitted to ApJ

We present a new timing solution for PSR B1534+12, based on coherently-dedispersed observations at Jodrell Bank and, recently, Arecibo. The new data have resulted in improved measurements of the post-Keplerian timing parameters, including the orbital period derivative. At present, the poorly-known distance to the pulsar limits the precision of the measurement of the intrinsic orbital period derivative, and hence the strength of the test of general relativity that results from this binary system. By assuming that general relativity is the correct theory of gravity, we may invert the test and find an improved value of the pulsar distance., Comment: Acknowledgments expanded. 9 pages. Talk given at "Gravitational Waves and Experimental Gravity," XXXIV Rencontres de Moriond, January 23-30, 1999

We report on observations of the field of GRB~970508 made in early August 1998, 454 days after outburst, with the STIS CCD camera onboard the Hubble Space Telescope. The images, taken in open filter (50CCD) mode, clearly reveal the presence of a galaxy which was obscured in earlier (June 1997) HST images by emission from the optical transient (OT). The galaxy is regular in shape: after correcting for the HST/STIS PSF, it is well-fitted by an exponential disk with a scale length of 0."046 +/- 0."006 and an ellipticity of 0.70 +/- 0.07. All observations are marginally consistent with a continuous decline in OT emission as t^{-1.3} beginning two days after outburst; however, we find no direct evidence in the image for emission from the OT, and the surface brightness profile of the galaxy is most regular if we assume that the OT emission is negligible, suggesting that the OT may have faded more rapidly at late times than is predicted by the power-law decay. Due to the wide bandwidth of the STIS clear mode, the estimated magnitude of the galaxy is dependent on the galaxy spectrum that is assumed. Using colors obtained from late-time ground-based observations to constrain the spectrum, we find V = 25.4 +/- 0.15, a few tenths of a magnitude brighter than earlier ground-based estimates that were obtained by observing the total light of the galaxy and the OT and then subtracting the estimated OT brightness assuming it fades as a single power-law. This again suggests that the OT may have faded faster at late time than the power-law predicts. The position of the OT agrees with that of the isophotal center of the galaxy to 0."01. This remarkable agreement raises the possibility that the GRB may have been associated with either an active galactic nucleus or a nuclear starburst., Comment: Submitted to the Astrophysical Journal (Letters). Thirteen pages, three encapsulated figures. Abstract slightly abridged

The millisecond pulsar PSR B1620-26, in the globular cluster M4, has a white dwarf companion in a half-year orbit. Anomalously large variations in the pulsar's apparent spin-down rate have suggested the presence of a second companion in a much wider orbit. Using timing observations made on more than seven hundred days spanning eleven years, we confirm this anomalous timing behavior. We explicitly demonstrate, for the first time, that a timing model consisting of the sum of two non-interacting Keplerian orbits can account for the observed signal. Both circular and elliptical orbits are allowed, although highly eccentric orbits require improbable orbital geometries. The motion of the pulsar in the inner orbit is very nearly a Keplerian ellipse, but the tidal effects of the outer companion cause variations in the orbital elements. We have measured the change in the projected semi-major axis of the orbit, which is dominated by precession-driven changes in the orbital inclination. This measurement, along with limits on the rate of change of other orbital elements, can be used to significantly restrict the properties of the outer orbit. We find that the second companion most likely has a mass m~0.01 Msun --- it is almost certainly below the hydrogen burning limit (m<0.036 Msun, 95% confidence) --- and has a current distance from the binary of ~35 AU and orbital period of order one hundred years. Circular (and near-circular) orbits are allowed only if the pulsar magnetic field is ~3x10^9 G, an order of magnitude higher than a typical millisecond pulsar field strength. In this case, the companion has mass m~1.2x10^-3 Msun and orbital period ~62 years., Comment: 12 pages, 6 figures, 3 tables. Very minor clarifications and rewording. Accepted for publication in the Astrophys. J

We report on HST and Palomar optical images of the field of GRB 990123, obtained on 8 and 9 February 1999. We find that the optical transient (OT) associated with GRB 990123 is located on an irregular galaxy, with magnitude V=24.20 +/- 0.15. The strong metal absorption lines seen in the spectrum of the OT, along with the low probability of a chance superposition, lead us to conclude that this galaxy is the host of the GRB. The OT is projected within the ~1'' visible stellar field of the host, nearer the edge than the center. We cannot, on this basis, rule out the galactic nucleus as the site of the GRB, since the unusual morphology of the host may be the result of an ongoing galactic merger, but our demonstration that this host galaxy has extremely blue optical to infrared colors more strongly supports an association between GRBs and star formation. We find that the OT magnitude on 1999 Feb 9.05, V = 25.45 +/- 0.15, is about 1.5 mag fainter than expected from extrapolation of the decay rate found in earlier observations. A detailed analysis of the OT light curve suggests that its fading has gone through three distinct phases: an early rapid decline (f_{nu} \propto t^{-1.6} for t < 0.1 days), a slower intermediate decline power-law decay (f_{nu} \propto t^{-1.1} for 0.1 < t < 2 days), and then a more rapid decay (at least as steep as (f_{\nu} \propto t^{-1.8} for t > 2 days). The break to steeper slope at late times may provide evidence that the optical emission from this GRB was highly beamed., Comment: Accepted for publication in Astrophysical Journal (Letters). Fourteen pages. Three encapsulated figures

We present a large sample of high-precision, coherently-dedispersed polarization profiles of millisecond pulsars (MSPs) at frequencies between 410 and 1414 MHz. These data include the first polarimetric observations of several of the pulsars, and the first low-frequency polarization profiles for others. Our observations support previous suggestions that the pulse shapes and polarimetry of MSPs are more complex than those of their slower relatives. An immediate conclusion is that polarimetry-based classification schemes proposed for young pulsars are of only limited use when applied to millisecond pulsars., Comment: 28 pages, 10 figures. Text matches version that appeared in ApJS. Full paper with high-resolution figures available at ftp://ftp.jb.man.ac.uk/pub/psr/papers/msppolpton.ps.gz

We report on observations of the fading optical counterpart of the gamma-ray burst GRB 970228, made with the Hubble Space Telescope STIS CCD approximately six months after outburst and with the HST/NICMOS and Keck/NIRC approximately one year after outburst. The unresolved counterpart is detected by STIS at V=28.0 +/- 0.25, consistent with a continued power-law decline with exponent -1.14 +/- 0.05. The counterpart is located within, but near the edge of, a faint extended source with diameter ~0."8 and integrated magnitude V=25.8 +/- 0.25. A reanalysis of HST and NTT observations performed shortly after the burst shows no evidence of proper motion of the point source or fading of the extended emission. Only the extended source is visible in the NICMOS images with a magnitude of H=23.3 +/- 0.1. The Keck observations find K = 22.8 +/- 0.3. Several distinct and independent means of deriving the foreground extinction in the direction of GRB 970228 all agree with A_V = 0.75 +/- 0.2. After adjusting for Galactic extinction, we find that the size of the observed extended emission is consistent with that of galaxies of comparable magnitude found in the Hubble Deep Field (HDF) and other deep HST images. Only 2% of the sky is covered by galaxies of similar or greater surface brightness; therefore the extended source is almost certainly the host galaxy. Additionally, we find that the extinction-corrected V - H and V - K colors of the host are as blue as any galaxy of comparable or brighter magnitude in the HDF. Taken in concert with recent observations of GRB 970508, GRB 971214, and GRB 980703 our work suggests that all four GRBs with spectroscopic identification or deep multicolor broad-band imaging of the host lie in rapidly star-forming galaxies., Comment: 24 pages, Latex, 4 PostScript figures, to appear in the May 10 issue of The Astrophysical Journal (Note: displayed abstract is abridged)

We have tested several models of GRB luminosity and redshift distribution functions for compatibility with the BATSE 4B number versus peak flux relation. Our results disagree with recent claims that current GRB observations can be used to strongly constrain the cosmic star formation history. Instead, we find that relaxing the assumption that GRBs are standard candles renders a very broad range of models consistent with the BATSE number-flux relation. We explicitly construct two sample distributions, one tracing the star formation history and one with a constant comoving density. We show that both distributions are compatible with the observed fluxes and redshifts of the bursts GRB970508, GRB971214, and GRB980703, and we discuss the measurements required to distinguish the two models., Comment: 12 pages, 2 postscript figures, uses AAS LaTex macros v4.0. To be published in Astrophysical Journal Letters, accepted August 20, 1998. Revised for publication

There are now about fifty known radio pulsars in binary systems, including at least five in double neutron star binaries. In some cases, the stellar masses can be directly determined from measurements of relativistic orbital effects. In others, only an indirect or statistical estimate of the masses is possible. We review the general problem of mass measurement in radio pulsar binaries, and critically discuss all current estimates of the masses of radio pulsars and their companions. We find that significant constraints exist on the masses of twenty-one radio pulsars, and on five neutron star companions of radio pulsars. All the measurements are consistent with a remarkably narrow underlying gaussian mass distribution, m = 1.35 +- 0.04 solar masses. There is no evidence that extensive mass accretion (delta m >~ 0.1 solar mass) has occurred in these systems. We also show that the observed inclinations of millisecond pulsar binaries are consistent with a random distribution, and thus find no evidence for either alignment or counteralignment of millisecond pulsar magnetic fields., Comment: 35 pages, 5 figures, submitted to the Astrophys. J. Also available at http://pulsar.princeton.edu/pulsar/publications.shtml

We have made timing observations of binary pulsar PSR B1534+12 with radio telescopes at Arecibo, Green Bank, and Jodrell Bank. By combining our new observations with data collected up to seven years earlier, we obtain a significantly improved solution for the astrometric, spin, and orbital parameters of the system. For the first time in any binary pulsar system, no fewer than five relativistic or "post-Keplerian" orbital parameters are measurable with useful accuracies in a theory-independent way. We find the orbital period of the system to be decreasing at a rate close to that expected from gravitational radiation damping, according to general relativity, although the precision of this test is limited to about 15% by the otherwise poorly known distance to the pulsar. The remaining post-Keplerian parameters are all consistent with one another and all but one of them have fractional accuracies better than 1%. By assuming that general relativity is the correct theory of gravity, at least to the accuracy demanded by this experiment, we find the masses of the pulsar and companion star each to be 1.339+-0.003 Msun and the system's distance to be d = 1.1+-0.2 kpc, marginally larger than the d ~ 0.7 kpc estimated from the dispersion measure. The increased distance reduces estimates of the projected rate of coalescence of double neutron-star systems in the universe, a quantity of considerable interest for experiments with terrestrial gravitational wave detectors such as LIGO., Comment: 17 pages, 4 figures, submitted to the ApJ

We report on Hubble Space Telescope observations of the optical transient discovered in the error box of GRB970508. The object was imaged on 1997 June 2 with the Space Telescope Imaging Spectrograph and Near-Infrared Camera and Multi-Object Spectrometer. The observations reveal a point-like source with R = 23.1 +- 0.2 and H = 20.6 +- 0.3, in agreement with the power-law temporal decay seen in previous ground-based monitoring. Unlike the case of GRB970228, no nebulosity is detected surrounding the OT of GRB970508, although Mg I absorption and [O II] emission seen in Keck spectra at a redshift of z = 0.835 suggest the presence of a dense, star-forming medium. The HST observations set very conservative upper limits of R ~24.5 and H ~22.2 on the brightness of any underlying extended source. If this subtends a substantial fraction of an arcsecond, then the R band limit is ~25.5. Subsequent photometry suggests a flattening of the light curve at later epochs. Assuming the OT decline follows a pure power-law and ascribing the flattening to the presence of an underlying component of constant flux, we find that this must have R = 25.4, consistent with the upper limits determined by HST. At z = 0.8, this would correspond to an absolute magnitude in the U band of ~-18, similar to that of the Large Magellanic Cloud. We propose a scenario in which the host galaxy of the GRB is of Magellanic type, possibly being a ``satellite'' of one of the bright galaxies located at few arcseconds from the optical transient., Comment: 5 pages, 1 figure, to appear in the Proceedings of the Fourth Gamma Ray Bursts Symposium, eds. C. A. Meegan, R. Preece, and T. Koshut

The Chandrasekhar mass limit sets the scale for the late evolutionary stages of massive stars, including the formation of neutron stars in core collapse supernovae. Because its value depends on the gravitational constant G, the masses of these neutron stars retain a record of past values of G. Using Bayesian statistical techniques, I show that measurements of the masses of young and old neutron stars in pulsar binaries limit \dot G/G=(-0.6 \pm 2.0) x 10^-12 yr^-1 (68% confidence) or \dot G/G=(-0.6 \pm 4.2) x 10^-12 yr^-1 (95% confidence)., Comment: RevTeX, 4 pages + 1 figure, to be published in Phys. Rev. Lett

Previous timing data for PSR~B1620$-$26 were consistent with a second companion mass $m_2$ anywhere in the range $\sim10^{-3}-1\,M_\odot$, i.e., from a Jupiter-type planet to a star. We present the latest timing parameters for the system, including a significant change in the projected semi-major axis of the inner binary, a marginal detection of the fourth time derivative of the pulse frequency, and the pulsar proper motion (which is in agreement with published values for the proper motion of M4), and use them to further constrain the mass $m_2$ and the orbital parameters. Using the observed value of $\stackrel{\ldots.}{f}$, we obtain a one-parameter family of solutions, all with $m_2 \lo 10^{-2}\,M_\odot$, i.e., excluding stellar masses. Varying $\stackrel{\ldots.}{f}$ within its formal $1\sigma$ error bar does not affect the mass range significantly. However, if we vary $\stackrel{\ldots.}{f}$ within a $4\sigma$ error bar, we find that stellar-mass solutions are still possible. We also calculate the predicted rate of change of the projected semi-major axis of the inner binary and show that it agrees with the measured value., Comment: uses PASP LaTeX macros, 6 pages, to appear in IAU Colloquium 160, Pulsars: Problems and Progress, also available from http://ensor.mit.edu/~rasio/

The observation by the BATSE instrument on the Compton Gamma Ray Observatory that gamma-ray bursts (GRBs) are distributed isotropically around the Earth but nonuniformly in distance has led to the widespread conclusion that GRBs are most likely to be at cosmological distances, making them the most luminous sources known in the Universe. If bursts arise from events that occur in normal galaxies, such as neutron star binary inspirals, then they will also occur in our Galaxy about every hundred thousand to million years. The gamma-ray flux at the Earth due to a Galactic GRB would far exceed that from even the largest solar flares. The absorption of this radiation in the atmosphere would substantially increase the stratospheric nitric oxide concentration through photodissociation of N$_2$, greatly reducing the ozone concentration for several years through NO$_x$ catalysis, with important biospheric effects due to increased solar ultraviolet flux. A nearby GRB may also leave traces in anomalous radionuclide abundances., Comment: uuencoded, gzip-ed postscript; 6 pages; submitted to ApJ Letters

The measurement or constraint of the masses of neutron stars and their binary companions tests theories of neutron star structure and of pulsar formation and evolution. We have measured the rate of the general relativistic advance of the longitude of periastron for the pulsar PSR B1802$-$07: $\dot\omega=0\fdg060\pm0\fdg009\,\mbox{yr}^{-1}$, which implies a total system mass, pulsar plus companion star, of $M=1.7\pm0.4\,\Msun$. We also present a much improved measurement of the rate of periastron advance for PSR B2303+46: $\dot\omega = 0\fdg0099\pm0\fdg0002\,\mbox{yr}^{-1}$, implying $M=2.53\pm0.08\,\Msun$ for this system. We discuss the available constraints on distribution of mass between the pulsars and their companions, and we compare the pulsar masses with other determinations of neutron star masses., Comment: 10 pages plus 3 figures, latex, accepted ApJ Letters, sets-001

The Gamma-ray Large Area Space Telescope (GLAST) will be launched less than a year from now, and its Large Area Telescope (LAT) is expected to discover scores to hundreds of gamma-ray pulsars. This poster discusses which of the over 1700 known pulsars, mostly visible only at radio Erequencies, are likely to emit greater than l00 MeV gamma rays with intensities detectable by the LAT. The main figure of merit used to select gamma-ray pulsar candidates is sqrt(E-dot)/d^2, where E-dot is the energy loss due to rotational spindown, and d is the distance to the pulsar. The figure of merit incorporates spin-down flux at earth (proportional to E-dot/d^2) times efficiency, assumed proportional to 1/sqrt(E-dot). A few individual objects are cited to illustrate the issues. Since large E-dot pulsars also tend to have large timing noise and occasional glitches, their ephemerides can become inaccurate in weeks to months. To detect and study the gamma-ray emission the photons must be accurately tagged with the pulse phase. With hours to days between gamma-ray photon arrival times from a pulsar and months to years of LAT exposure needed for good detections, GLAST will need timing measurements throughout the continuous gamma-ray observations. The poster will describe efforts to coordinate pulsar timing of the candidate gamma-ray pulsars.