Source author record

J. W. T. Hessels

J. W. T. Hessels appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph.IM astro-ph.SR astro-ph.GA gr-qc astro-ph.CO

Catalog footprint

What is connected

105works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

The FRB 20121102A November rain in 2018 observed with the Arecibo Telescope

We present 849 new bursts from FRB 20121102A detected with the 305-m Arecibo Telescope. Observations were conducted as part of our regular campaign to monitor activity and evolution of burst properties. The 10 reported observations were carried out between 1150 and 1730 MHz and fall in the active period around November 2018. All bursts were dedispersed at the same dispersion measure and are consistent with a single value of 562.4(1) pc/cm$^3$. The rate varies between 0 bursts and 218(16) bursts per hour, the highest rate observed to date. The times between consecutive bursts show a bimodal distribution. We find that a Poisson process with varying rate best describes arrival times with separations >0.1 s. Clustering on timescales of 22 ms reflects a characteristic timescale of the source and possibly the emission mechanism. We analyse the spectro-temporal structure of the bursts by fitting 2D Gaussians with a temporal drift to each sub-burst in the dynamic spectra. We find a linear relationship between the sub-burst's drift and its duration. At the same time, the drifts are consistent with coming from the sad-trombone effect. This has not been predicted by current models. The energy distribution shows an excess of high energy bursts and is insufficiently modelled by a single power-law even within single observations. We find long-term changes in the energy distribution, the average spectrum, and the sad-trombone drift, compared to earlier and later published observations. Despite the large burst rate, we find no strict short-term periodicity.

preprint2022arXiv

Arecibo observations of a burst storm from FRB 20121102A in 2016

FRB 20121102A is the first known fast radio burst (FRB) from which repeat bursts were detected, and one of the best-studied FRB sources in the literature. Here we report on the analysis of 478 bursts (333 previously unreported) from FRB 20121102A using the 305-m Arecibo telescope - detected during approximately 59 hours of observations between December 2015 and October 2016. The majority of bursts are from a burst storm around September 2016. This is the earliest available sample of a large number of FRB 20121102A bursts, and it thus provides an anchor point for long-term studies of the source's evolving properties. We observe that the bursts separate into two groups in the width-bandwidth-energy parameter space, which we refer to as the low-energy bursts (LEBs) and high-energy bursts (HEBs). The LEBs are typically longer duration and narrower bandwidth than the HEBs, reminiscent of the spectro-temporal differences observed between the bursts of repeating and non-repeating FRBs. We fit the cumulative burst rate-energy distribution with a broken power-law and find that it flattens out toward higher energies. The sample shows a diverse zoo of burst morphologies. Notably, burst emission seems to be more common at the top than the bottom of our 1150 - 1730 MHz observing band. We also observe that bursts from the same day appear to be more similar to each other than to those of other days, but this observation requires confirmation. The wait times and burst rates that we measure are consistent with previous studies. We discuss these results, primarily in the context of magnetar models.

preprint2022arXiv

Discoveries and Timing of Pulsars in NGC 6440

Using the MeerKAT radio telescope, a series of observations have been conducted to time the known pulsars and search for new pulsars in the globular cluster NGC 6440. As a result, two pulsars have been discovered, NGC 6440G and NGC 6440H, one of which is isolated and the other a non-eclipsing (at frequencies above 962 MHz) "Black Widow", with a very low mass companion (M$_{\rm c}$ > 0.006 M$_{\odot}$). It joins the other binary pulsars discovered so far in this cluster which all have low companion masses (M$_{\rm c}$ < 0.30 M$_{\odot}$). We present the results of long-term timing solutions obtained using data from both Green Bank and MeerKAT telescopes for these two new pulsars and an analysis of the pulsars NGC 6440C and NGC 6440D. For the isolated pulsar NGC 6440C, we searched for planets using a Markov Chain Monte Carlo technique. We find evidence for significant unmodelled variations but they cannot be well modelled as planets nor as part of a power-law red-noise process. Studies of the eclipses of the "Redback" pulsar NGC 6440D at two different frequency bands reveal a frequency dependence with longer and asymmetric eclipses at lower frequencies (962-1283 MHz).

preprint2022arXiv

Fast radio bursts at the dawn of the 2020s

Since the discovery of the first fast radio burst (FRB) in 2007, and their confirmation as an abundant extragalactic population in 2013, the study of these sources has expanded at an incredible rate. In our 2019 review on the subject we presented a growing, but still mysterious, population of FRBs -- 60 unique sources, 2 repeating FRBs, and only 1 identified host galaxy. However, in only a few short years new observations and discoveries have given us a wealth of information about these sources. The total FRB population now stands at over 600 published sources, 24 repeaters, and 19 host galaxies. Higher time resolution data, sustained monitoring, and precision localisations have given us insight into repeaters, host galaxies, burst morphology, source activity, progenitor models, and the use of FRBs as cosmological probes. The recent detection of a bright FRB-like burst from the Galactic magnetar SGR~1935+2154 provides an important link between FRBs and magnetars. There also continue to be surprising discoveries, like periodic modulation of activity from repeaters and the localisation of one FRB source to a relatively nearby globular cluster associated with the M81 galaxy. In this review, we summarise the exciting observational results from the past few years. We also highlight their impact on our understanding of the FRB population and proposed progenitor models. We build on the introduction to FRBs in our earlier review, update our readers on recent results, and discuss interesting avenues for exploration as the field enters a new regime where hundreds to thousands of new FRBs will be discovered and reported each year.

preprint2022arXiv

FRB 121102: drastic changes in the burst polarization contrasts with the stability of the persistent emission

We study milliarcsecond-scale properties of the persistent radio counterpart to FRB 121102 and investigate the spectro-polarimetric properties of a bright burst. For the former, we use European VLBI Network (EVN) observations in 2017 at 1.7 and 4.8 GHz. For the latter, we re-analyse the 1.7-GHz data from the 100-m Effelseberg telescope taken in 2016. These observations predate other polarimetric studies of FRB 121102, and yield the highest burst Faraday rotation measure (RM) to date, RM = 1.27*10^5 rad m^-2, consistent with the decreasing RM trend. The fractional polarization of the burst emission is 15% at 1.7 GHz. This can be reconciled with the high fractional polarization at higher frequencies if the Faraday width of the burst environment is 150 rad m^-2 - a bare 0.1% of the total Faraday rotation. The width may originate from minor non-uniformities in the Faraday screen, or from effects in the emitting region itself. The upper limit on the persistent source size is 1 pc, barely consistent with a young supernova (SN) scenario. The flux variability limit of <10% is not in favor of the young SN scenario, and challenges other interpretations as well. The fractional polarization of the faint persistent source is constrained at <25% at 4.8 GHz ruling out a common origin with the highly polarized individual bursts.

preprint2022arXiv

PRECISE localizations of repeating Fast Radio Bursts

Fast Radio Bursts (FRBs) are extremely luminous and brief signals (with duration of milliseconds or even shorter) of extragalactic origin. Despite the fact that hundreds of FRBs have been discovered to date, their nature still remains unclear. Precise localizations of FRBs can unveil their host galaxies and local environments -- and thus shed light on the physical processes that led to the burst production. However, this has only been achieved for a few FRBs to date. The European VLBI Network (EVN) is currently the only instrument capable of localizing FRBs down to the milliarcsecond level. This level of precision was critical to associate the first localized FRB, 20121102A, to a star-forming region in a low-metallicity dwarf galaxy and physically related it to a compact persistent radio source. Analogously, a second repeating FRB, 20180916B, was found to just outside the edge of a prominent star-forming region of a nearby spiral galaxy. The PRECISE project (Pinpointing REpeating ChIme Sources with EVN dishes), starting from 2019, has observed hundreds of hours per year with a subset of EVN telescopes with the goal of localizing repeating FRBs discovered by the CHIME/FRB Collaboration. The ultimate goal of PRECISE is to disentangling the environments where FRBs can be produced. Here we present the state of the art of the FRB field, the PRECISE project, and the localizations achieved until now, which have unveiled a variety of environments where FRBs can be found that challenges the current models.

preprint2022arXiv

Searching for pulsars associated with polarised point sources using LOFAR: Initial discoveries from the TULIPP project

Discovering radio pulsars, particularly millisecond pulsars (MSPs), is important for a range of astrophysical applications, such as testing theories of gravity or probing the magneto-ionic interstellar medium. We aim to discover pulsars that may have been missed in previous pulsar searches by leveraging known pulsar observables (primarily polarisation) in the sensitive, low-frequency radio images from the Low-Frequency Array (LOFAR) Two-metre Sky Survey (LoTSS), and have commenced the Targeted search, using LoTSS images, for polarised pulsars (TULIPP) survey. For this survey, we identified linearly and circularly polarised point sources with flux densities brighter than 2 mJy in LoTSS images at a centre frequency of 144 MHz with a 48 MHz bandwidth. Over 40 known pulsars, half of which are MSPs, were detected as polarised sources in the LoTSS images and excluded from the survey. We have obtained beam-formed LOFAR observations of 30 candidates, which were searched for pulsations using coherent de-dispersion. Here, we present the results of the first year of the TULIPP survey. We discovered two pulsars, PSRs J1049+5822 and J1602+3901, with rotational periods of P=0.73 s and 3.7 ms, respectively. We also detected a further five known pulsars (two slowly-rotating pulsars and three MSPs) for which accurate sky positions were not available to allow a unique cross-match with LoTSS sources. This targeted survey presents a relatively efficient method by which pulsars, particularly MSPs, may be discovered using the flexible observing modes of sensitive radio telescopes such as the Square Kilometre Array and its pathfinders/precursors, particularly since wide-area all-sky surveys using coherent de-dispersion are currently computationally infeasible.

preprint2022arXiv

X-ray burst and persistent emission properties of the magnetar SGR 1830-0645 in outburst

We report on NICER X-ray monitoring of the magnetar SGR 1830-0645 covering 223 days following its October 2020 outburst, as well as Chandra and radio observations. We present the most accurate spin ephemerides of the source so far: $ν=0.096008680(2)$~Hz, $\dotν=-6.2(1)\times10^{-14}$~Hz~s$^{-1}$, and a significant second and third frequency derivative terms indicative of non-negligible timing noise. The phase-averaged 0.8--7~keV spectrum is well fit with a double-blackbody (BB) model throughout the campaign. The BB temperatures remain constant at 0.46 and 1.2 keV. The areas and flux of each component decreased by a factor of 6, initially through a steep decay trend lasting about 46 days followed by a shallow long-term one. The pulse shape in the same energy range is initially complex, exhibiting three distinct peaks, yet with clear continuous evolution throughout the outburst towards a simpler, single-pulse shape. The rms pulsed fraction is high and increases from about 40% to 50%. We find no dependence of pulse shape or fraction on energy. These results suggest that multiple hotspots, possibly possessing temperature gradients, emerged at outburst-onset, and shrank as the outburst decayed. We detect 84 faint bursts with \nicer, having a strong preference for occurring close to the surface emission pulse maximum the first time this phenomenon is detected in such a large burst sample. This likely implies a very low altitude for the burst emission region, and a triggering mechanism connected to the surface active zone. Finally, our radio observations at several epochs and multiple frequencies reveal no evidence of pulsed or burst-like radio emission.

preprint2021arXiv

A repeating fast radio burst source in a globular cluster

Fast radio bursts (FRBs) are exceptionally luminous flashes of unknown physical origin, reaching us from other galaxies (Petroff et al. 2019). Most FRBs have only ever been seen once, while others flash repeatedly, though sporadically (Spitler et al. 2016, CHIME/FRB Collaboration et al. 2021). Many models invoke magnetically powered neutron stars (magnetars) as the engines producing FRB emission (Margalit & Metzger 2018, CHIME/FRB Collaboration et al. 2020). Recently, CHIME/FRB announced the discovery (Bhardwaj et al. 2021) of the repeating FRB 20200120E, coming from the direction of the nearby grand design spiral galaxy M81. Four potential counterparts at other observing wavelengths were identified (Bhardwaj et al. 2021) but no definitive association with these sources, or M81, could be made. Here we report an extremely precise localisation of FRB 20200120E, which allows us to associate it with a globular cluster (GC) in the M81 galactic system and to place it ~2pc offset from the optical center of light of the GC. This confirms (Bhardwaj et al. 2021) that FRB 20200120E is 40 times closer than any other known extragalactic FRB. Because such GCs host old stellar populations, this association strongly challenges FRB models that invoke young magnetars formed in a core-collapse supernova as powering FRB emission. We propose, instead, that FRB 20200120E is a highly magnetised neutron star formed via either accretion-induced collapse of a white dwarf or via merger of compact stars in a binary system (Margalit et al. 2019). Alternative scenarios involving compact binary systems, efficiently formed inside globular clusters, could also be responsible for the observed bursts.

preprint2021arXiv

Highly polarised microstructure from the repeating FRB 20180916B

Fast radio bursts (FRBs) are bright, coherent, short-duration radio transients of as-yet unknown extragalactic origin. FRBs exhibit a wide variety of spectral, temporal and polarimetric properties, which can unveil clues into their emission physics and propagation effects in the local medium. Here we present the high-time-resolution (down to 1 $μ$s) polarimetric properties of four 1.7-GHz bursts from the repeating FRB 20180916B, which were detected in voltage data during observations with the European VLBI Network (EVN). We observe a range of emission timescales spanning three orders of magnitude, with the shortest component width reaching 3-4 $μ$s (below which we are limited by scattering). This is the shortest timescale measured in any FRB, to date. We demonstrate that all four bursts are highly linearly polarised ($\gtrsim 80\%$), show no evidence for significant circular polarisation ($\lesssim 15\%$), and exhibit a constant polarisation position angle (PPA) during and between bursts. On short timescales ($\lesssim 100$ $μ$s), however, there appear to be subtle (few degree) PPA variations across the burst profiles. These observational results are most naturally explained in an FRB model where the emission is magnetospheric in origin, as opposed to models where the emission originates at larger distances in a relativistic shock.

preprint2021arXiv

Milliarcsecond localisation of the repeating FRB 20201124A

