• ### A Luminous and Highly-variable Gamma-ray Flare Following the 2017 Periastron of PSR B1259-63/LS 2883(1805.03537)

May 9, 2018 astro-ph.HE
Three periastron passages of the PSR B1259$-$63/LS 2883 binary system, consisting of a 48 ms rotation-powered pulsar and a $\sim30$ M$_{\odot}$ Be star, have been observed by the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope, in 2010, 2014, and 2017. During the most-recent periastron passage, sustained low-level gamma-ray emission was observed over a $\sim3$-week long interval immediately after periastron, which was followed by an interval of no emission. Sporadic flares were detected starting 40 days post-periastron and lasted approximately 50 days, during which the emission displayed significant spectral curvature, variability on timescales as short as 1.5 minutes, and peak flux levels well in excess of the pulsar spin-down power. By contrast, during the 2010 and 2014 periastron passages, significant gamma-ray emission was not observed with the LAT until 30 and 32 days post-periastron, respectively. The previous flares did not exhibit spectral curvature, showed no short term variability, and did not exceed the pulsar spin-down power. The high flux and short timescales observed in 2017 suggest significant beaming of the emission is required and constrain the size of the emission region. The flares occur long enough after periastron that the neutron star should already have passed through the extended disk-like outflow, thus constraining options for target material and seed photon sources for inverse Compton models.
• ### Einstein@Home discovers a radio-quiet gamma-ray millisecond pulsar(1803.06855)

March 19, 2018 astro-ph.HE
Millisecond pulsars (MSPs) are old neutron stars that spin hundreds of times per second and appear to pulsate as their emission beams cross our line of sight. To date, radio pulsations have been detected from all rotation-powered MSPs. In an attempt to discover radio-quiet gamma-ray MSPs, we used the aggregated power from the computers of tens of thousands of volunteers participating in the Einstein@Home distributed computing project to search for pulsations from unidentified gamma-ray sources in Fermi Large Area Telescope data. This survey discovered two isolated MSPs, one of which is the only known rotation-powered MSP to remain undetected in radio observations. These gamma-ray MSPs were discovered in completely blind searches without prior constraints from other observations, raising hopes for detecting MSPs from a predicted Galactic bulge population.
• ### The Einstein@Home Gamma-Ray Pulsar Survey II. Source Selection, Spectral Analysis and Multi-wavelength Follow-up(1712.05395)

Dec. 14, 2017 astro-ph.HE
We report on the analysis of 13 gamma-ray pulsars discovered in the Einstein@Home blind search survey using Fermi Large Area Telescope (LAT) Pass 8 data. The 13 new gamma-ray pulsars were discovered by searching 118 unassociated LAT sources from the third LAT source catalog (3FGL), selected using the Gaussian Mixture Model (GMM) machine learning algorithm on the basis of their gamma-ray emission properties being suggestive of pulsar magnetospheric emission. The new gamma-ray pulsars have pulse profiles and spectral properties similar to those of previously-detected young gamma-ray pulsars. Follow-up radio observations have revealed faint radio pulsations from two of the newly-discovered pulsars, and enabled us to derive upper limits on the radio emission from the others, demonstrating that they are likely radio-quiet gamma-ray pulsars. We also present results from modeling the gamma-ray pulse profiles and radio profiles, if available, using different geometric emission models of pulsars. The high discovery rate of this survey, despite the increasing difficulty of blind pulsar searches in gamma rays, suggests that new systematic surveys such as presented in this article should be continued when new LAT source catalogs become available.
• ### Extreme Scattering Events Towards Two Young Pulsars(1712.00426)

Dec. 1, 2017 astro-ph.HE
We have measured the scintillation properties of 151 young, energetic pulsars with the Parkes radio telescope and have identified two extreme scattering events (ESEs). Towards PSR J1057-5226 we discovered a three-year span of strengthened scattering during which the variability in flux density and the scintillation bandwidth decreased markedly. The transverse size of the scattering region is ~23 au, and strong flux density enhancement before and after the ESE may arise from refractive focusing. Long observations reveal scintillation arcs characteristic of interference between rays scattered at large angles, and the clearest arcs appear during the ESE. The arcs suggest scattering by a screen 100-200 pc from the earth, perhaps ionized filamentary structure associated with the boundary of the local bubble(s). Towards PSR J1740-3015 we observed a "double dip" in the measured flux density similar to ESEs observed towards compact extragalactic radio sources. The observed shape is consistent with that produced by a many-au scale diverging plasma lens with electron density ~500cm$^{-3}$ . The continuing ESE is at least 1500 d long, making it the longest detected event to date. These detections, with materially different observational signatures, indicate that well-calibrated pulsar monitoring is a keen tool for ESE detection and ISM diagnostics. They illustrate the strong r\^ole au-scale non-Kolmogorov density fluctuations and the local ISM structure play in such events and are key to understanding both their intrinsic physics and their impact on other phenomena, particularly fast radio bursts.
• ### Magnetar-like X-ray bursts suppress pulsar radio emission(1710.03718)

