Search

The interstellar medium of the Milky Way contains turbulent plasma with structures driven by energetic processes that fuel star formation and shape the evolution of our Galaxy. Radio waves from pulsars are scattered off the small (au-scale and below) structures, resulting in frequency-dependent interference patterns that are modulated in time because of the relative motions of the pulsar, Earth, and plasma. Power spectral analyses of these patterns show parabolic arcs with curvatures that encode the locations and kinematics of individual structures. Here we report the discovery of at least 25 distinct plasma structures in the direction of the brilliant millisecond pulsar, PSR J0437$-$4715, in observations obtained with the MeerKAT radio telescope. Four arcs reveal structures within 5000 au of the pulsar, from a series of shocks induced as the pulsar and its wind interact with the ambient insterstellar medium. The measured radial distance and velocity of the main shock allows us to solve the shock geometry and space velocity of the pulsar in three dimensions, while the velocity of another structure unexpectedly indicates a back flow from the direction of the shock or pulsar-wind tail. The remaining 21 arcs represent a surprising abundance of structures sustained by turbulence within the Local Bubble -- a region of the interstellar medium thought to be depleted of gas by a series of supernova explosions about 14 Myr ago., Comment: 46 pages, 10 figures, 1 table, submitted to Nature Astronomy

We present GTC-OSIRIS phase-resolved optical spectroscopy of three compact binary MSPs, or `spiders': PSR J1048+2339, PSR J1810+1744, and (for the first time) PSR J1908+2105. For the companion in each system, the temperature is traced throughout its orbit, and radial velocities are measured. The radial velocities are found to vary with the absorption features used when measuring them, resulting in a lower radial velocity curve semi-amplitude measured from the day side of two of the systems when compared to the night: for J1048 ($K_\mathrm{day} = 344 \pm 4$ km s$^{-1}$, $K_\mathrm{night} = 372 \pm 3$ km s$^{-1}$) and, tentatively, for J1810 ($K_\mathrm{day} = 448 \pm 19$ km s$^{-1}$, $K_\mathrm{night} = 491 \pm 32$ km s$^{-1}$). With existing inclination constraints, this gives the neutron star (NS) and companion masses $M_\mathrm{NS} = 1.50 - 2.04$ $M_\odot$ and $M_2 = 0.32 - 0.40$ $M_\odot$ for J1048, and $M_\mathrm{NS} > 1.7$ $M_\odot$ and $M_2 = 0.05 - 0.08$ $M_\odot$ for J1810. For J1908, we find an upper limit of $K_2 < 32$ km s$^{-1}$, which constrains its mass ratio $q = M_2 / M_\mathrm{NS} > 0.55$ and inclination $i < 6.0^\circ$, revealing the previously misunderstood system to be the highest mass ratio, lowest inclination redback yet. This raises questions for the origins of its substantial radio eclipses. Additionally, we find evidence of asymmetric heating in J1048 and J1810, and signs of metal enrichment in J1908. We also explore the impact of inclination on spectroscopic temperatures, and demonstrate that the temperature measured at quadrature ($\phi = 0.25, 0.75$) is essentially independent of inclination, and thus can provide additional constraints on photometric modelling., Comment: Submitted to MNRAS. 18 pages, 19 figures

We report on the results of an image-based search for pulsar candidates toward the Galactic bulge. We used mosaic images from the MeerKAT radio telescope, that were taken as part of a 173 deg**2 survey of the bulge and Galactic center of our Galaxy at L band (856-1712 MHz) in all four Stokes I, Q, U and V. The image root-mean-square noise levels of 12-17 uJy/ba represent a significant increase in sensitivity over past image-based pulsar searches. Our primary search criterion was circular polarization, but we used other criteria including linear polarization, in-band spectral index, compactness, variability and multi-wavelength counterparts to select pulsar candidates. We first demonstrate the efficacy of this technique by searching for polarized emission from known pulsars, and comparing our results with measurements from the literature. Our search resulted in a sample of 75 polarized pulsar candidates. Bright stars or young stellar objects were associated with 28 of these sources, including a small sample of highly polarized dwarf stars with pulsar-like steep spectra. Comparing the properties of this sample with the known pulsars, we identified 30 compelling candidates for pulsation follow-up, including two sources with both strong circular and linear polarization. The remaining 17 sources are either pulsars or stars, but we cannot rule out an extragalactic origin or image artifacts among the brighter, flat spectrum objects., Comment: 28 pages, 10 figures, 6 tables. ApJ, in press

In addition to being the most magnetic objects in the known universe, magnetars are the only objects observed to generate fast-radio-burst-like emissions. The formation mechanism of magnetars is still highly debated, and may potentially be probed with the magnetar velocity distribution. We carried out a 3-year-long astrometric campaign on Swift J1818.0-1607 -- the fastest-spinning magnetar, using the Very Long Baseline Array. After applying the phase-calibrating 1D interpolation strategy, we obtained a small proper motion of 8.5 $\mathrm{mas~yr^{-1}}$ magnitude, and a parallax of $0.12\pm0.02$ mas (uncertainties at $1\,\sigma$ confidence throughout the Letter) for Swift J1818.0-1607. The latter is the second magnetar parallax, and is among the smallest neutron star parallaxes ever determined. From the parallax, we derived the distance $9.4^{+2.0}_{-1.6}$ kpc, which locates Swift J1818.0-1607 at the far side of the Galactic central region. Combined with the distance, the small proper motion leads to a transverse peculiar velocity $v_\perp=48^{+50}_{-16}$ $\mathrm{km~s^{-1}}$ -- a new lower limit to magnetar $v_\perp$. Incorporating previous $v_\perp$ estimates of seven other magnetars, we acquired $v_\perp=149^{+132}_{-68}$ $\mathrm{km~s^{-1}}$ for the sample of astrometrically studied magnetars, corresponding to the three-dimensional space velocity $\sim190^{+168}_{-87}$ $\mathrm{km~s^{-1}}$, smaller than the average level of young pulsars. Additionally, we found that the magnetar velocity sample does not follow the unimodal young pulsar velocity distribution reported by Hobbs et al. at $>2\,\sigma$ confidence, while loosely agreeing with more recent bimodal young pulsar velocity distributions derived from relatively small samples of quality astrometric determinations., Comment: 13 pages, 4 figures, 4 tables, accepted for publication in ApJL

We conducted a search for HI emission of the gas-rich galaxies in the Vela region ($260^{\circ} \leq \ell \leq 290^{\circ}, -2^{\circ} \leq b \leq 1^{\circ}$) to explore the Vela Supercluster (VSCL) at $V_\mathrm{hel} \sim 18000$ km s$^{-1}$, largely obscured by Galactic dust. Within the mostly RFI-free band ($250 < V_\mathrm{hel} < 25000$ km s$^{-1}$) of MeerKAT, the analysis focuses on $157$ hexagonally distributed pointings extracted from the SARAO MeerKAT Galactic Plane Survey located in the Vela region (Vela$-$SMGPS). These were combined into 10 contiguous mosaics, covering a ${\sim}90$ deg$^2$ area. Among the $843$ HI detected sources, 39 were previously discovered in the Parkes HIZOA survey ($V_\mathrm{hel} < 12000$ km s$^{-1}$; rms $\sim 6$ mJy beam$^{-1}$). With the improved rms level of the Vela$-$SMGPS, i.e., $0.29 - 0.56$ mJy beam$^{-1}$, our study unveils nearly 12 times more detections (471 candidates) in that same velocity range. We furthermore could identify $187$ galaxy candidates with an HI mass limit reaching $\log (M_{\rm HI}/\rm M_{\odot}) = 9.44$ in the VSCL velocity range $V_\mathrm{hel} \sim 19500 \pm 3500$ km s$^{-1}$. We find indications of two wall-like overdensities that confirm the original suspicion that these walls intersect at low latitudes around longitudes of $\ell \sim 272^{\circ} - 278^{\circ}$. We also find a strong signature most likely associated with the Hydra/Antlia extension and evidence of a previously unknown narrow filament at $V_\mathrm{hel} \sim 12000$ km s$^{-1}$. This paper demonstrates the efficiency of systematic HI surveys with the SKA precursor MeerKAT, even in the most obscured part of the Zone of Avoidance (ZOA)., Comment: Accepted for publication in MNRAS, this version contains the full catalogue of detections

The Local Void is one of the nearest large voids, located at a distance of 23 Mpc. It lies largely behind the Galactic Bulge and is therefore extremely difficult to observe. We use HI 21 cm emission observations from the SARAO MeerKAT Galactic Plane Survey (SMGPS) to study the Local Void and its surroundings over the Galactic longitude range 329$^{\circ}< \ell <$ 55$^{\circ}$, Galactic latitude $|b| <$ 1.5$^{\circ}$, and redshift $cz <$ 7500 km/s. We have detected 291 galaxies to median rms sensitivity of 0.44 mJy per beam per 44 km/s channel. We find 17 galaxies deep inside the Void, 96 at the border of the Void, while the remaining 178 galaxies are in average density environments. The extent of the Void is ~ 58 Mpc. It is severely under-dense for the longitude range 350$^{\circ}< \ell <$ 35$^{\circ}$ up to redshift $z <$ 4500 km/s. The galaxies in the Void tend to have \HI masses that are lower (by approximately 0.25 dex) than their average density counterparts. We find several potential candidates for small groups of galaxies, of which two groups (with 3 members and 5 members) in the Void show signs of filamentary substructure within the Void., Comment: Accepted for publication in MNRAS; 16 pages, 14 figures, Supplementary data are available online at MNRAS

This paper presents the first HI results extracted from the SARAO MeerKAT Galactic Plane Survey (SMGPS) $-$ a narrow strip ($b \sim 3^\circ$) along the southern Milky Way. The primary goal consisted in tracing the Great Attractor (GA) Wall across the innermost Zone of Avoidance. We reduced a segment spanning the longitude range $302^\circ \leq \ell \leq 332^\circ$ for the redshift range $z \leq 0.08$. The superb SMGPS sensitivity (rms = 0.3-0.5 mJy beam$^{-1}$ per 44 kms$^{-1}$ channel) and angular resolution ($\sim$ 31" $\times$ 26") lead to a detection limit of log$(M_{\rm HI}/$M$_\odot) \geq$ 8.5 at the GA distance ($V_{\rm hel} \sim 3500 - 6500$ kms$^{-1}$). A total of 477 galaxy candidates were identified over the full redshift range. A comparison of the few HI detections with counterparts in the literature (mostly HIZOA) found the HI fluxes and other HI parameters to be highly consistent. The continuation of the GA Wall is confirmed through a prominent overdensity of $N = 214$ detections in the GA distance range. At higher latitudes, the wall moves to higher redshifts, supportive of a possible link with the Ophiuchus cluster located behind the Galactic Bulge. This deep interferometric HI survey demonstrates the power of the SMGPS in improving our insight of large-scale structures at these extremely low latitudes, despite the high obscuration and continuum background., Comment: 5 pages, 7 figures, 2 appendices of 12 pages. Journal reference: MNRAS Letters, accepted

Radio emission from magnetars provides a unique probe of the relativistic, magnetized plasma within the near-field environment of these ultra-magnetic neutron stars. The transmitted waves can undergo birefringent and dispersive propagation effects that result in frequency-dependent conversions of linear to circularly polarized radiation and vice-versa, thus necessitating classification when relating the measured polarization to the intrinsic properties of neutron star and fast radio burst (FRB) emission sites. We report the detection of such behavior in 0.7-4 GHz observations of the P = 5.54 s radio magnetar XTE J1810$-$197 following its 2018 outburst. The phenomenon is restricted to a narrow range of pulse phase centered around the magnetic meridian. Its temporal evolution is closely coupled to large-scale variations in magnetic topology that originate from either plastic motion of an active region on the magnetar surface or free precession of the neutron star crust. Our model of the effect deviates from simple theoretical expectations for radio waves propagating through a magnetized plasma. Birefringent self-coupling between the transmitted wave modes, line-of-sight variations in the magnetic field direction and differences in particle charge or energy distributions above the magnetic pole are explored as possible explanations. We discuss potential links between the immediate magneto-ionic environments of magnetars and those of FRB progenitors., Comment: 22 pages (inc. extended data and supplementary materials), 12 figures, 3 tables. Published in Nature Astronomy

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results., Comment: 142 pages. Accepted by the Astrophysical Journal Supplement

We report the radio and high-energy properties of a new outburst from the radio-loud magnetar 1E 1547.0$-$5408. Following the detection of a short burst from the source with Swift-BAT on 2022 April 7, observations by NICER detected an increased flux peaking at $(6.0 \pm 0.4) \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ in the soft X-ray band, falling to the baseline level of $1.7\times10^{-11}$ erg s$^{-1}$ cm$^{-2}$ over a 17-day period. Joint spectroscopic measurements by NICER and NuSTAR indicated no change in the hard non-thermal tail despite the prominent increase in soft X-rays. Observations at radio wavelengths with Murriyang, the 64-m Parkes radio telescope, revealed that the persistent radio emission from the magnetar disappeared at least 22 days prior to the initial Swift-BAT detection and was re-detected two weeks later. Such behavior is unprecedented in a radio-loud magnetar, and may point to an unnoticed slow rise in the high-energy activity prior to the detected short-bursts. Finally, our combined radio and X-ray timing revealed the outburst coincided with a spin-up glitch, where the spin-frequency and spin-down rate increased by $0.2 \pm 0.1$ $\mu$Hz and $(-2.4 \pm 0.1) \times 10^{-12}$ s$^{-2}$ respectively. A linear increase in spin-down rate of $(-2.0 \pm 0.1) \times 10^{-19}$ s$^{-3}$ was also observed over 147 d of post-outburst timing. Our results suggest that the outburst may have been associated with a reconfiguration of the quasi-polar field lines, likely signalling a changing twist, accompanied by spatially broader heating of the surface and a brief quenching of the radio signal, yet without any measurable impact on the hard X-ray properties., Comment: Accepted to ApJ, 18 pages, 7 figures. Fixed typo in abstract and discussion

We present the first 2.5 years of data from the MeerKAT Pulsar Timing Array (MPTA), part of MeerTime, a MeerKAT Large Survey Project. The MPTA aims to precisely measure pulse arrival times from an ensemble of 88 pulsars visible from the Southern Hemisphere, with the goal of contributing to the search, detection and study of nanohertz-frequency gravitational waves as part of the International Pulsar Timing Array. This project makes use of the MeerKAT telescope, and operates with a typical observing cadence of two weeks using the L-band receiver that records data from 856-1712 MHz. We provide a comprehensive description of the observing system, software, and pipelines used and developed for the MeerTime project. The data products made available as part of this data release are from the 78 pulsars that had at least $30$ observations between the start of the MeerTime programme in February 2019 and October 2021. These include both sub-banded and band-averaged arrival times, as well as the initial timing ephemerides, noise models, and the frequency-dependent standard templates (portraits) used to derive pulse arrival times. After accounting for detected noise processes in the data, the frequency-averaged residuals of $67$ of the pulsars achieved a root-mean-square residual precision of $< 1 \mu \rm{s}$. We also present a novel recovery of the clock correction waveform solely from pulsar timing residuals, and an exploration into preliminary findings of interest to the international pulsar timing community. The arrival times, standards and full Stokes parameter calibrated pulsar timing archives are publicly available., Comment: 18 pages, 6 figures

We report the discovery of PSR J1555-2908, a 1.79 ms radio and gamma-ray pulsar in a 5.6 hr binary system with a minimum companion mass of 0.052 $M_\odot$. This fast and energetic ($\dot E = 3 \times 10^{35}$ erg/s) millisecond pulsar was first detected as a gamma-ray point source in Fermi LAT sky survey observations. Guided by a steep spectrum radio point source in the Fermi error region, we performed a search at 820 MHz with the Green Bank Telescope that first discovered the pulsations. The initial radio pulse timing observations provided enough information to seed a search for gamma-ray pulsations in the LAT data, from which we derive a timing solution valid for the full Fermi mission. In addition to the radio and gamma-ray pulsation discovery and timing, we searched for X-ray pulsations using NICER but no significant pulsations were detected. We also obtained time-series r-band photometry that indicates strong heating of the companion star by the pulsar wind. Material blown off the heated companion eclipses the 820 MHz radio pulse during inferior conjunction of the companion for ~10% of the orbit, which is twice the angle subtended by its Roche lobe in an edge-on system., Comment: 15 pages, 6 figures, accepted by ApJ

Radio-loud magnetars display a wide variety of radio-pulse phenomenology seldom seen among the population of rotation-powered pulsars. Spectropolarimetry of the radio pulses from these objects has the potential to place constraints on their magnetic topology and unveil clues about the magnetar radio emission mechanism. Here we report on eight observations of the magnetar Swift J1818.0$-$1607 taken with the Parkes Ultra-Wideband Low receiver covering a wide frequency range from 0.7 to 4 GHz over a period of 5 months. The magnetar exhibits significant temporal profile evolution over this period, including the emergence of a new profile component with an inverted spectrum, two distinct types of radio emission mode switching, detected during two separate observations, and the appearance and disappearance of multiple polarization modes. These various phenomena are likely a result of ongoing reconfiguration of the plasma content and electric currents within the magnetosphere. Geometric fits to the linearly polarized position angle indicate we are viewing the magnetar at an angle of $\sim$99$^{\circ}$ from the spin axis, and its magnetic and rotation axes are misaligned by $\sim$112$^{\circ}$. While conducting these fits, we found the position angle swing had reversed direction on MJD 59062 compared to observations taken 15 days earlier and 12 days later. We speculate this phenomena may be evidence the radio emission from this magnetar originates from magnetic field lines associated with two co-located magnetic poles that are connected by a coronal loop., Comment: 13 pages, 10 figures, 4 tables. Accepted to MNRAS

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have evaluated the impact of differences in data processing or flux determination methods used by different groups, and attempted to mitigate these differences to provide a more uniform dataset. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.92\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.584\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.168\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al. (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle., Comment: Accepted by ApJ. Revised X-ray and up-to-date afterglow measurements are available in ASCII format on GitHub (https://github.com/kmooley/GW170817/) and mirrored on tauceti.caltech.edu (http://www.tauceti.caltech.edu/kunal/gw170817/). The online dataset will be continuously updated as new measurements become available

We present the final sample of 98 detached double white dwarf (WD) binaries found in the Extremely Low Mass (ELM) Survey, a spectroscopic survey targeting <0.3 Msun He-core WDs completed in the Sloan Digital Sky Survey footprint. Over the course of the survey we observed ancillary low mass WD candidates like GD278, which we show is a P=0.19 d double WD binary, as well as candidates that turn out to be field blue straggler/subdwarf A-type stars with luminosities too large to be WDs given their Gaia parallaxes. Here, we define a clean sample of ELM WDs that is complete within our target selection and magnitude range 15

The globular cluster Omega Centauri is the most massive and luminous cluster in the Galaxy. The $\gamma$-ray source FL8Y J1326.7$-$4729 is coincident with the core of the cluster, leading to speculation that hitherto unknown radio pulsars or annihilating dark matter may be present in the cluster core. Here we report on the discovery of five millisecond pulsars in Omega Centauri following observations with the Parkes radio telescope. Four of these pulsars are isolated with spin periods of 4.1, 4.2, 4.6 and 6.8 ms. The fifth has a spin period of 4.8 ms and is in an eclipsing binary system with an orbital period of 2.1 hours. Deep radio continuum images of the cluster centre with the Australian Telescope Compact Array reveal a small population of compact radio sources making it likely that other pulsars await discovery. We consider it highly likely that the millisecond pulsars are the source of the $\gamma$-ray emission. The long-term timing of these pulsars opens up opportunities to explore the dynamics and interstellar medium of the cluster., Comment: 8 pages, 2 figures, accepted for publication in the Astrophysical Journal Letters

We collect new and archival optical observations of nine "black-widow" millisecond pulsar binaries. New measurements include direct imaging with the Keck, Gemini-S, MDM, and LCO 2~m telescopes. This is supplemented by synthesized colors from Keck long-slit spectra. Four black-widow optical companions are presented here for the first time. Together these data provide multicolor photometry covering a large fraction of the orbital phase. We fit these light curves with a direct (photon) heating model using a version of the ICARUS light-curve modeling code. The fits provide distance and fill-factor estimates, inclinations, and heating powers. We compare the heating powers with the observed GeV luminosities, noting that the ratio is sensitive to pulsar distance and to the gamma-ray beaming. We make a specific correction for "outer-gap" model beams, but even then some sources are substantially discrepant, suggesting imperfect beaming corrections and/or errors in the fit distance. The fits prefer large metal abundance for half of the targets, a reasonable result for these wind-stripped secondaries. The companion radii indicate substantial Roche-lobe filling, $f_c \approx 0.7-1$ except for PSR J0952$-$0607, which with $f_c< 0.5$ has a companion density $\rho \approx 10\,{\rm g\,cm^{-3}}$, suggesting unusual evolution. We note that the direct-heating fits imply large heating powers and rather small inclinations, and we speculate that unmodeled effects can introduce such bias., Comment: To appear in the Astrophysical Journal Revised version, fixed duplicate panel in Fig 6

To assist with the commissioning (Jiang et al. 2019) of the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed a pulsar search, with the primary goal of developing and testing the pulsar data acquisition and processing pipelines. We tested and used three pipelines, two (P1 and P2 hereafter) searched for the periodic signature of pulsars whereas the other one was used to search for bright single pulses (P3 hereafter). A pulsar candidate was discovered in the observation on the 22nd August, 2017, and later confirmed by the Parkes radio telescope on the 10th September, 2017. This pulsar, named PSR J1900-0134, was the first pulsar discovered by FAST. The pulsar has a pulse period of 1.8 s and a dispersion measure (DM) of 188\,pc\,cm$^{-3}$., Comment: Sci. China-Phys. Mech. Astron. 62, 959508 (2019)

The anomalous X-ray pulsar XTE J1810$-$197 was the first magnetar found to emit pulsed radio emission. After spending almost a decade in a quiescent, radio-silent state, the magnetar was reported to have undergone a radio outburst in December, 2018. We observed radio pulsations from XTE J1810$-$197 during this early phase of its radio revival using the Ultra-Wideband Low receiver system of the Parkes radio telescope, obtaining wideband (704 MHz to 4032 MHz) polarization pulse profiles, single pulses and flux density measurements. Dramatic changes in polarization and rapid variations of the position angle of linear polarization across the main pulse and in time have been observed. The pulse profile exhibits similar structures throughout our three observations (over a week time scale), displaying a small amount of profile evolution in terms of polarization and pulse width across the wideband. We measured a flat radio spectrum across the band with a positive spectral index, in addition to small levels of flux and spectral index variability across our observing span. The observed wideband polarization properties are significantly different compared to those taken after the 2003 outburst, and therefore provide new information about the origin of radio emission., Comment: 7 pages, 4 figures, accepted by ApJ Letters

Identifying as many gamma-ray pulsars as possible in the Fermi Large Area Telescope (LAT) data helps test pulsar emission models by comparing predicted and observed properties for a large, varied sample with as little selection bias as possible. It also improves extrapolations from the observed population to estimate the contribution of unresolved pulsars to the diffuse gamma-ray emission. We use a recently developed method to determine the probability that a given gamma-ray photon comes from a known position in the sky, convolving the photon's energy with the LAT's energy-dependent point-spread-function (PSF), without the need for an accurate spatial and spectral model of the gamma-ray sky around the pulsar. The method is simple and fast and, importantly, provides probabilities, or _weights_, for gamma rays from pulsars too faint for phase-integrated detection. We applied the method to over a thousand pulsars for which we obtained rotation ephemerides from radio observations, and discovered gamma-ray pulsations from 16 pulsars, 12 young and 4 recycled. PSR J2208+4056 has spindown power Edot = 8e32 erg/s, about three times lower than the previous observed gamma-ray emission "deathline". PSRs J2208+4056 and J1816-0755 have radio interpulses, constraining their geometry and perhaps enhancing their gamma-ray luminosity. We discuss whether the deathline is an artifact of selection bias due to the pulsar distance., Comment: Accepted by the Astrophysical Journal

We report the discovery of a previously unidentified pulsar as part of a radio campaign to identify neutron star companions to low-mass white dwarfs (LMWDs) using the Robert C.\ Byrd Green Bank Telescope (GBT). PSR J0802-0955, which is coincident with the position of a WD with a mass of 0.2 solar masses, has a pulse period of 571 ms. Because of its relatively long pulse period, the lack of radial velocity (RV) variations in the radio data, and GBT's large beam size at the observing frequency of 340 MHz, we conclude that PSR J0802-0955 is unassociated with the LMWD at roughly the same position and distance., Comment: Accepted for publication in Research Notes of the AAS

We have constructed timing solutions for 81 gamma-ray pulsars covering more than five years of Fermi data. The sample includes 37 radio-quiet or radio-faint pulsars which cannot be timed with other telescopes. These timing solutions and the corresponding pulse times of arrival are prerequisites for further study, e.g. phase-resolved spectroscopy or searches for mode switches. Many gamma-ray pulsars are strongly affected by timing noise, and we present a new method for characterizing the noise process and mitigating its effects on other facets of the timing model. We present an analysis of timing noise over the population using a new metric for characterizing its strength and spectral shape, namely its time-domain correlation. The dependence of the strength on spin frequency and spin-down rate is in good agreement with previous studies. We find that noise process power spectra $S(f)$ for unrecycled pulsars are steep, with strong correlations over our entire data set and spectral indices $S(f)\propto f^{-\alpha}$ of 5 to 9. One possible explanation for these results is the occurrence of unmodelled, episodic 'microglitches'. Finally, we show that our treatment of timing noise results in robust parameter estimation, and in particular we measure a precise timing position for each pulsar. We extensively validate our results with multi-wavelength astrometry, and using our updated position, we firmly identify the X-ray counterpart of PSR J1418-6058., Comment: 13 pages, 11 figures, 2 tables; accepted in ApJ

PSR J2129-0429 is a "redback" eclipsing millisecond pulsar binary with an unusually long 15.2 hour orbit. It was discovered by the Green Bank Telescope in a targeted search of unidentified Fermi gamma-ray sources. The pulsar companion is optically bright (mean $m_R = 16.6$ mag), allowing us to construct the longest baseline photometric dataset available for such a system. We present ten years of archival and new photometry of the companion from LINEAR, CRTS, PTF, the Palomar 60-inch, and LCOGT. Radial velocity spectroscopy using the Double-Beam Spectrograph on the Palomar 200-inch indicates that the pulsar is massive: $1.74\pm0.18 M_\odot$. The G-type pulsar companion has mass $0.44\pm0.04 M_\odot$, one of the heaviest known redback companions. It is currently 95\% Roche-lobe filling and only mildly irradiated by the pulsar. We identify a clear 13.1 mmag yr$^{-1}$ secular decline in the mean magnitude of the companion as well as smaller-scale variations in the optical lightcurve shape. This behavior may indicate that the companion is cooling. Binary evolution calculations indicate that PSR J2129-0429 has an orbital period almost exactly at the bifurcation period between systems that converge into tighter orbits as black widows and redbacks and those that diverge into wider pulsar--white dwarf binaries. Its eventual fate may depend on whether it undergoes future episodes of mass transfer and increased irradiation., Comment: 13 pages, 8 figures. Submitted to ApJ

PSR J2032+4127 is a gamma-ray and radio-emitting pulsar which has been regarded as a young luminous isolated neutron star. However, its recent spin-down rate has extraordinarily increased by a factor of two. We present evidence that this is due to its motion as a member of a highly-eccentric binary system with a 15-solar-mass Be star, MT91~213. Timing observations show that, not only are the positions of the two stars coincident within 0.4 arcsec, but timing models of binary motion of the pulsar fit the data much better than a model of a young isolated pulsar. MT91~213, and hence the pulsar, lie in the Cyg~OB2 stellar association, which is at a distance of only 1.4-1.7 kpc. The pulsar is currently on the near side of, and accelerating towards, the Be star, with an orbital period of 20-30 years. The next periastron is well-constrained to occur in early 2018, providing an opportunity to observe enhanced high-energy emission as seen in other Be-star binary systems., Comment: 8 pages, 5 figures. Accepted for publication in MNRAS

XSS J12270-4859 is an X-ray binary associated with the Fermi LAT gamma-ray source 1FGL J1227.9-4852. In 2012 December, this source underwent a transition where the X-ray and optical luminosity dropped and the spectral signatures of an accretion disc disappeared. We report the discovery of a 1.69 millisecond pulsar (MSP), PSR J1227-4853, at a dispersion measure of 43.4 pc cm$^{-3}$ associated with this source, using the GMRT at 607 MHz. This demonstrates that, post-transition, the system hosts an active radio MSP. This is the third system after PSR J1023+0038 and PSR J1824-2452I showing evidence of state switching between radio MSP and low-mass X-ray binary (LMXB) states. We report timing observations of PSR J1227-4853 with the GMRT and Parkes, which give a precise determination of the rotational and orbital parameters of the system. The companion mass measurement of 0.17 to 0.46 M$_{sun}$ suggests that this is a redback system. PSR J1227-4853 is eclipsed for about 40% of its orbit at 607 MHz; with additional short-duration eclipses at all orbital phases. We also find that the pulsar is very energetic, with a spin-down luminosity of ~ 10$^{35}$ erg s$^{-1}$. We report simultaneous imaging and timing observations with the GMRT, which suggests that eclipses are caused by absorption, rather than dispersion smearing or scattering., Comment: 16 pages, 3 figures, accepted for publication in ApJ Letters

We present optical spectroscopy, astrometry, radio, and X-ray observations of the runaway binary LP 400-22. We refine the orbital parameters of the system based on our new radial velocity observations. Our parallax data indicate that LP 400-22 is significantly more distant (3 sigma lower limit of 840 pc) than initially predicted. LP 400-22 has a tangential velocity in excess of 830 km/s; it is unbound to the Galaxy. Our radio and X-ray observations fail to detect a recycled millisecond pulsar companion, indicating that LP 400-22 is a double white dwarf system. This essentially rules out a supernova runaway ejection mechanism. Based on its orbit, a Galactic center origin is also unlikely. However, its orbit intersects the locations of several globular clusters; dynamical interactions between LP 400-22 and other binary stars or a central black hole in a dense cluster could explain the origin of this unusual binary., Comment: MNRAS, in press

We describe the first X-ray observations of binary millisecond pulsars PSRs J0023+0923, J1810+1744, J2215+5135, and J2256-1024. All four are Fermi gamma-ray sources and three are 'black-widow' pulsars, with companions of mass < 0.1 solar masses. Data were taken using the Chandra X-Ray Observatory and covered a full binary orbit for each pulsar. Two pulsars, PSRs J2215+5135 and J2256-1024, show significant orbital variability and X-ray flux minima at the times of eclipses observed at radio wavelengths. This phenomenon is consistent with intrabinary shock emission characteristic of black-widow pulsars. The other two pulsars, PSRs J0023+0923 and J1810+1744, do not demonstrate significant variability, but are fainter than the other two sources. Spectral fits yield power-law indices that range from 1.4 to 2.3 and blackbody temperatures in the hundreds of eV. The spectrum for PSR J2215+5135 shows a significant hard X-ray component (41% of counts are above 2 keV), which is additional evidence for the presence of intrabinary shock emission., Comment: 3 pages, 2 figures, 1 table, 2012 IAU General Assembly

We report on Chandra X-ray and Parkes radio observations of IGR J11014-6103, which is a possible pulsar wind nebula with a complex X-ray morphology and a likely radio counterpart. With the superb angular resolution of Chandra, we find evidence that a portion of the extended emission may be related to a bow shock due to the putative pulsar moving through the interstellar medium. The inferred direction of motion is consistent with IGR J11014-6103 having been born in the event that produced the supernova remnant (SNR) MSH 11-61A. If this association is correct, then previous constraints on the expansion of MSH 11-61A imply a transverse velocity for IGR J11014-6103 of 2,400-2,900 km/s, depending on the SNR model used. This would surpass the kick velocities of any known pulsars and rival or surpass the velocities of any compact objects that are associated with SNRs. While it is important to confirm the nature of the source, our radio pulsation search did not yield a detection., Comment: Accepted for ApJ Letters, 5 pages, 5 figures

We have analysed the long- and short-term time dependence of the pulse arrival times and the pulse detection rates for eight Rotating Radio Transient (RRAT) sources from the Parkes Multi-beam Pulsar Survey (PMPS). We find significant periodicities in the individual pulse arrival times from six RRATs. These periodicities range from 30 minutes to 2100 days and from one to 16 independent (i.e. non-harmonically related) periodicities are detected for each RRAT. In addition, we find that pulse emission is a random (i.e. Poisson) process on short (hour-long) time scales but that most of the objects exhibit longer term (months-years) non-random behaviour. We find that PSRs J1819-1458 and J1317-5759 emit more doublets (two consecutive pulses) and triplets (three consecutive pulses) than is expected in random pulse distributions. No evidence for such an excess is found for the other RRATs. There are several different models for RRAT emission depending on both extrinsic and intrinsic factors which are consistent with these properties., Comment: Accepted by MNRAS on 2011 July 2. Contains 11 pages, 4 figures, 4 tables

We have used the Green Bank Telescope at 350MHz to search 50 faint, unidentified Fermi Gamma-ray sources for radio pulsations. So far, these searches have resulted in the discovery of 10 millisecond pulsars, which are plausible counterparts to these unidentified Fermi sources. Here we briefly describe this survey and the characteristics of the newly discovered MSPs., Comment: 4 pages, 2 figures, to appear in AIP Conference Proceedings of Pulsar Conference 2010 "Radio Pulsars: a key to unlock the secrets of the Universe", Sardinia, October 2010

We have detected an Halpha bow shock nebula around PSR J1741-2054, a pulsar discovered through its GeV gamma-ray pulsations. The pulsar is only ~1.5" behind the leading edge of the shock. Optical spectroscopy shows that the nebula is non-radiative, dominated by Balmer emission. The Halpha images and spectra suggest that the pulsar wind momentum is equatorially concentrated and implies a pulsar space velocity ~150km/s, directed 15+/-10deg out of the plane of the sky. The complex Halpha profile indicates that different portions of the post-shock flow dominate line emission as gas moves along the nebula and provide an opportunity to study the structure of this unusual slow non-radiative shock under a variety of conditions. CXO ACIS observations reveal an X-ray PWN within this nebula, with a compact ~2.5" equatorial structure and a trail extending several arcmin behind. Together these data support a close (<0.5kpc) distance, a spin geometry viewed edge-on and highly efficient gamma-ray production for this unusual, energetic pulsar., Comment: 7 pages, 8 figures, To appear in the Astrophysical Journal

We present the evolution of the radio emission from the 2.8-s pulsar of the double pulsar system PSR J0737-3039A/B. We provide an update on the Burgay et al. (2005) analysis by describing the changes in the pulse profile and flux density over five years of observations, culminating in the B pulsar's radio disappearance in 2008 March. Over this time, the flux density decreases by 0.177 mJy/yr at the brightest orbital phases and the pulse profile evolves from a single to a double peak, with a separation rate of 2.6 deg/yr. The pulse profile changes are most likely caused by relativistic spin precession, but can not be easily explained with a circular hollow-cone beam as in the model of Clifton & Weisberg (2008). Relativistic spin precession, coupled with an elliptical beam, can model the pulse profile evolution well. This particular beam shape predicts geometrical parameters for the two bright orbital phases which are consistent and similar to those derived by Breton et al. (2008). However, the observed decrease in flux over time and B's eventual disappearance cannot be easily explained by the model and may be due to the changing influence of A on B., Comment: 20 pages, 18 figures, Accepted by ApJ on 2 August 2010

We measure the radial velocity curve of the eclipsing detached white dwarf binary NLTT 11748. The primary exhibits velocity variations with a semi-amplitude of 273 km/s and an orbital period of 5.641 hr. We do not detect any spectral features from the secondary star, or any spectral changes during the secondary eclipse. We use our composite spectrum to constrain the temperature and surface gravity of the primary to be Teff = 8690 (140) K and log g = 6.54 (0.05), which correspond to a mass of 0.18 Msol. For an inclination angle of 89.9 deg. derived from the eclipse modeling, the mass function requires a 0.76 Msol companion. The merger time for the system is 7.2 Gyr. However, due to the extreme mass ratio of 0.24, the binary will most likely create an AM CVn system instead of a merger., Comment: ApJ Letters, in press

PSR J1802-2124 is a 12.6-ms pulsar in a 16.8-hour binary orbit with a relatively massive white dwarf (WD) companion. These properties make it a member of the intermediate-mass class of binary pulsar (IMBP) systems. We have been timing this pulsar since its discovery in 2002. Concentrated observations at the Green Bank Telescope, augmented with data from the Parkes and Nancay observatories, have allowed us to determine the general relativistic Shapiro delay. This has yielded pulsar and white dwarf mass measurements of 1.24(11) and 0.78(4) solar masses (68% confidence), respectively. The low mass of the pulsar, the high mass of the WD companion, the short orbital period, and the pulsar spin period may be explained by the system having gone through a common-envelope phase in its evolution. We argue that selection effects may contribute to the relatively small number of known IMBPs., Comment: 9 pages, 4 figures, 3 tables, accepted for publication in the Astrophysical Journal

SDSS J091709.55+463821.8 (hereafter J0917+4638) is the lowest surface gravity white dwarf (WD) currently known, with log g = 5.55 +/- 0.05 (M ~ 0.17 M_sun; Kilic et al. 2007a,b). Such low-mass white dwarfs (LMWDs) are believed to originate in binaries that evolve into WD/WD or WD/neutron star (NS) systems. An optical search for J0917+4638's companion showed that it must be a compact object with a mass >= 0.28 M_sun (Kilic 2007b). Here we report on Green Bank Telescope 820 MHz and XMM-Newton X-ray observations of J0917+4638 intended to uncover a potential NS companion to the LMWD. No convincing pulsar signal is detected in our radio data. Our X-ray observation also failed to detect X-ray emission from J0917+4638's companion, while we would have detected any of the millisecond radio pulsars in 47 Tuc. We conclude that the companion is almost certainly another WD., Comment: 4 pages, 1 table; to appear in the Astrophysical Journal Letters

We have conducted a search for pulsar companions to 15 low-mass white dwarfs (LMWDs; M < 0.4 M_Sun) at 820 MHz with the NRAO Green Bank Telescope (GBT). These LMWDs were spectroscopically identified in the Sloan Digital Sky Survey (SDSS), and do not show the photometric excess or spectroscopic signature associated with a companion in their discovery data. However, LMWDs are believed to evolve in binary systems and to have either a more massive WD or a neutron star as a companion. Indeed, evolutionary models of low-mass X-ray binaries, the precursors of millisecond pulsars (MSPs), produce significant numbers of LMWDs (e.g., Benvenuto & De Vito 2005), suggesting that the SDSS LMWDs may have neutron star companions. No convincing pulsar signal is detected in our data. This is consistent with the findings of van Leeuwen et al. (2007), who conducted a GBT search for radio pulsations at 340 MHz from unseen companions to eight SDSS WDs (five are still considered LMWDs; the three others are now classified as "ordinary" WDs). We discuss the constraints our non-detections place on the probability P_MSP that the companion to a given LMWD is a radio pulsar in the context of the luminosity and acceleration limits of our search; we find that P_MSP < 10 +4 -2 %., Comment: 5 pages, 1 figure, 1 table; accepted for publication in the Astrophysical Journal

We present timing data spanning 6.4 yr for the young and energetic PSR J0205+6449, in the supernova remnant 3C 58. Data were obtained with the Rossi X-ray Timing Explorer, the Jodrell Bank Observatory and the Green Bank Telescope. We present phase-coherent timing analyses showing timing noise and two spin-up glitches with fractional frequency increases of ~3.4E-7 near MJD 52555, and ~3.8E-6 between MJDs 52777 and 53062. These glitches are unusually large if the pulsar was created in the historical supernova in 1181 as has been suggested. For the X-ray timing we developed a new unbinned maximum-likelihood method for determining pulse arrival times which performs significantly better than the traditional binned techniques. In addition, we present an X-ray pulse profile analysis of four years of RXTE data showing that the pulsar is detected up to ~40 keV. We also present the first measurement of the phase offset between the radio and X-ray pulse for this source, showing that the radio pulse leads the X-ray pulse by phi=0.10+/-0.01 in phase. We compile all known measurements of the phase offsets between radio and X-ray and radio and gamma-ray pulses for X-ray and gamma-ray pulsars. We show that there is no relationship between pulse period and phase offset, supported by our measurement of the phase offset for PSR J0205+6449., Comment: 19 pages, 12 figures. Published in the Astrophysical Journal. Includes additional data analysis and two new figures