Very long baseline interferometric (VLBI) localisations of repeating fast radio bursts (FRBs) have demonstrated a diversity of local environments: from nearby star-forming regions to globular clusters. Here we report the VLBI localisation of FRB 20201124A using an ad-hoc array of dishes that also participate in the European VLBI Network (EVN). In our campaign, we detected 18 total bursts from FRB 20201124A at two separate epochs. By combining the visibilities from both observing epochs, we were able to localise FRB 20201124A with a 1-$σ$ error of 4.5 milliarcseconds (mas). We use the relatively large burst sample to investigate astrometric accuracy, and find that for $\gtrsim20$ baselines ($\gtrsim7$ dishes) that we can robustly reach milliarcsecond precision even using single-burst data sets. Sub-arcsecond precision is still possible for single bursts, even when only $\sim$ six baselines (four dishes) are available. We explore two methods for determining the individual burst positions: the peaks of the dirty maps and a Gaussian fit to the cross fringe pattern on the dirty maps. We found the latter to be more reliable due to the lower mean and standard deviation in the offsets from the FRB position. Our VLBI work places FRB 20201124A 705$\pm$26 mas (1-$σ$ errors) from the optical centre of the host galaxy, and consistent with originating from within the recently-discovered extended radio structure associated with star-formation in the host galaxy. Future high-resolution optical observations, e.g. with Hubble Space Telescope, can determine the proximity of our FRB 20201124A VLBI position to nearby knots of star formation.

preprint2021arXiv

Rotation Measure Evolution of the Repeating Fast Radio Burst Source FRB 121102

The repeating fast radio burst source FRB 121102 has been shown to have an exceptionally high and variable Faraday rotation measure (RM), which must be imparted within its host galaxy and likely by or within its local environment. In the redshifted ($z=0.193$) source reference frame, the RM decreased from $1.46\times10^5$~rad~m$^{-2}$ to $1.33\times10^5$~rad~m$^{-2}$ between January and August 2017, showing day-timescale variations of $\sim200$~rad~m$^{-2}$. Here we present sixteen FRB 121102 RMs from burst detections with the Arecibo 305-m radio telescope, the Effelsberg 100-m, and the Karl G. Jansky Very Large Array, providing a record of FRB 121102's RM over a 2.5-year timespan. Our observations show a decreasing trend in RM, although the trend is not linear, dropping by an average of 15\% year$^{-1}$ and is $\sim9.7\times10^4$~rad~m$^{-2}$ at the most recent epoch of August 2019. Erratic, short-term RM variations of $\sim10^3$~rad~m$^{-2}$ week$^{-1}$ were also observed between MJDs 58215--58247. A decades-old neutron star embedded within a still-compact supernova remnant or a neutron star near a massive black hole and its accretion torus have been proposed to explain the high RMs. We compare the observed RMs to theoretical models describing the RM evolution for FRBs originating within a supernova remnant. FRB 121102's age is unknown, and we find that the models agree for source ages of $\sim6-17$~years at the time of the first available RM measurements in 2017. We also draw comparisons to the decreasing RM of the Galactic center magnetar, PSR J1745--2900.

preprint2021arXiv

Study of 72 pulsars discovered in the PALFA survey: Timing analysis, glitch activity, emission variability, and a pulsar in an eccentric binary

We present new discoveries and results from long-term timing of 72 pulsars discovered in the Arecibo PALFA survey, including precise determination of astrometric and spin parameters, and flux density and scatter broadening measurements at 1.4 GHz. Notable discoveries include two young pulsars (characteristic ages $\sim$30 kyr) with no apparent supernova remnant associations, three mode changing, 12 nulling and two intermittent pulsars. We detected eight glitches in five pulsars. Among them is PSR J1939+2609, an apparently old pulsar (characteristic age $\sim$1 Gy), and PSR J1954+2529, which likely belongs to a newly-emerging class of binary pulsars. The latter is the only pulsar among the 72 that is clearly not isolated: a non-recycled neutron star with a 931-ms spin period in an eccentric ($e\,=\,0.114$) wide ($P_b\,=\,82.7\,$d) orbit with a companion of undetermined nature having a minimum mass of $\sim0.6\,M_{\odot}$. Since operations at Arecibo ceased in 2020 August, we give a final tally of PALFA sky coverage, and compare its 207 pulsar discoveries to the known population. On average, they are 50% more distant than other Galactic plane radio pulsars; PALFA millisecond pulsars (MSP) have twice the dispersion measure per unit spin period than the known population of MSP in the Plane. The four intermittent pulsars discovered by PALFA more than double the population of such objects, which should help to improve our understanding of pulsar magnetosphere physics. The statistics for these, RRATS, and nulling pulsars suggest that there are many more of these objects in the Galaxy than was previously thought.

preprint2020arXiv

A repeating fast radio burst source localised to a nearby spiral galaxy

Fast radio bursts (FRBs) are brief, bright, extragalactic radio flashes. Their physical origin remains unknown, but dozens of possible models have been postulated. Some FRB sources exhibit repeat bursts. Though over a hundred FRB sources have been discovered to date, only four have been localised and associated with a host galaxy, with just one of the four known to repeat. The properties of the host galaxies, and the local environments of FRBs, provide important clues about their physical origins. However, the first known repeating FRB has been localised to a low-metallicity, irregular dwarf galaxy, and the apparently non-repeating sources to higher-metallicity, massive elliptical or star-forming galaxies, suggesting that perhaps the repeating and apparently non-repeating sources could have distinct physical origins. Here we report the precise localisation of a second repeating FRB source, FRB 180916.J0158+65, to a star-forming region in a nearby (redshift $z = 0.0337 \pm 0.0002$) massive spiral galaxy, whose properties and proximity distinguish it from all known hosts. The lack of both a comparably luminous persistent radio counterpart and a high Faraday rotation measure further distinguish the local environment of FRB 180916.J0158+65 from that of the one previously localised repeating FRB source, FRB 121102. This demonstrates that repeating FRBs have a wide range of luminosities, and originate from diverse host galaxies and local environments.

preprint2020arXiv

Asymmetric mass ratios for bright double neutron-star mergers

The discovery of a radioactively powered kilonova associated with the binary neutron star merger GW170817 was the first - and still only - confirmed electromagnetic counterpart to a gravitational-wave event. However, observations of late-time electromagnetic emission are in tension with the expectations from standard neutron-star merger models. Although the large measured ejecta mass is potentially explained by a progenitor system that is asymmetric in terms of the stellar component masses, i.e. with a mass ratio $q$ of 0.7-0.8, the known Galactic population of merging double neutron star (DNS) systems (i.e. those that will coalesce within billions of years or less) has, until now, only consisted of nearly equal-mass ($q > 0.9$) binaries. PSR J1913+1102 is a DNS system in a 5-hour, low-eccentricity ($e = 0.09$) orbit, implying an orbital separation of 1.8 solar radii, with the two neutron stars predicted to coalesce in 470 million years due to gravitational-wave emission. Here we report that the masses of the two neutron stars, as measured by a dedicated pulsar timing campaign, are $1.62 \pm 0.03$ and $1.27 \pm 0.03$ solar masses for the pulsar and companion neutron star, respectively; with a measured mass ratio $q = 0.78 \pm 0.03$, it is the most asymmetric DNS among known merging systems. Based on this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30% (90% confidence) of the total population of merging DNS binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event.

preprint2020arXiv

Detection of Repeating FRB 180916.J0158+65 Down to Frequencies of 300 MHz

We report on the detection of seven bursts from the periodically active, repeating fast radio burst (FRB) source FRB 180916.J0158+65 in the 300-400-MHz frequency range with the Green Bank Telescope (GBT). Emission in multiple bursts is visible down to the bottom of the GBT band, suggesting that the cutoff frequency (if it exists) for FRB emission is lower than 300 MHz. Observations were conducted during predicted periods of activity of the source, and had simultaneous coverage with the Low Frequency Array (LOFAR) and the FRB backend on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. We find that one of the GBT-detected bursts has potentially associated emission in the CHIME band (400-800 MHz) but we detect no bursts in the LOFAR band (110-190 MHz), placing a limit of $α> -1.0$ on the spectral index of broadband emission from the source. We also find that emission from the source is severely band-limited with burst bandwidths as low as $\sim$40 MHz. In addition, we place the strictest constraint on observable scattering of the source, $<$ 1.7 ms, at 350 MHz, suggesting that the circumburst environment does not have strong scattering properties. Additionally, knowing that the circumburst environment is optically thin to free-free absorption at 300 MHz, we find evidence against the association of a hyper-compact HII region or a young supernova remnant (age $<$ 50 yr) with the source.

preprint2020arXiv

Periodic activity from a fast radio burst source

Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.

preprint2020arXiv

The LOFAR Tied-Array All-Sky Survey: Timing of 21 pulsars including the first binary pulsar discovered with LOFAR

We report on the multi-frequency timing observations of 21 pulsars discovered in the LOFAR Tied-Array All-Sky Survey (LOTAAS). The timing data were taken at central frequencies of 149 MHz (LOFAR) as well as 334 and 1532 MHz (Lovell Telecope). The sample of pulsars includes 20 isolated pulsars and the first binary pulsar discovered by the survey, PSR J1658$+$3630. We modelled the timing properties of the pulsars, which showed that they have, on average, larger characteristic ages. We present the pulse profiles of the pulsars across the three observing bands, where PSR J1643$+$1338 showed profile evolution that appears not to be well-described by the radius-to-frequency-mapping model. Furthermore, we modelled the spectra of the pulsars across the same observing bands, using a simple power law, and found an average spectral index of $-1.9 \pm 0.5$. Amongst the pulsars studied here, PSR J1657$+$3304 showed large flux density variations of a factor of 10 over 300 days, as well as mode changing and nulling on timescales of a few minutes. We modelled the rotational and orbital properties of PSR J1658$+$3630, which has a spin period of 33 ms in a binary orbit of 3.0 days with a companion of minimum mass of 0.87$M_{\odot}$, likely a Carbon-Oxygen or Oxygen-Neon-Magnesium type white dwarf. PSR J1658$+$3630 has a dispersion measure of 3.0 pc cm$^{-3}$, making it possibly one of the closest binary pulsars known.

preprint2019arXiv

A LOFAR census of non-recycled pulsars: extending below 80 MHz

We present the results from the low-frequency (40--78 MHz) extension of the first LOFAR pulsar census of non-recycled pulsars. We have used the Low-Band Antennas of the LOFAR core stations to observe 87 pulsars out of 158 that have been detected previously with the High-Band Antennas. Forty-three pulsars have been detected and we present here their flux densities and flux-calibrated profiles. Seventeen of these pulsars have not been, to our knowledge, detected before at such low frequencies. We re-calculate the spectral indices using the new low-frequency flux density measurements from the LOFAR census and discuss the prospects of studying pulsars at the very low frequencies with the current and upcoming facilities, such as NenuFAR.

preprint2019arXiv

Detection and timing of gamma-ray pulsations from the $707$ Hz pulsar J0952$-$0607

The Low-Frequency Array radio telescope discovered the $707$ Hz binary millisecond pulsar (MSP) J0952$-$0607 in a targeted radio pulsation search of an unidentified $\textit{Fermi}$ gamma-ray source. This source shows a weak energy flux of $F_γ= 2.6 \times 10^{-12}\,\text{erg}\,\text{cm}^{-2}\,\text{s}^{-1}$ in the energy range between $100\,\text{MeV}$ and $100\,\text{GeV}$. Here we report the detection of pulsed gamma-ray emission from PSR$\,$J0952$-$0607 in a very sensitive gamma-ray pulsation search. The pulsar's rotational, binary, and astrometric properties are measured over seven years of $\textit{Fermi}$-Large Area Telescope data. For this we take into account the uncertainty on the shape of the gamma-ray pulse profile. We present an updated radio-timing solution now spanning more than two years and show results from optical modeling of the black-widow-type companion based on new multi-band photometric data taken with HiPERCAM on the Gran Telescopio Canarias on La Palma and ULTRACAM on the New Technology Telescope at ESO La Silla. PSR$\,$J0952$-$0607 is now the fastest-spinning pulsar for which the intrinsic spin-down rate has been reliably constrained ($\dot{P}_\text{int} \lesssim 4.6 \times 10^{-21}\,\text{s}\,\text{s}^{-1}$). The inferred surface magnetic field strength of $B_\text{surf} \lesssim 8.2 \times 10^{7}\,\text{G}$ is among the ten lowest of all known pulsars. This discovery is another example of an extremely fast spinning black-widow pulsar hiding within an unidentified $\textit{Fermi} gamma-ray source. In the future such systems might help to pin down the maximum spin frequency and the minimum surface magnetic field strength of MSPs.

preprint2019arXiv

Quasi-simultaneous radio and X-ray observations of Aql X-1: probing low luminosities

Aql X-1 is one of the best-studied neutron star low-mass X-ray binaries. It was previously targeted using quasi-simultaneous radio and X-ray observations during at least 7 different accretion outbursts. Such observations allow us to probe the interplay between accretion inflow (X-ray) and jet outflow (radio). Thus far, these combined observations have only covered one order of magnitude in radio and X-ray luminosity range; this means that any potential radio - X-ray luminosity correlation, $L_R \propto L_X^β$, is not well constrained ($β\approx$ 0.4-0.9, based on various studies) or understood. Here we present quasi-simultaneous Very Large Array and Swift-XRT observations of Aql X-1's 2016 outburst, with which we probe one order of magnitude fainter in radio and X-ray luminosity compared to previous studies ($6 \times 10^{34} < L_X < 3 \times 10^{35}$ erg s$^{-1}$, i.e., the intermediate to low-luminosity regime between outburst peak and quiescence). The resulting radio non-detections indicate that Aql X-1's radio emission decays more rapidly at low X-ray luminosities than previously assumed - at least during the 2016 outburst. Assuming similar behaviour between outbursts, and combining all available data, this can be modelled as a steep $β=1.17^{+0.30}_{-0.21}$ power-law index or as a sharp radio cut-off at $L_X \lesssim 5 \times 10^{35}$ erg s$^{-1}$ (given our deep radio upper limits at X-ray luminosities below this value). We discuss these results in the context of other similar studies.

preprint2019arXiv

Radio and X-ray monitoring of the accreting millisecond X-ray pulsar IGR J17591-2342 in outburst