Oct. 10, 2017 astro-ph.HE
Rotation-powered pulsars and magnetars are two different observational manifestations of neutron stars: rotation powered pulsars are rapidly spinning objects that are mostly observed as pulsating radio sources, while magnetars, neutron stars with the highest known magnetic fields, often emit short-duration X-ray bursts. Here we report simultaneous observations of the high-magnetic-field radio pulsar PSR J1119-6127 at X-ray, with XMM-Newton & NuSTAR, and at radio energies with Parkes radio telescope, during a period of magnetar-like bursts. The rotationally powered radio emission shuts off coincident with the occurrence of multiple X-ray bursts, and recovers on a time scale of ~70 seconds. These observations of related radio and X-ray phenomena further solidify the connection between radio pulsars and magnetars, and suggest that the pair plasma produced in bursts can disrupt the acceleration mechanism of radio emitting particles.
• ### Spectral properties of 441 radio pulsars(1709.08864)

Sept. 26, 2017 astro-ph.HE
We present a study of the spectral properties of 441 pulsars observed with the Parkes radio telescope near the centre frequencies of 728, 1382 and 3100 MHz. The observations at 728 and 3100 MHz were conducted simultaneously using the dual-band 10-50cm receiver. These high-sensitivity, multi-frequency observations provide a systematic and uniform sample of pulsar flux densities. We combine our measurements with spectral data from the literature in order to derive the spectral properties of these pulsars. Using techniques from robust regression and information theory we classify the observed spectra in an objective, robust and unbiased way into five morphological classes: simple or broken power law, power law with either low or high-frequency cut-off and log-parabolic spectrum. While about $79 \%$ of the pulsars that could be classified have simple power law spectra, we find significant deviations in 73 pulsars, 35 of which have curved spectra, 25 with a spectral break and 10 with a low-frequency turn-over. We identify 11 gigahertz-peaked spectrum (GPS) pulsars, with 3 newly identified in this work and 8 confirmations of known GPS pulsars; 3 others show tentative evidence of GPS, but require further low-frequency measurements to support this classification. The weighted mean spectral index of all pulsars with simple power law spectra is $-1.60 \pm 0.03$. The observed spectral indices are well described by a shifted log-normal distribution. The strongest correlations of spectral index are with spin-down luminosity, magnetic field at the light-cylinder and spin-down rate. We also investigate the physical origin of the observed spectral features and determine emission altitudes for three pulsars.
• We present results from the first directed search for nontensorial gravitational waves. While general relativity allows for tensorial (plus and cross) modes only, a generic metric theory may, in principle, predict waves with up to six different polarizations. This analysis is sensitive to continuous signals of scalar, vector or tensor polarizations, and does not rely on any specific theory of gravity. After searching data from the first observation run of the advanced LIGO detectors for signals at twice the rotational frequency of 200 known pulsars, we find no evidence of gravitational waves of any polarization. We report the first upper limits for scalar and vector strains, finding values comparable in magnitude to previously-published limits for tensor strain. Our results may be translated into constraints on specific alternative theories of gravity.
• ### Locating the intense interstellar scattering towards the inner Galaxy(1707.03842)

