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.

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

In recent years, the number of pulsars with secure mass measurements has increased to a level that allows us to probe the underlying neutron star (NS) mass distribution in detail. We critically review the radio pulsar mass measurements. For the first time, we are able to analyze a sizable population of NSs with a flexible modeling approach that can effectively accommodate a skewed underlying distribution and asymmetric measurement errors. We find that NSs that have evolved through different evolutionary paths reflect distinctive signatures through dissimilar distribution peak and mass cutoff values. NSs in double neutron star and neutron star-white dwarf systems show consistent respective peaks at 1.33 Msun and 1.55 Msun suggesting significant mass accretion (delta m~0.22 Msun) has occurred during the spin-up phase. The width of the mass distribution implied by double NS systems is indicative of a tight initial mass function while the inferred mass range is significantly wider for NSs that have gone through recycling. We find a mass cutoff at ~2.1 Msun for NSs with white dwarf companions which establishes a firm lower bound for the maximum NS mass. This rules out the majority of strange quark and soft equation of state models as viable configurations for NS matter. The lack of truncation close to the maximum mass cutoff along with the skewed nature of the inferred mass distribution both enforce the suggestion that the 2.1 Msun limit is set by evolutionary constraints rather than nuclear physics or general relativity, and the existence of rare super-massive NSs is possible.

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.

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

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.

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.

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.

In recent years, the number of pulsars with secure mass measurements has increased to a level that allows us to probe the underlying neutron star mass distribution in detail. We critically review radio pulsar mass measurements and present a detailed examination through which we are able to put stringent constraints on the underlying neutron star mass distribution. For the first time, we are able to analyze a sizable population of neutron star-white dwarf systems in addition to double neutron star systems with a technique that accounts for systematically different measurement errors. We find that neutron stars that have evolved through different evolutionary paths reflect distinctive signatures through dissimilar distribution peak and mass cutoff values. Neutron stars in double neutron star and neutron star-white dwarf systems show consistent respective peaks at 1.35 Msun and 1.50 Msun which suggest significant mass accretion (Delta m~0.15 Msun) has occurred during the spin up phase. The width of the mass distribution implied by double neutron star systems is indicative of a tight initial mass function while the inferred mass range is significantly wider for neutron stars that have gone through recycling. We find a mass cutoff at 2 Msun for neutron stars with white dwarf companions which establishes a firm lower bound for the maximum neutron star mass. This rules out the majority of strange quark and soft equation of state models as viable configurations for neutron star matter. The lack of truncation close to the maximum mass cutoff suggests that the 2 Msun limit is set by evolutionary constraints rather than nuclear physics or general relativity, and the existence of rare super-massive neutron stars is possible.

In the absence of constraints from the binary companion or supernova remnant, the standard method for estimating pulsar ages is to infer an age from the rate of spin-down. While the generic spin-down age may give realistic estimates for normal pulsars, it can fail for pulsars with very short periods. Details of the spin-up process during the low mass X-ray binary phase pose additional constraints on the period (P) and spin-down rates (Pdot) that may consequently affect the age estimate. Here, we propose a new recipe to estimate millisecond pulsar (MSP) ages that parametrically incorporates constraints arising from binary evolution and limiting physics. We show that the standard method can be improved by this approach to achieve age estimates closer to the true age whilst the standard spin-down age may over- or under-estimate the age of the pulsar by more than a factor of ~10 in the millisecond regime. We use this approach to analyze the population on a broader scale. For instance, in order to understand the dominant energy loss mechanism after the onset of radio emission, we test for a range of plausible braking indices. We find that a braking index of n=3 is consistent with the observed MSP population. We demonstrate the existence and quantify the potential contributions of two main sources of age corruption: the previously known "age bias" due to secular acceleration and "age contamination" driven by sub-Eddington progenitor accretion rates. We explicitly show that descendants of LMXBs that have accreted at very low rates will exhibit ages that appear older than the age of the Galaxy. We further elaborate on this technique, the implications and potential solutions it offers regarding MSP evolution, the underlying age distribution and the post-accretion energy loss mechanism.