IGR J17591$-$2342 is a new accreting millisecond X-ray pulsar (AMXP) that was recently discovered in outburst in 2018. Early observations revealed that the source's radio emission is brighter than that of any other known neutron star low-mass X-ray binary (NS-LMXB) at comparable X-ray luminosity, and assuming its likely $\gtrsim 6$ kpc distance. It is comparably radio bright to black hole LMXBs at similar X-ray luminosities. In this work, we present the results of our extensive radio and X-ray monitoring campaign of the 2018 outburst of IGR J17591$-$2342. In total we collected 10 quasi-simultaneous radio (VLA, ATCA) and X-ray (Swift-XRT) observations, which make IGR J17591$-$2342 one of the best-sampled NS-LMXBs. We use these to fit a power-law correlation index $β= 0.37^{+0.42}_{-0.40}$ between observed radio and X-ray luminosities ( $L_\mathrm{R}\propto L_\mathrm{X}^β$). However, our monitoring revealed a large scatter in IGR J17591$-$2342's radio luminosity (at a similar X-ray luminosity, $L_\mathrm{X} \sim 10^{36}$ erg s$^{-1}$, and spectral state), with $L_\mathrm{R} \sim 4 \times 10^{29}$ erg s$^{-1}$ during the first three reported observations, and up to a factor of 4 lower $L_\mathrm{R}$ during later radio observations. Nonetheless, the average radio luminosity of IGR J17591$-$2342 is still one of the highest among NS-LMXBs, and we discuss possible reasons for the wide range of radio luminosities observed in such systems during outburst. We found no evidence for radio pulsations from IGR J17591$-$2342 in our Green Bank Telescope observations performed shortly after the source returned to quiescence. Nonetheless, we cannot rule out that IGR J17591$-$2342 becomes a radio millisecond pulsar during quiescence.

preprint2018arXiv

Constraining very-high-energy and optical emission from FRB 121102 with the MAGIC telescopes

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$σ$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $σ$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

preprint2016arXiv

A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

preprint2016arXiv

A Repeating Fast Radio Burst

Fast Radio Bursts are millisecond-duration astronomical radio pulses of unknown physical origin that appear to come from extragalactic distances. Previous follow-up observations have failed to find additional bursts at the same dispersion measures (i.e. integrated column density of free electrons between source and telescope) and sky position as the original detections. The apparent non-repeating nature of the fast radio bursts has led several authors to hypothesise that they originate in cataclysmic astrophysical events. Here we report the detection of ten additional bursts from the direction of FRB121102, using the 305-m Arecibo telescope. These new bursts have dispersion measures and sky positions consistent with the original burst. This unambiguously identifies FRB121102 as repeating and demonstrates that its source survives the energetic events that cause the bursts. Additionally, the bursts from FRB121102 show a wide range of spectral shapes that appear to be predominantly intrinsic to the source and which vary on timescales of minutes or shorter. While there may be multiple physical origins for the population of fast radio bursts, the repeat bursts with high dispersion measure and variable spectra specifically seen from FRB121102 support models that propose an origin in a young, highly magnetised, extragalactic neutron star.

preprint2016arXiv

Disk reflection and a possible disk wind during a soft X-ray state in the neutron star low-mass X-ray binary 1RXS J180408.9-342058

1RXS J180408.9-342058 is a transient neutron star low-mass X-ray binary that exhibited a bright accretion outburst in 2015. We present Nustar, Swift, and Chandra observations obtained around the peak of this outburst. The source was in a soft X-ray spectral state and displayed an X-ray luminosity of Lx~(2-3)E37 (D/5.8 kpc)^2 erg cm-2 (0.5-10 keV). The Nustar data reveal a broad Fe-K emission line that we model as relativistically broadened reflection to constrain the accretion geometry. We found that the accretion disk is viewed at an inclination of i~27-35 degrees and extended close to the neutron star, down to Rin~5-7.5 gravitational radii (~11-17 km). This inner disk radius suggests that the neutron star magnetic field strength is B<2E8 G. We find a narrow absorption line in the Chandra/HEG data at an energy of ~7.64 keV with a significance of ~4.8 sigma. This feature could correspond to blue-shifted Fe xxvi and arise from an accretion disk wind, which would imply an outflow velocity of v~0.086c (~25800 km s-1). However, this would be extreme for an X-ray binary and it is unclear if a disk wind should be visible at the low inclination angle that we infer from our reflection analysis. Finally, we discuss how the X-ray and optical properties of 1RXS J180408.9-342058 are consistent with a relatively small (Porb<3 hr) binary orbit.

preprint2016arXiv

Enabling pulsar and fast transient searches using coherent dedispersion

We present an implementation of the coherent dedispersion algorithm capable of dedispersing high-time-resolution radio observations to many different dispersion measures (DMs). This approach allows the removal of the dispersive effects of the interstellar medium and enables searches for pulsed emission from pulsars and other millisecond-duration transients at low observing frequencies and/or high DMs where time broadening of the signal due to dispersive smearing would otherwise severely reduce the sensitivity. The implementation, called 'cdmt', for Coherent Dispersion Measure Trials, exploits the parallel processing capability of general-purpose graphics processing units to accelerate the computations. We describe the coherent dedispersion algorithm and detail how cdmt implements the algorithm to efficiently compute many coherent DM trials. We present the concept of a semi-coherent dedispersion search, where coherently dedispersed trials at coarsely separated DMs are subsequently incoherently dedispersed at finer steps in DM. The software is used in an ongoing LOFAR pilot survey to test the feasibility of performing semi-coherent dedispersion searches for millisecond pulsars at 135MHz. This pilot survey has led to the discovery of a radio millisecond pulsar -- the first at these low frequencies. This is the first time that such a broad and comprehensive search in DM-space has been done using coherent dedispersion, and we argue that future low-frequency pulsar searches using this approach are both scientifically compelling and feasible. Finally, we compare the performance of cdmt with other available alternatives.

preprint2016arXiv

From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release

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.

preprint2016arXiv

High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array

We report on the high-precision timing of 42 radio millisecond pulsars (MSPs) observed by the European Pulsar Timing Array (EPTA). This EPTA Data Release 1.0 extends up to mid-2014 and baselines range from 7-18 years. It forms the basis for the stochastic gravitational-wave background, anisotropic background, and continuous-wave limits recently presented by the EPTA elsewhere. The Bayesian timing analysis performed with TempoNest yields the detection of several new parameters: seven parallaxes, nine proper motions and, in the case of six binary pulsars, an apparent change of the semi-major axis. We find the NE2001 Galactic electron density model to be a better match to our parallax distances (after correction from the Lutz-Kelker bias) than the M2 and M3 models by Schnitzeler (2012). However, we measure an average uncertainty of 80\% (fractional) for NE2001, three times larger than what is typically assumed in the literature. We revisit the transverse velocity distribution for a set of 19 isolated and 57 binary MSPs and find no statistical difference between these two populations. We detect Shapiro delay in the timing residuals of PSRs J1600$-$3053 and J1918$-$0642, implying pulsar and companion masses $m_p=1.22_{-0.35}^{+0.5} \text{M}_{\odot}$, $m_c = 0.21_{-0.04}^{+0.06} \text{M}_{\odot }$ and $m_p=1.25_{-0.4}^{+0.6} \text{M}_{\odot}$, $m_c = 0.23_{-0.05}^{+0.07} \text{M}_{\odot }$, respectively. Finally, we use the measurement of the orbital period derivative to set a stringent constraint on the distance to PSRs J1012$+$5307 and J1909$-$3744, and set limits on the longitude of ascending node through the search of the annual-orbital parallax for PSRs J1600$-$3053 and J1909$-$3744.

preprint2016arXiv

Low-radio-frequency eclipses of the redback pulsar J2215+5135 observed in the image plane with LOFAR

The eclipses of certain types of binary millisecond pulsars (i.e. `black widows' and `redbacks') are often studied using high-time-resolution, `beamformed' radio observations. However, they may also be detected in images generated from interferometric data. As part of a larger imaging project to characterize the variable and transient sky at radio frequencies <200 MHz, we have blindly detected the redback system PSR J2215+5135 as a variable source of interest with the Low-Frequency Array (LOFAR). Using observations with cadences of 2 weeks - 6 months, we find preliminary evidence that the eclipse duration is frequency dependent ($\propto ν^{-0.4}$), such that the pulsar is eclipsed for longer at lower frequencies, in broad agreement with beamformed studies of other similar sources. Furthermore, the detection of the eclipses in imaging data suggests an eclipsing medium that absorbs the pulsed emission, rather than scattering it. Our study is also a demonstration of the prospects of finding pulsars in wide-field imaging surveys with the current generation of low-frequency radio telescopes.

preprint2016arXiv

New methods to constrain the radio transient rate: results from a survey of four fields with LOFAR

We report on the results of a search for radio transients between 115 and 190\,MHz with the LOw-Frequency ARray (LOFAR). Four fields have been monitored with cadences between 15 minutes and several months. A total of 151 images were obtained, giving a total survey area of 2275 deg$^2$. We analysed our data using standard LOFAR tools and searched for radio transients using the LOFAR Transients Pipeline (TraP). No credible radio transient candidate has been detected; however, we are able to set upper limits on the surface density of radio transient sources at low radio frequencies. We also show that low-frequency radio surveys are more sensitive to steep-spectrum coherent transient sources than GHz radio surveys. We used two new statistical methods to determine the upper limits on the transient surface density. One is free of assumptions on the flux distribution of the sources, while the other assumes a power-law distribution in flux and sets more stringent constraints on the transient surface density. Both of these methods provide better constraints than the approach used in previous works. The best value for the upper limit we can set for the transient surface density, using the method assuming a power-law flux distribution, is 1.3$\cdot$10$^{-3}$ deg$^{-2}$ for transients brighter than 0.3 Jy with a time-scale of 15 min, at a frequency of 150 MHz. We also calculated for the first time upper limits for the transient surface density for transients of different time-scales. We find that the results can differ by orders of magnitude from previously reported, simplified estimates.

preprint2016arXiv

Search for transient gravitational waves in coincidence with short duration radio transients during 2007-2013

We present an archival search for transient gravitational-wave bursts in coincidence with 27 single pulse triggers from Green Bank Telescope pulsar surveys, using the LIGO, Virgo and GEO interferometer network. We also discuss a check for gravitational-wave signals in coincidence with Parkes Fast Radio Bursts using similar methods. Data analyzed in these searches were collected between 2007 and 2013. Possible sources of emission of both short-duration radio signals and transient gravitational-wave emission include starquakes on neutron stars, binary coalescence of neutron stars, and cosmic string cusps. While no evidence for gravitational-wave emission in coincidence with these radio transients was found, the current analysis serves as a prototype for similar future searches using more sensitive second-generation interferometers.

preprint2016arXiv

The International Pulsar Timing Array: First Data Release

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.

preprint2016arXiv

The noise properties of 42 millisecond pulsars from the European Pulsar Timing Array and their impact on gravitational wave searches

The sensitivity of Pulsar Timing Arrays to gravitational waves depends on the noise present in the individual pulsar timing data. Noise may be either intrinsic or extrinsic to the pulsar. Intrinsic sources of noise will include rotational instabilities, for example. Extrinsic sources of noise include contributions from physical processes which are not sufficiently well modelled, for example, dispersion and scattering effects, analysis errors and instrumental instabilities. We present the results from a noise analysis for 42 millisecond pulsars (MSPs) observed with the European Pulsar Timing Array. For characterising the low-frequency, stochastic and achromatic noise component, or "timing noise", we employ two methods, based on Bayesian and frequentist statistics. For 25 MSPs, we achieve statistically significant measurements of their timing noise parameters and find that the two methods give consistent results. For the remaining 17 MSPs, we place upper limits on the timing noise amplitude at the 95% confidence level. We additionally place an upper limit on the contribution to the pulsar noise budget from errors in the reference terrestrial time standards (below 1%), and we find evidence for a noise component which is present only in the data of one of the four used telescopes. Finally, we estimate that the timing noise of individual pulsars reduces the sensitivity of this data set to an isotropic, stochastic GW background by a factor of >9.1 and by a factor of >2.3 for continuous GWs from resolvable, inspiralling supermassive black-hole binaries with circular orbits.

preprint2016arXiv

The repeating Fast Radio Burst FRB 121102: Multi-wavelength observations and additional bursts

We report on radio and X-ray observations of the only known repeating Fast Radio Burst (FRB) source, FRB 121102. We have detected six additional radio bursts from this source: five with the Green Bank Telescope at 2 GHz, and one at 1.4 GHz at the Arecibo Observatory for a total of 17 bursts from this source. All have dispersion measures consistent with a single value ($\sim559$ pc cm$^{-3}$) that is three times the predicted maximum Galactic value. The 2-GHz bursts have highly variable spectra like those at 1.4 GHz, indicating that the frequency structure seen across the individual 1.4 and 2-GHz bandpasses is part of a wideband process. X-ray observations of the FRB 121102 field with the Swift and Chandra observatories show at least one possible counterpart; however, the probability of chance superposition is high. A radio imaging observation of the field with the Jansky Very Large Array at 1.6 GHz yields a 5$σ$ upper limit of 0.3 mJy on any point-source continuum emission. This upper limit, combined with archival WISE 22-$μ$m and IPHAS H$α$ surveys, rules out the presence of an intervening Galactic HII region. We update our estimate of the FRB detection rate in the PALFA survey to be 1.1$^{+3.7}_{-1.0} \times 10^4$ FRBs sky$^{-1}$ day$^{-1}$ (95% confidence) for peak flux density at 1.4 GHz above 300 mJy. We find that the intrinsic widths of the 12 FRB 121102 bursts from Arecibo are, on average, significantly longer than the intrinsic widths of the 13 single-component FRBs detected with the Parkes telescope.

preprint2015arXiv

A LOFAR census of millisecond pulsars

We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOFAR in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths, and flux densities and compare these with higher observing frequencies. The flux-calibrated, multifrequency LOFAR pulse profiles are publicly available via the EPN Database of Pulsar Profiles. We also present average values of dispersion measures (DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered profiles, and the rest are only weakly detected. A qualitative comparison of LOFAR profiles with those at higher radio frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As previously demonstrated this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterising and correcting these systematic effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect gravitational waves.

preprint2015arXiv

Arecibo Pulsar Survey Using ALFA. IV. Mock Spectrometer Data Analysis, Survey Sensitivity, and the Discovery of 41 Pulsars

The on-going PALFA survey at the Arecibo Observatory began in 2004 and is searching for radio pulsars in the Galactic plane at 1.4 GHz. Observations since 2009 have been made with new wider-bandwidth spectrometers than were previously employed in this survey. A new data reduction pipeline has been in place since mid-2011 which consists of standard methods using dedispersion, searches for accelerated periodic sources, and search for single pulses, as well as new interference-excision strategies and candidate selection heuristics. This pipeline has been used to discover 41 pulsars, including 8 millisecond pulsars (MSPs; P < 10 ms), bringing the PALFA survey's discovery totals to 145 pulsars, including 17 MSPs, and one Fast Radio Burst (FRB). The pipeline presented here has also re-detected 188 previously known pulsars including 60 found in PALFA data by re-analyzing observations previously searched by other pipelines. A comprehensive description of the survey sensitivity, including the effect of interference and red noise, has been determined using synthetic pulsar signals with various parameters and amplitudes injected into real survey observations and subsequently recovered with the data reduction pipeline. We have confirmed that the PALFA survey achieves the sensitivity to MSPs predicted by theoretical models. However, we also find that compared to theoretical survey sensitivity models commonly used there is a degradation in sensitivity to pulsars with periods P >= 100 ms that gradually becomes up to a factor of ~10 worse for P > 4 s at DM < 150 pc/cc. This degradation of sensitivity at long periods is largely due to red noise. We find that 35 +- 3% of pulsars are missed despite being bright enough to be detected in the absence of red noise. This reduced sensitivity could have implications on estimates of the number of long-period pulsars in the Galaxy.

preprint2015arXiv

Discovery and Follow-up of Rotating Radio Transients with the Green Bank and LOFAR Telescopes

We have discovered 21 Rotating Radio Transients (RRATs) in data from the Green Bank Telescope (GBT) 350-MHz Drift-scan and the Green Bank North Celestial Cap pulsar surveys using a new candidate sifting algorithm. RRATs are pulsars with sporadic emission that are detected through their bright single pulses rather than Fourier domain searches. We have developed {\tt RRATtrap}, a single-pulse sifting algorithm that can be integrated into pulsar survey data analysis pipelines in order to find RRATs and Fast Radio Bursts. We have conducted follow-up observations of our newly discovered sources at several radio frequencies using the GBT and Low Frequency Array (LOFAR), yielding improved positions and measurements of their periods, dispersion measures, and burst rates, as well as phase-coherent timing solutions for four of them. The new RRATs have dispersion measures (DMs) ranging from 15 to 97 pc cm$^{-3}$, periods of 240 ms to 3.4 s, and estimated burst rates of 20 to 400 pulses hr$^{-1}$ at 350 MHz. We use this new sample of RRATs to perform statistical comparisons between RRATs and canonical pulsars in order to shed light on the relationship between the two populations. We find that the DM and spatial distributions of the RRATs agree with those of the pulsars found in the same survey. We find evidence that slower pulsars (i.e. $P>200$ ms) are preferentially more likely to emit bright single pulses than are faster pulsars ($P<200$ ms), although this conclusion is tentative. Our results are consistent with the proposed link between RRATs, transient pulsars, and canonical pulsars as sources in various parts of the pulse activity spectrum.

preprint2015arXiv

Einstein@Home Discovery of a PALFA Millisecond Pulsar in an Eccentric Binary Orbit

We report the discovery of the millisecond pulsar (MSP) PSR J1950+2414 ($P=4.3$ ms) in a binary system with an eccentric ($e=0.08$) 22-day orbit in Pulsar ALFA survey observations with the Arecibo telescope. Its companion star has a median mass of 0.3 $M_\odot$ and is most likely a white dwarf. Fully recycled MSPs like this one are thought to be old neutron stars spun-up by mass transfer from a companion star. This process should circularize the orbit, as is observed for the vast majority of binary MSPs, which predominantly have orbital eccentricities $e < 0.001$. However, four recently discovered binary MSPs have orbits with $0.027 < e < 0.44$; PSR J1950+2414 is the fifth such system to be discovered. The upper limits for its intrinsic spin period derivative and inferred surface magnetic field strength are comparable to those of the general MSP population. The large eccentricities are incompatible with the predictions of the standard recycling scenario: something unusual happened during their evolution. Proposed scenarios are a) initial evolution of the pulsar in a triple system which became dynamically unstable, b) origin in an exchange encounter in an environment with high stellar density, c) rotationally delayed accretion-induced collapse of a super-Chandrasekhar white dwarf, and d) dynamical interaction of the binary with a circumbinary disk. We compare the properties of all five known eccentric MSPs with the predictions of these formation channels. Future measurements of the masses and proper motion might allow us to firmly exclude some of the proposed formation scenarios.