July 12, 2017 astro-ph.GA
We use VLBA+VLA observations to measure the sizes of the scatter-broadened images of 6 of the most heavily scattered known pulsars: 3 within the Galactic Centre (GC) and 3 elsewhere in the inner Galactic plane. By combining the measured sizes with temporal pulse broadening data from the literature and using the thin-screen approximation, we locate the scattering medium along the line of sight to these 6 pulsars. At least two scattering screens are needed to explain the observations of the GC sample. We show that the screen inferred by previous observations of SGR J1745-2900 and Sgr A*, which must be located far from the GC, falls off in strength on scales < 0.2 degree. A second scattering component closer to (< 2 kpc) or even (tentatively) within (< 700 pc) the GC produces most or all of the temporal broadening observed in the other GC pulsars. Outside the GC, the scattering locations for all three pulsars are ~2 kpc from Earth, consistent with the distance of the Carina-Sagittarius or Scutum spiral arm. For each object the 3D scattering origin coincides with a known HII region (and in one case also a supernova remnant), suggesting that such objects preferentially cause the intense interstellar scattering seen towards the Galactic plane. We show that the HII regions should contribute > 25% of the total dispersion measure (DM) towards these pulsars, and calculate reduced DM distances. Those distances for other pulsars lying behind HII regions may be similarly overestimated.
• We present the result of searches for gravitational waves from 200 pulsars using data from the first observing run of the Advanced LIGO detectors. We find no significant evidence for a gravitational-wave signal from any of these pulsars, but we are able to set the most constraining upper limits yet on their gravitational-wave amplitudes and ellipticities. For eight of these pulsars, our upper limits give bounds that are improvements over the indirect spin-down limit values. For another 32, we are within a factor of 10 of the spin-down limit, and it is likely that some of these will be reachable in future runs of the advanced detector. Taken as a whole, these new results improve on previous limits by more than a factor of two.
• ### Comparison of Pulsar Positions from Timing and Very Long Baseline Astrometry(1704.01011)

Pulsar positions can be measured with high precision using both pulsar timing methods and very-long-baseline interferometry (VLBI). Pulsar timing positions are referenced to a solar-system ephemeris, whereas VLBI positions are referenced to distant quasars. Here we compare pulsar positions from published VLBI measurements with those obtained from pulsar timing data from the Nanshan and Parkes radio telescopes in order to relate the two reference frames. We find that the timing positions differ significantly from the VLBI positions (and also differ between different ephemerides). A statistically significant change in the obliquity of the ecliptic of $2.16\pm0.33$\,mas is found for the JPL ephemeris DE405, but no significant rotation is found in subsequent JPL ephemerides. The accuracy with which we can relate the two frames is limited by the current uncertainties in the VLBI reference source positions and in matching the pulsars to their reference source. Not only do the timing positions depend on the ephemeris used in computing them, but also different segments of the timing data lead to varying position estimates. These variations are mostly common to all ephemerides, but slight changes are seen at the 10$\mu$as level between ephemerides.
• ### A transient, flat spectrum radio pulsar near the Galactic centre(1703.02544)

March 7, 2017 astro-ph.HE
• ### Robust Estimation of Scattering in Pulsar Timing Analysis(1703.02108)

March 6, 2017 astro-ph.IM, astro-ph.HE
We present a robust approach to incorporating models for the time-variable broadening of the pulse profile due to scattering in the ionized interstellar medium into profile-domain pulsar timing analysis. We use this approach to simultaneously estimate temporal variations in both the dispersion measure (DM) and scattering, together with a model for the pulse profile that includes smooth evolution as a function of frequency, and the pulsar's timing model. We show that fixing the scattering timescales when forming time-of-arrival estimates, as has been suggested in the context of traditional pulsar timing analysis, can significantly underestimate the uncertainties in both DM, and the arrival time of the pulse, leading to bias in the timing parameters. We apply our method using a new, publicly available, GPU accelerated code, both to simulations, and observations of the millisecond pulsar PSR J1643$-$1224. This pulsar is known to exhibit significant scattering variability compared to typical millisecond pulsars, and we find including low-frequency ($< 1$ GHz) data without a model for these scattering variations leads to significant periodic structure in the DM, and also biases the astrometric parameters at the $4\sigma$ level, for example, changing proper motion in right ascension by $0.50 \pm 0.12$. If low frequency observations are to be included when significant scattering variations are present, we conclude it is necessary to not just model those variations, but also to sample the parameters that describe the variations simultaneously with all other parameters in the model, a task for which profile domain pulsar timing is ideally suited.
• ### Wide-band Profile Domain Pulsar Timing Analysis(1612.05258)

We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the latter case we introduce a new parameterisation of the posterior that is extremely efficient in the low signal-to-noise regime and can be readily applied to a wide range of scientific problems. We apply this methodology to a series of simulations, and to between seven and nine yr of observations for PSRs J1713$+$0747, J1744$-$1134, and J1909$-$3744 with frequency coverage that spans 700-3600MHz. We use a smooth model for profile evolution across the full frequency range, and compare smooth and piecewise models for the temporal variations in DM. We find the profile domain framework consistently results in improved timing precision compared to the standard analysis paradigm by as much as 40% for timing parameters. Incorporating smoothness in the DM variations into the model further improves timing precision by as much as 30%. For PSR J1713+0747 we also detect pulse shape variation uncorrelated between epochs, which we attribute to variation intrinsic to the pulsar at a level consistent with previously published analyses. Not accounting for this shape variation biases the measured arrival times at the level of $\sim$30ns, the same order of magnitude as the expected shift due to gravitational-waves in the pulsar timing band.
• ### The magnetic field and turbulence of the cosmic web measured using a brilliant fast radio burst(1611.05758)