The discovery of millisecond pulsations from neutron stars in low mass X-ray binary (LMXB) systems has substantiated the theoretical prediction that links millisecond radio pulsars (MSRPs) and LMXBs. Since then, the process that produces millisecond radio pulsars from LMXBs, followed by spin-down due to dipole radiation has been conceived as the 'standard evolution' of millisecond pulsars. However, the question whether all the observed millisecond radio pulsars could be produced by LMXBs has not been quantitatively addressed until now. The standard evolutionary process produces millisecond pulsars with periods (P) and spin-down rates (Pdot) that are not entirely independent. The possible P-Pdot values that millisecond radio pulsars can attain are jointly constrained. In order to test whether the observed millisecond radio pulsars are the unequivocal descendants of millisecond X-ray pulsars (MSXP), we have produced a probability map that represents the expected distribution of millisecond radio pulsars for the standard model. We show with more than 95 % confidence that the fastest spinning millisecond radio pulsars with high magnetic fields, e.g. PSR B1937+21, cannot be produced by the observed millisecond X-ray pulsars within the framework of the standard model.

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.

We calculate the joint period-spindown (P-Pdot) distributions of millisecond radio pulsars (MSRP) for the standard evolutionary model in order to test whether the observed MSRPs are the unequivocal descendants of millisecond X-ray pulsars (MSXP). The P-Pdot densities implied by the standard evolutionary model compared with observations suggest that there is a statistically significant overabundance of young/high magnetic field MSRPs. Taking biases due to observational selection effects into account, it is unlikely that MSRPs have evolved from a single coherent progenitor population that loses energy via magnetic dipole radiation after the onset of radio emission. By producing the P-Pdot probability map, we show with more than 95% confidence that the fastest spinning millisecond pulsars with high magnetic fields, e.g. PSR B1937+21, cannot be produced by the observed MSXPs within the framework of the standard model.

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/ .

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.

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.

We present near-infrared (nIR) and optical observations of the afterglow of GRB 030115. Discovered in an infrared search at Kitt Peak 5 hours after the burst trigger, this afterglow is amongst the faintest observed in the R-band at an early epoch, and exhibits very red colors, with $R-K\approx 6$. The magnitude of the optical afterglow of GRB 030115 is fainter than many upper limits for other bursts, suggesting that without early nIR observations it would have been classified as a ``dark'' burst. Both the color and optical magnitude of the afterglow are likely due to dust extinction and indicate that at least some optical afterglows are very faint due to dust along the line of sight. Multicolor {\it Hubble Space Telescope} observations were also taken of the host galaxy and the surrounding field. Photometric redshifts imply that the host, and a substantial number of faint galaxies in the field are at $z \sim 2.5$. The overdensity of galaxies is sufficiently great that GRB 030115 may have occurred in a rich high-redshift cluster. The host galaxy shows extremely red colors (R-K=5) and is the first GRB host to be classified as an Extremely Red Object (ERO). Some of the galaxies surrounding the host also show very red colors, while the majority of the cluster are much bluer, indicating ongoing unobscured star formation. As it is thought that much of high redshift star formation occurs in highly obscured environments it may be that GRB 030115 represent a transition object, between the relatively unobscured afterglows seen to date and a population which are very heavily extinguished, even in the nIR.

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.

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.

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.

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.

