-
We present the longest-term timing study so far of three Rotating Radio
Transients (RRATs) - J1819-1458, J1840-1419 and J1913+1330 - performed using
the Lovell, Parkes and Green Bank telescopes over the past decade. We study
long-term and short- term variations of the pulse emission rate from these
RRATs and report a marginal indication of a long-term increase in pulse
detection rate over time for PSR J1819-1458 and J1913+1330. For PSR J1913+1330,
we also observe a two orders of magnitude variation in the observed pulse
detection rates across individual epochs, which may constrain the models
explaining the origin of RRAT pulses. PSR J1913+1330 is also observed to
exhibit a weak persistent emission mode. We investigate the post-glitch timing
properties of J1819-1458 (the only RRAT for which glitches are observed) and
discuss the implications for possible glitch models. Its post-glitch
over-recovery of the frequency derivative is magnetar-like and similar
behaviour is only observed for two other pulsars, both of which have relatively
high magnetic field strengths. Following the over-recovery we also observe that
some fraction of the pre-glitch frequency derivative is gradually recovered.
-
The MeerKAT telescope represents an outstanding opportunity for radio pulsar
timing science with its unique combination of a large collecting area and
aperture efficiency (effective area $\sim$7500 m$^2$), system temperature
($T<20$K), high slew speeds (1-2 deg/s), large bandwidths (770 MHz at 20cm
wavelengths), southern hemisphere location (latitude $\sim -30^\circ$) and
ability to form up to four sub-arrays. The MeerTime project is a five-year
program on the MeerKAT array by an international consortium that will regularly
time over 1000 radio pulsars to perform tests of relativistic gravity, search
for the gravitational wave signature induced by supermassive black hole
binaries in the timing residuals of millisecond pulsars, explore the interiors
of neutron stars through a pulsar glitch monitoring programme, explore the
origin and evolution of binary pulsars, monitor the swarms of pulsars that
inhabit globular clusters and monitor radio magnetars. In addition to these
primary programmes, over 1000 pulsars will have their arrival times monitored
and the data made immediately public. The MeerTime pulsar backend comprises two
server-class machines each of which possess four Graphics Processing Units. Up
to four pulsars can be coherently dedispersed simultaneously up to dispersion
measures of over 1000 pc cm$^{-3}$. All data will be provided in psrfits
format. The MeerTime backend will be capable of producing coherently
dedispersed filterbank data for timing multiple pulsars in the cores of
globular clusters that is useful for pulsar searches of tied array beams. All
MeerTime data will ultimately be made available for public use, and any
published results will include the arrival times and profiles used in the
results.
-
We have investigated the mode-changing properties of PSR B0329+54 using 31
epochs of simultaneous 13 cm/3 cm single-pulse observations obtained with
Shanghai Tian Ma 65 m telescope. The pulsar was found in the abnormal emission
mode 17 times, accounting for ~13% of the 41.6 hours total observation time.
Single pulse analyses indicate that mode changes took place simultaneously at
13 cm/3 cm within a few rotational periods. We detected occasional bright and
narrow pulses whose peak flux densities were 10 times higher than that of the
integrated profile in both bands. At 3 cm, about 0.66% and 0.27% of single
pulses were bright in the normal mode and abnormal mode respectively, but at 13
cm the occurrence rate was only about 0.007%. We divided the pulsar radiation
window into three components (C1, C2 and C3) corresponding to the main peaks of
the integrated profile. The bright pulses preferentially occurred at pulse
phases corresponding to the peaks of C2 and C3. Fluctuation spectra showed that
C2 had excess red noise in the normal mode, but broad quasi-periodic features
with central frequencies around 0.12 cycles/period in the abnormal mode. At 3
cm, C3 had a stronger quasi-periodic modulation centered around 0.06
cycles/period in the abnormal mode. Although there were some asymmetries in the
two-dimensional fluctuation spectra, we found no clear evidence for systematic
subpulse drifting. Consistent with previous low-frequency observations, we
found a very low nulling probability for B0329+54 with upper limits of 0.13%
and 1.68% at 13 cm/3 cm respectively.
-
We have observed a complex and continuous change in the integrated pulse
profile of PSR B2217+47, manifested as additional components trailing the main
peak. These transient components are detected over 6 years at $150$ MHz using
the LOw Frequency ARray (LOFAR), but they are not seen in contemporaneous
Lovell observations at $1.5$ GHz. We argue that propagation effects in the
ionized interstellar medium (IISM) are the most likely cause. The putative
structures in the IISM causing the profile variation are roughly half-way
between the pulsar and the Earth and have transverse radii $R \sim 30$ AU. We
consider different models for the structures. Under the assumption of spherical
symmetry, their implied average electron density is $\overline{n}_e \sim 100$
cm$^{-3}$. Since PSR B2217+47 is more than an order of magnitude brighter than
the average pulsar population visible to LOFAR, similar profile variations
would not have been identified in most pulsars, suggesting that subtle profile
variations in low-frequency profiles might be more common than we have observed
to date. Systematic studies of these variations at low frequencies can provide
a new tool to investigate the proprieties of the IISM and the limits to the
precision of pulsar timing.
-
We have observed and analysed the eclipses of the black widow pulsar
J1810+1744 at low radio frequencies. Using LOw-Frequency ARray (LOFAR) and
Westerbork Synthesis Radio Telescope observations between 2011--2015 we have
measured variations in flux density, dispersion measure and scattering around
eclipses. High-time-resolution, simultaneous beamformed and interferometric
imaging LOFAR observations show concurrent disappearance of pulsations and
total flux from the source during the eclipses, with a $3\sigma$ upper limit of
36 mJy ($<10\%$ of the pulsar's averaged out-of-eclipse flux density). The
dispersion measure variations are highly asymmetric, suggesting a tail of
material swept back due to orbital motion. The egress deviations are variable
on timescales shorter than the 3.6 hr orbital period and are indicative of a
clumpy medium. Additional pulse broadening detected during egress is typically
$<20\%$ of the pulsar's spin period, showing no evidence of scattering the
pulses beyond detectability in the beamformed data. The eclipses, lasting
$\sim13\%$ of the orbit at 149 MHz, are shown to be frequency-dependent with
total duration scaling as $\propto\nu^{-0.41\pm0.03}$. The results are
discussed in the context of the physical parameters of the system, and an
examination of eclipse mechanisms reveals cyclotron-synchrotron absorption as
the most likely primary cause, although non-linear scattering mechanisms cannot
be quantitatively ruled out. The inferred mass loss rate is a similar
order-of-magnitude to the mean rate required to fully evaporate the companion
in a Hubble time.
-
An understanding of pulsar timing noise offers the potential to improve the
timing precision of a large number of pulsars as well as facilitating our
understanding of pulsar magnetospheres. For some sources, timing noise is
attributable to a pulsar switching between two different spin-down rates
$(\dot{\nu})$. Such transitions may be common but difficult to resolve using
current techniques. In this work, we use simulations of $\dot{\nu}$-variable
pulsars to investigate the likelihood of resolving individual $\dot{\nu}$
transitions. We inject step-changes in the value of $\dot{\nu}$ with a wide
range of amplitudes and switching timescales. We then attempt to redetect these
transitions using standard pulsar timing techniques. The pulse arrival-time
precision and the observing cadence are varied. Limits on $\dot{\nu}$
detectability based on the effects such transitions have on the timing
residuals are derived. With the typical cadences and timing precision of
current timing programs, we find we are insensitive to a large region of
$\Delta \dot{\nu}$ parameter space which encompasses small, short timescale
switches. We find, where the rotation and emission states are correlated, that
using changes to the pulse shape to estimate $\dot{\nu}$ transition epochs, can
improve detectability in certain scenarios. The effects of cadence on $\Delta
\dot{\nu}$ detectability are discussed and we make comparisons with a known
population of intermittent and mode-switching pulsars. We conclude that for
short timescale, small switches, cadence should not be compromised when new
generations of ultra-sensitive radio telescopes are online.
-
We present results from the first directed search for nontensorial
gravitational waves. While general relativity allows for tensorial (plus and
cross) modes only, a generic metric theory may, in principle, predict waves
with up to six different polarizations. This analysis is sensitive to
continuous signals of scalar, vector or tensor polarizations, and does not rely
on any specific theory of gravity. After searching data from the first
observation run of the advanced LIGO detectors for signals at twice the
rotational frequency of 200 known pulsars, we find no evidence of gravitational
waves of any polarization. We report the first upper limits for scalar and
vector strains, finding values comparable in magnitude to previously-published
limits for tensor strain. Our results may be translated into constraints on
specific alternative theories of gravity.
-
We present the result of searches for gravitational waves from 200 pulsars
using data from the first observing run of the Advanced LIGO detectors. We find
no significant evidence for a gravitational-wave signal from any of these
pulsars, but we are able to set the most constraining upper limits yet on their
gravitational-wave amplitudes and ellipticities. For eight of these pulsars,
our upper limits give bounds that are improvements over the indirect spin-down
limit values. For another 32, we are within a factor of 10 of the spin-down
limit, and it is likely that some of these will be reachable in future runs of
the advanced detector. Taken as a whole, these new results improve on previous
limits by more than a factor of two.
-
Single pulses preserve information about the pulsar radio emission and
propagation in the pulsar magnetosphere, and understanding the behaviour of
their variability is essential for estimating the fundamental limit on the
achievable pulsar timing precision. Here we report the findings of our analysis
of single pulses from PSR J1713+0747 with data collected by the Large European
Array for Pulsars (LEAP). We present statistical studies of the pulse
properties that include distributions of their energy, phase and width. Two
modes of systematic sub-pulse drifting have been detected, with a periodicity
of 7 and 3 pulse periods. The two modes appear at different ranges of pulse
longitude but overlap under the main peak of the integrated profile. No
evidence for pulse micro-structure is seen with a time resolution down to 140
ns. In addition, we show that the fractional polarisation of single pulses
increases with their pulse peak flux density. By mapping the probability
density of linear polarisation position angle with pulse longitude, we reveal
the existence of two orthogonal polarisation modes. Finally, we find that the
resulting phase jitter of integrated profiles caused by single pulse
variability can be described by a Gaussian probability distribution only when
at least 100 pulses are used for integration. Pulses of different flux
densities and widths contribute approximately equally to the phase jitter, and
no improvement on timing precision is achieved by using a sub-set of pulses
with a specific range of flux density or width.
-
We present first results from a LOFAR census of non-recycled pulsars. The
census includes almost all such pulsars known (194 sources) at declinations
${\rm Dec}> 8^\circ$ and Galactic latitudes $|{\rm Gb}| > 3^\circ$, regardless
of their expected flux densities and scattering times. Each pulsar was observed
for $\geq 20$ minutes in the contiguous frequency range of 110--188 MHz.
Full-Stokes data were recorded. We present the dispersion measures, flux
densities, and calibrated total intensity profiles for the 158 pulsars detected
in the sample. The median uncertainty in census dispersion measures ($1.5
\times 10^{-3}$ pc cm$^{-3}$) is ten times smaller, on average, than in the
ATNF pulsar catalogue. We combined census flux densities with those in the
literature and fitted the resulting broadband spectra with single or broken
power-law functions. For 48 census pulsars such fits are being published for
the first time. Typically, the choice between single and broken power-laws, as
well as the location of the spectral break, were highly influenced by the
spectral coverage of the available flux density measurements. In particular,
the inclusion of measurements below 100 MHz appears essential for investigating
the low-frequency turnover in the spectra for most of the census pulsars. For
several pulsars, we compared the spectral indices from different works and
found the typical spread of values to be within 0.5--1.5, suggesting a
prevailing underestimation of spectral index errors in the literature. The
census observations yielded some unexpected individual source results, as we
describe in the paper. Lastly, we will provide this unique sample of wide-band,
low-frequency pulse profiles via the European Pulsar Network Database.
-
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/).
-
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.
-
We study the spin-down changes of PSR B1859$+$07 over a period of more than
28 years of radio observation. We identify that the time derivative of the
rotational frequency ($\nu$) varies quasi-periodically with a period of
$\sim$350 days, switching mainly between two spin-down states. The profile
shape of the pulsar is correlated with the $\dot \nu$ variation, producing two
slightly different profile shapes corresponding to high- and low-$\dot \nu$
states. In addition to these two normal emission states, we confirm the
occasional flare-state of the pulsar, in which the emission appears early in
spin phase compared to that of the common normal emission. The profile of the
flare-state is significantly different from that of the two normal emission
states. The correlation analysis further shows that the flare-state is not
directly linked with the $\dot \nu$ changes. With a simple emission beam model,
we estimate the emission altitude of the normal emission to be 240~km, and
explain the origin of the flare-state as an emission height variation from the
leading edge of the beam. We also argue that the emission of these states can
be explained with a partially active beam model. In this scenario, the trailing
portion of the radio beam is usually active and the normal emission is
produced. The flare-state occurs when the leading edge of the beam becomes
active while the trailing part is being blocked. This model estimates a fixed
emission altitude of 360~km. However, the cause of the flare-state (i.e. the
emission height variation, or the time-dependent activity across the radio
beam) is not easily explained.
-
In this paper, we report radio observations of the Galactic Center magnetar
PSR J1745-2900 at six epochs between June and October, 2014. These observations
were carried out using the new Shanghai Tian Ma Radio Telescope at a frequency
of 8.6 GHz. Both the flux density and integrated profile of PSR J1745-2900 show
dramatic changes from epoch to epoch showing that the pulsar was in its
"erratic" phase. On MJD 56836, the flux density of this magnetar was about 8.7
mJy, which was ten times large than that reported at the time of discovery,
enabling a single-pulse analysis. The emission is dominated by narrow "spiky"
pulses which follow a log-normal distribution in peak flux density. From 1913
pulses, we detected 53 pulses whose peak flux density is ten times greater than
that of the integrated profile. They are concentrated in pulse phase at the
peaks of the integrated profile. The pulse widths at the 50% level of these
bright pulses was between 0.2 to 0.9 deg, much narrower than that of integrated
profile (~12 deg). The observed pulse widths may be limited by interstellar
scattering. No clear correlation was found between the widths and peak flux
density of these pulses and no evidence was found for subpulse drifting.
Relatively strong spiky pulses are also detected in the other five epochs of
observation, showing the same properties as that detected in MJD 56836. These
strong spiky pulses cannot be classified as "giant" pulses but are more closely
related to normal pulse emission.
-
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\c{c}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 \nu^{\alpha}$), with
$\alpha\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.
-
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.
-
We present an analysis of approximately 200 hours of observations of the
pulsars J1634$-$5107, J1717$-$4054 and J1853$+$0505, taken over the course of
14.7 yr. We show that all of these objects exhibit long term nulls and
radio-emitting phases (i.e. minutes to many hours), as well as considerable
nulling fractions (NFs) in the range $\sim67\,\% - 90\,\%$. PSR J1717$-$4054 is
also found to exhibit short timescale nulls ($1 - 40~P$) and burst phases
($\lesssim 200~P$) during its radio-emitting phases. This behaviour acts to
modulate the NF, and therefore the detection rate of the source, over
timescales of minutes. Furthermore, PSR J1853$+$0505 is shown to exhibit a weak
emission state, in addition to its strong and null states, after sufficient
pulse integration. This further indicates that nulls may often only represent
transitions to weaker emission states which are below the sensitivity
thresholds of particular observing systems. In addition, we detected a
peak-to-peak variation of $33\pm1\,\%$ in the spin-down rate of PSR
J1717$-$4054, over timescales of hundreds of days. However, no long-term
correlation with emission variation was found.
-
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.
-
The SKA will discover tens of thousands of pulsars and provide unprecedented
data quality on these, as well as the currently known population, due to its
unrivalled sensitivity. Here, we outline the state of the art of our
understanding of magnetospheric radio emission from pulsars and how we will use
the SKA to solve the open problems in pulsar magnetospheric physics.
-
Providing a link between magnetars and radio pulsars, high magnetic field
neutron stars are ideal targets to investigate how bursting/magnetospheric
activity and braking torque variations are connected to rotational glitches.
The last spin-up glitch of the highly magnetised pulsar J1119$-$6127 back in
2007 was the first glitch in a rotationally powered radio pulsar to be
accompanied by radiative changes. Moreover, it was followed by an uncommon
glitch relaxation that resulted in a smaller spin-down rate relative to the
prediction of the pre-glitch timing model. Here, we present 4 years of new
radio timing observations and analyse the total of 16 years of timing data for
this source. The new data uncover an ongoing evolution of the spin-down rate,
thereby allowing us to exclude permanent changes in the external or internal
torque as a stand-alone cause of the peculiar features of the glitch recovery.
Furthermore, no additional variations of the radio pulse profile are detected,
strengthening the association of the previously observed transient emission
features with the glitching activity. A self-consistent measurement of the
braking index yields a value $n\simeq2.7$, indicating a trajectory in the
$P-\dot{P}$ plane inclined towards the magnetars. Such a potential evolutionary
link might be strengthened by a, possibly permanent, reduction of $\sim15\%$ in
$n$ at the epoch of the 2007 glitch.
-
We study the spin-down properties of PSR B0919+06 based on almost 30 years of
radio observations. We confirm that the time derivative of the rotational
frequency $\dot \nu$ is modulated quasi-periodically and show that it exhibits
a repeating double-peaked structure throughout the entire observation span. We
model the $\dot \nu$ variation of the pulsar assuming two spin-down rates with
sudden switches between them in time. Our results show that the double-peak
structure in $\dot \nu$ has a repetition time of about 630 days until MJD 52000
(April 2001) and 550 days since then. During this cycle, the pulsar spin varies
from the lower spin-down rate to the upper spin-down rate twice with different
amounts of time spent in each state, resulting in a further quasi-stable
secondary modulation of the two-state switching. This particular spin-down
state switching is broadly consistent with free precession of the pulsar,
however, a strong evidence linked with this mechanism is not clearly
established. We also confirm that the pulsar occasionally emits groups of
pulses which appear early in pulse phase, so-called "flares", and these events
significantly contribute to the pulse profile shape. We find the $\dot \nu$
modulation and the pulse shape variations are correlated throughout the
observations. However, the flare-state is not entirely responsible for this
correlation. In addition to the flare-state, we detect flare-like events from
the pulsar in single pulse observations. During these events, the shift in
pulse phase is small compared to that of the main flare-state and clearly
visible only in single pulse observations.
-
We present radio intensity and polarisation profiles of 28
$\gamma$-ray-detected pulsars with the aim of putting constraints on their
viewing geometries using data from the Parkes telescope. Constraints are formed
both from the goodness-of-fit of the position angles to the Rotating Vector
Model and from the beam opening angle considering aberration and retardation
effects. Uncertainties on the relevant parameters are systematically taken into
account in order to produce a more robust constraint, using a new approach.
Surprisingly, we find that the distribution of the magnetic inclination angle
($\alpha$) in this subset of pulsars peaks at low values, contrary to the
predictions of $\gamma$-ray models. We find a lack of correlation between these
and a set of $\alpha$ values which were derived using $\gamma$-ray light
curves, suggesting a problem in the interpretation of the data in one or both
of these domains. Finally, we also show that the $\alpha$ distribution of
pulsars with multiple radio components is no different to that of
single-component pulsars.
-
Several studies have shown the distribution of pulsars' magnetic inclination
angles to be skewed towards low values compared with the distribution expected
if the rotation and magnetic axes are placed randomly on the star. Here we
focus on a sample of 28 $\gamma$-ray-detected pulsars using data taken as part
of the Parkes telescope's \emph{FERMI} timing program. In doing so we find a
preference in the sample for low magnetic inclination angles, $\alpha$, in
stark contrast to both the expectation that the magnetic and rotation axes are
orientated randomly at the birth of the pulsar and to
$\gamma$-ray-emission-model-based expected biases. In this paper, after
exploring potential explanations, we conclude that there are two possible
causes of this preference, namely that low $\alpha$ values are intrinsic to the
sample, or that the emission regions extend outside what is traditionally
thought to be the open-field-line region in a way which is dependent on the
magnetic inclination. Each possibility is expected to have important
consequences, ranging from supernova physics to population studies of pulsars
and considerations of the radio beaming fraction. We also present a simple
conversion scheme between the observed and intrinsic magnetic inclinations
which is valid under the assumption that the observed skew is not intrinsic and
which can be applied to all existing measurements. We argue that extending the
active field-line region will help to resolve the existing tension between
emission geometries derived from radio polarisation measurements and those
required to model $\gamma$-ray light curves.
-
We present broadband, low-frequency (25-80 MHz and 110-190 MHz) LOFAR
observations of PSR B0943+10, with the goal of better illuminating the nature
of its enigmatic mode-switching behaviour. This pulsar shows two relatively
stable states: a Bright (B) and Quiet (Q) mode, each with different
characteristic brightness, profile morphology, and single-pulse properties. We
model the average profile evolution both in frequency and time from the onset
of each mode, and highlight the differences between the two modes. In both
modes, the profile evolution can be well explained by radius-to-frequency
mapping at altitudes within a few hundred kilometres of the stellar surface. If
both B and Q-mode emission originate at the same magnetic latitude, then we
find that the change of emission height between the modes is less than 6%. We
also find that, during B-mode, the average profile is gradually shifting
towards later spin phase and then resets its position at the next Q-to-B
transition. The observed B-mode profile delay is frequency-independent (at
least from 25-80 MHz) and asymptotically changes towards a stable value of
about 0.004 in spin phase by the end of mode instance, much too large to be due
to changing spin-down rate. Such a delay can be interpreted as a gradual
movement of the emission cone against the pulsar's direction of rotation, with
different field lines being illuminated over time. Another interesting
explanation is a possible variation of accelerating potential inside the polar
gap. This explanation connects the observed profile delay to the gradually
evolving subpulse drift rate, which depends on the gradient of the potential
across the field lines.
-
We have conducted two pilot surveys for radio pulsars and fast transients
with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the
first low-frequency fast-radio burst limit and the discovery of two new
pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a
large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times.
Each observation covered ~75 sq. deg using 7 independent fields formed by
incoherently summing the high-band antenna fields. The second pilot survey, the
LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold
increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR
"Superterp" stations, we formed 19 tied-array beams, together covering 4 sq.
deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of
dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy
for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected
65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR
or any digital aperture array. LOTAS also re-detected 27 previously known
pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky
surveys for pulsars and fast transients, and they set the stage for further
surveying efforts using LOFAR and the planned low-frequency component of the
Square Kilometer Array.
An understanding of pulsar timing noise offers the potential to improve the timing precision of a large number of pulsars as well as facilitating our understanding of pulsar magnetospheres. For some sources, timing noise is attributable to a pulsar switching between two different spin-down rates $(\dot{\nu})$. Such transitions may be common but difficult to resolve using current techniques. In this work, we use simulations of $\dot{\nu}$-variable pulsars to investigate the likelihood of resolving individual $\dot{\nu}$ transitions. We inject step-changes in the value of $\dot{\nu}$ with a wide range of amplitudes and switching timescales. We then attempt to redetect these transitions using standard pulsar timing techniques. The pulse arrival-time precision and the observing cadence are varied. Limits on $\dot{\nu}$ detectability based on the effects such transitions have on the timing residuals are derived. With the typical cadences and timing precision of current timing programs, we find we are insensitive to a large region of $\Delta \dot{\nu}$ parameter space which encompasses small, short timescale switches. We find, where the rotation and emission states are correlated, that using changes to the pulse shape to estimate $\dot{\nu}$ transition epochs, can improve detectability in certain scenarios. The effects of cadence on $\Delta \dot{\nu}$ detectability are discussed and we make comparisons with a known population of intermittent and mode-switching pulsars. We conclude that for short timescale, small switches, cadence should not be compromised when new generations of ultra-sensitive radio telescopes are online.
We present the result of searches for gravitational waves from 200 pulsars using data from the first observing run of the Advanced LIGO detectors. We find no significant evidence for a gravitational-wave signal from any of these pulsars, but we are able to set the most constraining upper limits yet on their gravitational-wave amplitudes and ellipticities. For eight of these pulsars, our upper limits give bounds that are improvements over the indirect spin-down limit values. For another 32, we are within a factor of 10 of the spin-down limit, and it is likely that some of these will be reachable in future runs of the advanced detector. Taken as a whole, these new results improve on previous limits by more than a factor of two.
Single pulses preserve information about the pulsar radio emission and propagation in the pulsar magnetosphere, and understanding the behaviour of their variability is essential for estimating the fundamental limit on the achievable pulsar timing precision. Here we report the findings of our analysis of single pulses from PSR J1713+0747 with data collected by the Large European Array for Pulsars (LEAP). We present statistical studies of the pulse properties that include distributions of their energy, phase and width. Two modes of systematic sub-pulse drifting have been detected, with a periodicity of 7 and 3 pulse periods. The two modes appear at different ranges of pulse longitude but overlap under the main peak of the integrated profile. No evidence for pulse micro-structure is seen with a time resolution down to 140 ns. In addition, we show that the fractional polarisation of single pulses increases with their pulse peak flux density. By mapping the probability density of linear polarisation position angle with pulse longitude, we reveal the existence of two orthogonal polarisation modes. Finally, we find that the resulting phase jitter of integrated profiles caused by single pulse variability can be described by a Gaussian probability distribution only when at least 100 pulses are used for integration. Pulses of different flux densities and widths contribute approximately equally to the phase jitter, and no improvement on timing precision is achieved by using a sub-set of pulses with a specific range of flux density or width.
We present first results from a LOFAR census of non-recycled pulsars. The census includes almost all such pulsars known (194 sources) at declinations ${\rm Dec}> 8^\circ$ and Galactic latitudes $|{\rm Gb}| > 3^\circ$, regardless of their expected flux densities and scattering times. Each pulsar was observed for $\geq 20$ minutes in the contiguous frequency range of 110--188 MHz. Full-Stokes data were recorded. We present the dispersion measures, flux densities, and calibrated total intensity profiles for the 158 pulsars detected in the sample. The median uncertainty in census dispersion measures ($1.5 \times 10^{-3}$ pc cm$^{-3}$) is ten times smaller, on average, than in the ATNF pulsar catalogue. We combined census flux densities with those in the literature and fitted the resulting broadband spectra with single or broken power-law functions. For 48 census pulsars such fits are being published for the first time. Typically, the choice between single and broken power-laws, as well as the location of the spectral break, were highly influenced by the spectral coverage of the available flux density measurements. In particular, the inclusion of measurements below 100 MHz appears essential for investigating the low-frequency turnover in the spectra for most of the census pulsars. For several pulsars, we compared the spectral indices from different works and found the typical spread of values to be within 0.5--1.5, suggesting a prevailing underestimation of spectral index errors in the literature. The census observations yielded some unexpected individual source results, as we describe in the paper. Lastly, we will provide this unique sample of wide-band, low-frequency pulse profiles via the European Pulsar Network Database.
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/).
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.
We study the spin-down changes of PSR B1859$+$07 over a period of more than 28 years of radio observation. We identify that the time derivative of the rotational frequency ($\nu$) varies quasi-periodically with a period of $\sim$350 days, switching mainly between two spin-down states. The profile shape of the pulsar is correlated with the $\dot \nu$ variation, producing two slightly different profile shapes corresponding to high- and low-$\dot \nu$ states. In addition to these two normal emission states, we confirm the occasional flare-state of the pulsar, in which the emission appears early in spin phase compared to that of the common normal emission. The profile of the flare-state is significantly different from that of the two normal emission states. The correlation analysis further shows that the flare-state is not directly linked with the $\dot \nu$ changes. With a simple emission beam model, we estimate the emission altitude of the normal emission to be 240~km, and explain the origin of the flare-state as an emission height variation from the leading edge of the beam. We also argue that the emission of these states can be explained with a partially active beam model. In this scenario, the trailing portion of the radio beam is usually active and the normal emission is produced. The flare-state occurs when the leading edge of the beam becomes active while the trailing part is being blocked. This model estimates a fixed emission altitude of 360~km. However, the cause of the flare-state (i.e. the emission height variation, or the time-dependent activity across the radio beam) is not easily explained.
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\c{c}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 \nu^{\alpha}$), with $\alpha\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.
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.
We present an analysis of approximately 200 hours of observations of the pulsars J1634$-$5107, J1717$-$4054 and J1853$+$0505, taken over the course of 14.7 yr. We show that all of these objects exhibit long term nulls and radio-emitting phases (i.e. minutes to many hours), as well as considerable nulling fractions (NFs) in the range $\sim67\,\% - 90\,\%$. PSR J1717$-$4054 is also found to exhibit short timescale nulls ($1 - 40~P$) and burst phases ($\lesssim 200~P$) during its radio-emitting phases. This behaviour acts to modulate the NF, and therefore the detection rate of the source, over timescales of minutes. Furthermore, PSR J1853$+$0505 is shown to exhibit a weak emission state, in addition to its strong and null states, after sufficient pulse integration. This further indicates that nulls may often only represent transitions to weaker emission states which are below the sensitivity thresholds of particular observing systems. In addition, we detected a peak-to-peak variation of $33\pm1\,\%$ in the spin-down rate of PSR J1717$-$4054, over timescales of hundreds of days. However, no long-term correlation with emission variation was found.
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.
The SKA will discover tens of thousands of pulsars and provide unprecedented data quality on these, as well as the currently known population, due to its unrivalled sensitivity. Here, we outline the state of the art of our understanding of magnetospheric radio emission from pulsars and how we will use the SKA to solve the open problems in pulsar magnetospheric physics.
Providing a link between magnetars and radio pulsars, high magnetic field neutron stars are ideal targets to investigate how bursting/magnetospheric activity and braking torque variations are connected to rotational glitches. The last spin-up glitch of the highly magnetised pulsar J1119$-$6127 back in 2007 was the first glitch in a rotationally powered radio pulsar to be accompanied by radiative changes. Moreover, it was followed by an uncommon glitch relaxation that resulted in a smaller spin-down rate relative to the prediction of the pre-glitch timing model. Here, we present 4 years of new radio timing observations and analyse the total of 16 years of timing data for this source. The new data uncover an ongoing evolution of the spin-down rate, thereby allowing us to exclude permanent changes in the external or internal torque as a stand-alone cause of the peculiar features of the glitch recovery. Furthermore, no additional variations of the radio pulse profile are detected, strengthening the association of the previously observed transient emission features with the glitching activity. A self-consistent measurement of the braking index yields a value $n\simeq2.7$, indicating a trajectory in the $P-\dot{P}$ plane inclined towards the magnetars. Such a potential evolutionary link might be strengthened by a, possibly permanent, reduction of $\sim15\%$ in $n$ at the epoch of the 2007 glitch.
We study the spin-down properties of PSR B0919+06 based on almost 30 years of radio observations. We confirm that the time derivative of the rotational frequency $\dot \nu$ is modulated quasi-periodically and show that it exhibits a repeating double-peaked structure throughout the entire observation span. We model the $\dot \nu$ variation of the pulsar assuming two spin-down rates with sudden switches between them in time. Our results show that the double-peak structure in $\dot \nu$ has a repetition time of about 630 days until MJD 52000 (April 2001) and 550 days since then. During this cycle, the pulsar spin varies from the lower spin-down rate to the upper spin-down rate twice with different amounts of time spent in each state, resulting in a further quasi-stable secondary modulation of the two-state switching. This particular spin-down state switching is broadly consistent with free precession of the pulsar, however, a strong evidence linked with this mechanism is not clearly established. We also confirm that the pulsar occasionally emits groups of pulses which appear early in pulse phase, so-called "flares", and these events significantly contribute to the pulse profile shape. We find the $\dot \nu$ modulation and the pulse shape variations are correlated throughout the observations. However, the flare-state is not entirely responsible for this correlation. In addition to the flare-state, we detect flare-like events from the pulsar in single pulse observations. During these events, the shift in pulse phase is small compared to that of the main flare-state and clearly visible only in single pulse observations.
We present broadband, low-frequency (25-80 MHz and 110-190 MHz) LOFAR observations of PSR B0943+10, with the goal of better illuminating the nature of its enigmatic mode-switching behaviour. This pulsar shows two relatively stable states: a Bright (B) and Quiet (Q) mode, each with different characteristic brightness, profile morphology, and single-pulse properties. We model the average profile evolution both in frequency and time from the onset of each mode, and highlight the differences between the two modes. In both modes, the profile evolution can be well explained by radius-to-frequency mapping at altitudes within a few hundred kilometres of the stellar surface. If both B and Q-mode emission originate at the same magnetic latitude, then we find that the change of emission height between the modes is less than 6%. We also find that, during B-mode, the average profile is gradually shifting towards later spin phase and then resets its position at the next Q-to-B transition. The observed B-mode profile delay is frequency-independent (at least from 25-80 MHz) and asymptotically changes towards a stable value of about 0.004 in spin phase by the end of mode instance, much too large to be due to changing spin-down rate. Such a delay can be interpreted as a gradual movement of the emission cone against the pulsar's direction of rotation, with different field lines being illuminated over time. Another interesting explanation is a possible variation of accelerating potential inside the polar gap. This explanation connects the observed profile delay to the gradually evolving subpulse drift rate, which depends on the gradient of the potential across the field lines.
We have conducted two pilot surveys for radio pulsars and fast transients with the Low-Frequency Array (LOFAR) around 140 MHz and here report on the first low-frequency fast-radio burst limit and the discovery of two new pulsars. The first survey, the LOFAR Pilot Pulsar Survey (LPPS), observed a large fraction of the northern sky, ~1.4 x 10^4 sq. deg, with 1-hr dwell times. Each observation covered ~75 sq. deg using 7 independent fields formed by incoherently summing the high-band antenna fields. The second pilot survey, the LOFAR Tied-Array Survey (LOTAS), spanned ~600 sq. deg, with roughly a 5-fold increase in sensitivity compared with LPPS. Using a coherent sum of the 6 LOFAR "Superterp" stations, we formed 19 tied-array beams, together covering 4 sq. deg per pointing. From LPPS we derive a limit on the occurrence, at 142 MHz, of dispersed radio bursts of < 150 /day/sky, for bursts brighter than S > 107 Jy for the narrowest searched burst duration of 0.66 ms. In LPPS, we re-detected 65 previously known pulsars. LOTAS discovered two pulsars, the first with LOFAR or any digital aperture array. LOTAS also re-detected 27 previously known pulsars. These pilot studies show that LOFAR can efficiently carry out all-sky surveys for pulsars and fast transients, and they set the stage for further surveying efforts using LOFAR and the planned low-frequency component of the Square Kilometer Array.