Nov. 17, 2016 astro-ph.CO, astro-ph.HE
Fast radio bursts (FRBs) are millisecond-duration events thought to originate beyond the Milky Way galaxy. Uncertainty surrounding the burst sources, and their propagation through intervening plasma, has limited their use as cosmological probes. We report on a mildly dispersed (dispersion measure 266.5+-0.1 pc cm^-3), exceptionally intense (120+-30 Jy), linearly polarized, scintillating burst (FRB 150807) that we directly localize to 9 arcmin^2. Based on a low Faraday rotation (12.0+-0.7 rad m^-2), we infer negligible magnetization in the circum-burst plasma and constrain the net magnetization of the cosmic web along this sightline to <21 nG, parallel to the line-of-sight. The burst scintillation suggests weak turbulence in the ionized intergalactic medium.
• ### The Braking Index of a Radio-quiet Gamma-ray Pulsar(1611.01292)

Nov. 7, 2016 astro-ph.HE
We report the discovery and timing measurements of PSR J1208-6238, a young and highly magnetized gamma-ray pulsar, with a spin period of 440 ms. The pulsar was discovered in gamma-ray photon data from the Fermi Large Area Telescope (LAT) during a blind-search survey of unidentified LAT sources, running on the distributed volunteer computing system Einstein@Home. No radio pulsations were detected in dedicated follow-up searches with the Parkes radio telescope, with a flux density upper limit at 1369 MHz of 30 $\mu$Jy. By timing this pulsar's gamma-ray pulsations, we measure its braking index over five years of LAT observations to be $n = 2.598 \pm 0.001 \pm 0.1$, where the first uncertainty is statistical and the second estimates the bias due to timing noise. Assuming its braking index has been similar since birth, the pulsar has an estimated age of around 2,700 yr, making it the youngest pulsar to be found in a blind search of gamma-ray data and the youngest known radio-quiet gamma-ray pulsar. Despite its young age the pulsar is not associated with any known supernova remnant or pulsar wind nebula. The pulsar's inferred dipolar surface magnetic field strength is $3.8 \times 10^{13}$ G, almost 90% of the quantum-critical level. We investigate some potential physical causes of the braking index deviating from the simple dipole model but find that LAT data covering a longer time interval will be necessary to distinguish between these.
• ### Fermi Large Area Telescope Observations of the Monoceros Loop Supernova Remnant(1608.06380)

Aug. 23, 2016 astro-ph.HE
We present an analysis of the gamma-ray measurements by the Large Area Telescope onboard the \textit{Fermi Gamma-ray Space Telescope} in the region of the supernova remnant~(SNR) Monoceros Loop~(G205.5$+$0.5). The brightest gamma-ray peak is spatially correlated with the Rosette Nebula, which is a molecular cloud complex adjacent to the southeast edge of the SNR. After subtraction of this emission by spatial modeling, the gamma-ray emission from the SNR emerges, which is extended and fit by a Gaussian spatial template. The gamma-ray spectra are significantly better reproduced by a curved shape than a simple power law. The luminosities between 0.2--300~GeV are $\sim$~$4 \times 10^{34}$~erg~s$^{-1}$ for the SNR and $\sim$~$3 \times 10^{34}$~erg~s$^{-1}$ for the Rosette Nebula, respectively. We argue that the gamma rays likely originate from the interactions of particles accelerated in the SNR. The decay of neutral pions produced in nucleon-nucleon interactions of accelerated hadrons with interstellar gas provides a reasonable explanation for the gamma-ray emission of both the Rosette Nebula and the Monoceros SNR.
• ### The disturbance of a millisecond pulsar magnetosphere(1608.02163)