(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.

X-ray Flashes (XRFs) are, like Gamma-Ray Bursts (GRBs), thought to signal the collapse of massive stars in distant galaxies. Many models posit that the isotropic equivalent energies of XRFs are lower than those for GRBs, such that they are visible from a reduced range of distances when compared with GRBs. Here we present the results of two epoch Hubble Space Telescope imaging of two XRFs. These images taken approximately 45 and 200 days post burst reveal no evidence for an associated supernova in either case. Supernovae such as SN 1998bw would have been visible out to z ~1.5 in each case, while fainter supernovae such as SN 2002ap would have been visible to z ~ 1. If the XRFs lie at such large distances, their energies would not fit the observed correlation between the GRB peak energy and isotropic energy release, in which soft bursts are less energetic. We conclude that, should these XRFs reside at low redshifts ($z<0.6$), either their line of sight is heavily extinguished, or they are associated with extremely faint supernovae, or, unlike GRBs, these XRFs do not have temporally coincident supernovae.

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.

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.

The combination of ever more precise radio timing data and serendipitous {\it HST} data has confirmed that the outer companion to PSR B1620$-$26 is a planet. Here we summarize the observational situation, including preliminary new timing solutions and the implications of the measured system parameters. We detail the proposed formation scenarios, discussing the advantages and problems of each for explaining the origin of the triple, and we speculate on some of the implications for planet formation in the early universe. Future data on this system will provide additional constraints on fundamental modes of planet formation. We predict that many more exchanged planets will be discovered orbiting recycled pulsars in globular clusters as the sensitivity and duration of radio timing increases. Strong observational tests of some of the alternative formation models should be possible with additional data.

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.

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.

The pulsar B1620-26 has two companions, one of stellar mass and one of planetary mass. We detected the stellar companion using Hubble Space Telescope observations. The color and magnitude of the stellar companion indicate that it is an undermassive white dwarf ($0.34 \pm 0.04 M_{\odot}$) of age $480\pm 140$ Myr. This places a constraint on the recent history of this triple system and supports a scenario in which the current configuration arose through a dynamical exchange interaction in the cluster core. This implies that planets may be relatively common in low-metallicity globular clusters and that planet formation is more widespread and happened earlier than previously believed.

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.

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.

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.

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").

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.

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.

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.

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.

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.

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.

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.

We present results of long-term timing of eclipsing binaries PSR B1744-24A and PSR B1957+20 at Arecibo, the VLA, and Green Bank. Both pulsars exhibit irregularities in pulsar rotation and orbital motion. Increases and decreases of the orbital period of PSR B1957+20 are of order Delta-P_b/P_b ~ 10^-7, varying on a time scale of a few years. Over a decade of observations, the orbital period of PSR B1744-24A has only decreased, with time scale |P_b/P_b-dot| ~ 200 Myr. When the effects of orbital motion are removed from the timing data, long-term trends remain in the pulse phase residuals, with amplitudes of order 30 and 500 us, respectively, for B1957+20 and B1744-24A. Such large "timing noise" is not seen in other spun-up pulsars (isolated or binary), leading us to conclude that it is a consequence of mass flow in the system. Possible causes include variations in the rotation of the pulsars and movement of the binary systems along the line of sight (perhaps due to gravitational interactions with outflowing matter).

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.

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.

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.

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.

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.

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.

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.

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.

We report on observations of the fading optical counterpart of the gamma-ray burst GRB 970228, made on 4~September~1997 using the STIS CCD on the Hubble Space Telescope. The unresolved counterpart is detected at V=28 +/- 0.25, consistent with a continued power-law decline with exponent -1.14 +/- 0.05. No proper motion is detected, in contradiction of an earlier claim. The counterpart is located within, but near the edge of, a faint extended source with diameter ~0."8 and integrated magnitude 25.7 +/- 0.25. Comparison with WFPC2 data taken one month after the initial burst and NTT data taken on March 13 shows no evidence for fading of the extended emission. After adjusting for the probable Galactic extinction in the direction of GRB 970228 of A_v=0.7, we find that the observed nebula is consistent with the sizes of galaxies of comparable magnitude found in the Hubble Deep Field and other deep HST images, and that only 2% of the sky is covered by galaxies of similar or greater surface brightness. We therefore conclude that the extended source observed about GRB 970228 is most likely a galaxy at moderate redshift, and is almost certainly the host of the gamma-ray burst.

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.

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.

We report on Hubble Space Telescope (HST) observations of the optical transient (OT) discovered in the error box of the gamma-ray burst GRB970508. The object was imaged on 1997 June 2 with the Space Telescope Imaging Spectrograph (STIS) and Near-Infrared Camera and Multi-Object Spectrometer (NICMOS). 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 ground-based monitoring. Unlike the case of GRB970228, no nebulosity is detected surrounding the OT of GRB970508. We 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. In combination with Keck spectra that show Mg I absorption and [O II] emission at a redshift of z = 0.835, our observations suggest that the OT is located in a star-forming galaxy with total luminosity one order of magnitude lower than the knee of the galaxy luminosity function, L*. Such galaxies are now thought to harbor the majority of star formation at z ~ 1; therefore, these observations may provide support for a link between GRBs and star formation.

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).

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.

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.

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.