preprint2015arXiv

Imaging Jupiter's radiation belts down to 127 MHz with LOFAR

Context. Observing Jupiter's synchrotron emission from the Earth remains today the sole method to scrutinize the distribution and dynamical behavior of the ultra energetic electrons magnetically trapped around the planet (because in-situ particle data are limited in the inner magnetosphere). Aims. We perform the first resolved and low-frequency imaging of the synchrotron emission with LOFAR at 127 MHz. The radiation comes from low energy electrons (~1-30 MeV) which map a broad region of Jupiter's inner magnetosphere. Methods (see article for complete abstract) Results. The first resolved images of Jupiter's radiation belts at 127-172 MHz are obtained along with total integrated flux densities. They are compared with previous observations at higher frequencies and show a larger extent of the synchrotron emission source (>=4 $R_J$). The asymmetry and the dynamic of east-west emission peaks are measured and the presence of a hot spot at lambda_III=230 ° $\pm$ 25 °. Spectral flux density measurements are on the low side of previous (unresolved) ones, suggesting a low-frequency turnover and/or time variations of the emission spectrum. Conclusions. LOFAR is a powerful and flexible planetary imager. The observations at 127 MHz depict an extended emission up to ~4-5 planetary radii. The similarities with high frequency results reinforce the conclusion that: i) the magnetic field morphology primarily shapes the brightness distribution of the emission and ii) the radiating electrons are likely radially and latitudinally distributed inside about 2 $R_J$. Nonetheless, the larger extent of the brightness combined with the overall lower flux density, yields new information on Jupiter's electron distribution, that may shed light on the origin and mode of transport of these particles.

preprint2015arXiv

Limits on anisotropy in the nanohertz stochastic gravitational-wave background

The paucity of observed supermassive black hole binaries (SMBHBs) may imply that the gravitational wave background (GWB) from this population is anisotropic, rendering existing analyses sub-optimal. We present the first constraints on the angular distribution of a nanohertz stochastic GWB from circular, inspiral-driven SMBHBs using the $2015$ European Pulsar Timing Array data [Desvignes et al. (in prep.)]. Our analysis of the GWB in the $\sim 2 - 90$ nHz band shows consistency with isotropy, with the strain amplitude in $l>0$ spherical harmonic multipoles $\lesssim 40\%$ of the monopole value. We expect that these more general techniques will become standard tools to probe the angular distribution of source populations.

preprint2015arXiv

Limits on Fast Radio Bursts at 145 MHz with ARTEMIS, a real-time software backend

Fast Radio Bursts (FRBs), are millisecond radio signals that exhibit dispersion larger than what the Galactic electron density can account for. We have conducted a 1446 hour survey for Fast Radio Bursts (FRBs) at 145~MHz, covering a total of 4193 sq. deg on the sky. We used the UK station of the LOFAR radio telescope -- the Rawlings Array -- , accompanied for a majority of the time by the LOFAR station at Nançay, observing the same fields at the same frequency. Our real-time search backend, ARTEMIS, utilizes graphics processing units to search for pulses with dispersion measures up to 320 cm$^{-3}$ pc. Previous derived FRB rates from surveys around 1.4~GHz, and favoured FRB interpretations, motivated this survey, despite all previous detections occurring at higher dispersion measures. We detected no new FRBs above a signal-to-noise threshold of 10, leading to the most stringent upper limit yet on the FRB event rate at these frequencies: 29 sky$^{-1}$ day$^{-1}$ for 5~ms-duration pulses above 62~Jy. The non-detection could be due to scatter-broadening, limitations on the volume and time searched, or the shape of FRB flux density spectra. Assuming the latter and that FRBs are standard candles, the non-detection is compatible with the published FRB sky rate, if their spectra follow a power law with frequency ($\propto ν^α$), with $α\gtrsim+0.1$, demonstrating a marked difference from pulsar spectra. Our results suggest that surveys at higher frequencies, including the low frequency component of the Square Kilometre Array, will have better chances to detect, estimate rates and understand the origin and properties of FRBs.

preprint2015arXiv

LOFAR discovery of a quiet emission mode in PSR B0823+26

PSR B0823+26, a 0.53-s radio pulsar, displays a host of emission phenomena over timescales of seconds to (at least) hours, including nulling, subpulse drifting, and mode-changing. Studying pulsars like PSR B0823+26 provides further insight into the relationship between these various emission phenomena and what they might teach us about pulsar magnetospheres. Here we report on the LOFAR discovery that PSR B0823+26 has a weak and sporadically emitting 'quiet' (Q) emission mode that is over 100 times weaker (on average) and has a nulling fraction forty-times greater than that of the more regularly-emitting 'bright' (B) mode. Previously, the pulsar has been undetected in the Q-mode, and was assumed to be nulling continuously. PSR B0823+26 shows a further decrease in average flux just before the transition into the B-mode, and perhaps truly turns off completely at these times. Furthermore, simultaneous observations taken with the LOFAR, Westerbork, Lovell, and Effelsberg telescopes between 110 MHz and 2.7 GHz demonstrate that the transition between the Q-mode and B-mode occurs within one single rotation of the neutron star, and that it is concurrent across the range of frequencies observed.

preprint2015arXiv

LOFAR MSSS: Detection of a low-frequency radio transient in 400 hrs of monitoring of the North Celestial Pole

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

preprint2015arXiv

LOFAR tied-array imaging and spectroscopy of solar S bursts

Context. The Sun is an active source of radio emission that is often associated with energetic phenomena ranging from nanoflares to coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), numerous millisecond duration radio bursts have been reported, such as radio spikes or solar S bursts (where S stands for short). To date, these have neither been studied extensively nor imaged because of the instrumental limitations of previous radio telescopes. Aims. Here, Low Frequency Array (LOFAR) observations were used to study the spectral and spatial characteristics of a multitude of S bursts, as well as their origin and possible emission mechanisms. Methods. We used 170 simultaneous tied-array beams for spectroscopy and imaging of S bursts. Since S bursts have short timescales and fine frequency structures, high cadence (~50 ms) tied-array images were used instead of standard interferometric imaging, that is currently limited to one image per second. Results. On 9 July 2013, over 3000 S bursts were observed over a time period of ~8 hours. S bursts were found to appear as groups of short-lived (<1 s) and narrow-bandwidth (~2.5 MHz) features, the majority drifting at ~3.5 MHz/s and a wide range of circular polarisation degrees (2-8 times more polarised than the accompanying Type III bursts). Extrapolation of the photospheric magnetic field using the potential field source surface (PFSS) model suggests that S bursts are associated with a trans-equatorial loop system that connects an active region in the southern hemisphere to a bipolar region of plage in the northern hemisphere. Conclusions. We have identified polarised, short-lived solar radio bursts that have never been imaged before. They are observed at a height and frequency range where plasma emission is the dominant emission mechanism, however they possess some of the characteristics of electron-cyclotron maser emission.

preprint2015arXiv

Probing Atmospheric Electric Fields in Thunderstorms through Radio Emission from Cosmic-Ray-Induced Air Showers

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields.

preprint2015arXiv

Properties and Evolution of the Redback Millisecond Pulsar Binary PSR J2129-0429

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.

preprint2015arXiv

Pulsar polarisation below 200 MHz: Average profiles and propagation effects

We present the highest-quality polarisation profiles to date of 16 non-recycled pulsars and four millisecond pulsars, observed below 200 MHz with the LOFAR high-band antennas. Based on the observed profiles, we perform an initial investigation of expected observational effects resulting from the propagation of polarised emission in the pulsar magnetosphere and the interstellar medium. The predictions of magnetospheric birefringence in pulsars have been tested using spectra of the pulse width and fractional polarisation from multifrequency data. The derived spectra offer only partial support for the expected effects of birefringence on the polarisation properties, with only about half of our sample being consistent with the model's predictions. It is noted that for some pulsars these measurements are contaminated by the effects of interstellar scattering. For a number of pulsars in our sample, we have observed significant variations in the amount of Faraday rotation as a function of pulse phase, which is possibly an artefact of scattering. These variations are typically two orders of magnitude smaller than that observed at 1400 MHz by Noutsos et al. (2009), for a different sample of southern pulsars. In this paper we present a possible explanation for the difference in magnitude of this effect between the two frequencies, based on scattering. Finally, we have estimated the magnetospheric emission heights of low-frequency radiation from four pulsars, based on the phase lags between the flux-density and the PA profiles, and the theoretical framework of Blaskiewicz, Cordes & Wasserman (1991). These estimates yielded heights of a few hundred km; at least for PSR B1133+16, this is consistent with emission heights derived based on radius-to-frequency mapping, but is up to a few times larger than the recent upper limit based on pulsar timing.

preprint2015arXiv

Pulsar Wind Nebulae in the SKA era

Neutron stars lose the bulk of their rotational energy in the form of a pulsar wind: an ultra-relativistic outflow of predominantly electrons and positrons. This pulsar wind significantly impacts the environment and possible binary companion of the neutron star, and studying the resultant pulsar wind nebulae is critical for understanding the formation of neutron stars and millisecond pulsars, the physics of the neutron star magnetosphere, the acceleration of leptons up to PeV energies, and how these particles impact the interstellar medium. With the SKA1 and the SKA2, it could be possible to study literally hundreds of PWNe in detail, critical for understanding the many open questions in the topics listed above.

preprint2015arXiv

Radio imaging observations of PSR J1023+0038 in an LMXB state