Aug. 6, 2016 astro-ph.SR, astro-ph.HE
Pulsar timing has enabled some of the strongest tests of fundamental physics. Central to the technique is the assumption that the detected radio pulses can be used to accurately measure the rotation of the pulsar. Here we report on a broad-band variation in the pulse profile of the millisecond pulsar J1643-1224. A new component of emission suddenly appears in the pulse profile, decays over 4 months, and results in a permanently modified pulse shape. Profile variations such as these may be the origin of timing noise observed in other millisecond pulsars. The sensitivity of pulsar-timing observations to gravitational radiation can be increased by accounting for this variability.
• ### Characterising the rotational irregularities of the Vela pulsar from 21 yr of phase-coherent timing(1604.02661)

April 10, 2016 astro-ph.SR, astro-ph.HE
Pulsars show two classes of rotational irregularities that can be used to understand neutron-star interiors and magnetospheres: glitches and timing noise. Here we present an analysis of the Vela pulsar spanning nearly 21 yr of observation and including 8 glitches. We identify the relative pulse number of all of the observations between glitches, with the only pulse-number ambiguities existing over glitch events. We use the phase coherence of the timing solution to simultaneously model the timing noise and glitches in a Bayesian framework, allowing us to select preferred models for both. We find the glitches can be described using only permanent and transient changes in spin frequency, i.e., no step changes in frequency derivative. For all of the glitches, we only need two exponentially decaying changes in spin frequency to model the transient components. In contrast to previous studies, we find that the dominant transient components decay on a common $\approx$ 1300 d time scale, and that a larger fraction ( $\gtrsim 25\%$) of glitch amplitudes are associated with these transient components. We also detect shorter-duration transient components of $\approx$ 25 d, as previously observed, but are limited in sensitivity to events with shorter durations by the cadence of our observations. The timing noise is well described by a steep power-law process that is independent of the glitches and subdominant to the glitch recovery. The braking index is constrained to be $<$ 8 with 95% confidence. This methodology can be used to robustly measure the properties of glitches and timing noise in other pulsars.
• ### Discovery of a millisecond pulsar in the 5.4 day binary 3FGL J1417.5-4402: observing the late phase of pulsar recycling(1603.00142)

March 1, 2016 astro-ph.SR, astro-ph.HE
In a search of the unidentified Fermi gamma-ray source 3FGL J1417.5-4402 with the Parkes radio telescope, we discovered PSR J1417-4402, a 2.66 ms pulsar having the same 5.4 day orbital period as the optical and X-ray binary identified by Strader et al. The existence of radio pulsations implies that the neutron star is currently not accreting. Substantial outflows from the companion render the radio pulsar undetectable for more than half of the orbit, and may contribute to the observed Halpha emission. Our initial pulsar observations, together with the optically inferred orbit and inclination, imply a mass ratio of 0.171+/-0.002, a companion mass of M_2=0.33+/-0.03 Msun, and a neutron star mass in the range 1.77<M_1<2.13 Msun. However, there remains a discrepancy between the distance of 4.4 kpc inferred from the optical properties of the companion and the smaller radio dispersion measure distance of 1.6 kpc. The smaller distance would reduce the inferred Roche-lobe filling factor, increase the inferred inclination angle, and decrease the masses. As a wide binary, PSR J1417-4402 differs from the radio-eclipsing black widow and redback pulsars being discovered in large numbers by Fermi. It is probably a system that began mass transfer onto the neutron star after the companion star left the main sequence. The companion should end its evolution as a He white dwarf in a 6-20 day orbit, i.e., as a typical binary millisecond pulsar companion.
• ### From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release(1602.05570)

We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.
• ### The International Pulsar Timing Array: First Data Release(1602.03640)

Feb. 11, 2016 astro-ph.IM
The highly stable spin of neutron stars can be exploited for a variety of (astro-)physical investigations. In particular arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such "Pulsar Timing Arrays" (PTAs) have been set up around the world over the past decades and collectively form the "International" PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available online) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits placed by individual PTAs by a factor of ~2 and provides a 2-sigma limit on the dimensionless amplitude of a stochastic GWB of 1.7x10^{-15} at a frequency of 1 yr^{-1}. This is 1.7 times less constraining than the limit placed by (Shannon et al. 2015), due mostly to the more recent, high-quality data they used.
• ### Six New Millisecond Pulsars from Arecibo Searches of Fermi Gamma-Ray Sources(1601.05343)

Jan. 20, 2016 astro-ph.SR, astro-ph.HE
We have discovered six radio millisecond pulsars (MSPs) in a search with the Arecibo telescope of 34 unidentified gamma-ray sources from the Fermi Large Area Telescope (LAT) 4-year point source catalog. Among the 34 sources, we also detected two MSPs previously discovered elsewhere. Each source was observed at a center frequency of 327 MHz, typically at three epochs with individual integration times of 15 minutes. The new MSP spin periods range from 1.99 to 4.66 ms. Five of the six pulsars are in interacting compact binaries (period < 8.1 hr), while the sixth is a more typical neutron star-white dwarf binary with an 83-day orbital period. This is a higher proportion of interacting binaries than for equivalent Fermi-LAT searches elsewhere. The reason is that Arecibo's large gain afforded us the opportunity to limit integration times to 15 minutes, which significantly increased our sensitivity to these highly accelerated systems. Seventeen of the remaining 26 gamma-ray sources are still categorized as strong MSP candidates, and will be re-searched.
• ### Multiwavelength Observations of the Redback Millisecond Pulsar J1048+2339(1601.03681)

Jan. 14, 2016 astro-ph.SR, astro-ph.HE
We report on radio timing and multiwavelength observations of the 4.66 ms redback pulsar J1048+2339, which was discovered in an Arecibo search targeting the Fermi-LAT source 3FGLJ1048.6+2338. Two years of timing allowed us to derive precise astrometric and orbital parameters for the pulsar. PSR J1048+2339 is in a 6-hour binary, and exhibits radio eclipses over half the orbital period and rapid orbital period variations. The companion has a minimum mass of 0.3 solar masses, and we have identified a $V \sim 20$ variable optical counterpart in data from several surveys. The phasing of its $\sim 1$~mag modulation at the orbital period suggests highly efficient and asymmetric heating by the pulsar wind, which may be due to an intrabinary shock that is distorted near the companion, or to the companion's magnetic field channeling the pulsar wind to specific locations on its surface. We also present gamma-ray spectral analysis of the source and preliminary results from searches for gamma-ray pulsations using the radio ephemeris.
• ### A pilot ASKAP survey of radio transient events in the region around the intermittent pulsar PSR J1107-5907(1512.02702)

We use observations from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to search for transient radio sources in the field around the intermittent pulsar PSR J1107-5907. The pulsar is thought to switch between an "off" state in which no emission is detectable, a weak state and a strong state. We ran three independent transient detection pipelines on two-minute snapshot images from a 13 hour BETA observation in order to 1) study the emission from the pulsar, 2) search for other transient emission from elsewhere in the image and 3) to compare the results from the different transient detection pipelines. The pulsar was easily detected as a transient source and, over the course of the observations, it switched into the strong state three times giving a typical timescale between the strong emission states of 3.7 hours. After the first switch it remained in the strong state for almost 40 minutes. The other strong states lasted less than 4 minutes. The second state change was confirmed using observations with the Parkes radio telescope. No other transient events were found and we place constraints on the surface density of such events on these timescales. The high sensitivity Parkes observations enabled us to detect individual bright pulses during the weak state and to study the strong state over a wide observing band. We conclude by showing that future transient surveys with ASKAP will have the potential to probe the intermittent pulsar population.
• ### Emission-rotation correlation in pulsars: new discoveries with optimal techniques(1511.05481)

Nov. 18, 2015 astro-ph.HE
Pulsars are known to display short-term variability. Recently, examples of longer-term emission variability have emerged that are often correlated with changes in the rotational properties of the pulsar. To further illuminate this relationship, we have developed techniques to identify emission and rotation variability in pulsar data, and determine correlations between the two. Individual observations may be too noisy to identify subtle changes in the pulse profile. We use Gaussian process (GP) regression to model noisy observations and produce a continuous map of pulse profile variability. Generally, multiple observing epochs are required to obtain the pulsar spin frequency derivative ($\dot{\nu}$). GP regression is, therefore, also used to obtain $\dot{\nu}$, under the hypothesis that pulsar timing noise is primarily caused by unmodelled changes in $\dot{\nu}$. Our techniques distinguish between two types of variability: changes in the total flux density versus changes in the pulse shape. We have applied these techniques to 168 pulsars observed by the Parkes radio telescope, and see that although variations in flux density are ubiquitous, substantial changes in the shape of the pulse profile are rare. We reproduce previously published results and present examples of profile shape changing in seven pulsars; in particular, a clear new example of correlated changes in profile shape and rotation is found in PSR~J1602$-$5100. In the shape changing pulsars, a more complex picture than the previously proposed two state model emerges. We conclude that our simple assumption that all timing noise can be interpreted as $\dot{\nu}$ variability is insufficient to explain our dataset.