The transitional millisecond pulsar binary system PSR J1023+0038 re-entered an accreting state in 2013 June, in which it bears many similarities to low-mass X-ray binaries (LMXBs) in quiescence or near-quiescence. At a distance of just 1.37 kpc, PSR J1023+0038 offers an unsurpassed ability to study low-level accretion onto a highly-magnetized compact object. We have monitored PSR J1023+0038 intensively using radio imaging with the Karl G. Jansky Very Large Array, the European VLBI Network and LOFAR, seeing rapidly variable, flat spectrum emission that persists over a period of 6 months. The flat spectrum and variability are indicative of synchrotron emission originating in an outflow from the system, most likely in the form of a compact, partially self-absorbed jet, as is seen in LMXBs at higher accretion rates. The radio brightness, however, greatly exceeds extrapolations made from observations of more vigorously accreting neutron star LMXB systems. We postulate that PSR J1023+0038 is undergoing radiatively inefficient 'propeller-mode' accretion, with the jet carrying away a dominant fraction of the liberated accretion luminosity. We confirm that the enhanced gamma-ray emission seen in PSR J1023+0038 since it re-entered an accreting state has been maintained; the increased gamma-ray emission in this state can also potentially be associated with propeller-mode accretion. Similar accretion modes can be invoked to explain the radio and X-ray properties of the other two known transitional millisecond pulsar systems XSS J12270-4859 and PSR J1824-2452I (M28I), suggesting that radiatively inefficient accretion may be an ubiquitous phenomenon amongst (at least one class of) neutron star binaries at low accretion rates.

preprint2015arXiv

Radio Timing and Optical Photometry of the Black Widow System PSR J1953+1846A in the Globular Cluster M71

We report on the determination of the astrometric, spin and orbital parameters for PSR J1953+1846A, a "black widow" binary millisecond pulsar in the globular cluster M71. By using the accurate position and orbital parameters obtained from radio timing, we identified the optical companion in ACS/Hubble Space Telescope images. It turns out to be a faint (m_F606W>=24, m_F814W>=23) and variable star located at only ~0.06" from the pulsar timing position. The light curve shows a maximum at the pulsar inferior conjunction and a minimum at the pulsar superior conjunction, thus confirming the association with the system. The shape of the optical modulation suggests that the companion star is heated, likely by the pulsar wind. The comparison with the X-ray light curve possibly suggests the presence of an intra-binary shock due to the interaction between the pulsar wind and the material released by the companion. This is the second identification (after COM-M5C) of an optical companion to a black widow pulsar in a globular cluster. Interestingly, the two companions show a similar light curve and share the same position in the color magnitude diagram.

preprint2015arXiv

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

preprint2015arXiv

Timing of five millisecond pulsars discovered in the PALFA survey

We present the discovery of five millisecond pulsars (MSPs) from the PALFA Galactic plane survey using Arecibo. Four of these (PSRs J0557+1551, J1850+0244, J1902+0300, and J1943+2210) are binary pulsars whose companions are likely white dwarfs, and one (PSR J1905+0453) is isolated. Phase-coherent timing solutions, ranging from $\sim$1 to $\sim$3 years in length, and based on observations from the Jodrell Bank and Arecibo telescopes, provide precise determinations of spin, orbital, and astrometric parameters. All five pulsars have large dispersion measures ($>100$ pc cm$^{-3}$, within the top 20% of all known Galactic field MSPs) and are faint (1.4 GHz flux density < 0.1 mJy, within the faintest 5% of all known Galactic field MSPs), illustrating PALFA's ability to find increasingly faint, distant MSPs in the Galactic plane. In particular, PSR J1850+0244 has a dispersion measure of 540 pc cm$^{-3}$, the highest of all known MSPs. Such distant, faint MSPs are important input for accurately modeling the total Galactic MSP population.

preprint2015arXiv

Wide-Band, Low-Frequency Pulse Profiles of 100 Radio Pulsars with LOFAR

LOFAR offers the unique capability of observing pulsars across the 10-240 MHz frequency range with a fractional bandwidth of roughly 50%. This spectral range is well-suited for studying the frequency evolution of pulse profile morphology caused by both intrinsic and extrinsic effects: such as changing emission altitude in the pulsar magnetosphere or scatter broadening by the interstellar medium, respectively. The magnitude of most of these effects increases rapidly towards low frequencies. LOFAR can thus address a number of open questions about the nature of radio pulsar emission and its propagation through the interstellar medium. We present the average pulse profiles of 100 pulsars observed in the two LOFAR frequency bands: High Band (120-167 MHz, 100 profiles) and Low Band (15-62 MHz, 26 profiles). We compare them with Westerbork Synthesis Radio Telescope (WSRT) and Lovell Telescope observations at higher frequencies (350 and1400 MHz) in order to study the profile evolution. The profiles are aligned in absolute phase by folding with a new set of timing solutions from the Lovell Telescope, which we present along with precise dispersion measures obtained with LOFAR. We find that the profile evolution with decreasing radio frequency does not follow a specific trend but, depending on the geometry of the pulsar, new components can enter into, or be hidden from, view. Nonetheless, in general our observations confirm the widening of pulsar profiles at low frequencies, as expected from radius-to-frequency mapping or birefringence theories. We offer this catalog of low-frequency pulsar profiles in a user friendly way via the EPN Database of Pulsar Profiles (http://www.epta.eu.org/epndb/).

preprint2014arXiv

A 24-Hour Global Campaign To Assess Precision Timing of the Millisecond Pulsar J1713+0747

The radio millisecond pulsar J1713+0747 is regarded as one of the highest-precision clocks in the sky, and is regularly timed for the purpose of detecting gravitational waves. The International Pulsar Timing Array collaboration undertook a 24-hour global observation of PSR J1713+0747 in an effort to better quantify sources of timing noise in this pulsar, particularly on intermediate (1 - 24 hr) timescales. We observed the pulsar continuously over 24 hr with the Arecibo, Effelsberg, GMRT, Green Bank, LOFAR, Lovell, Nancay, Parkes, and WSRT radio telescopes. The combined pulse times-of-arrival presented here provide an estimate of what sources of timing noise, excluding DM variations, would be present as compared to an idealized root-N improvement in timing precision, where N is the number of pulses analyzed. In the case of this particular pulsar, we find that intrinsic pulse phase jitter dominates arrival time precision when the S/N of single pulses exceeds unity, as measured using the eight telescopes that observed at L-band/1.4 GHz. We present first results of specific phenomena probed on the unusually long timescale (for a single continuous observing session) of tens of hours, in particular interstellar scintillation, and discuss the degree to which scintillation and profile evolution affect precision timing. This paper presents the data set as a basis for future, deeper studies.

preprint2014arXiv

A Cosmic Census of Radio Pulsars with the SKA

The Square Kilometre Array (SKA) will make ground breaking discoveries in pulsar science. In this chapter we outline the SKA surveys for new pulsars, as well as how we will perform the necessary follow-up timing observations. The SKA's wide field-of-view, high sensitivity, multi-beaming and sub-arraying capabilities, coupled with advanced pulsar search backends, will result in the discovery of a large population of pulsars. These will enable the SKA's pulsar science goals (tests of General Relativity with pulsar binary systems, investigating black hole theorems with pulsar-black hole binaries, and direct detection of gravitational waves in a pulsar timing array). Using SKA1-MID and SKA1-LOW we will survey the Milky Way to unprecedented depth, increasing the number of known pulsars by more than an order of magnitude. SKA2 will potentially find all the Galactic radio-emitting pulsars in the SKA sky which are beamed in our direction. This will give a clear picture of the birth properties of pulsars and of the gravitational potential, magnetic field structure and interstellar matter content of the Galaxy. Targeted searches will enable detection of exotic systems, such as the ~1000 pulsars we infer to be closely orbiting Sgr A*, the supermassive black hole in the Galactic Centre. In addition, the SKA's sensitivity will be sufficient to detect pulsars in local group galaxies. To derive the spin characteristics of the discoveries we will perform live searches, and use sub-arraying and dynamic scheduling to time pulsars as soon as they are discovered, while simultaneously continuing survey observations. The large projected number of discoveries suggests that we will uncover currently unknown rare systems that can be exploited to push the boundaries of our understanding of astrophysics and provide tools for testing physics, as has been done by the pulsar community in the past.

preprint2014arXiv

A millisecond pulsar in a stellar triple system

Gravitationally bound three-body systems have been studied for hundreds of years and are common in our Galaxy. They show complex orbital interactions, which can constrain the compositions, masses, and interior structures of the bodies and test theories of gravity, if sufficiently precise measurements are available. A triple system containing a radio pulsar could provide such measurements, but the only previously known such system, B1620-26 (with a millisecond pulsar, a white dwarf, and a planetary-mass object in an orbit of several decades), shows only weak interactions. Here we report precision timing and multi-wavelength observations of PSR J0337+1715, a millisecond pulsar in a hierarchical triple system with two other stars. Strong gravitational interactions are apparent and provide the masses of the pulsar (1.4378(13) Msun, where Msun is the solar mass and the parentheses contain the uncertainty in the final decimal places) and the two white dwarf companions (0.19751(15) Msun and 0.4101(3) Msun), as well as the inclinations of the orbits (both approximately 39.2 degrees). The unexpectedly coplanar and nearly circular orbits indicate a complex and exotic evolutionary past that differs from those of known stellar systems. The gravitational field of the outer white dwarf strongly accelerates the inner binary containing the neutron star, and the system will thus provide an ideal laboratory in which to test the strong equivalence principle of general relativity.

preprint2014arXiv

A state change in the low-mass X-ray binary XSS J12270-4859

Millisecond radio pulsars acquire their rapid rotation rates through mass and angular momentum transfer in a low-mass X-ray binary system. Recent studies of PSR J1824-2452I and PSR J1023+0038 have observationally demonstrated this link, and they have also shown that such systems can repeatedly transition back-and-forth between the radio millisecond pulsar and low-mass X-ray binary states. This also suggests that a fraction of such systems are not newly born radio millisecond pulsars but are rather suspended in a back-and-forth state switching phase, perhaps for giga-years. XSS J12270-4859 has been previously suggested to be a low-mass X-ray binary, and until recently the only such system to be seen at MeV-GeV energies. We present radio, optical and X-ray observations that offer compelling evidence that XSS J12270-4859 is a low-mass X-ray binary which transitioned to a radio millisecond pulsar state between 2012 November 14 and 2012 December 21. Though radio pulsations remain to be detected, we use optical and X-ray photometry/spectroscopy to show that the system has undergone a sudden dimming and no longer shows evidence for an accretion disk. The optical observations constrain the orbital period to 6.913+-0.002 hr.

preprint2014arXiv

Arecibo Pulsar Survey Using ALFA. III. Precursor Survey and Population Synthesis

The Pulsar Arecibo L-band Feed Array (PALFA) Survey uses the ALFA 7-beam receiver to search both inner and outer Galactic sectors visible from Arecibo ($32^{\circ}\lesssim \ell \lesssim 77^{\circ}$ and $168^{\circ}\lesssim \ell \lesssim 214^{\circ}$) close to the Galactic plane ($|b|\lesssim5^{\circ}$) for pulsars. In this paper we detail a precursor survey of this region with PALFA, which observed a subset of the full region (slightly more restrictive in $\ell$ and $|b|\lesssim1^{\circ}$) and detected 45 pulsars. For both Galactic millisecond and normal pulsar populations, we compare the survey's detections with simulations to model these populations and, in particular, to estimate the number of observable pulsars in the Galaxy. We place 95\% confidence intervals of 82,000 to 143,000 on the number of detectable normal pulsars and 9,000 to 100,000 on the number of detectable millisecond pulsars in the Galactic disk. These are consistent with previous estimates. Given the most likely population size in each case (107,000 and 15,000 for normal and millisecond pulsars, respectively) we extend survey detection simulations to predict that, when complete, the full PALFA survey should have detected $1,000\substack{+330 \\ -230}$ normal pulsars and $30\substack{+200 \\ -20}$ millisecond pulsars. Identical estimation techniques predict that $490\substack{+160 \\ -115}$ normal pulsars and $12\substack{+70 \\ -5}$ millisecond pulsars would be detected by the beginning of 2014; at the time, the PALFA survey had detected 283 normal pulsars and 31 millisecond pulsars, respectively.

preprint2014arXiv

Discovery of Carbon Radio Recombination Lines in absorption towards Cygnus~A

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $α$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $α$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$σ$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$σ$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

preprint2014arXiv

Fast Radio Burst Discovered in the Arecibo Pulsar ALFA Survey

Recent work has exploited pulsar survey data to identify temporally isolated, millisecond-duration radio bursts with large dispersion measures (DMs). These bursts have been interpreted as arising from a population of extragalactic sources, in which case they would provide unprecedented opportunities for probing the intergalactic medium; they may also be linked to new source classes. Until now, however, all so-called fast radio bursts (FRBs) have been detected with the Parkes radio telescope and its 13-beam receiver, casting some concern about the astrophysical nature of these signals. Here we present FRB 121102, the first FRB discovery from a geographic location other than Parkes. FRB 121102 was found in the Galactic anti-center region in the 1.4-GHz Pulsar ALFA survey with the Arecibo Observatory with a DM = 557.4 $\pm$ 3 pc cm$^{-3}$, pulse width of $3\; \pm 0.5$ ms, and no evidence of interstellar scattering. The observed delay of the signal arrival time with frequency agrees precisely with the expectation of dispersion through an ionized medium. Despite its low Galactic latitude ($b = -0.2^{\circ}$), the burst has three times the maximum Galactic DM expected along this particular line-of-sight, suggesting an extragalactic origin. A peculiar aspect of the signal is an inverted spectrum; we interpret this as a consequence of being detected in a sidelobe of the ALFA receiver. FRB 121102's brightness, duration, and the inferred event rate are all consistent with the properties of the previously detected Parkes bursts.

preprint2014arXiv

Gravitational waves from known pulsars: results from the initial detector era

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.

preprint2014arXiv

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

preprint2014arXiv

LOFAR tied-array imaging of Type III solar radio bursts

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

preprint2014arXiv

Pulsars in Globular Clusters with the SKA

Globular clusters are highly efficient radio pulsar factories. These pulsars can be used as precision probes of the clusters' structure, gas content, magnetic field, and formation history; some of them are also highly interesting in their own right because they probe exotic stellar evolution scenarios as well as the physics of dense matter, accretion, and gravity. Deep searches with SKA1-MID and SKA1-LOW will plausibly double to triple the known population. Such searches will only require one to a few tied-array beams, and can be done during early commissioning of the telescope - before an all-sky pulsar survey using hundreds to thousands of tied-array beams is feasible. With SKA2 it will be possible to observe most of the active radio pulsars within a large fraction of the Galactic globular clusters, an estimated population of 600 - 3700 observable pulsars (those beamed towards us). This rivals the total population of millisecond pulsars that can be found in the Galactic field; fully characterizing it will provide the best-possible physical laboratories as well as a rich dynamical history of the Galactic globular cluster system.

preprint2014arXiv

Radio timing and optical photometry of the black widow system PSR J1518+0204C in the globular cluster M5

We report on the determination of astrometric, spin and orbital parameters for PSR J1518+0204C, a "black widow" binary millisecond pulsar in the globular cluster M5. The accurate position and orbital parameters obtained from radio timing allowed us to search for the optical companion. By using WFC3/HST images we identified a very faint variable star (m_F390W > 24.8, m_F606W > 24.3, m_F814W > 23.1) located at only 0.25" from the pulsar's timing position. Due to its strong variability, this star is visible only in a sub-sample of images. However, the light curve obtained folding the available data with the orbital parameters of the pulsar shows a maximum at the pulsar inferior conjunction and a possible minimum at the pulsar superior conjunction. Furthermore, the shape of the optical modulation indicates a heating process possibly due to the pulsar wind. This is the first identification of an optical companion to a black widow pulsar in the dense stellar environment of a globular cluster.

preprint2014arXiv

Swings between rotation and accretion power in a millisecond binary pulsar

It is thought that neutron stars in low-mass binary systems can accrete matter and angular momentum from the companion star and be spun-up to millisecond rotational periods. During the accretion stage, the system is called a low-mass X-ray binary, and bright X-ray emission is observed. When the rate of mass transfer decreases in the later evolutionary stages, these binaries host a radio millisecond pulsar whose emission is powered by the neutron star's rotating magnetic field. This evolutionary model is supported by the detection of millisecond X-ray pulsations from several accreting neutron stars and also by the evidence for a past accretion disc in a rotation-powered millisecond pulsar. It has been proposed that a rotation-powered pulsar may temporarily switch on during periods of low mass inflow in some such systems. Only indirect evidence for this transition has hitherto been observed. Here we report observations of accretion-powered, millisecond X-ray pulsations from a neutron star previously seen as a rotation-powered radio pulsar. Within a few days after a month-long X-ray outburst, radio pulses were again detected. This not only shows the evolutionary link between accretion and rotation-powered millisecond pulsars, but also that some systems can swing between the two states on very short timescales.

preprint2014arXiv

The Green Bank Northern Celestial Cap Pulsar Survey - I: Survey Description, Data Analysis, and Initial Results

We describe an ongoing search for pulsars and dispersed pulses of radio emission, such as those from rotating radio transients (RRATs) and fast radio bursts (FRBs), at 350 MHz using the Green Bank Telescope. With the Green Bank Ultimate Pulsar Processing Instrument, we record 100 MHz of bandwidth divided into 4,096 channels every 81.92 $μs$. This survey will cover the entire sky visible to the Green Bank Telescope ($δ> -40^\circ$, or 82% of the sky) and outside of the Galactic Plane will be sensitive enough to detect slow pulsars and low dispersion measure ($<$30 $\mathrm{pc\,cm^{-3}}$) millisecond pulsars (MSPs) with a 0.08 duty cycle down to 1.1 mJy. For pulsars with a spectral index of $-$1.6, we will be 2.5 times more sensitive than previous and ongoing surveys over much of our survey region. Here we describe the survey, the data analysis pipeline, initial discovery parameters for 62 pulsars, and timing solutions for 5 new pulsars. PSR J0214$+$5222 is an MSP in a long-period (512 days) orbit and has an optical counterpart identified in archival data. PSR J0636$+$5129 is an MSP in a very short-period (96 minutes) orbit with a very low mass companion (8 $M_\mathrm{J}$). PSR J0645$+$5158 is an isolated MSP with a timing residual RMS of 500 ns and has been added to pulsar timing array experiments. PSR J1434$+$7257 is an isolated, intermediate-period pulsar that has been partially recycled. PSR J1816$+$4510 is an eclipsing MSP in a short-period orbit (8.7 hours) and may have recently completed its spin-up phase.

preprint2013arXiv

A New Accretion Disk Around the Missing Link Binary System PSR J1023+0038

PSR J1023+0038 is an exceptional system for understanding how slowly rotating neutron stars are spun up to millisecond rotational periods through accretion from a companion star. Observed as a radio pulsar from 2007-2013, optical data showed that the system had an accretion disk in 2000/2001. Starting at the end of 2013 June, the radio pulsar has become undetectable, suggesting a return to the previous accretion-disk state, where the system more closely resembles an X-ray binary. In this Letter we report the first targeted X-ray observations ever performed of the active phase and complement them with UV/Optical and radio observations collected in 2013 October. We find strong evidence that indeed an accretion disk has recently formed in the system and we report the detection of fast X-ray changes spanning about two orders of magnitude in luminosity. No radio pulsations are seen during low flux states in the X-ray light-curve or at any other times.

preprint2013arXiv

A state change in the missing link binary pulsar system PSR J1023+0038

We present radio, X-ray, and gamma-ray observations which reveal that the binary millisecond pulsar / low-mass X-ray binary transition system PSR J1023+0038 has undergone a transformation in state. Whereas until recently the system harbored a bright millisecond radio pulsar, the radio pulsations at frequencies between 300 to 5000 MHz have now become undetectable. Concurrent with this radio disappearance, the gamma-ray flux of the system has quintupled. We conclude that, though the radio pulsar is currently not detectable, the pulsar mechanism is still active and the pulsar wind, as well as a newly formed accretion disk, are together providing the necessary conditions to create the gamma-ray increase. The system is the first example of a transient, compact, low-mass gamma-ray binary and will continue to provide an exceptional test bed for better understanding the formation of millisecond pulsars as well as accretion onto neutron stars in general.

preprint2013arXiv

Calibrating High-Precision Faraday Rotation Measurements for LOFAR and the Next Generation of Low-Frequency Radio Telescopes

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

preprint2013arXiv

Detecting cosmic rays with the LOFAR radio telescope

The low frequency array (LOFAR), is the first radio telescope designed with the capability to measure radio emission from cosmic-ray induced air showers in parallel with interferometric observations. In the first $\sim 2\,\mathrm{years}$ of observing, 405 cosmic-ray events in the energy range of $10^{16} - 10^{18}\,\mathrm{eV}$ have been detected in the band from $30 - 80\,\mathrm{MHz}$. Each of these air showers is registered with up to $\sim1000$ independent antennas resulting in measurements of the radio emission with unprecedented detail. This article describes the dataset, as well as the analysis pipeline, and serves as a reference for future papers based on these data. All steps necessary to achieve a full reconstruction of the electric field at every antenna position are explained, including removal of radio frequency interference, correcting for the antenna response and identification of the pulsed signal.

preprint2013arXiv

Differential Frequency-dependent Delay from the Pulsar Magnetosphere

Some radio pulsars show clear drifting subpulses, in which subpulses are seen to drift in pulse longitude in a systematic pattern. Here we examine how the drifting subpulses of PSR B0809+74 evolve with time and observing frequency. We show that the subpulse period (P3) is constant on timescales of days, months and years, and between 14-5100 MHz. Despite this, the shapes of the driftbands change radically with frequency. Previous studies have concluded that, while the subpulses appear to move through the pulse window approximately linearly at low frequencies (< 500 MHz), a discrete step of 180 degrees in subpulse phase is observed at higher frequencies (> 820 MHz) near to the peak of the average pulse profile. We use LOFAR, GMRT, GBT, WSRT and Effelsberg 100-m data to explore the frequency-dependence of this phase step. We show that the size of the subpulse phase step increases gradually, and is observable even at low frequencies. We attribute the subpulse phase step to the presence of two separate driftbands, whose relative arrival times vary with frequency - one driftband arriving 30 pulses earlier at 20 MHz than it does at 1380 MHz, whilst the other arrives simultaneously at all frequencies. The drifting pattern which is observed here cannot be explained by either the rotating carousel model or the surface oscillation model, and could provide new insight into the physical processes happening within the pulsar magnetosphere.

preprint2013arXiv

Discovery of the Optical Counterparts to Four Energetic Fermi Millisecond Pulsars

In the last few years, over 43 millisecond radio pulsars have been discovered by targeted searches of unidentified gamma-ray sources found by the Fermi Gamma-Ray Space Telescope. A large fraction of these millisecond pulsars are in compact binaries with low-mass companions. These systems often show eclipses of the pulsar signal and are commonly known as black widows and redbacks because the pulsar is gradually destroying its companion. In this paper, we report on the optical discovery of four strongly irradiated millisecond pulsar companions. All four sources show modulations of their color and luminosity at the known orbital periods from radio timing. Light curve modelling of our exploratory data shows that the equilibrium temperature reached on the companion's dayside with respect to their nightside is consistent with about 10-30% of the available spin-down energy from the pulsar being reprocessed to increase the companion's dayside temperature. This value compares well with the range observed in other irradiated pulsar binaries and offers insights about the energetics of the pulsar wind and the production of gamma-ray emission. In addition, this provides a simple way of estimating the brightness of irradiated pulsar companions given the pulsar spin-down luminosity. Our analysis also suggests that two of the four new irradiated pulsar companions are only partially filling their Roche lobe. Some of these sources are relatively bright and represent good targets for spectroscopic follow-up. These measurements could enable, among other things, mass determination of the neutron stars in these systems.

preprint2013arXiv

IGR J19308+0530: Roche lobe overflow on to a compact object from a donor 1.8 times as massive

We present phase-resolved spectroscopy and photometry of the optical counterpart to the X-ray binary IGR J19308+0530. Ellipsoidal modulations in the light curve show that the F-type companion star in the system is Roche lobe filling. The optical spectra are dominated by absorption features from the donor star, with ~10-20 per cent disc contribution to the optical continuum. We measure an orbital period of 14.662+-0.001 h, a radial velocity semiamplitude for the companion star of K2 = 91.4+-1.4 kms-1 and a rotational broadening of vsini = 108.9+-0.6 kms-1. From K2 and vsini, given that the donor star is filling its Roche lobe, we derive a mass ratio of q = M2/M1=1.78+-0.04, which is typically considered to be too large for stable Roche lobe overflow. Our observations support an inclination of ~50 degrees. The accretor in IGR J19308+0530 is most likely a white dwarf, although a neutron star cannot entirely be excluded.

preprint2013arXiv

Initial deep LOFAR observations of Epoch of Reionization windows: I. The North Celestial Pole

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

preprint2013arXiv

LOFAR detections of low-frequency radio recombination lines towards Cassiopeia A

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

preprint2013arXiv

LOFAR: The LOw-Frequency ARray

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

preprint2013arXiv

PEACE: Pulsar Evaluation Algorithm for Candidate Extraction -- A software package for post-analysis processing of pulsar survey candidates

Modern radio pulsar surveys produce a large volume of prospective candidates, the majority of which are polluted by human-created radio frequency interference or other forms of noise. Typically, large numbers of candidates need to be visually inspected in order to determine if they are real pulsars. This process can be labor intensive. In this paper, we introduce an algorithm called PEACE (Pulsar Evaluation Algorithm for Candidate Extraction) which improves the efficiency of identifying pulsar signals. The algorithm ranks the candidates based on a score function. Unlike popular machine-learning based algorithms, no prior training data sets are required. This algorithm has been applied to data from several large-scale radio pulsar surveys. Using the human-based ranking results generated by students in the Arecibo Remote Command enter programme, the statistical performance of PEACE was evaluated. It was found that PEACE ranked 68% of the student-identified pulsars within the top 0.17% of sorted candidates, 95% within the top 0.34%, and 100% within the top 3.7%. This clearly demonstrates that PEACE significantly increases the pulsar identification rate by a factor of about 50 to 1000. To date, PEACE has been directly responsible for the discovery of 47 new pulsars, 5 of which are millisecond pulsars that may be useful for pulsar timing based gravitational-wave detection projects.

preprint2013arXiv

Searching for pulsars using image pattern recognition

In this paper, we present a novel artificial intelligence (AI) program that identifies pulsars from recent surveys using image pattern recognition with deep neural nets---the PICS (Pulsar Image-based Classification System) AI. The AI mimics human experts and distinguishes pulsars from noise and interferences by looking for patterns from candidate. The information from each pulsar candidate is synthesized in four diagnostic plots, which consist of up to thousands pixel of image data. The AI takes these data from each candidate as its input and uses thousands of such candidates to train its ~9000 neurons. Different from other pulsar selection programs which use pre-designed patterns, the PICS AI teaches itself the salient features of different pulsars from a set of human-labeled candidates through machine learning. The deep neural networks in this AI system grant it superior ability in recognizing various types of pulsars as well as their harmonic signals. The trained AI's performance has been validated with a large set of candidates different from the training set. In this completely independent test, PICS ranked 264 out of 277 pulsar-related candidates, including all 56 previously known pulsars, to the top 961 (1%) of 90008 test candidates, missing only 13 harmonics. The first non-pulsar candidate appears at rank 187, following 45 pulsars and 141 harmonics. In other words, 100% of the pulsars were ranked in the top 1% of all candidates, while 80% were ranked higher than any noise or interference. The performance of this system can be improved over time as more training data are accumulated. This AI system has been integrated into the PALFA survey pipeline and has discovered six new pulsars to date.

preprint2013arXiv

Studying Galactic interstellar turbulence through fluctuations in synchrotron emission: First LOFAR Galactic foreground detection

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

preprint2013arXiv

Synchronous X-ray and Radio Mode Switches: a Rapid Global Transformation of the Pulsar Magnetosphere

Pulsars emit low-frequency radio waves through to high-energy gamma-rays that are generated anywhere from the surface out to the edges of the magnetosphere. Detecting correlated mode changes in the multi-wavelength emission is therefore key to understanding the physical relationship between these emission sites. Through simultaneous observations, we have detected synchronous switching in the radio and X-ray emission properties of PSR B0943+10. When the pulsar is in a sustained radio 'bright' mode, the X-rays show only an un-pulsed, non-thermal component. Conversely, when the pulsar is in a radio 'quiet' mode, the X-ray luminosity more than doubles and a 100%-pulsed thermal component is observed along with the non-thermal component. This indicates rapid, global changes to the conditions in the magnetosphere, which challenge all proposed pulsar emission theories.

preprint2013arXiv

The Einstein@Home search for radio pulsars and PSR J2007+2722 discovery

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 193 countries, to search for new neutron stars using data from electromagnetic and gravitational-wave detectors. This paper presents a detailed description of the search for new radio pulsars using Pulsar ALFA survey data from the Arecibo Observatory. The enormous computing power allows this search to cover a new region of parameter space; it can detect pulsars in binary systems with orbital periods as short as 11 minutes. We also describe the first Einstein@Home discovery, the 40.8 Hz isolated pulsar PSR J2007+2722, and provide a full timing model. PSR J2007+2722's pulse profile is remarkably wide with emission over almost the entire spin period. This neutron star is most likely a disrupted recycled pulsar, about as old as its characteristic spin-down age of 404 Myr. However there is a small chance that it was born recently, with a low magnetic field. If so, upper limits on the X-ray flux suggest but can not prove that PSR J2007+2722 is at least ~ 100 kyr old. In the future, we expect that the massive computing power provided by volunteers should enable many additional radio pulsar discoveries.

preprint2013arXiv

The missing GeV γ-ray binary: Searching for HESS J0632+057 with Fermi-LAT

The very high energy (VHE; >100 GeV) source HESS J0632+057 has been recently confirmed as a γ-ray binary, a subclass of the high mass X-ray binary (HMXB) population, through the detection of an orbital period of 321 days. We performed a deep search for the emission of HESS J0632+057 in the GeV energy range using data from the Fermi Large Area Telescope (LAT). The analysis was challenging due to the source being located in close proximity to the bright γ-ray pulsar PSR J0633+0632 and lying in a crowded region of the Galactic plane where there is prominent diffuse emission. We formulated a Bayesian block algorithm adapted to work with weighted photon counts, in order to define the off-pulse phases of PSR J0633+0632. A detailed spectral-spatial model of a 5 deg circular region centred on the known location of HESS J0632+057 was generated to accurately model the LAT data. No significant emission from the location of HESS J0632+057 was detected in the 0.1-100 GeV energy range integrating over ~3.5 years of data; with a 95% flux upper limit of F_{0.1-100 GeV} < 3 x 10-8 ph cm-2 s-1. A search for emission over different phases of the orbit also yielded no significant detection. A search for source emission on shorter timescales (days--months) did not yield any significant detections. We also report the results of a search for radio pulsations using the 100-m Green Bank Telescope (GBT). No periodic signals or individual dispersed bursts of a likely astronomical origin were detected. We estimated the flux density limit of < 90/40 μJy at 2/9 GHz. The LAT flux upper limits combined with the detection of HESS J0632+057 in the 136-400 TeV energy band by the MAGIC collaboration imply that the VHE spectrum must turn over at energies <136 GeV placing constraints on any theoretical models invoked to explain the γ-ray emission.

preprint2013arXiv

Timing and Interstellar Scattering of Thirty-five Distant Pulsars Discovered in the PALFA Survey

We have made extensive observations of 35 distant slow (non-recycled) pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Timing observations of these pulsars over several years at Arecibo Observatory and Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation properties. Despite being a relatively distant population, these pulsars have properties that mirror those of the previously known pulsar population. Many of the sources exhibit timing noise, and one underwent a small glitch. We have used multifrequency data to measure the interstellar scattering properties of these pulsars. We find scattering to be higher than predicted along some lines of sight, particularly in the Cygnus region. Lastly, we present XMM-Newton and Chandra observations of the youngest and most energetic of the pulsars, J1856+0245, which has previously been associated with the GeV-TeV pulsar wind nebula HESS J1857+026.

preprint2012arXiv

A search for radio pulsars and fast transients in M31 using the WSRT

We present the results of the most sensitive and comprehensive survey yet undertaken for radio pulsars and fast transients in the Andromeda galaxy (M31) and its satellites, using the Westerbork Synthesis Radio Telescope (WSRT) at a central frequency of 328 MHz. We used the WSRT in a special configuration called 8gr8 (eight-grate) mode, which provides a large instantaneous field-of-view, about 5 square degrees per pointing, with good sensitivity, long dwell times (up to 8 hours per pointing), and good spatial resolution (a few arc minutes) for locating sources. We have searched for both periodicities and single pulses in our data, aiming to detect bright, persistent radio pulsars and rotating radio transients (RRATs) of either Galactic or extragalactic origin. Our searches did not reveal any confirmed periodic signals or bright single bursts from (potentially) cosmological distances. However, we do report the detection of several single pulse events, some repeating at the same dispersion measure, which could potentially originate from neutron stars in M31. One in particular was seen multiple times, including a burst of six pulses in 2000 seconds, at a dispersion measure of 54.7 pc cm^-3, which potentially places the origin of this source outside of our Galaxy. Our results are compared to a range of hypothetical populations of pulsars and RRATs in M31 and allow us to constrain the luminosity function of pulsars in M31. They also show that, unless the pulsar population in M31 is much dimmer than in our Galaxy, there is no need to invoke any violation of the inverse square law of the distance for pulsar fluxes.

preprint2012arXiv

Discovery of the Optical/Ultraviolet/Gamma-ray Counterpart to the Eclipsing Millisecond Pulsar J1816+4510

The energetic, eclipsing millisecond pulsar J1816+4510 was recently discovered in a low-frequency radio survey with the Green Bank Telescope. With an orbital period of 8.7 hr and minimum companion mass of 0.16 Msun it appears to belong to an increasingly important class of pulsars that are ablating their low-mass companions. We report the discovery of the gamma-ray counterpart to this pulsar, and present a likely optical/ultraviolet counterpart as well. Using the radio ephemeris we detect pulsations in the unclassified gamma-ray source 2FGL J1816.5+4511, implying an efficiency of ~25% in converting the pulsar's spin-down luminosity into gamma-rays and adding PSR J1816+4510 to the large number of millisecond pulsars detected by Fermi. The likely optical/UV counterpart was identified through position coincidence (<0.1") and unusual colors. Assuming that it is the companion, with R=18.27+/-0.03 mag and effective temperature >15,000 K it would be among the brightest and hottest of low-mass pulsar companions, and appears qualitatively different from other eclipsing pulsar systems. In particular, current data suggest that it is a factor of two larger than most white dwarfs of its mass, but a factor of four smaller than its Roche lobe. We discuss possible reasons for its high temperature and odd size, and suggest that it recently underwent a violent episode of mass loss. Regardless of origin, its brightness and the relative unimportance of irradiation make it an ideal target for a mass, and hence a neutron star mass, determination.

preprint2012arXiv

Fermi-LAT constraints on the Pulsar Wind Nebula nature of HESS J1857+026

Since its launch, the Fermi satellite has firmly identified 5 pulsar wind nebulae plus a large number of candidates, all powered by young and energetic pulsars. HESS J1857+026 is a spatially extended gamma-ray source detected by H.E.S.S. and classified as a possible pulsar wind nebula candidate powered by PSR J1856+0245. We search for gamma-ray pulsations from PSR J1856+0245 and explore the characteristics of its associated pulsar wind nebula. Using a rotational ephemeris obtained from the Lovell telescope at Jodrell Bank Observatory at 1.5 GHz, we phase-fold 36 months of gamma-ray data acquired by the Large Area Telescope (LAT) aboard Fermi. We also perform a complete gamma-ray spectral and morphological analysis. No gamma-ray pulsations were detected from PSR J1856+0245. However, significant emission is detected at a position coincident with the TeV source HESS J1857+026. The gamma-ray spectrum is well described by a simple power-law with a spectral index of 1.53 \pm 0.11_{\rm stat} \pm 0.55_{\rm syst}$ and an energy flux of $G(0.1$--100 GeV$)=(2.71 \pm 0.52_{\rm stat} \pm 1.51_{\rm syst}) \times 10^{-11}$ ergs cm$^{-2}$ s$^{-1}$. The $γ$-ray luminosity is $L_{PWN}^γ (0.1$--100 GeV$)=(2.5 \pm 0.5_{stat} \pm 1.5_{syst}) \times 10^{35} (\frac{d}{9 kpc})^2$ ergs s$^{-1}$, assuming a distance of 9 kpc. This implies a $γ-$ray efficiency of $\sim$ 5% for $\dot{E}=4.6 \times 10^{36}$ erg $s^{-1}$, in the range expected for pulsar wind nebulae. Detailed multi-wavelength modeling provides new constraints on its pulsar wind nebula nature.

preprint2012arXiv

First LOFAR observations at very low frequencies of cluster-scale non-thermal emission: the case of Abell 2256

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

preprint2012arXiv

Four Highly Dispersed Millisecond Pulsars Discovered in the Arecibo PALFA Galactic Plane Survey

We present the discovery and phase-coherent timing of four highly dispersed millisecond pulsars (MSPs) from the Arecibo PALFA Galactic plane survey: PSRs J1844+0115, J1850+0124, J1900+0308, and J1944+2236. Three of the four pulsars are in binary systems with low-mass companions, which are most likely white dwarfs, and which have orbital periods on the order of days. The fourth pulsar is isolated. All four pulsars have large dispersion measures (DM > 100 pc cm-3), are distant (> 3.4 kpc), faint at 1.4 GHz (< 0.2 mJy), and are fully recycled (with spin periods P between 3.5 and 4.9 ms). The three binaries also have very small orbital eccentricities, as expected for tidally circularized, fully recycled systems with low-mass companions. These four pulsars have DM/P ratios that are among the highest values for field MSPs in the Galaxy. These discoveries bring the total number of confirmed MSPs from the PALFA survey to fifteen. The discovery of these MSPs illustrates the power of PALFA for finding weak, distant MSPs at low-Galactic latitudes. This is important for accurate estimates of the Galactic MSP population and for the number of MSPs that the Square Kilometer Array can be expected to detect.

preprint2012arXiv

Radio Searches of Fermi LAT Sources and Blind Search Pulsars: The Fermi Pulsar Search Consortium

We present a summary of the Fermi Pulsar Search Consortium (PSC), an international collaboration of radio astronomers and members of the Large Area Telescope (LAT) collaboration, whose goal is to organize radio follow-up observations of Fermi pulsars and pulsar candidates among the LAT gamma-ray source population. The PSC includes pulsar observers with expertise using the world's largest radio telescopes that together cover the full sky. We have performed very deep observations of all 35 pulsars discovered in blind frequency searches of the LAT data, resulting in the discovery of radio pulsations from four of them. We have also searched over 300 LAT gamma-ray sources that do not have strong associations with known gamma-ray emitting source classes and have pulsar-like spectra and variability characteristics. These searches have led to the discovery of a total of 43 new radio millisecond pulsars (MSPs) and four normal pulsars. These discoveries greatly increase the known population of MSPs in the Galactic disk, more than double the known population of so-called `black widow' pulsars, and contain many promising candidates for inclusion in pulsar timing arrays.

preprint2012arXiv

The PALFA Survey: Going to great depths to find radio pulsars

The on-going PALFA survey is searching the Galactic plane (|b| < 5 deg., 32 < l < 77 deg. and 168 < l < 214 deg.) for radio pulsars at 1.4 GHz using ALFA, the 7-beam receiver installed at the Arecibo Observatory. By the end of August 2012, the PALFA survey has discovered 100 pulsars, including 17 millisecond pulsars (P < 30 ms). Many of these discoveries are among the pulsars with the largest DM/P ratios, proving that the PALFA survey is capable of probing the Galactic plane for millisecond pulsars to a much greater depth than any previous survey. This is due to the survey's high sensitivity, relatively high observing frequency, and its high time and frequency resolution. Recently the rate of discoveries has increased, due to a new more sensitive spectrometer, two updated complementary search pipelines, the development of online collaborative tools, and access to new computing resources. Looking forward, focus has shifted to the application of artificial intelligence systems to identify pulsar-like candidates, and the development of an improved full-resolution pipeline incorporating more sophisticated radio interference rejection. The new pipeline will be used in a complete second analysis of data already taken, and will be applied to future survey observations. An overview of recent developments, and highlights of exciting discoveries will be presented.

preprint2012arXiv

Two Millisecond Pulsars Discovered by the PALFA Survey and a Shapiro Delay Measurement

We present two millisecond pulsar discoveries from the PALFA survey of the Galactic plane with the Arecibo telescope. PSR J1955+2527 is an isolated pulsar with a period of 4.87 ms, and PSR J1949+3106 has a period of 13.14 ms and is in a 1.9-day binary system with a massive companion. Their timing solutions, based on 4 years of timing measurements with the Arecibo, Green Bank, Nançay and Jodrell Bank telescopes, allow precise determination of spin and astrometric parameters, including precise determinations of their proper motions. For PSR J1949+3106, we can clearly detect the Shapiro delay. From this we measure the pulsar mass to be 1.47(+0.43/-0.31) solar masses, the companion mass to be 0.85(+0.14/-0.11) solar masses and the orbital inclination to be i = 79.9(+1.6/-1.9) degrees, where uncertainties correspond to +/- 1-σ confidence levels. With continued timing, we expect to also be able to detect the advance of periastron for the J1949+3106 system. This effect, combined with the Shapiro delay, will eventually provide very precise mass measurements for this system and a test of general relativity.

preprint2012arXiv

Wide-band Simultaneous Observations of Pulsars: Disentangling Dispersion Measure and Profile Variations

Dispersion in the interstellar medium is a well known phenomenon that follows a simple relationship, which has been used to predict the time delay of dispersed radio pulses since the late 1960s. We performed wide-band simultaneous observations of four pulsars with LOFAR (at 40-190 MHz), the 76-m Lovell Telescope (at 1400 MHz) and the Effelsberg 100-m Telescope (at 8000 MHz) to test the accuracy of the dispersion law over a broad frequency range. In this paper we present the results of these observations which show that the dispersion law is accurate to better than 1 part in 100000 across our observing band. We use this fact to constrain some of the properties of the ISM along the line-of-sight and use the lack of any aberration or retardation effects to determine upper limits on emission heights in the pulsar magnetosphere. We also discuss the effect of pulse profile evolution on our observations, and the implications that it could have for precision pulsar timing projects such as the detection of gravitational waves with pulsar timing arrays.

preprint2011arXiv

Arecibo PALFA Survey and Einstein@Home: Binary Pulsar Discovery by Volunteer Computing

We report the discovery of the 20.7 ms binary pulsar J1952+2630, made using the distributed computing project Einstein@Home in Pulsar ALFA survey observations with the Arecibo telescope. Follow-up observations with the Arecibo telescope confirm the binary nature of the system. We obtain a circular orbital solution with an orbital period of 9.4 hr, a projected orbital radius of 2.8 lt-s, and a mass function of f = 0.15 solar masses by analysis of spin period measurements. No evidence of orbital eccentricity is apparent; we set a 2-sigma upper limit e < 1.7e-3. The orbital parameters suggest a massive white dwarf companion with a minimum mass of 0.95 solar masses, assuming a pulsar mass of 1.4 solar masses. Most likely, this pulsar belongs to the rare class of intermediate mass binary pulsars. Future timing observations will aim to determine the parameters of this system further, measure relativistic effects, and elucidate the nature of the companion star.

preprint2011arXiv

Discovery of Nine Gamma-Ray Pulsars in Fermi-LAT Data Using a New Blind Search Method

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorporates several advances, including a metric-based gridding of the search parameter space (frequency, frequency derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs J1803-2149 and J2111+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6e35 erg/s and ages below 100 kyr). The seven remaining pulsars, PSRs J0106+4855, J0622+3749, J1620-4927, J1746-3239, J2028+3332, J2030+4415, J2139+4716, are older and less energetic; two of them are located at higher Galactic latitudes (|b| > 10 deg). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface magnetic field (2e11 G) of all LAT blind-search pulsars. PSR J2139+4716 has the lowest spin-down power (3e33 erg/s) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong to the increasing population of radio-quiet gamma-ray pulsars.

preprint2011arXiv

Observing pulsars and fast transients with LOFAR

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

preprint2011arXiv

Placing limits on the stochastic gravitational-wave background using European Pulsar Timing Array data

Direct detection of low-frequency gravitational waves ($10^{-9} - 10^{-8}$ Hz) is the main goal of pulsar timing array (PTA) projects. One of the main targets for the PTAs is to measure the stochastic background of gravitational waves (GWB) whose characteristic strain is expected to approximately follow a power-law of the form $h_c(f)=A (f/\hbox{yr}^{-1})^α$, where $f$ is the gravitational-wave frequency. In this paper we use the current data from the European PTA to determine an upper limit on the GWB amplitude $A$ as a function of the unknown spectral slope $α$ with a Bayesian algorithm, by modelling the GWB as a random Gaussian process. For the case $α=-2/3$, which is expected if the GWB is produced by supermassive black-hole binaries, we obtain a 95% confidence upper limit on $A$ of $6\times 10^{-15}$, which is 1.8 times lower than the 95% confidence GWB limit obtained by the Parkes PTA in 2006. Our approach to the data analysis incorporates the multi-telescope nature of the European PTA and thus can serve as a useful template for future intercontinental PTA collaborations.

preprint2010arXiv

Discovery of a new Soft Gamma Repeater, SGR J1833-0832

On 2010 March 19, the Swift/Burst Alert Telescope triggered on a short burst with temporal and spectral characteristics similar to those of Soft Gamma Repeater (SGR) bursts. The source location, however, did not coincide with any known SGR. Subsequent observations of the source error box with the Swift/X-ray Telescope and the Rossi X-ray Timing Explorer (RXTE) led to the discovery of a new X-ray source, with a spin period of 7.56 s, confirming SGR J1833-0832 as a new magnetar. Based on our detailed temporal and spectral analyses, we show that the new SGR is rapidly spinning down (4 x 10^{-12} s/s) and find an inferred dipole magnetic field of 1.8 x 10^{14} G. We also show that the X-ray flux of SGR J1833-0832 remained constant for approximately 20 days following the burst and then started to decline. We derived an accurate location of the source with the Chandra X-ray Observatory and we searched for a counterpart in deep optical and infrared observations of SGR J1833-0832, and for radio pulsed emission with the Westerbork Radio Synthesis Telescope. Finally, we compare the spectral and temporal properties of the source to the other magnetar candidates.

preprint2010arXiv

On the nature and evolution of the unique binary pulsar J1903+0327

(abridged) PSR J1903+0327, a millisecond pulsar in an eccentric (e = 0.44) 95-day orbit with a (~ 1Msun) companion poses a challenge to our understanding of stellar evolution in binary and multiple-star systems. Here we describe optical and radio observations which rule out most of the scenarios proposed to explain formation of this system. Radio timing measurements of three post-Keplerian effects yield the most precise measurement of the mass of a millisecond pulsar to date: 1.667 +/- 0.021 solar masses (99.7% confidence limit) (...). Optical spectroscopy of a proposed main sequence counterpart show that its orbital motion mirrors the pulsar's 95-day orbit; being therefore its binary companion (...) The optical detection also provides a measurement of the systemic radial velocity of the binary; this and the proper motion measured from pulsar timing allow the determination of the systemic 3-D velocity in the Galaxy. We find that the system is always within 270 pc of the plane of the Galaxy, but always more than 3 kpc away from the Galactic centre. Thus an exchange interaction in a dense stellar environment (like a globular cluster or the Galactic centre) is not likely to be the origin of this system. We suggest that after the supernova that formed it, the neutron star was in a tight orbit with a main-sequence star, the present companion was a tertiary farther out. The neutron star then accreted matter from its evolving inner MS companion, forming a millisecond pulsar. The former donor star then disappears, either due to a chaotic 3-body interaction with the outer star (caused by the expansion of the inner orbit that necessarily results from mass transfer), or in the case of a very compact inner system, due to ablation/accretion by the newly formed millisecond pulsar.

preprint2010arXiv

Pulsar Discovery by Global Volunteer Computing

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 192 countries to "mine" large data sets. It has now found a 40.8 Hz isolated pulsar in radio survey data from the Arecibo Observatory taken in February 2007. Additional timing observations indicate that this pulsar is likely a disrupted recycled pulsar. PSR J2007+2722's pulse profile is remarkably wide with emission over almost the entire spin period; the pulsar likely has closely aligned magnetic and spin axes. The massive computing power provided by volunteers should enable many more such discoveries.

preprint2010arXiv

Searches for gravitational waves from known pulsars with S5 LIGO data

We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358.

preprint2010arXiv

The LOFAR Known Pulsar Data Pipeline

Transient radio phenomena and pulsars are one of six LOFAR Key Science Projects (KSPs). As part of the Transients KSP, the Pulsar Working Group (PWG) has been developing the LOFAR Pulsar Data Pipelines to both study known pulsars as well as search for new ones. The pipelines are being developed for the Blue Gene/P (BG/P) supercomputer and a large Linux cluster in order to utilize enormous amounts of computational capabilities (50Tflops) to process data streams of up to 23TB/hour. The LOFAR pipeline output will be using the Hierarchical Data Format 5 (HDF5) to efficiently store large amounts of numerical data, and to manage complex data encompassing a variety of data types, across distributed storage and processing architectures. We present the LOFAR Known Pulsar Data Pipeline overview, the pulsar beam-formed data format, the status of the pipeline processing as well as our future plans for developing the LOFAR Pulsar Search Pipeline. These LOFAR pipelines and software tools are being developed as the next generation toolset for pulsar processing in Radio Astronomy.

preprint2009arXiv

Arecibo Pulsar Survey Using ALFA: Probing Radio Pulsar Intermittency and Transients

We present radio transient search algorithms, results, and statistics from the ongoing Arecibo Pulsar ALFA (PALFA) Survey of the Galactic plane. We have discovered seven objects by detecting isolated dispersed pulses and one of the new discoveries has a duty cycle of 0.01%, the smallest known. The impact of selection effects on the detectability and classification of intermittent sources is discussed, and the relative efficiencies of periodicity vs. single pulse searches are compared for various pulsar classes. We find that scintillation, off-axis detection and few rotation periods within an observation may misrepresent normal periodic pulsars as intermittent sources. Finally, we derive constraints on transient pulse rate and flux density from the PALFA survey parameters and results.

preprint2009arXiv

The Observed Spin Distributions of Millisecond Radio and X-ray Pulsars

We consider the currently observed spin distributions of various types of neutron stars, including isolated and binary radio millisecond pulsars in the Galactic plane and globular cluster system as well as neutron stars in low-mass X-ray binary systems where the spin rate is known either through coherent pulsations or burst oscillations. We find that the spin distributions of isolated and binary radio millisecond pulsars are statistically different, at least for those residing in globular clusters, with the binary pulsars being on average faster spinning. This result is likely to hold despite observational biases still affecting the observed spin distribution. A possible explanation for this is that the isolated radio millisecond pulsars are on average older than those in binary systems.

J. W. T. Hessels

What is connected

Connect this record

See the researcher in context

Building this map preview

105 published item(s)

The FRB 20121102A November rain in 2018 observed with the Arecibo Telescope

Arecibo observations of a burst storm from FRB 20121102A in 2016

Discoveries and Timing of Pulsars in NGC 6440

Fast radio bursts at the dawn of the 2020s

FRB 121102: drastic changes in the burst polarization contrasts with the stability of the persistent emission

PRECISE localizations of repeating Fast Radio Bursts

Searching for pulsars associated with polarised point sources using LOFAR: Initial discoveries from the TULIPP project

X-ray burst and persistent emission properties of the magnetar SGR 1830-0645 in outburst

A repeating fast radio burst source in a globular cluster

Highly polarised microstructure from the repeating FRB 20180916B

Milliarcsecond localisation of the repeating FRB 20201124A

Rotation Measure Evolution of the Repeating Fast Radio Burst Source FRB 121102

Study of 72 pulsars discovered in the PALFA survey: Timing analysis, glitch activity, emission variability, and a pulsar in an eccentric binary

A repeating fast radio burst source localised to a nearby spiral galaxy

Asymmetric mass ratios for bright double neutron-star mergers

Detection of Repeating FRB 180916.J0158+65 Down to Frequencies of 300 MHz

Periodic activity from a fast radio burst source

The LOFAR Tied-Array All-Sky Survey: Timing of 21 pulsars including the first binary pulsar discovered with LOFAR

A LOFAR census of non-recycled pulsars: extending below 80 MHz

Detection and timing of gamma-ray pulsations from the $707$ Hz pulsar J0952$-$0607

Quasi-simultaneous radio and X-ray observations of Aql X-1: probing low luminosities

Radio and X-ray monitoring of the accreting millisecond X-ray pulsar IGR J17591-2342 in outburst

Constraining very-high-energy and optical emission from FRB 121102 with the MAGIC telescopes

A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations

A Repeating Fast Radio Burst

Disk reflection and a possible disk wind during a soft X-ray state in the neutron star low-mass X-ray binary 1RXS J180408.9-342058

Enabling pulsar and fast transient searches using coherent dedispersion

From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release

High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array

Low-radio-frequency eclipses of the redback pulsar J2215+5135 observed in the image plane with LOFAR

New methods to constrain the radio transient rate: results from a survey of four fields with LOFAR

Search for transient gravitational waves in coincidence with short duration radio transients during 2007-2013

The International Pulsar Timing Array: First Data Release

The noise properties of 42 millisecond pulsars from the European Pulsar Timing Array and their impact on gravitational wave searches

The repeating Fast Radio Burst FRB 121102: Multi-wavelength observations and additional bursts

A LOFAR census of millisecond pulsars

Arecibo Pulsar Survey Using ALFA. IV. Mock Spectrometer Data Analysis, Survey Sensitivity, and the Discovery of 41 Pulsars

Discovery and Follow-up of Rotating Radio Transients with the Green Bank and LOFAR Telescopes

Einstein@Home Discovery of a PALFA Millisecond Pulsar in an Eccentric Binary Orbit

Imaging Jupiter's radiation belts down to 127 MHz with LOFAR

Limits on anisotropy in the nanohertz stochastic gravitational-wave background

Limits on Fast Radio Bursts at 145 MHz with ARTEMIS, a real-time software backend

LOFAR discovery of a quiet emission mode in PSR B0823+26

LOFAR MSSS: Detection of a low-frequency radio transient in 400 hrs of monitoring of the North Celestial Pole

LOFAR tied-array imaging and spectroscopy of solar S bursts

Probing Atmospheric Electric Fields in Thunderstorms through Radio Emission from Cosmic-Ray-Induced Air Showers

Properties and Evolution of the Redback Millisecond Pulsar Binary PSR J2129-0429

Pulsar polarisation below 200 MHz: Average profiles and propagation effects

Pulsar Wind Nebulae in the SKA era

Radio imaging observations of PSR J1023+0038 in an LMXB state

Radio Timing and Optical Photometry of the Black Widow System PSR J1953+1846A in the Globular Cluster M71

The LOFAR Multifrequency Snapshot Sky Survey (MSSS) I. Survey description and first results

Timing of five millisecond pulsars discovered in the PALFA survey

Wide-Band, Low-Frequency Pulse Profiles of 100 Radio Pulsars with LOFAR

A 24-Hour Global Campaign To Assess Precision Timing of the Millisecond Pulsar J1713+0747

A Cosmic Census of Radio Pulsars with the SKA

A millisecond pulsar in a stellar triple system

A state change in the low-mass X-ray binary XSS J12270-4859

Arecibo Pulsar Survey Using ALFA. III. Precursor Survey and Population Synthesis

Discovery of Carbon Radio Recombination Lines in absorption towards Cygnus~A

Fast Radio Burst Discovered in the Arecibo Pulsar ALFA Survey

Gravitational waves from known pulsars: results from the initial detector era

Initial LOFAR observations of Epoch of Reionization windows: II. Diffuse polarized emission in the ELAIS-N1 field

LOFAR tied-array imaging of Type III solar radio bursts

Pulsars in Globular Clusters with the SKA

Radio timing and optical photometry of the black widow system PSR J1518+0204C in the globular cluster M5

Swings between rotation and accretion power in a millisecond binary pulsar

The Green Bank Northern Celestial Cap Pulsar Survey - I: Survey Description, Data Analysis, and Initial Results

A New Accretion Disk Around the Missing Link Binary System PSR J1023+0038

A state change in the missing link binary pulsar system PSR J1023+0038

Calibrating High-Precision Faraday Rotation Measurements for LOFAR and the Next Generation of Low-Frequency Radio Telescopes

Detecting cosmic rays with the LOFAR radio telescope

Differential Frequency-dependent Delay from the Pulsar Magnetosphere

Discovery of the Optical Counterparts to Four Energetic Fermi Millisecond Pulsars