Source author record

M. Kramer

M. Kramer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

259works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

259 published item(s)

preprint2026arXiv

Radio timing constraints on the orbital orientation and component masses of PSR J1455$-$3330

PSR J1455$-$3330 is a $\sim$7.98 ms pulsar in a $\sim$76.17 day nearly circular orbit with a white dwarf companion. In this work, we combine the available Lovell, Nançay decimetric Radio Telescope, Green Bank, and MeerKAT pulsar timing data spanning $\sim$ 30 years to measure the kinematic and relativistic effects of PSR J1455$-$3330 to constrain its 3D orbital geometry and component masses. We detect a relativistic Shapiro delay signal. We measure a significant orthometric amplitude $h_3 = 0.307^{+0.022}_{-0.026}$ $μ$s and an orthometric ratio $ς= 0.551^{+0.057}_{-0.054}$. We measure the change in projected semi-major axis $\dot{x} = -202.1^{+2.5}_{-2.7} \times10^{-16} \, \rm s\,s^{-1}$ with high significance, parallax, $\varpi$ = 1.11(6) mas, parallax derived distance 0.90(5) kpc, and a precise total proper motion magnitude of 12.432(2) mas yr$^{-1}$. A self-consistent analysis of all kinematic and relativistic effects, assuming general relativity, yields two solutions: (1) a pulsar mass $M_{\rm p} = 1.39^{+0.38}_{-0.18}\, \rm M_{\odot}$, a companion mass $M_{\rm c} = 0.293^{+0.056}_{-0.026}$ $\rm M_{\odot}$, an orbital inclination, $i = 63(2)^{\circ}$, and longitude of the ascending node, $Ω= 212(12)^{\circ}$ or (2) a pulsar mass $M_{\rm p} = 1.53^{+1.10}_{-0.22} \, \rm M_{\odot}$, a companion mass $M_{\rm c} = 0.309^{+0.163}_{-0.026}\, \rm M_{\odot}$, an orbital inclination, $i = 123(4)^{\circ}$, and longitude of the ascending node, $Ω= 334(12)^{\circ}$. All uncertainties represent the 68.27$\%$ credibility region. These results strongly favour a helium-dominated white dwarf companion.

preprint2026arXiv

Reconstruction of atmospheric neutrinos in DUNE's horizontal-drift far-detector module

This paper reports on the capabilities in reconstructing and identifying atmospheric neutrino interactions in one of the Deep Underground Neutrino Experiment's (DUNE) far detector modules, a liquid argon time projection chamber (LArTPC) with horizontal drift (FD-HD) of ionization electrons. The reconstruction is based upon the workflow developed for DUNE's long-baseline oscillation analysis, with some necessary machine-learning models' retraining and the addition of features relevant only to atmospheric neutrinos such as the neutrino direction reconstruction. Where relevant, the impact of the detection of the charged particles of the hadronic system is emphasized, and comparisons are carried out between the case when lepton-only information is considered in the reconstruction (as is the case for many neutrino oscillation experiments), versus when all particles identified in the LArTPC were included. Three neutrino direction reconstruction methods have been developed and studied for the atmospheric analyses: using lepton-only information, using all reconstructed particles, and using only correlations from reconstructed hits. The results indicate that incorporating more than just lepton information significantly improves the resolution of both neutrino direction and energy reconstruction. The angle reconstruction algorithms developed in this work result in no strong dependence on particle direction for reconstruction efficiencies or neutrino flavor identification. This comprehensive review of the reconstruction of atmospheric neutrinos in DUNE's FD-HD LArTPC is the first step towards developing a first neutrino oscillation sensitivity analysis, which will ready DUNE for its first measurements.

preprint2025arXiv

Populations of evolved massive binary stars in the Small Magellanic Cloud I: Predictions from detailed evolution models

Context. The majority of massive stars are born with a close binary companion. How this affects their evolution and fate is still largely uncertain, especially at low metallicity. Aims. We derive synthetic populations of massive post-interaction binary products and compare them with corresponding observed populations in the Small Magellanic Cloud (SMC). Methods. We analyse 53298 detailed binary evolutionary models computed with MESA. Our models include the physics of rotation, mass and angular momentum transfer, magnetic internal angular momentum transport, and tidal spin-orbit coupling. They cover initial primary masses of 5-100Msun, initial mass ratios of 0.3-0.95, and all initial periods for which interaction is expected. They are evolved through the first mass transfer and the donor star death, a possible ensuing Be/X-ray binary phase, and they end when the mass gainer leaves the main sequence. Results.In our fiducial synthetic population, 8% of the OB stars in the SMC are post-mass transfer systems, and 7% are merger products. In many of our models, the mass gainers are spun up and form Oe/Be stars. While our model underpredicts the number of Be/X-ray binaries in the SMC, it reproduces the main features of their orbital period distribution and the observed number of SMC binary WR stars. We expect $\sim$50 OB+BH binaries below and $\sim$170 above 20d orbital period. The latter might produce merging double BHs. However, their progenitors, the predicted long-period WR+OB binaries, are not observed. Conclusions. While the comparison with the observed SMC stars supports many physics assumptions in our high-mass binary models, a better match of the large number of observed OBe stars and Be/X-ray binaries likely requires a lower merger rate and/or a higher mass transfer efficiency during the first mass transfer. The fate of the initially wide O star binaries remains uncertain.

preprint2025arXiv

Populations of evolved massive binary stars in the Small Magellanic Cloud II: Predictions from rapid binary evolution

Massive star evolution plays a crucial role in astrophysics but bares large uncertainties. This problem becomes more severe by the majority of massive stars being born in close binary systems, whose evolution is affected by the interaction of their components. We want to constrain major uncertainties in massive binary star evolution, in particular the efficiency and the stability of the first mass transfer phase. We use the rapid population synthesis code ComBinE to generate synthetic populations of post-interaction binaries, assuming constant mass-transfer efficiency. We employ a new merger criterion that adjusts self-consistently to any prescribed mass-transfer efficiency. We tailor our synthetic populations to be comparable to the expected binary populations in the Small Magellanic Cloud (SMC). We find that the observed populations of evolved massive binaries can not be reproduced with a single mass-transfer efficiency. Instead, a rather high efficiency (>50%) is needed to reproduce the number of Be stars and Be/X-ray binaries in the SMC, while a low efficiency (~10%) leads to a better agreement with the observed number of Wolf-Rayet stars. We construct a corresponding mass-dependent mass-transfer efficiency recipe to produce our fiducial synthetic SMC post-interaction binary population. It reproduces the observed number and properties of the Be/X-ray and WR-binaries rather well, and is not in stark disagreement with the observed OBe star population. It further predicts two large, yet unobserved populations of OB+BH binaries, that is ~100 OB+BH systems with rather small orbital periods (<20 days) and ~40 longer period OBe+BH systems.

preprint2023arXiv

MeerKAT discovery of 13 new pulsars in Omega Centauri

The most massive globular cluster in our Galaxy, Omega Centauri, is an interesting target for pulsar searches, because of its multiple stellar populations and the intriguing possibility that it was once the nucleus of a galaxy that was absorbed into the Milky Way. The recent discoveries of pulsars in this globular cluster and their association with known X-ray sources was a hint that, given the large number of known X-ray sources, there is a much larger undiscovered pulsar population. We used the superior sensitivity of the MeerKAT radio telescope to search for pulsars in Omega Centauri. In this paper, we present some of the first results of this survey, including the discovery of 13 new pulsars; the total number of known pulsars in this cluster currently stands at 18. At least half of them are in binary systems and preliminary orbital constraints suggest that most of the binaries have light companions. We also discuss the ratio between isolated and binaries pulsars and how they were formed in this cluster.

preprint2022arXiv

A Gamma-ray Pulsar Timing Array Constrains the Nanohertz Gravitational Wave Background

After large galaxies merge, their central supermassive black holes are expected to form binary systems whose orbital motion generates a gravitational wave background (GWB) at nanohertz frequencies. Searches for this background utilize pulsar timing arrays, which perform long-term monitoring of millisecond pulsars (MSPs) at radio wavelengths. We use 12.5 years of Fermi Large Area Telescope data to form a gamma-ray pulsar timing array. Results from 35 bright gamma-ray pulsars place a 95\% credible limit on the GWB characteristic strain of $1.0\times10^{-14}$ at 1 yr$^{-1}$, which scales as the observing time span $t_{\mathrm{obs}}^{-13/6}$. This direct measurement provides an independent probe of the GWB while offering a check on radio noise models.

preprint2022arXiv

Arecibo and FAST Timing Follow-up of twelve Millisecond Pulsars Discovered in Commensal Radio Astronomy FAST Survey

We report the phase-connected timing ephemeris, polarization pulse profiles, Faraday rotation measurements, and Rotating-Vector-Model (RVM) fitting results of twelve millisecond pulsars (MSPs) discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in the Commensal radio Astronomy FAST survey (CRAFTS). The timing campaigns were carried out with FAST and Arecibo over three years. Eleven of the twelve pulsars are in neutron star - white dwarf binary systems, with orbital periods between 2.4 and 100 d. Ten of them have spin periods, companion masses, and orbital eccentricities that are consistent with the theoretical expectations for MSP - Helium white dwarf (He WD) systems. The last binary pulsar (PSR J1912$-$0952) has a significantly smaller spin frequency and a smaller companion mass, the latter could be caused by a low orbital inclination for the system. Its orbital period of 29 days is well within the range of orbital periods where some MSP - He WD systems have shown anomalous eccentricities, however, the eccentricity of PSR J1912$-$0952 is typical of what one finds for the remaining MSP - He WD systems.

preprint2022arXiv

Detection of quasi-periodic micro-structure in three millisecond pulsars with the Large European Array for Pulsars

We report on the detection of quasi-periodic micro-structure in three millisecond pulsars (MSPs), PSRs J1022+1001, J2145-0750 and J1744-1134, using high time resolution data acquired with the Large European Array for Pulsars at a radio frequency of 1.4 GHz. The occurrence rate of quasi-periodic micro-structure is consistent among pulses with different peak flux densities. Using an auto-correlation analysis, we measure the periodicity and width of the micro-structure in these three pulsars. The detected micro-structure from PSRs J1022+1001 and J1744-1134 is often highly linearly polarised. In PSR J1022+1001, the linear polarisation position angles of micro-structure pulses are in general flat with a small degree of variation. Using these results, we further examine the frequency and rotational period dependency of micro-structure properties established in previous work, along with the angular beaming and temporal modulation models that explains the appearance of micro-structure. We also discuss a possible link of micro-structure to the properties of some of the recently discovered fast radio bursts which exhibit a very similar emission morphology.

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

First discoveries and localisations of Fast Radio Bursts with MeerTRAP: a real-time, commensal MeerKAT survey

We report on the discovery and localization of fast radio bursts (FRBs) from the MeerTRAP project, a commensal fast radio transient-detection programme at MeerKAT in South Africa. Our hybrid approach combines a coherent search with an average field-of-view of 0.4 $\rm deg^{2}$ with an incoherent search utilizing a field-of-view of $\sim$1.27 $\rm deg^{2}$ (both at 1284~MHz). Here, we present results on the first three FRBs: FRB 20200413A (DM=1990.05 pc cm$^{-3}$), FRB 20200915A (DM=740.65 pc cm$^{-3}$), and FRB 20201123A (DM=433.55 pc cm$^{-3}$). FRB 20200413A was discovered only in the incoherent beam. FRB 20200915A (also discovered only in the incoherent beam) shows speckled emission in the dynamic spectrum which cannot be explained by interstellar scintillation in our Galaxy or plasma lensing, and might be intrinsic to the source. FRB 20201123A shows a faint post-cursor burst about 200 ms after the main burst and warrants further follow-up to confirm whether it is a repeating FRB. FRB 20201123A also exhibits significant temporal broadening consistent with scattering by a turbulent medium. The broadening exceeds that predicted for medium along the sightline through our Galaxy. We associate this scattering with the turbulent medium in the environment of the FRB in the host galaxy. Within the approximately $1'$ localization region of FRB 20201123A, we identify one luminous galaxy ($r \approx 15.67$; J173438.35$-$504550.4) that dominates the posterior probability for a host association. The galaxy's measured properties are consistent with other FRB hosts with secure associations.

preprint2022arXiv

First measurement of high-energy reactor antineutrinos at Daya Bay

This Letter reports the first measurement of high-energy reactor antineutrinos at Daya Bay, with nearly 9000 inverse beta decay candidates in the prompt energy region of 8-12~MeV observed over 1958 days of data collection. A multivariate analysis is used to separate 2500 signal events from background statistically. The hypothesis of no reactor antineutrinos with neutrino energy above 10~MeV is rejected with a significance of 6.2 standard deviations. A 29\% antineutrino flux deficit in the prompt energy region of 8-11~MeV is observed compared to a recent model prediction. We provide the unfolded antineutrino spectrum above 7 MeV as a data-based reference for other experiments. This result provides the first direct observation of the production of antineutrinos from several high-$Q_β$ isotopes in commercial reactors.

preprint2022arXiv

Four pulsar discoveries in NGC 6624 by TRAPUM using MeerKAT

We report 4 new pulsars discovered in the core-collapsed globular cluster (GC) NGC 6624 by the TRAPUM Large Survey Project with the MeerKAT telescope. All of the new pulsars found are isolated. PSR J1823$-$3021I and PSR J1823$-$3021K are millisecond pulsars with period of respectively 4.319 ms and 2.768 ms. PSR J1823$-$3021J is mildly recycled with a period of 20.899 ms, and PSR J1823$-$3022 is a long period pulsar with a period of 2.497 s. The pulsars J1823$-$3021I, J1823$-$3021J, and J1823$-$3021K have position and dispersion measure (DM) compatible with being members of the GC and are therefore associated with NGC 6624. Pulsar J1823$-$3022 is the only pulsar bright enough to be re-detected in archival observations of the cluster. This allowed the determination of a timing solution that spans over two decades. It is not possible at the moment to claim the association of pulsar J1823$-$3022 with the GC given the long period and large offset in position ($\sim 3$ arcminutes) and DM (with a fractional difference of 11 percent compared the average of the pulsars in NGC 6624). The discoveries made use of the beamforming capability of the TRAPUM backend to generate multiple beams in the same field of view which allows sensitive searches to be performed over a few half-light radii from the cluster center and can simultaneously localise the discoveries. The discoveries reflect the properties expected for pulsars in core-collapsed GCs.

preprint2022arXiv

Joint Determination of Reactor Antineutrino Spectra from $^{235}$U and $^{239}$Pu Fission by Daya Bay and PROSPECT

A joint determination of the reactor antineutrino spectra resulting from the fission of $^{235}$U and $^{239}$Pu has been carried out by the Daya Bay and PROSPECT collaborations. This Letter reports the level of consistency of $^{235}$U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant $^{235}$U and $^{239}$Pu isotopes and improves the uncertainty of the $^{235}$U spectral shape to about 3\%. The ${}^{235}$U and $^{239}$Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the $^{235}$U and $^{239}$Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

preprint2022arXiv

Long term radio and X-ray evolution of the magnetar Swift J1818.0-1607

We report on the the long term monitoring campaign of the seemingly youngest magnetar Swift~J1818.0-1607 at radio and X-ray wavelengths over a span of one year. We obtained a coherent timing solution for the magnetar over the same time span. The frequency derivative of the magnetar shows systematic variation with the values oscillating about a mean value of $-$1.37$\times$10$^{-11}$ Hz s$^{-1}$. The magnitude of the variation in the frequency derivative reduces with time before converging on the mean value. We were able to identify four states in the spin-frequency derivative that were quantified by the amount of modulation about the mean value and the transition between these states seem to be correlated with the change in the radio emission of the magnetar while no correlation is seen in the average radio profile variability on a shorter timescale (days). The 0.5--12 keV X-ray flux shows a monotonic decrease that can be attributed to thermal emission from a hot-spot on the surface of the neutron star that is reducing in size. Such decrease is consistent with what is seen in other magnetars. The potential correlation between the radio emission mode and the behaviour of the spin-down rate hints to a global change in the magnetopshere of the magnetar akin to the correlation seen in a subset of mode-changing radio pulsars and suggests a physical link between the two sub-populations.

preprint2022arXiv

Long-term rotational and emission variability of 17 radio pulsars

With the ever-increasing sensitivity and timing baselines of modern radio telescopes, a growing number of pulsars are being shown to exhibit transitions in their rotational and radio emission properties. In many of these cases, the two are correlated with pulsars assuming a unique spin-down rate ($\dotν$) for each of their specific emission states. In this work we revisit 17 radio pulsars previously shown to exhibit spin-down rate variations. Using a Gaussian process regression (GPR) method to model the timing residuals and the evolution of the profile shape, we confirm the transitions already observed and reveal new transitions in 8 years of extended monitoring with greater time resolution and enhanced observing bandwidth. We confirm that 7 of these sources show emission-correlated $\dotν$ transitions ($Δ\dotν$) and we characterise this correlation for one additional pulsar, PSR B1642$-$03. We demonstrate that GPR is able to reveal extremely subtle profile variations given sufficient data quality. We also corroborate the dependence of $Δ\dotν$ amplitude on $\dotν$ and pulsar characteristic age. Linking $Δ\dotν$ to changes in the global magnetospheric charge density $Δρ$, we speculate that $\dotν$ transitions associated with large $Δρ$ values may be exhibiting detectable profile changes with improved data quality, in cases where they have not previously been observed.

preprint2022arXiv

MeerTRAP: Twelve Galactic fast transients detected in a real-time, commensal MeerKAT survey

MeerTRAP is a real-time untargeted search project using the MeerKAT telescope to find single pulses from fast radio transients and pulsars. It is performed commensally with the MeerKAT large survey projects (LSPs), using data from up to 64 of MeerKAT's 13.96~m dishes to form hundreds of coherent beams on sky, each of which is processed in real time to search for millisecond-duration pulses. We present the first twelve Galactic sources discovered by MeerTRAP, with DMs in the range of 33--381~pc~cm$^{-3}$. One source may be Galactic or extragalactic depending on the Galactic electron density model assumed. Follow-up observations performed with the MeerKAT, Lovell, and Parkes radio telescopes have detected repeat pulses from seven of the twelve sources. Pulse periods have been determined for four sources. Another four sources could be localised to the arcsecond-level using a novel implementation of the tied-array beam localisation method.

preprint2022arXiv

News and views regarding PSR J1757-1854, a highly-relativistic binary pulsar

We provide an update on the ongoing monitoring and study of the highly-relativistic double neutron star binary system PSR J1757-1854, a 21.5-ms pulsar in a highly eccentric, 4.4-hour orbit. The extreme nature of this pulsar's orbit allows it to probe a parameter space largely unexplored by other relativistic binary pulsars. For example, it displays one of the highest gravitational wave (GW) luminosities of any known binary pulsar, as well as the highest rate or orbital decay due to GW damping. PSR J1757-1854 is also notable in that it is an excellent candidate for exploring new tests of General Relativity and other gravitational theories, with possible measurements of both Lense-Thirring precession and relativistic orbital deformation (through the post-Keplerian parameter $δ_θ$) anticipated within the next 3-5 years. Here we present a summary of the latest interim results from the ongoing monitoring of this pulsar as part of an international, multi-telescope campaign. This includes an update of the pulsar's long-term timing and post-Keplerian parameters, new constraints on the pulsar's proper motion and corresponding Shklovskii kinematic correction, and new limits on the pulsar's geodetic precession as determined by monitoring for secular changes in the pulse profile. We also highlight prospects for future work, including an updated timeline on new relativistic tests following the introduction of MeerKAT observations.

preprint2022arXiv

The Case for Space Environmentalism

The shell bound by the Karman line at a height of 80 to 100km above the Earth's surface, and Geosynchronous Orbit, at 36,000km, is defined as the orbital space surrounding the Earth. It is within this region, and especially in Low Earth Orbit (LEO), where environmental issues are becoming urgent because of the rapid growth of the anthropogenic space object population, including satellite "mega-constellations". In this Perspective, we summarise the case that the orbital space around the Earth should be considered an additional ecosystem, and so subject to the same care and concerns and the same broad regulations as, for example, the oceans and the atmosphere. We rely on the orbital space environment by looking through it as well as by working within it. Hence, we should consider damage to professional astronomy, public stargazing and the cultural importance of the sky, as well as the sustainability of commercial, civic and military activity in space. Damage to the orbital space environment has problematic features in common with other types of environmental issue. First, the observed and predicted damage is incremental and complex, with many contributors. Second, whether or not space is formally and legally seen as a global commons, the growing commercial exploitation of what may appear a "free" resource is in fact externalising the true costs.

preprint2022arXiv

The eccentric millisecond pulsar, PSR J0955$-$6150 I: Pulse profile analysis, mass measurements and constraints on binary evolution

PSR J0955$-$6150 is a member of a class of eccentric MSP+He WD systems (eMSPs), whose binary evolution is poorly understood and believed to be different to that of traditional MSP+He WD systems. Measuring the masses of the stars in this system is important for testing hypotheses for the formation of eMSPs. We have carried out observations of this pulsar with the Parkes and MeerKAT radio telescopes. Our observations reveal a strong frequency evolution of this pulsar's intensity, with a spectral index ($α$) of $-3.13(2)$. The sensitivity of MeerKAT has resulted in a $>10$-fold improvement in the timing precision compared to older Parkes observations. Combined with the 8-year timing baseline, it has allowed precise measurements of a proper motion and three orbital "post-Keplerian" parameters: the rate of advance of periastron, $\dotω = 0.00152(1) \, {\rm deg} \, yr^{-1}$ and the orthometric Shapiro delay parameters, $h_3 = 0.89(7) \, μ$s and $ς= 0.88(2)$. Assuming general relativity, we obtain $M_{p} = 1.71(2) \, M_{\odot}$ for the mass of the pulsar and $M_{c} = 0.254(2) \, M_{\odot}$ for the mass of the companion; the orbital inclination is 83.2(4) degrees. We find that the spin axis has a misalignment relative to the orbital angular momentum of $> 4.8$ degrees at 99% CI. While the value of $M_{\rm p}$ falls within the wide range observed in eMSPs, $M_{\rm c}$ is significantly smaller than expected, allowing several formation hypotheses being ruled out. $M_{\rm c}$ is also significantly different from the expected value for an ideal low mass X-ray binary evolution scenario. The putative misalignment between the spin axis of the pulsar and the orbital angular momentum suggests that the unknown process that created the orbital eccentricity of the binary was also capable of changing its orbital orientation, an important evidence for understanding the origin of eMSPs.

preprint2022arXiv

The High Time Resolution Universe Pulsar Survey -- XVII. PSR J1325-6253, a low eccentricity double neutron star system from an ultra-stripped supernova

The observable population of double neutron star (DNS) systems in the Milky Way allow us to understand the nature of supernovae and binary stellar evolution. Until now, all DNS systems in wide orbits ($ P_{\textrm{orb}}>$ 1~day) have been found to have orbital eccentricities, $e > 0.1$. In this paper, we report the discovery of pulsar PSR J1325$-$6253: a DNS system in a 1.81 day orbit with a surprisingly low eccentricity of just $e = 0.064$. Through 1.4 yr of dedicated timing with the Parkes radio telescope we have been able to measure its rate of advance of periastron, $\dotω=0.138 \pm 0.002$ $\rm deg$ $\rm yr^{-1}$. If this induced $\dotω$ is solely due to general relativity then the total mass of the system is, $M_{\rm sys} = 2.57 \pm 0.06$ M$_{\odot}$. Assuming an edge-on orbit the minimum companion mass is constrained to be $M_\mathrm{c,min}>0.98$ M$_{\odot}$ which implies the pulsar mass is $M_\mathrm{p,max}<1.59 $ M$_{\odot}$. Its location in the $P$-$\dot{P}$ diagram suggests that, like other DNS systems, PSR J1325$-$6253 is a recycled pulsar and if its mass is similar to the known examples ($>1.3$ M$_\odot$), then the companion neutron star is probably less than $\sim1.25$ M$_\odot$ and the system is inclined at about $50^{\circ}$-$60^{\circ}$. The low eccentricity along with the wide orbit of the system strongly favours a formation scenario involving an ultra-stripped supernova explosion.

preprint2022arXiv

The International Pulsar Timing Array second data release: Search for an isotropic Gravitational Wave Background

We searched for an isotropic stochastic gravitational wave background in the second data release of the International Pulsar Timing Array, a global collaboration synthesizing decadal-length pulsar-timing campaigns in North America, Europe, and Australia. In our reference search for a power law strain spectrum of the form $h_c = A(f/1\,\mathrm{yr}^{-1})^α$, we found strong evidence for a spectrally-similar low-frequency stochastic process of amplitude $A = 3.8^{+6.3}_{-2.5}\times10^{-15}$ and spectral index $α= -0.5 \pm 0.5$, where the uncertainties represent 95\% credible regions, using information from the auto- and cross-correlation terms between the pulsars in the array. For a spectral index of $α= -2/3$, as expected from a population of inspiralling supermassive black hole binaries, the recovered amplitude is $A = 2.8^{+1.2}_{-0.8}\times10^{-15}$. Nonetheless, no significant evidence of the Hellings-Downs correlations that would indicate a gravitational-wave origin was found. We also analyzed the constituent data from the individual pulsar timing arrays in a consistent way, and clearly demonstrate that the combined international data set is more sensitive. Furthermore, we demonstrate that this combined data set produces comparable constraints to recent single-array data sets which have more data than the constituent parts of the combination. Future international data releases will deliver increased sensitivity to gravitational wave radiation, and significantly increase the detection probability.

preprint2022arXiv

The MeerTime Pulsar Timing Array -- A Census of Emission Properties and Timing Potential

MeerTime is a five-year Large Survey Project to time pulsars with MeerKAT, the 64-dish South African precursor to the Square Kilometre Array. The science goals for the programme include timing millisecond pulsars (MSPs) to high precision (< 1 $μ$s) to study the Galactic MSP population and to contribute to global efforts to detect nanohertz gravitational waves with the International Pulsar Timing Array (IPTA). In order to plan for the remainder of the programme and to use the allocated time most efficiently, we have conducted an initial census with the MeerKAT "L-band" receiver of 189 MSPs visible to MeerKAT and here present their dispersion measures, polarization profiles, polarization fractions, rotation measures, flux density measurements, spectral indices, and timing potential. As all of these observations are taken with the same instrument (which uses coherent dedispersion, interferometric polarization calibration techniques, and a uniform flux scale), they present an excellent resource for population studies. We used wideband pulse portraits as timing standards for each MSP and demonstrated that the MeerTime Pulsar Timing Array (MPTA) can already contribute significantly to the IPTA as it currently achieves better than 1 $μ$s timing accuracy on 89 MSPs (observed with fortnightly cadence). By the conclusion of the initial five-year MeerTime programme in July 2024, the MPTA will be extremely significant in global efforts to detect the gravitational wave background with a contribution to the detection statistic comparable to other long-standing timing programmes.

preprint2022arXiv

TRAPUM discovery of thirteen new pulsars in NGC 1851 using MeerKAT

We report the discovery of 13 new pulsars in the globular cluster NGC 1851 by the TRAPUM Large Survey Project using the MeerKAT radio telescope. The discoveries consist of six isolated millisecond pulsars (MSPs) and seven binary pulsars, of which six are MSPs and one is mildly recycled. For all the pulsars, we present the basic kinematic, astrometric, and orbital parameters, where applicable, as well as their polarimetric properties, when these are measurable. Two of the binary MSPs (PSR J0514-4002D and PSR J0514-4002E) are in wide and extremely eccentric (e > 0.7) orbits with a heavy white dwarf and a neutron star as their companion, respectively. With these discoveries, NGC 1851 is now tied with M28 as the cluster with the third largest number of known pulsars (14). Its pulsar population shows remarkable similarities with that of M28, Terzan 5 and other clusters with comparable structural parameters. The newly-found pulsars are all located in the innermost regions of NGC 1851 and will likely enable, among other things, detailed studies of the cluster structure and dynamics.

preprint2022arXiv

Upgraded GMRT survey for pulsars in globular clusters. I: Discovery of a millisecond binary pulsar in NGC 6652

Globular clusters contain a unique pulsar population, with many exotic systems that can form only in their dense stellar environments. The leap in sensitivity of the upgraded Giant Metrewave Radio Telescope (uGMRT) in India, especially at low radio frequencies ($<$ 1 GHz) has motivated a new search for radio pulsars in a group of eight Southern globular clusters. We discovered PSR J1835$-$3259B, a 1.83-ms pulsar in NGC 6652; this is in a near-circular wide orbit of 28.7 hr with a low-mass ($ \sim 0.2 \, M_{\rm \odot}$) companion, likely a Helium white dwarf. We derived a 10-year timing solution for this system. We also present measurements of scattering, flux densities and spectral indices for some of the previously known pulsars in these GCs. A significant fraction of the pulsars in these clusters have steep spectral indices. Additionally, we detected eight radio point sources not associated with any known pulsar positions in the radio images. There are four newly identified sources, three in NGC 6652 and one in NGC 6539, and one previously identified source each in NGC 1851, NGC 6440, NGC 6544, and Terzan 5. Surprisingly, our images show that our newly discovered pulsar, PSR J1835$-$3259B, is the brightest pulsar in all GCs we have imaged; like other pulsars with broad profiles (Terzan 5 C and O), its flux density in the radio images is much larger than in its pulsations. This indicates that their pulsed emission is only a fraction of their total emission. The detection of radio sources outside the core radii but well within the tidal radii of these clusters show that future GC surveys should complement the search analysis by using the imaging capability of interferometers, and preferentially synthesize large number of search beams in order to obtain a larger field of view.

preprint2021arXiv

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst.

preprint2021arXiv

Measurements of pulse jitter and single-pulse variability in millisecond pulsars using MeerKAT

Using the state-of-the-art SKA precursor, the MeerKAT radio telescope, we explore the limits to precision pulsar timing of millisecond pulsars achievable due to pulse stochasticity (jitter). We report new jitter measurements in 15 of the 29 pulsars in our sample and find that the levels of jitter can vary dramatically between them. For some, like the 2.2~ms pulsar PSR J2241--5236, we measure an implied jitter of just $\sim$ 4~ns/hr, while others like the 3.9~ms PSR J0636--3044 are limited to $\sim$ 100 ns/hr. While it is well known that jitter plays a central role to limiting the precision measurements of arrival times for high signal-to-noise ratio observations, its role in the measurement of dispersion measure (DM) has not been reported, particularly in broad-band observations. Using the exceptional sensitivity of MeerKAT, we explored this on the bright millisecond pulsar PSR J0437--4715 by exploring the DM of literally every pulse. We found that the derived single pulse DMs vary by typically 0.0085 cm$^{-3}$ pc from the mean, and that the best DM estimate is limited by the differential pulse jitter across the band. We postulate that all millisecond pulsars will have their own limit on DM precision which can only be overcome with longer integrations. Using high-time resolution filterbank data of 9 $μ$s, we also present a statistical analysis of single pulse phenomenology. Finally, we discuss optimization strategies for the MeerKAT pulsar timing program and its role in the context of the International Pulsar Timing Array (IPTA).

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.

preprint2020arXiv

Detection of the magnetar XTE J1810-197 at 150 and 260 GHz with the NIKA2 Kinetic Inductance Detector camera

The investigation of pulsars between millimetre and optical wavelengths is challenging due to the faintness of the pulsar signals and the relative low sensitivity of the available facilities compared to 100-m class telescopes operating in the centimetre band. The Kinetic Inductance Detector (KID) technology offers large instantaneous bandwidths and a high sensitivity that can help to substantially increase the ability of existing observatories at short wavelengths to detect pulsars and transient emission. To investigate the feasibility of detecting pulsars with KIDs, we observed the anomalous X-ray pulsar XTE J1810-197 with the New IRAM KIDs Array-2 (NIKA2) camera installed at the IRAM 30-m Telescope in Spain. We detected the pulsations from the pulsar with NIKA2 at its two operating frequency bands, 150 and 260 GHz ($λ$=2.0 and 1.15 mm, respectively). This is the first time that a pulsar is detected with a receiver based on KID technology in the millimetre band. In addition, this is the first report of short millimetre emission from XTE J1810-197 after its reactivation in December 2018, and it is the first time that the source is detected at 260 GHz, which gives us new insights into the radio emission process of the star.

preprint2020arXiv

Dispersion measure variability for 36 millisecond pulsars at 150MHz with LOFAR

Radio pulses from pulsars are affected by plasma dispersion, which results in a frequency-dependent propagation delay. Variations in the magnitude of this effect lead to an additional source of red noise in pulsar timing experiments, including pulsar timing arrays that aim to detect nanohertz gravitational waves. We aim to quantify the time-variable dispersion with much improved precision and characterise the spectrum of these variations. We use the pulsar timing technique to obtain highly precise dispersion measure (DM) time series. Our dataset consists of observations of 36 millisecond pulsars, which were observed for up to 7.1 years with the LOFAR telescope at a centre frequency of ~150 MHz. Seventeen of these sources were observed with a weekly cadence, while the rest were observed at monthly cadence. We achieve a median DM precision of the order of 10^-5 cm^-3 pc for a significant fraction of our sources. We detect significant variations of the DM in all pulsars with a median DM uncertainty of less than 2x10^-4 cm^-3 pc. The noise contribution to pulsar timing experiments at higher frequencies is calculated to be at a level of 0.1-10 us at 1.4 GHz over a timespan of a few years, which is in many cases larger than the typical timing precision of 1 us or better that PTAs aim for. We found no evidence for a dependence of DM on radio frequency for any of the sources in our sample. The DM time series we obtained using LOFAR could in principle be used to correct higher-frequency data for the variations of the dispersive delay. However, there is currently the practical restriction that pulsars tend to provide either highly precise times of arrival (ToAs) at 1.4 GHz or a high DM precision at low frequencies, but not both, due to spectral properties. Combining the higher-frequency ToAs with those from LOFAR to measure the infinite-frequency ToA and DM would improve the result.

preprint2020arXiv

Giant pulses from J1823-3021A observed with the MeerKAT telescope

The millisecond pulsar J1823-3021A is a very active giant pulse emitter in the globular cluster NGC 6624. New observations with the MeerKAT radio telescope have revealed 14350 giant pulses over 5 hours of integration time, with an average wait time of about 1 second between giant pulses. The giant pulses occur in phases compatible with the ordinary radio emission, follow a power-law distribution with an index of -2.63 $\pm$ 0.02 and contribute 4 percent of the total integrated flux. The spectral index of the giant pulses follows a Gaussian distribution centered around -1.9 with a standard deviation of 0.6 and is on average flatter than the integrated emission, which has a spectral index of -2.81 $\pm$ 0.02. The waiting times between the GPs are accurately described by a Poissonian distribution, suggesting that the time of occurrence of a GP is independent from the times of occurrence of other GPs. 76 GPs show multiple peaks within the same rotation, a rate that is also compatible with the mutual independence of the GP times of occurrence. We studied the polarization properties of the giant pulses finding, on average, linear polarization only at the 1 percent level and circular polarization at the 3 percent level, similar to the polarization percentages of the total integrated emission. In 4 cases it was possible to measure the RM of the GPs which are highly variable and, in two cases, is inconsistent with the mean RM of the total integrated pulsar signal.

preprint2020arXiv

Improved Constraints on Sterile Neutrino Mixing from Disappearance Searches in the MINOS, MINOS+, Daya Bay, and Bugey-3 Experiments

Searches for electron antineutrino, muon neutrino, and muon antineutrino disappearance driven by sterile neutrino mixing have been carried out by the Daya Bay and MINOS+ collaborations. This Letter presents the combined results of these searches, along with exclusion results from the Bugey-3 reactor experiment, framed in a minimally extended four-neutrino scenario. Significantly improved constraints on the $θ_{μe}$ mixing angle are derived that constitute the most stringent limits to date over five orders of magnitude in the sterile mass-squared splitting $Δm^2_{41}$, excluding the 90% C.L. sterile-neutrino parameter space allowed by the LSND and MiniBooNE observations at 90% CL$_s$ for $Δm^2_{41}<5\,$eV$^2$.Furthermore, the LSND and MiniBooNE 99% C.L. allowed regions are excluded at 99% CL$_s$ for $Δm^2_{41}$ $<$ 1.2 eV$^2$.

preprint2020arXiv

Optimization of the JUNO liquid scintillator composition using a Daya Bay antineutrino detector

To maximize the light yield of the liquid scintillator (LS) for the Jiangmen Underground Neutrino Observatory (JUNO), a 20 t LS sample was produced in a pilot plant at Daya Bay. The optical properties of the new LS in various compositions were studied by replacing the gadolinium-loaded LS in one antineutrino detector. The concentrations of the fluor, PPO, and the wavelength shifter, bis-MSB, were increased in 12 steps from 0.5 g/L and <0.01 mg/L to 4 g/L and 13 mg/L, respectively. The numbers of total detected photoelectrons suggest that, with the optically purified solvent, the bis-MSB concentration does not need to be more than 4 mg/L. To bridge the one order of magnitude in the detector size difference between Daya Bay and JUNO, the Daya Bay data were used to tune the parameters of a newly developed optical model. Then, the model and tuned parameters were used in the JUNO simulation. This enabled to determine the optimal composition for the JUNO LS: purified solvent LAB with 2.5 g/L PPO, and 1 to 4 mg/L bis-MSB.

preprint2020arXiv

Radio Afterglows of Very High Energy Gamma-Ray Bursts 190829A and 180720B

We present high cadence multi-frequency radio observations of the long Gamma-Ray Burst (GRB) 190829A, which was detected at photon energies above 100 GeV by the High Energy Stereoscopic System (H.E.S.S.). Observations with the Meer Karoo Array Telescope (MeerKAT, 1.3 GHz), and Arcminute Microkelvin Imager - Large Array (AMI-LA, 15.5 GHz) began one day post-burst and lasted nearly 200 days. We used complementary data from Swift X-Ray Telescope (XRT), which ran to 100 days post-burst. We detected a likely forward shock component with both MeerKAT and XRT up to over 100 days post-burst. Conversely, the AMI-LA light curve appears to be dominated by reverse shock emission until around 70 days post-burst when the afterglow flux drops below the level of the host galaxy. We also present previously unpublished observations of the other H.E.S.S.-detected GRB, GRB 180720B from AMI-LA, which shows likely forward shock emission that fades in less than 10 days. We present a comparison between the radio emission from the three GRBs with detected very high energy (VHE) gamma-ray emission and a sensitivity-limited radio afterglow sample. GRB 190829A has the lowest isotropic radio luminosity of any GRB in our sample, but the distribution of luminosities is otherwise consistent, as expected, with the VHE GRBs being drawn from the same parent distribution as the other radio-detected long GRBs.

preprint2020arXiv

Shapiro delay in PSR~J1811$-$2405

This paper presents the first detection of Shapiro delay from the binary millisecond pulsar PSR J1811-2405. We report a 11-$σ$ measurement of the orthometric amplitude, h$_3$=6.8(6)x10$^{-7}$, and a 16-$σ$ measurement of the orthometric ratio, $ς$= 0.81(5). Given the relatively high orbital inclination, i= 79(2) deg, of this binary system, we obtain constraints on the companion mass of m$_c$= 0.31$\substack{+0.08 \\ -0.06}M_\odot$. The pulsar mass is currently less well constrained, with a value of 2.0$\substack{+0.8 \\ -0.5}M_\odot$. The companion mass and the orbital period are in agreement with the prediction made by previous numerical calculations of the evolution of compact binary systems. From a study of the polarization, we find that the orbital inclination angle is ~100 deg and that PSR J1811-2405 is an orthogonal rotator. In addition, the $μ$s-level timing precision together with its narrow profile make PSR J1811-2405 a good candidate for inclusion in the pulsar timing arrays being used to detect nHz Gravitational waves.

preprint2020arXiv

Simultaneous multi-telescope observations of FRB 121102

We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and resulted in 4 simultaneous detections. The observations with MeerKAT's wide band receiver, which extends down to relatively low frequencies (900-1670 MHz usable L-band range), have allowed us to get a detailed look at the complex frequency structure, intensity variations and frequency-dependent sub-pulse drifting. The drift rates we measure for the full-band and sub-banded data are consistent with those published between 600-6500 MHz with a slope of -0.147 +/- 0.014 ms^-1. Two of the detected bursts exhibit fainter 'precursors' separated from the brighter main pulse by ~28 ms and ~34 ms. A follow-up multi-telescope campaign on the 6th and 8th October 2019 to better understand these frequency drifts and structures over a wide and continuous band was undertaken. No detections resulted, indicating that the source was 'inactive' over a broad frequency range during this time.

preprint2020arXiv

The Galactic population and properties of young, highly-energetic pulsars

The population of young, non-recycled pulsars with spin down energies Edot >10^35 erg/s is sampled predominantly at gamma-ray and radio wavelengths. A total of 137 such pulsars are known, with partial overlap between the sources detectable in radio and gamma-rays. We use a very small set of assumptions in an attempt to test whether the observed pulsar sample can be explained by a single underlying population of neutron stars. For radio emission we assume a canonical conal beam with a fixed emission height of 300~km across all spin periods and a luminosity law which depends on Edot^{0.25}. For gamma-ray emission we assume the outer-gap model and a luminosity law which depends on Edot^{0.5}. We synthesise a population of fast-spinning pulsars with a birth rate of one per 100 years. We find that this simple model can reproduce most characteristics of the observed population with two caveats. The first is a deficit of gamma-ray pulsars at the highest Edot which we surmise to be an observational selection effect due to the difficulties of finding gamma-ray pulsars in the presence of glitches without prior knowledge from radio frequencies. The second is a deficit of radio pulsars with interpulse emission, which may be related to radio emission physics. We discuss the implications of these findings.

preprint2020arXiv

The High Time Resolution Universe Pulsar Survey -- XVI. Discovery and timing of 40 pulsars from the southern Galactic plane

We present the results of processing an additional 44% of the High Time Resolution Universe South Low Latitude (HTRU-S LowLat) pulsar survey, the most sensitive blind pulsar survey of the southern Galactic plane to date. Our partially-coherent segmented acceleration search pipeline is designed to enable the discovery of pulsars in short, highly-accelerated orbits, while our 72-min integration lengths will allow us to discover pulsars at the lower end of the pulsar luminosity distribution. We report the discovery of 40 pulsars, including three millisecond pulsar-white dwarf binary systems (PSRs J1537-5312, J1547-5709 and J1618-4624), a black-widow binary system (PSR J1745-23) and a candidate black-widow binary system (PSR J1727-2951), a glitching pulsar (PSR J1706-4434), an eclipsing binary pulsar with a 1.5-yr orbital period (PSR J1653-45), and a pair of long spin-period binary pulsars which display either nulling or intermittent behaviour (PSRs J1812-15 and J1831-04). We show that the total population of 100 pulsars discovered in the HTRU-S LowLat survey to date represents both an older and lower-luminosity population, and indicates that we have yet to reach the bottom of the luminosity distribution function. We present evaluations of the performance of our search technique and of the overall yield of the survey, considering the 94% of the survey which we have processed to date. We show that our pulsar yield falls below earlier predictions by approximately 25% (especially in the case of millisecond pulsars), and discuss explanations for this discrepancy as well as future adaptations in RFI mitigation and searching techniques which may address these shortfalls.

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.

preprint2020arXiv

The MeerKAT Telescope as a Pulsar Facility: System verification and early science results from MeerTime

We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

preprint2020arXiv

The SUrvey for Pulsars and Extragalactic Radio Bursts IV: Discovery and polarimetry of a 12.1-second radio pulsar

We report the discovery of PSR~J2251$-$3711, a radio pulsar with a spin period of 12.1 seconds, the second longest currently known. Its timing parameters imply a characteristic age of 15 Myr, a surface magnetic field of $1.3 \times 10^{13}$~G and a spin-down luminosity of $2.9 \times 10^{29}~\mathrm{erg~s}^{-1}$. Its dispersion measure of 12.12(1)~$\mathrm{pc}~\mathrm{cm}^{-3}$ leads to distance estimates of 0.5 and 1.3 kpc according to the NE2001 and YMW16 Galactic free electron density models, respectively. Some of its single pulses show an uninterrupted 180 degree sweep of the phase-resolved polarization position angle, with an S-shape reminiscent of the rotating vector model prediction. However, the fact that this sweep occurs at different phases from one pulse to another is remarkable and without straightforward explanation. Although PSR~J2251$-$3711 lies in the region of the $P-\dot{P}$ parameter space occupied by the X-ray Isolated Neutron Stars (XINS), there is no evidence for an X-ray counterpart in our Swift XRT observation; this places a 99\%-confidence upper bound on its unabsorbed bolometric thermal luminosity of $1.1 \times 10^{31}~(d / 1~\mathrm{kpc})^2~\mathrm{erg/s}$ for an assumed temperature of 85 eV, where $d$ is the distance to the pulsar. Further observations are needed to determine whether it is a rotation-powered pulsar with a true age of at least several Myr, or a much younger object such as an XINS or a recently cooled magnetar. Extreme specimens like PSR J2251$-$3711 help bridge populations in the so-called neutron star zoo in an attempt to understand their origins and evolution.

preprint2020arXiv

The SUrvey for Pulsars and Extragalactic Radio Bursts V: Recent Discoveries and Full Timing Solutions

The SUrvey for Pulsars and Extragalactic Radio Bursts ran from 2014 April to 2019 August, covering a large fraction of the southern hemisphere at mid- to high-galactic latitudes, and consisting of 9-minute pointings taken with the 20-cm multibeam receiver on the Parkes Radio Telescope. Data up to 2017 September 21 have been searched using standard Fourier techniques, single-pulse searches, and Fast Folding Algorithm searches. We present 19 new discoveries, bringing the total to 27 discoveries in the programme, and we report the results of follow-up timing observations at Parkes for 26 of these pulsars, including the millisecond pulsar PSR J1421-4409; the faint, highly-modulated, slow pulsar PSR J1646-1910; and the nulling pulsar PSR J1337-4441. We present new timing solutions for 23 pulsars, and we report flux densities, modulation indices, and polarization properties.

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

A pulsar-based timescale from the International Pulsar Timing Array

We have constructed a new timescale, TT(IPTA16), based on observations of radio pulsars presented in the first data release from the International Pulsar Timing Array (IPTA). We used two analysis techniques with independent estimates of the noise models for the pulsar observations and different algorithms for obtaining the pulsar timescale. The two analyses agree within the estimated uncertainties and both agree with TT(BIPM17), a post-corrected timescale produced by the Bureau International des Poids et Mesures (BIPM). We show that both methods could detect significant errors in TT(BIPM17) if they were present. We estimate the stability of the atomic clocks from which TT(BIPM17) is derived using observations of four rubidium fountain clocks at the US Naval Observatory. Comparing the power spectrum of TT(IPTA16) with that of these fountain clocks suggests that pulsar-based timescales are unlikely to contribute to the stability of the best timescales over the next decade, but they will remain a valuable independent check on atomic timescales. We also find that the stability of the pulsar-based timescale is likely to be limited by our knowledge of solar-system dynamics, and that errors in TT(BIPM17) will not be a limiting factor for the primary goal of the IPTA, which is to search for the signatures of nano-Hertz gravitational waves.

preprint2019arXiv

An ultra-wide bandwidth (704 to 4032 MHz) receiver for the Parkes radio telescope

We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver including its astronomical objectives, as well as the feed, receiver, digitiser and signal-processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration and timing stability.

preprint2016arXiv

21-year timing of the black-widow pulsar J2051-0827

Timing results for the black-widow pulsar J2051-0827 are presented, using a 21-year dataset from four European Pulsar Timing Array telescopes and the Parkes radio telescope. This dataset, which is the longest published to date for a black-widow system, allows for an improved analysis that addresses previously unknown biases. While secular variations, as identified in previous analyses, are recovered, short-term variations are detected for the first time. Concurrently, a significant decrease of approx. 2.5x10-3 cm-3 pc in the dispersion measure associated with PSR J2051-0827 is measured for the first time and improvements are also made to estimates of the proper motion. Finally, PSR J2051-0827 is shown to have entered a relatively stable state suggesting the possibility of its eventual inclusion in pulsar timing arrays.

preprint2016arXiv

A Fast Radio Burst Host Galaxy

In recent years, millisecond duration radio signals originating from distant galaxies appear to have been discovered in the so-called Fast Radio Bursts. These signals are dispersed according to a precise physical law and this dispersion is a key observable quantity which, in tandem with a redshift measurement, can be used for fundamental physical investigations. While every fast radio burst has a dispersion measurement, none before now have had a redshift measurement, due to the difficulty in pinpointing their celestial coordinates. Here we present the discovery of a fast radio burst and the identification of a fading radio transient lasting $\sim 6$ days after the event, which we use to identify the host galaxy; we measure the galaxy's redshift to be $z=0.492\pm0.008$. The dispersion measure and redshift, in combination, provide a direct measurement of the cosmic density of ionised baryons in the intergalactic medium of $Ω_{\mathrm{IGM}}=4.9 \pm 1.3\%$, in agreement with the expectation from WMAP, and including all of the so-called "missing baryons". The $\sim6$-day transient is largely consistent with a short gamma-ray burst radio afterglow, and its existence and timescale do not support progenitor models such as giant pulses from pulsars, and supernovae. This contrasts with the interpretation of another recently discovered fast radio burst, suggesting there are at least two classes of bursts.

preprint2016arXiv

A glitch in the millisecond pulsar J0613-0200

We present evidence for a small glitch in the spin evolution of the millisecond pulsar J0613$-$0200, using the EPTA Data Release 1.0, combined with Jodrell Bank analogue filterbank TOAs recorded with the Lovell telescope and Effelsberg Pulsar Observing System TOAs. A spin frequency step of 0.82(3) nHz and frequency derivative step of ${-1.6(39) \times 10^{-19}\,\text{Hz} \ \text{s}^{-1}}$ are measured at the epoch of MJD 50888(30). After PSR B1821$-$24A, this is only the second glitch ever observed in a millisecond pulsar, with a fractional size in frequency of ${Δν/ν=2.5(1) \times 10^{-12}}$, which is several times smaller than the previous smallest glitch. PSR J0613$-$0200 is used in gravitational wave searches with pulsar timing arrays, and is to date only the second such pulsar to have experienced a glitch in a combined 886 pulsar-years of observations. We find that accurately modelling the glitch does not impact the timing precision for pulsar timing array applications. We estimate that for the current set of millisecond pulsars included in the International Pulsar Timing Array, there is a probability of $\sim 50$% that another glitch will be observed in a timing array pulsar within 10 years.

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 LOFAR census of non-recycled pulsars: average profiles, dispersion measures, flux densities, and spectra

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.

preprint2016arXiv

A Massive Millisecond Pulsar in an Eccentric Binary

The recent discovery of a population of eccentric (e ~ 0.1) millisecond pulsar (MSP) binaries with low-mass white dwarf companions in the Galactic field represents a challenge to evolutionary models that explain MSP formation as recycling: all such models predict that the orbits become highly circularised during a long period of accretion. The members of this new population exhibit remarkably similar properties (orbital periods, eccentricities, companion masses, spin periods) and several models have been put forward that suggest a common formation channel. In this work we present the results of an extensive timing campaign focusing on one member of this new population, PSR J1946+3417. Through measurement of the both the advance of periastron and Shapiro delay for this system, we determine the mass of the pulsar, companion and the inclination of the orbit to be 1.828(22) Msun, 0.2656(19) Msun and 76.4(6) , under the assumption that general relativity is the true description of gravity. Notably, this is the third highest mass measured for any pulsar. Using these masses and the astrometric properties of PSR J1946+3417 we examine three proposed formation channels for eccentric MSP binaries. While our results are consistent with eccentricity growth driven by a circumbinary disk or neutron star to strange star phase transition, we rule out rotationally delayed accretion-induced collapse as the mechanism responsible for the configuration of the PSR J1946+3417 system.

preprint2016arXiv

A millisecond pulsar in an extremely wide binary system

We report on 22 yrs of radio timing observations of the millisecond pulsar J1024$-$0719 by the telescopes participating in the European Pulsar Timing Array (EPTA). These observations reveal a significant second derivative of the pulsar spin frequency and confirm the discrepancy between the parallax and Shklovskii distances that has been reported earlier. We also present optical astrometry, photometry and spectroscopy of 2MASS J10243869$-$0719190. We find that it is a low-metallicity main-sequence star (K7V spectral type, $\mathrm{[M/H]}=-1.0$, $T_\mathrm{eff}=4050\pm50$ K) and that its position, proper motion and distance are consistent with those of PSR J1024$-$0719. We conclude that PSR J1024$-$0719 and 2MASS J10243869$-$0719190 form a common proper motion pair and are gravitationally bound. The gravitational interaction between the main-sequence star and the pulsar accounts for the spin frequency derivatives, which in turn resolves the distance discrepancy. Our observations suggest that the pulsar and main-sequence star are in an extremely wide ($P_\mathrm{b}>200$ yr) orbit. Combining the radial velocity of the companion and proper motion of the pulsar, we find that the binary system has a high spatial velocity of $384\pm45$ km s$^{-1}$ with respect to the local standard of rest and has a Galactic orbit consistent with halo objects. Since the observed main-sequence companion star cannot have recycled the pulsar to millisecond spin periods, an exotic formation scenario is required. We demonstrate that this extremely wide-orbit binary could have evolved from a triple system that underwent an asymmetric supernova explosion, though find that significant fine-tuning during the explosion is required. Finally, we discuss the implications of the long period orbit on the timing stability of PSR J1024$-$0719 in light of its inclusion in pulsar timing arrays.

preprint2016arXiv

Detection of the magnetar SGR J1745-2900 up to 291 GHz with evidence of polarized millimetre emission

In Torne et al. (2015), we showed detections of SGR J1745-2900 up to 225 GHz (1.33 mm); at that time the highest radio frequency detection of pulsar emission. In this work, we present the results of new observations of the same magnetar with detections up to 291 GHz (1.03 mm), together with evidence of linear polarization in its millimetre emission. SGR J1745-2900 continues to show variability and is, on average, a factor $\sim$4 brighter in the millimetre band than in our observations of July 2014. The new measured spectrum is slightly inverted, with $\left<α\right> = +0.4\pm0.2$ (for $S_ν \propto ν^α)$. However, the spectrum does not seem to be well described by a single power law, which might be due to the intrinsic variability of the source, or perhaps a turn-up somewhere between 8.35 and 87 GHz. These results may help us to improve our still incomplete model of pulsar emission and, in addition, they further support the search for and study of pulsars located at the Galactic Centre using millimetre wavelengths.

preprint2016arXiv

FRBCAT: The Fast Radio Burst Catalogue

Here we present a catalogue of known Fast Radio Burst (FRB) sources in the form of an online catalogue, FRBCAT. The catalogue includes information about the instrumentation used for the observations for each detected burst, the measured quantities from each observation, and model-dependent quantities derived from observed quantities. To aid in consistent comparisons of burst properties such as width and signal-to-noise ratios we have reprocessed all the bursts for which we have access to the raw data, with software which we make available. The originally derived properties are also listed for comparison. The catalogue is hosted online as a MySQL database which can also be downloaded in tabular or plain text format for off-line use. This database will be maintained for use by the community for studies of the FRB population as it grows.

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

Improved Search for a Light Sterile Neutrino with the Full Configuration of the Daya Bay Experiment

This Letter reports an improved search for light sterile neutrino mixing in the electron antineutrino disappearance channel with the full configuration of the Daya Bay Reactor Neutrino Experiment. With an additional 404 days of data collected in eight antineutrino detectors, this search benefits from 3.6 times the statistics available to the previous publication, as well as from improvements in energy calibration and background reduction. A relative comparison of the rate and energy spectrum of reactor antineutrinos in the three experimental halls yields no evidence of sterile neutrino mixing in the $2\times10^{-4} \lesssim |Δm^{2}_{41}| \lesssim 0.3$ eV$^{2}$ mass range. The resulting limits on $\sin^{2}2θ_{14}$ are improved by approximately a factor of 2 over previous results and constitute the most stringent constraints to date in the $|Δm^{2}_{41}| \lesssim 0.2$ eV$^{2}$ region.

preprint2016arXiv

Limits on Active to Sterile Neutrino Oscillations from Disappearance Searches in the MINOS, Daya Bay, and Bugey-3 Experiments

Searches for a light sterile neutrino have been performed independently by the MINOS and the Daya Bay experiments using the muon (anti)neutrino and electron antineutrino disappearance channels, respectively. In this Letter, results from both experiments are combined with those from the Bugey-3 reactor neutrino experiment to constrain oscillations into light sterile neutrinos. The three experiments are sensitive to complementary regions of parameter space, enabling the combined analysis to probe regions allowed by the LSND and MiniBooNE experiments in a minimally extended four-neutrino flavor framework. Stringent limits on $\sin^2 2θ_{μe}$ are set over 6 orders of magnitude in the sterile mass-squared splitting $Δm^2_{41}$. The sterile-neutrino mixing phase space allowed by the LSND and MiniBooNE experiments is excluded for $Δm^2_{41} < 0.8$ eV$^2$ at 95% CL$_s$.

preprint2016arXiv

Long-term observations of the pulsars in 47 Tucanae. I. A study of four elusive binary systems

For the past couple of decades, the Parkes radio telescope has been regularly observing the millisecond pulsars in 47 Tucanae (47 Tuc). This long-term timing program was designed to address a wide range of scientific issues related to these pulsars and the globular cluster where they are located. In this paper, the first of a series, we address one of these objectives: the characterization of four previously known binary pulsars for which no precise orbital parameters were known, namely 47 Tuc P, V, W and X (pulsars 47 Tuc R and Y are discussed elsewhere). We determined the previously unknown orbital parameters of 47 Tuc V and X and greatly improved those of 47 Tuc P and W. For pulsars W and X we obtained, for the first time, full coherent timing solutions across the whole data span, which allowed a much more detailed characterization of these systems. 47 Tuc W, a well-known tight eclipsing binary pulsar, exhibits a large orbital period variability, as expected for a system of its class. 47 Tuc X turns out to be in a wide, extremely circular, 10.9-day long binary orbit and its position is ~3.8 arcmin away from the cluster center, more than three times the distance of any other pulsar in 47 Tuc. These characteristics make 47 Tuc X a very different object with respect to the other pulsars of the cluster.

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 measurement of $θ_{13}$ via neutron capture on hydrogen at Daya Bay

This article reports an improved independent measurement of neutrino mixing angle $θ_{13}$ at the Daya Bay Reactor Neutrino Experiment. Electron antineutrinos were identified by inverse $β$-decays with the emitted neutron captured by hydrogen, yielding a data-set with principally distinct uncertainties from that with neutrons captured by gadolinium. With the final two of eight antineutrino detectors installed, this study used 621 days of data including the previously reported 217-day data set with six detectors. The dominant statistical uncertainty was reduced by 49%. Intensive studies of the cosmogenic muon-induced $^9$Li and fast neutron backgrounds and the neutron-capture energy selection efficiency, resulted in a reduction of the systematic uncertainty by 26%. The deficit in the detected number of antineutrinos at the far detectors relative to the expected number based on the near detectors yielded $\sin^22θ_{13} = 0.071 \pm 0.011$ in the three-neutrino-oscillation framework. The combination of this result with the gadolinium-capture result is also reported.

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

Optical and radio astrometry of the galaxy associated with FRB150418

A fading radio source, coincident in time and position with the fast radio burst FRB150418, has been associated with the galaxy WISE J071634.59-190039.2. Subsequent observations of this galaxy have revealed that it contains a persistent, but variable, radio source. We present e-MERLIN, VLBA, and ATCA radio observations and Subaru optical observations of WISE J071634.59-190039.2 and find that the persistent radio source is unresolved and must be compact (<0.01 kpc), and that its location is consistent with the optical centre of the galaxy. We conclude that it is likely that WISE J071634.59-190039.2 contains a weak radio AGN.

preprint2016arXiv

Prospects for High-Precision Pulsar Timing with the New Effelsberg PSRIX Backend

The PSRIX backend is the primary pulsar timing instrument of the Effelsberg 100-m radio telescope since early 2011. This new ROACH-based system enables bandwidths up to 500 MHz to be recorded, significantly more than what was possible with its predecessor, the Effelsberg-Berkeley Pulsar Processor (EBPP). We review the first four years of PSRIX timing data for 33 pulsars collected as part of the monthly European Pulsar Timing Array (EPTA) observations. We describe the automated data analysis pipeline, CoastGuard, that we developed to reduce these observations. We also introduce TOASTER, the EPTA timing database used to store timing results, processing information and observation metadata. Using these new tools, we measure the phase-averaged flux densities at 1.4 GHz of all 33 pulsars. For 7 of these pulsars, our flux density measurements are the first values ever reported. For the other 26 pulsars, we compare our flux density measurements with previously published values. By comparing PSRIX data with EBPP data, we find an improvement of ~2-5 times in signal-to-noise ratio achievable, which translates to an increase of ~2-5 times in pulse time-of-arrival (TOA) precision. We show that such an improvement in TOA precision will improve the sensitivity to the stochastic gravitational wave background. Finally, we showcase the flexibility of the new PSRIX backend by observing several millisecond-period pulsars (MSPs) at 5 and 9 GHz. Motivated by our detections, we discuss the potential for complementing existing pulsar timing array data sets with MSP monitoring campaigns at these frequencies.

preprint2016arXiv

PSR B0329+54: Statistics of Substructure Discovered within the Scattering Disk on RadioAstron Baselines of up to 235,000 km

We discovered fine-scale structure within the scattering disk of PSR B0329+54 in observations with the RadioAstron ground-space radio interferometer. Here, we describe this phenomenon, characterize it with averages and correlation functions, and interpret it as the result of decorrelation of the impulse-response function of interstellar scattering between the widely-separated antennas. This instrument included the 10-m Space Radio Telescope, the 110-m Green Bank Telescope, the 14x25-m Westerbork Synthesis Radio Telescope, and the 64-m Kalyazin Radio Telescope. The observations were performed at 324 MHz, on baselines of up to 235,000 km in November 2012 and January 2014. In the delay domain, on long baselines the interferometric visibility consists of many discrete spikes within a limited range of delays. On short baselines it consists of a sharp spike surrounded by lower spikes. The average envelope of correlations of the visibility function show two exponential scales, with characteristic delays of $τ_1=4.1\pm 0.3\ μ{\rm s}$ and $τ_2=23\pm 3\ μ{\rm s}$, indicating the presence of two scales of scattering in the interstellar medium. These two scales are present in the pulse-broadening function. The longer scale contains 0.38 times the scattered power of the shorter one. We suggest that the longer tail arises from highly-scattered paths, possibly from anisotropic scattering or from substructure at large angles.

preprint2016arXiv

PSR B0329+54: Substructure in the scatter-broadened image discovered with RadioAstron on baselines up to 330,000 km

We have resolved the scatter-broadened image of PSR B0329+54 and detected substructure within it. These results are not influenced by any extended structure of a source but instead are directly attributed to the interstellar medium. We obtained these results at 324 MHz with the ground-space interferometer RadioAstron which included the space radio telescope (SRT), ground-based Westerbork Synthesis Radio Telescope and 64-m Kalyazin Radio Telescope on baseline projections up to 330,000 km in 2013 November 22 and 2014 January 1 to 2. At short 15,000 to 35,000 km ground-space baseline projections the visibility amplitude decreases with baseline length providing a direct measurement of the size of the scattering disk of 4.8$\pm$0.8 mas. At longer baselines no visibility detections from the scattering disk would be expected. However, significant detections were obtained with visibility amplitudes of 3 to 5% of the maximum scattered around a mean and approximately constant up to 330,000 km. These visibilities reflect substructure from scattering in the interstellar medium and offer a new probe of ionized interstellar material. The size of the diffraction spot near Earth is 17,000$\pm$3,000 km. With the assumption of turbulent irregularities in the plasma of the interstellar medium, we estimate that the effective scattering screen is located 0.6$\pm$0.1 of the distance from Earth toward the pulsar.

preprint2016arXiv

Pulsar observations with European telescopes for testing gravity and detecting gravitational waves

A background of nanohertz gravitational waves from supermassive black hole binaries could soon be detected by pulsar timing arrays, which measure the times-of-arrival of radio pulses from millisecond pulsars with very high precision. The European Pulsar Timing Array uses five large European radio telescopes to monitor high-precision millisecond pulsars, imposing in this way strong constraints on a gravitational wave background. To achieve the necessary precision needed to detect gravitational waves, the Large European Array for Pulsars (LEAP) performs simultaneous observations of pulsars with all five telescopes, which allows us to coherently add the radio pulses, maximize the signal-to-noise of pulsar signals and increase the precision of times-of-arrival. We report on the progress made and results obtained by the LEAP collaboration, and in particular on the addition of the Sardinia Radio Telescope to the LEAP observations during its scientific validation phase. In addition, we discuss how LEAP can be used to monitor strong-gravity systems such as double neutron star systems and impose strong constraints on post-keplerian parameters.

preprint2016arXiv

Radio light curve of the galaxy possibly associated with FRB150418

We present observations made with the Australia Telescope Compact Array (ATCA), the Jansky Very Large Array (JVLA) and the Giant Metre-Wave Telescope of the radio source within the galaxy WISE~J071634.59-190039.2, claimed to be host of FRB~150418 by Keane et al. (2016). We have established a common flux density scale between the ATCA and JVLA observations, the main result of which is to increase the flux densities obtained by Keane et al. At a frequency of 5.5 GHz, the source has a mean flux density of 140uJy and is variable on short timescales with a modulation index of 0.36. Statistical analysis of the flux densities shows that the variations seen are consistent with refractive interstellar scintillation of the weak active galactic nucleus at the centre of the galaxy. It may therefore be the case that the FRB and the galaxy are not associated. However, taking into account the rarity of highly variable sources in the radio sky, and our lack of knowledge of the progenitors of FRBs as a class, the association between WISE~J071634.59-190039.2 and FRB~150418 remains a possibility.

preprint2016arXiv

Radio polarimetry of Galactic centre pulsars

To study the strength and structure of the magnetic field in the Galactic centre (GC) we measured Faraday rotation of the radio emission of pulsars which are seen towards the GC. Three of these pulsars have the largest rotation measures (RMs) observed in any Galactic object with the exception of Sgr A*. Their large dispersion measures, RMs and the large RM variation between these pulsars and other known objects in the GC implies that the pulsars lie in the GC and are not merely seen in projection towards the GC. The large RMs of these pulsars indicate large line-of-sight magnetic field components between ~ 16-33 microgauss; combined with recent model predictions for the strength of the magnetic field in the GC this implies that the large-scale magnetic field has a very small inclination angle with respect to the plane of the sky (~ 12 degrees). Foreground objects like the Radio Arc or possibly an ablated, ionized halo around the molecular cloud G0.11-0.11 could contribute to the large RMs of two of the pulsars. If these pulsars lie behind the Radio Arc or G0.11-0.11 then this proves that low-scattering corridors with lengths >~ 100 pc must exist in the GC. This also suggests that future, sensitive observations will be able to detect additional pulsars in the GC. Finally, we show that the GC component in our most accurate electron density model oversimplifies structure in the GC.

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

Study of the wave packet treatment of neutrino oscillation at Daya Bay

The disappearance of reactor $\barν_e$ observed by the Daya Bay experiment is examined in the framework of a model in which the neutrino is described by a wave packet with a relative intrinsic momentum dispersion $σ_\text{rel}$. Three pairs of nuclear reactors and eight antineutrino detectors, each with good energy resolution, distributed among three experimental halls, supply a high-statistics sample of $\barν_e$ acquired at nine different baselines. This provides a unique platform to test the effects which arise from the wave packet treatment of neutrino oscillation. The modified survival probability formula was used to fit Daya Bay data, providing the first experimental limits: $2.38 \cdot 10^{-17} < σ_{\rm rel} < 0.23$. Treating the dimensions of the reactor cores and detectors as constraints, the limits are improved: $10^{-14} \lesssim σ_{\rm rel} < 0.23$, and an upper limit of $σ_{\rm rel} <0.20$ is obtained. All limits correspond to a 95\% C.L. Furthermore, the effect due to the wave packet nature of neutrino oscillation is found to be insignificant for reactor antineutrinos detected by the Daya Bay experiment thus ensuring an unbiased measurement of the oscillation parameters $\sin^22θ_{13}$ and $Δm^2_{32}$ within the plane wave model.

preprint2016arXiv

The Braking Index of a Radio-quiet Gamma-ray Pulsar

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

preprint2016arXiv

The Detector System of The Daya Bay Reactor Neutrino Experiment

The Daya Bay experiment was the first to report simultaneous measurements of reactor antineutrinos at multiple baselines leading to the discovery of $\barν_e$ oscillations over km-baselines. Subsequent data has provided the world's most precise measurement of $\rm{sin}^22θ_{13}$ and the effective mass splitting $Δm_{ee}^2$. The experiment is located in Daya Bay, China where the cluster of six nuclear reactors is among the world's most prolific sources of electron antineutrinos. Multiple antineutrino detectors are deployed in three underground water pools at different distances from the reactor cores to search for deviations in the antineutrino rate and energy spectrum due to neutrino mixing. Instrumented with photomultiplier tubes (PMTs), the water pools serve as shielding against natural radioactivity from the surrounding rock and provide efficient muon tagging. Arrays of resistive plate chambers over the top of each pool provide additional muon detection. The antineutrino detectors were specifically designed for measurements of the antineutrino flux with minimal systematic uncertainty. Relative detector efficiencies between the near and far detectors are known to better than 0.2%. With the unblinding of the final two detectors' baselines and target masses, a complete description and comparison of the eight antineutrino detectors can now be presented. This paper describes the Daya Bay detector systems, consisting of eight antineutrino detectors in three instrumented water pools in three underground halls, and their operation through the first year of eight detector data-taking.

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

Three-dimensional Tomography of the Galactic and Extragalactic Magnetoionic Medium with the SKA

The magneto-ionic structures of the interstellar medium of the Milky Way and the intergalactic medium are still poorly understood, especially at distances larger than a few kiloparsecs from the Sun. The three-dimensional (3D) structure of the Galactic magnetic field and electron density distribution may be probed through observations of radio pulsars, primarily owing to their compact nature, high velocities, and highly-polarized short-duration radio pulses. Phase 1 of the SKA, i.e. SKA1, will increase the known pulsar population by an order of magnitude, and the full SKA, i.e. SKA2, will discover pulsars in the most distant regions of our Galaxy. SKA1-VLBI will produce model-independent distances to a large number of pulsars, and wide-band polarization observations by SKA1-LOW and SKA1-MID will yield high precision dispersion measure, scattering measure, and rotation measure estimates along thousands of lines of sight. When combined, these observations will enable detailed tomography of the large-scale magneto-ionic structure of both the Galactic disk and the Galactic halo. Turbulence in the interstellar medium can be studied through the variations of these observables and the dynamic spectra of pulsar flux densities. SKA1-LOW and SKA1-MID will monitor interstellar weather and produce sensitive dynamic and secondary spectra of pulsar scintillation, which can be used to make speckle images of the ISM, study turbulence on scales between ~10^8 and ~10^13 m, and probe pulsar emission regions on scales down to $\sim$10 km. In addition, extragalactic pulsars or fast radio bursts to be discovered by SKA1 and SKA2 can be used to probe the electron density distribution and magnetic fields in the intergalactic medium beyond the Milky Way.

preprint2016arXiv

Variability, polarimetry, and timing properties of single pulses from PSR J1713+0747 using the Large European Array for Pulsars

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.

preprint2015arXiv

A Bayesian method for pulsar template generation

Extracting Times of Arrival from pulsar radio signals depends on the knowledge of the pulsars pulse profile and how this template is generated. We examine pulsar template generation with Bayesian methods. We will contrast the classical generation mechanism of averaging intensity profiles with a new approach based on Bayesian inference. We introduce the Bayesian measurement model imposed and derive the algorithm to reconstruct a "statistical template" out of noisy data. The properties of these "statistical templates" are analysed with simulated and real measurement data from PSR B1133+16. We explain how to put this new form of template to use in analysing secondary parameters of interest and give various examples: We implement a nonlinear filter for determining ToAs of pulsars. Applying this method to data from PSR J1713+0747 we derive ToAs self consistently, meaning all epochs were timed and we used the same epochs for template generation. While the average template contains fluctuations and noise as unavoidable artifacts, we find that the "statistical template" derived by Bayesian inference quantifies fluctuations and remaining uncertainty. This is why the algorithm suggested turns out to reconstruct templates of statistical significance from ten to fifty single pulses. A moving data window of fifty pulses, taking out one single pulse at the beginning and adding one at the end of the window unravels the characteristics of the methods to be compared. It shows that the change induced in the classical reconstruction is dominated by random fluctuations for the average template, while statistically significant changes drive the dynamics of the proposed method's reconstruction. The analysis of phase shifts with simulated data reveals that the proposed nonlinear algorithm is able to reconstruct correct phase information along with an acceptable estimation of the remaining uncertainty.

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

A new measurement of antineutrino oscillation with the full detector configuration at Daya Bay

We report a new measurement of electron antineutrino disappearance using the fully-constructed Daya Bay Reactor Neutrino Experiment. The final two of eight antineutrino detectors were installed in the summer of 2012. Including the 404 days of data collected from October 2012 to November 2013 resulted in a total exposure of 6.9$\times$10$^5$ GW$_{\rm th}$-ton-days, a 3.6 times increase over our previous results. Improvements in energy calibration limited variations between detectors to 0.2%. Removal of six $^{241}$Am-$^{13}$C radioactive calibration sources reduced the background by a factor of two for the detectors in the experimental hall furthest from the reactors. Direct prediction of the antineutrino signal in the far detectors based on the measurements in the near detectors explicitly minimized the dependence of the measurement on models of reactor antineutrino emission. The uncertainties in our estimates of $\sin^{2}2θ_{13}$ and $|Δm^2_{ee}|$ were halved as a result of these improvements. Analysis of the relative antineutrino rates and energy spectra between detectors gave $\sin^{2}2θ_{13} = 0.084\pm0.005$ and $|Δm^{2}_{ee}|= (2.42\pm0.11) \times 10^{-3}$ eV$^2$ in the three-neutrino framework.

preprint2015arXiv

A survey of FRB fields: Limits on repeatability

Several theories exist to explain the source of the bright, millisecond duration pulses known as fast radio bursts (FRBs). If the progenitors of FRBs are non-cataclysmic, such as giant pulses from pulsars, pulsar-planet binaries, or magnetar flares, FRB emission may be seen to repeat. We have undertaken a survey of the fields of eight known FRBs from the High Time Resolution Universe survey to search for repeating pulses. Although no repeat pulses were detected the survey yielded the detection of a new FRB, described in Petroff et al. (2015a). From our observations we rule out periodic repeating sources with periods P $\leq$ 8.6 hours and rule out sources with periods 8.6 < P < 21 hours at the 90% confidence level. At P $\geq$ 21 hours our limits fall off as ~1/P. Dedicated and persistent observations of FRB source fields are needed to rule out repetition on longer timescales, a task well-suited to next generation wide-field transient detectors.

preprint2015arXiv

Cool white dwarf companions to four millisecond pulsars

We report on photometric and spectroscopic observations of white dwarf companions to four binary radio millisecond pulsars, leading to the discovery of companions to PSRs J0614-3329, J1231-1411 and J2017+0603. We place limits on the brightness of the companion to PSR J0613-0200. Optical spectroscopy of the companion to PSR J0614-3329 identifies it as a DA type white dwarf with a temperature of Teff=6460+-80 K, a surface gravity log g=7.0+-0.2 cgs and a mass of Mwd=0.24+-0.04 Msun. We find that the distance to PSR J0614-3329 is smaller than previously estimated, removing the need for the pulsar to have an unrealistically high gamma-ray efficiency. Comparing the photometry with predictions from white dwarf cooling models allows us to estimate temperatures and cooling ages of the companions to PSRs J0613-0200, J1231-1411 and J2017+0603. We find that the white dwarfs in these systems are cool Teff<4000 K and old >5 Gyr. Thin Hydrogen envelopes are required for these white dwarfs to cool to the observed temperatures, and we suggest that besides Hydrogen shell flashes, irradiation driven mass loss by the pulsar may have been important.

preprint2015arXiv

Five new Fast Radio Bursts from the HTRU high latitude survey: first evidence for two-component bursts

The detection of five new fast radio bursts (FRBs) found in the High Time Resolution Universe high latitude survey is presented. The rate implied is 6$^{+4}_{-3}\times~10^3$ (95%) FRBs sky$^{-1}$ day$^{-1}$ above a fluence of between 0.13 and 5.9 Jy ms for FRBs between 0.128 and 262 ms in duration. One of these FRBs has a clear two-component profile, each component is similar to the known population of single component FRBs and are separated by 2.4(4) ms. All the FRB components appear to be unresolved following deconvolution with a scattering tail and accounting for intra-channel smearing. The two-component FRB also has the highest dispersion measure (1629 pc cm$^{-3}$) of any FRB to-date. Many of the proposed models to explain FRBs use a single high energy event involving compact objects (such as neutron star mergers) and therefore cannot easily explain a two-component FRB. Models that are based on extreme versions of flaring, pulsing or orbital events however could produce multiple component profiles. The compatibility of these models and the FRB rate implied by these detections is discussed.

preprint2015arXiv

LEAP: the large European array for pulsars

The Large European Array for Pulsars (LEAP) is an experiment that harvests the collective power of Europe's largest radio telescopes in order to increase the sensitivity of high-precision pulsar timing. As part of the ongoing effort of the European Pulsar Timing Array (EPTA), LEAP aims to go beyond the sensitivity threshold needed to deliver the first direct detection of gravitational waves. The five telescopes presently included in LEAP are: the Effelsberg telescope, the Lovell telescope at Jodrell Bank, the Nan\c cay radio telescope, the Sardinia Radio Telescope and the Westerbork Synthesis Radio Telescope. Dual polarization, Nyquist-sampled time-series of the incoming radio waves are recorded and processed offline to form the coherent sum, resulting in a tied-array telescope with an effective aperture equivalent to a 195-m diameter circular dish. All observations are performed using a bandwidth of 128 MHz centered at a frequency of 1396 MHz. In this paper, we present the design of the LEAP experiment, the instrumentation, the storage and transfer of data, and the processing hardware and software. In particular, we present the software pipeline that was designed to process the Nyquist-sampled time-series, measure the phase and time delays between each individual telescope and a reference telescope and apply these delays to form the tied-array coherent addition. The pipeline includes polarization calibration and interference mitigation. We also present the first results from LEAP and demonstrate the resulting increase in sensitivity, which leads to an improvement in the pulse arrival times.

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 Sparse Image Reconstruction

Context. The LOw Frequency ARray (LOFAR) radio telescope is a giant digital phased array interferometer with multiple antennas distributed in Europe. It provides discrete sets of Fourier components of the sky brightness. Recovering the original brightness distribution with aperture synthesis forms an inverse problem that can be solved by various deconvolution and minimization methods Aims. Recent papers have established a clear link between the discrete nature of radio interferometry measurement and the "compressed sensing" (CS) theory, which supports sparse reconstruction methods to form an image from the measured visibilities. Empowered by proximal theory, CS offers a sound framework for efficient global minimization and sparse data representation using fast algorithms. Combined with instrumental direction-dependent effects (DDE) in the scope of a real instrument, we developed and validated a new method based on this framework Methods. We implemented a sparse reconstruction method in the standard LOFAR imaging tool and compared the photometric and resolution performance of this new imager with that of CLEAN-based methods (CLEAN and MS-CLEAN) with simulated and real LOFAR data Results. We show that i) sparse reconstruction performs as well as CLEAN in recovering the flux of point sources; ii) performs much better on extended objects (the root mean square error is reduced by a factor of up to 10); and iii) provides a solution with an effective angular resolution 2-3 times better than the CLEAN images. Conclusions. Sparse recovery gives a correct photometry on high dynamic and wide-field images and improved realistic structures of extended sources (of simulated and real LOFAR datasets). This sparse reconstruction method is compatible with modern interferometric imagers that handle DDE corrections (A- and W-projections) required for current and future instruments such as LOFAR and SKA

preprint2015arXiv

Long Term Study of the Double Pulsar J0737-3039 with XMM-Newton: pulsar timing

The relativistic double neutron star binary PSR J0737-3039 shows clear evidence of orbital phase-dependent wind-companion interaction, both in radio and X-rays. In this paper we present the results of timing analysis of PSR J0737-3039 performed during 2006 and 2011 XMM-Newton Large Programs that collected ~20,000 X-ray counts from the system. We detected pulsations from PSR J0737-3039A (PSR A) through the most accurate timing measurement obtained by XMM-Newton so far, the spin period error being of 2x10^-13 s. PSR A's pulse profile in X-rays is very stable despite significant relativistic spin precession that occurred within the time span of observations. This yields a constraint on the misalignment between the spin axis and the orbital momentum axis Delta_A ~6.6^{+1.3}_{-5.4} deg, consistent with estimates based on radio data. We confirmed pulsed emission from PSR J0737-3039B (PSR B) in X-rays even after its disappearance in radio. The unusual phenomenology of PSR B's X-ray emission includes orbital pulsed flux and profile variations as well as a loss of pulsar phase coherence on time scales of years. We hypothesize that this is due to the interaction of PSR A's wind with PSR B's magnetosphere and orbital-dependent penetration of the wind plasma onto PSR B closed field lines. Finally, the analysis of the full XMM-Newton dataset provided evidences of orbital flux variability (~7%) for the first time, involving a bow-shock scenario between PSR A's wind and PSR B's magnetosphere.

preprint2015arXiv

Long-term Observations of Three Nulling Pulsars

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.

preprint2015arXiv

Low-Frequency Spectral Turn-Overs in Millisecond Pulsars Studied from Imaging Observations

Measurements of pulsar flux densities are of great importance for understanding the pulsar emission mechanism and for predictions of pulsar survey yields and the pulsar population at large. Typically these flux densities are determined from phase-averaged "pulse profiles", but this method has limited applicability at low frequencies because the observed pulses can easily be spread out by interstellar effects like scattering or dispersion, leading to a non-pulsed continuum component that is necessarily ignored in this type of analysis. In particular for the class of the millisecond pulsars (MSPs) at frequencies below 200MHz, such interstellar effects can seriously compromise de- tectability and measured flux densities. In this paper we investigate MSP spectra based on a complementary approach, namely through investigation of archival con- tinuum imaging data. Even though these images lose sensitivity to pulsars since the on-pulse emission is averaged with off-pulse noise, they are insensitive to effects from scattering and provide a reliable way to determine the flux density and spectral indices of MSPs based on both pulsed and unpulsed components. Using the 74MHz VLSSr as well as the 325MHz WENSS and 1.4GHz NVSS catalogues, we investigate the imaging flux densities of MSPs and evaluate the likelihood of spectral turn-overs in this population. We determine three new MSP spectral indices and identify six new MSPs with likely spectral turn-overs.

preprint2015arXiv

Measurement of the Reactor Antineutrino Flux and Spectrum at Daya Bay

This Letter reports a measurement of the flux and energy spectrum of electron antineutrinos from six 2.9~GW$_{th}$ nuclear reactors with six detectors deployed in two near (effective baselines 512~m and 561~m) and one far (1,579~m) underground experimental halls in the Daya Bay experiment. Using 217 days of data, 296,721 and 41,589 inverse beta decay (IBD) candidates were detected in the near and far halls, respectively. The measured IBD yield is (1.55 $\pm$ 0.04) $\times$ 10$^{-18}$~cm$^2$/GW/day or (5.92 $\pm$ 0.14) $\times$ 10$^{-43}$~cm$^2$/fission. This flux measurement is consistent with previous short-baseline reactor antineutrino experiments and is $0.946\pm0.022$ ($0.991\pm0.023$) relative to the flux predicted with the Huber+Mueller (ILL+Vogel) fissile antineutrino model. The measured IBD positron energy spectrum deviates from both spectral predictions by more than 2$σ$ over the full energy range with a local significance of up to $\sim$4$σ$ between 4-6 MeV. A reactor antineutrino spectrum of IBD reactions is extracted from the measured positron energy spectrum for model-independent predictions.

preprint2015arXiv

Observing Radio Pulsars in the Galactic Centre with the Square Kilometre Array

The discovery and timing of radio pulsars within the Galactic centre is a fundamental aspect of the SKA Science Case, responding to the topic of "Strong Field Tests of Gravity with Pulsars and Black Holes" (Kramer et al. 2004; Cordes et al. 2004). Pulsars have in many ways proven to be excellent tools for testing the General theory of Relativity and alternative gravity theories (see Wex (2014) for a recent review). Timing a pulsar in orbit around a companion, provides a unique way of probing the relativistic dynamics and spacetime of such a system. The strictest tests of gravity, in strong field conditions, are expected to come from a pulsar orbiting a black hole. In this sense, a pulsar in a close orbit ($P_{\rm orb}$ < 1 yr) around our nearest supermassive black hole candidate, Sagittarius A* - at a distance of ~8.3 kpc in the Galactic centre (Gillessen et al. 2009a) - would be the ideal tool. Given the size of the orbit and the relativistic effects associated with it, even a slowly spinning pulsar would allow the black hole spacetime to be explored in great detail (Liu et al. 2012). For example, measurement of the frame dragging caused by the rotation of the supermassive black hole, would allow a test of the "cosmic censorship conjecture." The "no-hair theorem" can be tested by measuring the quadrupole moment of the black hole. These are two of the prime examples for the fundamental studies of gravity one could do with a pulsar around Sagittarius A*. As will be shown here, SKA1-MID and ultimately the SKA will provide the opportunity to begin to find and time the pulsars in this extreme environment.

preprint2015arXiv

PSR J1906+0722: An Elusive Gamma-ray Pulsar

We report the discovery of PSR J1906+0722, a gamma-ray pulsar detected as part of a blind survey of unidentified Fermi Large Area Telescope (LAT) sources being carried out on the volunteer distributed computing system, Einstein@Home. This newly discovered pulsar previously appeared as the most significant remaining unidentified gamma-ray source without a known association in the second Fermi-LAT source catalog (2FGL) and was among the top ten most significant unassociated sources in the recent third catalog (3FGL). PSR J1906+0722 is a young, energetic, isolated pulsar, with a spin frequency of $8.9$ Hz, a characteristic age of $49$ kyr, and spin-down power $1.0 \times 10^{36}$ erg s$^{-1}$. In 2009 August it suffered one of the largest glitches detected from a gamma-ray pulsar ($Δf / f \approx 4.5\times10^{-6}$). Remaining undetected in dedicated radio follow-up observations, the pulsar is likely radio-quiet. An off-pulse analysis of the gamma-ray flux from the location of PSR J1906+0722 revealed the presence of an additional nearby source, which may be emission from the interaction between a neighboring supernova remnant and a molecular cloud. We discuss possible effects which may have hindered the detection of PSR J1906+0722 in previous searches and describe the methods by which these effects were mitigated in this survey. We also demonstrate the use of advanced timing methods for estimating the positional, spin and glitch parameters of difficult-to-time pulsars such as this.

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

Simultaneous multifrequency radio observations of the Galactic Centre magnetar SGR J1745-2900

We report on simultaneous observations of the magnetar SGR J1745-2900 at frequencies $ν= 2.54$ to $225\,\rm{GHz}$ using the Nancay 94-m equivalent, Effelsberg 100-m, and IRAM 30-m radio telescopes. We detect SGR J1745-2900 up to 225 GHz, the highest radio frequency detection of pulsed emission from a neutron star to date. Strong single pulses are also observed from 4.85 up to 154 GHz. At the millimetre band we see significant flux density and spectral index variabilities on time scales of tens of minutes, plus variability between days at all frequencies. Additionally, SGR J1745-2900 was observed at a different epoch at frequencies 296 to 472 GHz using the APEX 12-m radio telescope, with no detections. Over the period MJD 56859.83-56862.93 the fitted spectrum yields a spectral index of $\left<α\right> = -0.4 \pm 0.1$ for a reference flux density $\left< S_{154} \right> = 1.1 \pm 0.2\rm{\,mJy}$ (with $S_ν \propto ν^α)$, a flat spectrum alike those of the other radio-loud magnetars. These results show that strongly magnetized neutron stars can be effective radio emitters at frequencies notably higher to what was previously known and that pulsar searches in the Galactic Centre are possible in the millimetre band.

preprint2015arXiv

Single pulse and profile variability study of PSR J1022+1001

Millisecond pulsars (MSPs) are known as highly stable celestial clocks. Nevertheless, recent studies have revealed the unstable nature of their integrated pulse profiles, which may limit the achievable pulsar timing precision. In this paper, we present a case study on the pulse profile variability of PSR J1022+1001. We have detected approximately 14,000 sub-pulses (components of single pulses) in 35-hr long observations, mostly located at the trailing component of the integrated profile. Their flux densities and fractional polarisation suggest that they represent the bright end of the energy distribution in ordinary emission mode and are not giant pulses. The occurrence of sub-pulses from the leading and trailing components of the integrated profile is shown to be correlated. For sub-pulses from the latter, a preferred pulse width of approximately 0.25 ms has been found. Using simultaneous observations from the Effelsberg 100-m telescope and the Westerbork Synthesis Radio Telescope, we have found that the integrated profile varies on a timescale of a few tens of minutes. We show that improper polarisation calibration and diffractive scintillation cannot be the sole reason for the observed instability. In addition, we demonstrate that timing residuals generated from averages of the detected sub-pulses are dominated by phase jitter, and place an upper limit of ~700 ns for jitter noise based on continuous 1-min integrations.

preprint2015arXiv

The High Time Resolution Universe Pulsar Survey XII : Galactic plane acceleration search and the discovery of 60 pulsars

We present initial results from the low-latitude Galactic plane region of the High Time Resolution Universe pulsar survey conducted at the Parkes 64-m radio telescope. We discuss the computational challenges arising from the processing of the terabyte-sized survey data. Two new radio interference mitigation techniques are introduced, as well as a partially-coherent segmented acceleration search algorithm which aims to increase our chances of discovering highly-relativistic short-orbit binary systems, covering a parameter space including potential pulsar-black hole binaries. We show that under a constant acceleration approximation, a ratio of data length over orbital period of ~0.1 results in the highest effectiveness for this search algorithm. From the 50 per cent of data processed thus far, we have re-detected 435 previously known pulsars and discovered a further 60 pulsars, two of which are fast-spinning pulsars with periods less than 30ms. PSR J1101-6424 is a millisecond pulsar whose heavy white dwarf (WD) companion and short spin period of 5.1ms indicate a rare example of full-recycling via Case A Roche lobe overflow. PSR J1757-27 appears to be an isolated recycled pulsar with a relatively long spin period of 17ms. In addition, PSR J1244-6359 is a mildly-recycled binary system with a heavy WD companion, PSR J1755-25 has a significant orbital eccentricity of 0.09, and PSR J1759-24 is likely to be a long-orbit eclipsing binary with orbital period of the order of tens of years. Comparison of our newly-discovered pulsar sample to the known population suggests that they belong to an older population. Furthermore, we demonstrate that our current pulsar detection yield is as expected from population synthesis.

preprint2015arXiv

The Parkes multibeam pulsar survey: VII. Timing of four millisecond pulsars and the underlying spin period distribution of the Galactic millisecond pulsar population

We present timing observations of four millisecond pulsars discovered in the Parkes 20-cm multibeam pulsar survey of the Galactic plane. PSRs J1552-4937 and J1843-1448 are isolated objects with spin periods of 6.28 and 5.47 ms respectively. PSR J1727-2946 is in a 40-day binary orbit and has a spin period of 27 ms. The 4.43-ms pulsar J1813-2621 is in a circular 8.16-day binary orbit around a low-mass companion star with a minimum companion mass of 0.2 solar masses. Combining these results with detections from five other Parkes multibeam surveys, gives a well-defined sample of 56 pulsars with spin periods below 20 ms. We develop a likelihood analysis to constrain the functional form which best describes the underlying distribution of spin periods for millisecond pulsars. The best results were obtained with a log-normal distribution. A gamma distribution is less favoured, but still compatible with the observations. Uniform, power-law and Gaussian distributions are found to be inconsistent with the data. Galactic millisecond pulsars being found by current surveys appear to be in agreement with a log-normal distribution which allows for the existence of pulsars with periods below 1.5 ms.

preprint2015arXiv

The peculiar radio galaxy 4C 35.06: a case for recurrent AGN activity?

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in the core of the galaxy cluster Abell 407. At higher spatial resolution (~4"), the source was known to have two inner radio lobes spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution (~30"), the structure appears to have a helical morphology. We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it moves in the dense environment in the cluster core. Using LOFAR, we can trace the relic plasma from that episode of activity out to greater distances from the core than ever before. Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad (FWHM of 288 km/s), similar to what is found in other clusters. Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms at low frequencies using morphological and spectral information promises to bring new important insights in this field.

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 real-time fast radio burst: polarization detection and multiwavelength follow-up

Fast radio bursts (FRBs) are one of the most tantalizing mysteries of the radio sky; their progenitors and origins remain unknown and until now no rapid multiwavelength follow-up of an FRB has been possible. New instrumentation has decreased the time between observation and discovery from years to seconds, and enables polarimetry to be performed on FRBs for the first time. We have discovered an FRB (FRB 140514) in real-time on 14 May, 2014 at 17:14:11.06 UTC at the Parkes radio telescope and triggered follow-up at other wavelengths within hours of the event. FRB 140514 was found with a dispersion measure (DM) of 562.7(6) cm$^{-3}$ pc, giving an upper limit on source redshift of $z \lesssim 0.5$. FRB 140514 was found to be 21$\pm$7% (3-$σ$) circularly polarized on the leading edge with a 1-$σ$ upper limit on linear polarization $<10%$. We conclude that this polarization is intrinsic to the FRB. If there was any intrinsic linear polarization, as might be expected from coherent emission, then it may have been depolarized by Faraday rotation caused by passing through strong magnetic fields and/or high density environments. FRB 140514 was discovered during a campaign to re-observe known FRB fields, and lies close to a previous discovery, FRB 110220; based on the difference in DMs of these bursts and time-on-sky arguments, we attribute the proximity to sampling bias and conclude that they are distinct objects. Follow-up conducted by 12 telescopes observing from X-ray to radio wavelengths was unable to identify a variable multiwavelength counterpart, allowing us to rule out models in which FRBs originate from nearby ($z < 0.3$) supernovae and long duration gamma-ray bursts.

preprint2014arXiv

An absence of fast radio bursts at intermediate galactic latitudes

Fast radio bursts (FRBs) are an emerging class of bright, highly dispersed radio pulses. Recent work by Thornton et al. (2013) has revealed a population of FRBs in the High Time Resolution Universe (HTRU) survey at high Galactic latitudes. A variety of progenitors have been proposed including cataclysmic events at cosmological distances, Galactic flare stars, and terrestrial radio frequency interference. Here we report on a search for FRBs at intermediate Galactic latitudes ($-15^{\circ}$ $< b <$ 15$^{\circ}$) in data taken as part of the HTRU survey. No FRBs were discovered in this region. Several effects such as dispersion, scattering, sky temperature and scintillation decrease the sensitivity by more than 3$σ$ in $\sim$20\% of survey pointings. Including all of these effects, we exclude the hypothesis that FRBs are uniformly distributed on the sky with 99\% confidence. This low probability implies that additional factors -- not accounted for by standard Galactic models -- must be included to ease the discrepancy between the detection rates at high and low Galactic latitudes. A revised rate estimate or another strong and heretofore unknown selection effect in Galactic latitude would provide closer agreement between the surveys' detection rates. The dearth of detections at low Galactic latitude disfavors a Galactic origin for these bursts.

preprint2014arXiv

Constraints on the Emission Geometries and Spin Evolution of Gamma-ray Millisecond Pulsars

Millisecond pulsars (MSPs) are a growing class of gamma-ray emitters. Pulsed gamma-ray signals have been detected from more than 40 MSPs with the Fermi Large Area Telescope (LAT). The wider radio beams and more compact magnetospheres of MSPs enable studies of emission geometries over a broader range of phase space than non-recycled radio-loud gamma-ray pulsars. We have modeled the gamma-ray light curves of 40 LAT-detected MSPs using geometric emission models assuming a vacuum retarded-dipole magnetic field. We modeled the radio profiles using a single-altitude hollow-cone beam, with a core component when indicated by polarimetry; however, for MSPs with gamma-ray and radio light curve peaks occurring at nearly the same rotational phase we assume that the radio emission is co-located with the gamma rays and caustic in nature. The best-fit parameters and confidence intervals are determined using a maximum likelihood technique. We divide the light curves into three model classes, with gamma-ray peaks trailing (Class I), aligned (Class II) or leading (Class III) the radio peaks. Outer gap and slot gap (two-pole caustic) models best fit roughly equal numbers of Class I and II, while Class III are exclusively fit with pair-starved polar cap models. Distinguishing between the model classes based on typical derived parameters is difficult. We explore the evolution of magnetic inclination angle with period and spin-down power, finding possible correlations. While the presence of significant off-peak emission can often be used as a discriminator between outer gap and slot gap models, a hybrid model may be needed.

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

Future mmVLBI Research with ALMA: A European vision

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

preprint2014arXiv

Gravitational wave astronomy with the SKA

On a time scale of years to decades, gravitational wave (GW) astronomy will become a reality. Low frequency (nanoHz) GWs are detectable through long-term timing observations of the most stable pulsars. Radio observatories worldwide are currently carrying out observing programmes to detect GWs, with data sets being shared through the International Pulsar Timing Array project. One of the most likely sources of low frequency GWs are supermassive black hole binaries (SMBHBs), detectable as a background due to a large number of binaries, or as continuous or burst emission from individual sources. No GW signal has yet been detected, but stringent constraints are already being placed on galaxy evolution models. The SKA will bring this research to fruition. In this chapter, we describe how timing observations using SKA1 will contribute to detecting GWs, or can confirm a detection if a first signal already has been identified when SKA1 commences observations. We describe how SKA observations will identify the source(s) of a GW signal, search for anisotropies in the background, improve models of galaxy evolution, test theories of gravity, and characterise the early inspiral phase of a SMBHB system. We describe the impact of the large number of millisecond pulsars to be discovered by the SKA; and the observing cadence, observation durations, and instrumentation required to reach the necessary sensitivity. We describe the noise processes that will influence the achievable precision with the SKA. We assume a long-term timing programme using the SKA1-MID array and consider the implications of modifications to the current design. We describe the possible benefits from observations using SKA1-LOW. Finally, we describe GW detection prospects with SKA1 and SKA2, and end with a description of the expectations of GW astronomy.

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

Independent Measurement of Theta13 via Neutron Capture on Hydrogen at Daya Bay

A new measurement of the $θ_{13}$ mixing angle has been obtained at the Daya Bay Reactor Neutrino Experiment via the detection of inverse beta decays tagged by neutron capture on hydrogen. The antineutrino events for hydrogen capture are distinct from those for gadolinium capture with largely different systematic uncertainties, allowing a determination independent of the gadolinium-capture result and an improvement on the precision of $θ_{13}$ measurement. With a 217-day antineutrino data set obtained with six antineutrino detectors and from six 2.9 GW$_{th}$ reactors, the rate deficit observed at the far hall is interpreted as $\sin^22θ_{13}=0.083\pm0.018$ in the three-flavor oscillation model. When combined with the gadolinium-capture result from Daya Bay, we obtain $\sin^22θ_{13}=0.089\pm0.008$ as the final result for the six-antineutrino-detector configuration of the Daya Bay experiment.

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

Interstellar Scintillation of the Double Pulsar J0737$-$3039

We report here a series of observations of the interstellar scintillation (ISS) of the double pulsar J0737$-$3039 over the course of 18 months. As in earlier work (Coles et al., 2005) the basic phenomenon is the variation in the ISS caused by the changing transverse velocities of each pulsar, the ionized interstellar medium (IISM), and the Earth. The transverse velocity of the binary system can be determined both by VLBI and timing observations. The orbital velocity and inclination is almost completely determined from timing observations, but the direction of the orbital angular momentum is not known. Since the Earth's velocity is known, and can be compared with the orbital velocity by its effect on the timescale of the ISS, we can determine the orientation $Ω$ of the pulsar orbit with respect to equatorial coordinates ($Ω= 65\pm2$ deg). We also resolve the ambiguity ($i= 88.7$ or $91.3$ deg) in the inclination of the orbit deduced from the measured Shapiro delay by our estimate $i=88.1\pm0.5$ deg. This relies on analysis of the ISS over both frequency and time and provides a model for the location, anisotropy, turbulence level and transverse phase gradient of the IISM. We find that the IISM can be well-modeled during each observation, typically of a few orbital periods, but its turbulence level and mean velocity vary significantly over the 18 months.

preprint2014arXiv

LOFAR observations of PSR B0943+10: profile evolution and discovery of a systematically changing profile delay in Bright mode

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.

preprint2014arXiv

LOFAR Observations of Swift J1644+57 and Implications for Short-Duration Transients

Aims: We have analyzed low frequency radio data of tidal disruption event (TDE) Swift J1644+57 to search for a counterpart. We consider how brief transient signals (on the order of seconds or minutes) originating from this location would appear in our data. We also consider how automatic radio frequency interference (RFI) flagging at radio telescope observatories might affect these and other transient observations in the future, particularly with brief transients of a few seconds duration. Methods: We observed the field in the low-frequency regime at 149 MHz with data obtained over several months with the Low Frequency Array (LOFAR). We also present simulations where a brief transient is injected into the data in order to see how it would appear in our measurement sets, and how it would be affected by RFI flagging. Finally, both based on simulation work and the weighted average of the observed background over the course of the individual observations, we present the possibility of brief radio transients in the data. Results: Our observations of Swift J1644+57 yielded no detection of the source and a peak flux density at this position of 24.7 $\pm$ 8.9 mJy. Our upper limit on the transient rate of the snapshot surface density in this field at sensitivities < 0.5 Jy is $ρ< 2.2 \times10^{-2}$ deg$^{-2}$. We also conclude that we did not observe any brief transient signals originating specifically from the Swift J1644+57 source itself, and searches for such transients are severely limited by automatic RFI flagging algorithms which flag transients of less than 2 minutes duration. As such, careful consideration of RFI flagging techniques must occur when searching for transient signals.

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

Measuring a Cherenkov ring in the radio emission from air showers at 110-190 MHz with LOFAR

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle.

preprint2014arXiv

Measuring pulse times of arrival from broadband pulsar observations

In recent years, instrumentation enabling pulsar observations with unprecedentedly high fractional bandwidth has been under development which can be used to substantially improve the precision of pulsar timing experiments. The traditional template-matching method used to calculate pulse times-of-arrival (ToAs), may not function effectively on these broadband data due to a variety of effects such as diffractive scintillation in the interstellar medium, profile variation as a function of frequency, dispersion measure (DM) evolution and so forth. In this paper, we describe the channelised Discrete Fourier Transform method that can greatly mitigate the influence of the aforementioned effects when measuring ToAs from broadband timing data. The method is tested on simulated data, and its potential in improving timing precision is shown. We further apply the method to PSR J1909$-$3744 data collected at the Nançay Radio Telescope with the Nançay Ultimate Pulsar Processing Instrument. We demonstrate a removal of systematics due to the scintillation effect as well as improvement on ToA measurement uncertainties. Our method also determines temporal variations in dispersion measure, which are consistent with multi-channel timing approaches used earlier.

preprint2014arXiv

Model-based asymptotically optimal dispersion measure correction for pulsar timing

In order to reach the sensitivity required to detect gravitational waves, pulsar timing array experiments need to mitigate as much noise as possible in timing data. A dominant amount of noise is likely due to variations in the dispersion measure. To correct for such variations, we develop a statistical method inspired by the maximum likelihood estimator and optimal filtering. Our method consists of two major steps. First, the spectral index and amplitude of dispersion measure variations are measured via a time-domain spectral analysis. Second, the linear optimal filter is constructed based on the model parameters found in the first step, and is used to extract the dispersion measure variation waveforms. Compared to current existing methods, this method has better time resolution for the study of short timescale dispersion variations, and generally produces smaller errors in waveform estimations. This method can process irregularly sampled data without any interpolation because of its time-domain nature. Furthermore, it offers the possibility to interpolate or extrapolate the waveform estimation to regions where no data is available. Examples using simulated data sets are included for demonstration.

preprint2014arXiv

On the Apparent Nulls and Extreme Variability of PSR J1107-5907

We present an analysis of the emission behaviour of PSR J1107-5907, a source known to exhibit separate modes of emission, using observations obtained over approximately 10 yr. We find that the object exhibits two distinct modes of emission; a strong mode with a broad profile and a weak mode with a narrow profile. During the strong mode of emission, the pulsar typically radiates very energetic emission over sequences of ~200-6000 pulses (~60 s-24 min), with apparent nulls over time-scales of up to a few pulses at a time. Emission during the weak mode is observed outside of these strong-mode sequences and manifests as occasional bursts of up to a few clearly detectable pulses at a time, as well as low-level underlying emission which is only detected through profile integration. This implies that the previously described null mode may in fact be representative of the bottom-end of the pulse intensity distribution for the source. This is supported by the dramatic pulse-to-pulse intensity modulation and rarity of exceptionally bright pulses observed during both modes of emission. Coupled with the fact that the source could be interpreted as a rotating radio transient (RRAT)-like object for the vast majority of the time, if placed at a further distance, we advance that this object likely represents a bridge between RRATs and extreme moding pulsars. Further to these emission properties, we also show that the source is consistent with being a near-aligned rotator and that it does not exhibit any measurable spin-down rate variation. These results suggest that nulls observed in other intermittent objects may in fact be representative of very weak emission without the need for complete cessation. As such, we argue that longer (> 1 h) observations of pulsars are required to discern their true modulation properties.

preprint2014arXiv

Pulsar-black hole binaries: prospects for new gravity tests with future radio telescopes

The anticipated discovery of a pulsar in orbit with a black hole is expected to provide a unique laboratory for black hole physics and gravity. In this context, the next generation of radio telescopes, like the Five-hundred-metre Aperture Spherical radio Telescope (FAST) and the Square Kilometre Array (SKA), with their unprecedented sensitivity, will play a key role. In this paper, we investigate the capability of future radio telescopes to probe the spacetime of a black hole and test gravity theories, by timing a pulsar orbiting a stellar-mass-black-hole (SBH). Based on mock data simulations, we show that a few years of timing observations of a sufficiently compact pulsar-SBH (PSR-SBH) system with future radio telescopes would allow precise measurements of the black hole mass and spin. A measurement precision of one per cent can be expected for the spin. Measuring the quadrupole moment of the black hole, needed to test GR's no-hair theorem, requires extreme system configurations with compact orbits and a large SBH mass. Additionally, we show that a PSR-SBH system can lead to greatly improved constraints on alternative gravity theories even if they predict black holes (practically) identical to GR's. This is demonstrated for a specific class of scalar-tensor theories. Finally, we investigate the requirements for searching for PSR-SBH systems. It is shown that the high sensitivity of the next generation of radio telescopes is key for discovering compact PSR-SBH systems, as it will allow for sufficiently short survey integration times.

preprint2014arXiv

RadioAstron Studies of the Nearby, Turbulent Interstellar Plasma With the Longest Space-Ground Interferometer Baseline

RadioAstron space-ground VLBI observations of the pulsar B0950+08, conducted with the 10-m space radio telescope in conjunction with the Arecibo 300-m telescope and Westerbork Synthesis Radio Telescope at a frequency of 324 MHz, were analyzed in order to investigate plasma inhomogeneities in the direction of this nearby pulsar. The observations were conducted at a spacecraft distance of 330,000 km, resulting in a projected baseline of 220,000 km, providing the greatest angular resolution ever achieved at meter wavelengths. Our analysis is based on fundamental behavior of structure and coherence functions. We find that the pulsar shows scintillation on two frequency scales, both much less than the observing frequency; but modulation is less than 100%. We infer that the scattering is weak, but a refracting wedge disperses the scintillation pattern. The refraction angle of this "cosmic prism" is measured as theta_0=1.1 - 4.4 mas, with the refraction direction being approximately perpendicular to the observer velocity. We show that the observed parameters of scintillation effects indicate that two plasma layers lie along the line of sight to the pulsar, at distances of 4.4 - 16.4 pc and 26 - 170 pc, and traveling in different directions relative to the line of sight. Spectra of turbulence for the two layers are found to follow a power law with the indices gamma_1 = gamma_2 = 3.00 +/- 0.08, significantly different from the index expected for a Kolmogorov spectrum of turbulence, gamma=11/3.

preprint2014arXiv

Realistic Modeling of the Pulse Profile of PSR J0737-3039A

The Double Pulsar, PSR J$0737$$-$$3039$A/B, is a unique system in which both neutron stars have been detected as radio pulsars. As shown in Ferdman et al., there is no evidence for pulse profile evolution of the A pulsar, and the geometry of the pulsar was fit well with a double-pole circular radio beam model. Assuming a more realistic polar cap model with a vacuum retarded dipole magnetosphere configuration including special relativistic effects, we create synthesized pulse profiles for A given the best-fit geometry from the simple circular beam model. By fitting synthesized pulse profiles to those observed from pulsar A, we constrain the geometry of the radio beam, namely the half-opening angle and the emission altitude, to be $30^\circ$ and $10$ neutron star radii, respectively. Combining the observational constraints of PSR J$0737$$-$$3039$A/B, we are able to construct the full three-dimensional orbital geometry of the Double Pulsar. The relative angle between the spin axes of the two pulsars ($Δ_S$) is estimated to be ($138^\circ \pm 5^\circ$) at the current epoch and will likely remain constant until tidal interactions become important in $\sim$$85$ Myr, at merger.

preprint2014arXiv

Search for a Light Sterile Neutrino at Daya Bay

A search for light sterile neutrino mixing was performed with the first 217 days of data from the Daya Bay Reactor Antineutrino Experiment. The experiment's unique configuration of multiple baselines from six 2.9~GW$_{\rm th}$ nuclear reactors to six antineutrino detectors deployed in two near (effective baselines 512~m and 561~m) and one far (1579~m) underground experimental halls makes it possible to test for oscillations to a fourth (sterile) neutrino in the $10^{\rm -3}~{\rm eV}^{2} < |Δm_{41}^{2}| < 0.3~{\rm eV}^{2}$ range. The relative spectral distortion due to electron antineutrino disappearance was found to be consistent with that of the three-flavor oscillation model. The derived limits on $\sin^22θ_{14}$ cover the $10^{-3}~{\rm eV}^{2} \lesssim |Δm^{2}_{41}| \lesssim 0.1~{\rm eV}^{2}$ region, which was largely unexplored.

preprint2014arXiv

Six Faint Gamma-ray Pulsars Seen with the Fermi Large Area Telescope -- Towards a Sample Blending into the Background

Context: GeV gamma-ray pulsations from over 140 pulsars have been characterized using the Fermi Large Area Telescope, enabling improved understanding of the emission regions within the neutron star magnetospheres, and the contributions of pulsars to high energy electrons and diffuse gamma rays in the Milky Way. The first gamma-ray pulsars to be detected were the most intense and/or those with narrow pulses. Aims: As the Fermi mission progresses, progressively fainter objects can be studied. In addition to more distant pulsars (thus probing a larger volume of the Galaxy), or ones in high background regions (thus improving the sampling uniformity across the Galactic plane), we detect pulsars with broader pulses or lower luminosity. Adding pulsars to our catalog with inclination angles that are rare in the observed sample, and/or with lower spindown power, will reduce the bias in the currently known gamma-ray pulsar population. Methods: We use rotation ephemerides derived from radio observations to phase-fold gamma rays recorded by the Fermi Large Area Telescope, to then determine the pulse profile properties. Spectral analysis provides the luminosities and, when the signal-to-noise ratio allows, the cutoff energies. We constrain the pulsar distances by different means in order to minimize the luminosity uncertainties. Results: We present six new gamma-ray pulsars with an eclectic mix of properties. Three are young, and three are recycled. They include the farthest, the lowest power, two of the highest duty-cycle pulsars seen, and only the fourth young gamma-ray pulsar with a radio interpulse. We discuss the biases existing in the current gamma-ray pulsar catalog, and steps to be taken to mitigate the bias.

preprint2014arXiv

Spectral measurement of electron antineutrino oscillation amplitude and frequency at Daya Bay

A measurement of the energy dependence of antineutrino disappearance at the Daya Bay Reactor Neutrino Experiment is reported. Electron antineutrinos ($\overlineν_{e}$) from six $2.9$ GW$_{\rm th}$ reactors were detected with six detectors deployed in two near (effective baselines 512 m and 561 m) and one far (1579 m) underground experimental halls. Using 217 days of data, 41589 (203809 and 92912) antineutrino candidates were detected in the far hall (near halls). An improved measurement of the oscillation amplitude $\sin^{2}2θ_{13} = 0.090^{+0.008}_{-0.009} $ and the first direct measurement of the $\overlineν_{e}$ mass-squared difference $|Δm^{2}_{ee}|= (2.59_{-0.20}^{+0.19}) \times 10^{-3}\ {\rm eV}^2 $ is obtained using the observed $\overlineν_{e}$ rates and energy spectra in a three-neutrino framework. This value of $|Δm^{2}_{ee}|$ is consistent with $|Δm^{2}_{μμ}|$ measured by muon neutrino disappearance, supporting the three-flavor oscillation model.

preprint2014arXiv

The Binary Companion of Young, Relativistic Pulsar J1906+0746

PSR J1906+0746 is a young pulsar in the relativistic binary with the second-shortest known orbital period, of 3.98 hours. We here present a timing study based on five years of observations, conducted with the 5 largest radio telescopes in the world, aimed at determining the companion nature. Through the measurement of three post-Keplerian orbital parameters we find the pulsar mass to be 1.291(11) M_sol, and the companion mass 1.322(11) M_sol respectively. These masses fit well in the observed collection of double neutron stars, but are also compatible with other white dwarfs around young pulsars such as J1906+0746. Neither radio pulsations nor dispersion-inducing outflows that could have further established the companion nature were detected. We derive an HI-absorption distance, which indicates that an optical confirmation of a white dwarf companion is very challenging. The pulsar is fading fast due to geodetic precession, limiting future timing improvements. We conclude that young pulsar J1906+0746 is likely part of a double neutron star, or is otherwise orbited by an older white dwarf, in an exotic system formed through two stages of mass transfer.

preprint2014arXiv

The High Time Resolution Universe Pulsar Survey X: Discovery of four millisecond pulsars and updated timing solutions of a further 12

We report on the discovery of four millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) pulsar survey being conducted at the Parkes 64-m radio telescope. All four MSPs are in binary systems and are likely to have white dwarf companions. In addition, we present updated timing solutions for 12 previously published HTRU MSPs, revealing new observational parameters such as five proper motion measurements and significant temporal dispersion measure variations in PSR J1017-7156. We discuss the case of PSR J1801-3210, which shows no significant period derivative after four years of timing data. Our best-fit solution shows a period derivative of the order of $10^{-23}$, an extremely small number compared to that of a typical MSP. However, it is likely that the pulsar lies beyond the Galactic Centre, and an unremarkable intrinsic period derivative is reduced to close to zero by the Galactic potential acceleration. Furthermore, we highlight the potential to employ PSR J1801-3210 in the strong equivalence principle test due to its wide and circular orbit. In a broader comparison with the known MSP population, we suggest a correlation between higher mass functions and the presence of eclipses in `very low-mass binary pulsars', implying that eclipses are observed in systems with high orbital inclinations. We also suggest that the distribution of the total mass of binary systems is inversely-related to the Galactic height distribution. Finally, we report on the first detection of PSRs J1543-5149 and J1811-2404 as gamma-ray pulsars.

preprint2014arXiv

The High Time Resolution Universe Survey - XI. Discovery of five recycled pulsars and the optical detectability of survey white dwarf companions

We present the discovery of a further five recycled pulsar systems in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods ranging from 2 ms to 66 ms, and four are in binary systems with orbital periods between 10.8 hours and 9.0 days. Three of these binary systems are particularly interesting; PSR J1227-6208 has a pulse period of 34.5 ms and the highest mass function of all pulsars with near-circular orbits. The circular orbit suggests that the companion is not another neutron star, so future timing experiments may reveal one of the heaviest white dwarfs ever found ($>$ 1.3 M$_\odot$). Timing observations of PSR J1431$-$4715 indicate that it is eclipsed by its companion which has a mass indicating it belongs to the redback class of eclipsing millisecond pulsars. PSR J1653-2054 has a companion with a minimum mass of only $0.08$ M$_\odot$, placing it among the class of pulsars with low-mass companions. Unlike the majority of such systems, however, no evidence of eclipses is seen at 1.4 GHz.

preprint2014arXiv

The LOFAR long baseline snapshot calibrator survey

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

preprint2014arXiv

The LOFAR Pilot Surveys for Pulsars and Fast Radio Transients

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.

preprint2013arXiv

A coherent acceleration search of the Parkes multi-beam pulsar survey - techniques and the discovery and timing of 16 pulsars

A fully coherent acceleration search algorithm has been applied to the Parkes multi- beam pulsar survey of the Galactic plane to search for previously undiscovered relativistic binary pulsars. The search has resulted in the discovery of 16 pulsars including a binary millisecond pulsar and an intermittent pulsar. Despite a number of promising candidates there have been no new discoveries of relativistic binary pulsars. Here we detail the acceleration search performed in our analysis and present coherent timing solutions for each of pulsars discovered. In light of the lack of discoveries of relativistic binary pulsars, we also discuss the technique of acceleration searching and its effectiveness in finding these systems.

preprint2013arXiv

A Population of Fast Radio Bursts at Cosmological Distances

Searches for transient astrophysical sources often reveal unexpected classes of objects that are useful physical laboratories. In a recent survey for pulsars and fast transients we have uncovered four millisecond-duration radio transients all more than 40° from the Galactic plane. The bursts' properties indicate that they are of celestial rather than terrestrial origin. Host galaxy and intergalactic medium models suggest that they have cosmological redshifts of 0.5 to 1, and distances of up to 3 gigaparsecs. No temporally coincident x- or gamma-ray signature was identified in association with the bursts. Characterization of the source population and identification of host galaxies offers an opportunity to determine the baryonic content of the Universe.

preprint2013arXiv

A strong magnetic field around the supermassive black hole at the centre of the Galaxy

The centre of our Milky Way harbours the closest candidate for a supermassive black hole. The source is thought to be powered by radiatively inefficient accretion of gas from its environment. This form of accretion is a standard mode of energy supply for most galactic nuclei. X-ray measurements have already resolved a tenuous hot gas component from which it can be fed. However, the magnetization of the gas, a crucial parameter determining the structure of the accretion flow, remains unknown. Strong magnetic fields can influence the dynamics of the accretion, remove angular momentum from the infalling gas, expel matter through relativistic jets and lead to the observed synchrotron emission. Here we report multi-frequency measurements with several radio telescopes of a newly discovered pulsar close to the Galactic Centre and show that its unusually large Faraday rotation indicates a dynamically relevant magnetic field near the black hole. If this field is accreted down to the event horizon it provides enough magnetic flux to explain the observed emission from the black hole, from radio to X-rays.

preprint2013arXiv

An Asteroid Belt Interpretation for the Timing Variations of the Millisecond Pulsar B1937+21

Pulsar timing observations have revealed companions to neutron stars that include other neutron stars, white dwarfs, main-sequence stars, and planets. We demonstrate that the correlated and apparently stochastic residual times of arrival from the millisecond pulsar B1937+21 are consistent with the signature of an asteroid belt having a total mass less than approximately 0.05 Earth masses. Unlike the solar system's asteroid belt, the best fit pulsar asteroid belt extends over a wide range of radii, consistent with the absence of any shepherding companions. We suggest that any pulsar that has undergone accretion-driven spin-up and subsequently evaporated its companion may harbor orbiting asteroid mass objects. The resulting timing variations may fundamentally limit the timing precision of some of the other millisecond pulsars. Observational tests of the asteroid belt model include identifying periodicities from individual asteroids, which are difficult; testing for statistical stationarity that become possible when observations are conducted over a longer observing span; and searching for reflected radio emission.

preprint2013arXiv

Assembly and Installation of the Daya Bay Antineutrino Detectors

The Daya Bay reactor antineutrino experiment is designed to make a precision measurement of the neutrino mixing angle theta13, and recently made the definitive discovery of its nonzero value. It utilizes a set of eight, functionally identical antineutrino detectors to measure the reactor flux and spectrum at baselines of 300 - 2000m from the Daya Bay and Ling Ao Nuclear Power Plants. The Daya Bay antineutrino detectors were built in an above-ground facility and deployed side-by-side at three underground experimental sites near and far from the nuclear reactors. This configuration allows the experiment to make a precision measurement of reactor antineutrino disappearance over km-long baselines and reduces relative systematic uncertainties between detectors and nuclear reactors. This paper describes the assembly and installation of the Daya Bay antineutrino detectors.

preprint2013arXiv

Broadband Pulsations from PSR B1821-24: Implications for Emission Models and the Pulsar Population of M28

We report a 5.4σ detection of pulsed gamma rays from PSR B1821-24 in the globular cluster M28 using ~44 months of Fermi Large Area Telescope (LAT) data that have been reprocessed with improved instrument calibration constants. We constructed a phase-coherent ephemeris, with post-fit residual RMS of 3 μs, using radio data spanning ~23.2 years, enabling measurements of the multi-wavelength light curve properties of PSR B1821-24 at the milliperiod level. We fold RXTE observations of PSR B1821-24 from 1996 to 2007 and discuss implications on the emission zones. The gamma-ray light curve consists of two peaks, separated by 0.41$\pm$0.02 in phase, with the first gamma-ray peak lagging the first radio peak by 0.05$\pm$0.02 in phase, consistent with the phase of giant radio pulses. We observe significant emission in the off-peak interval of PSR B1821-24 with a best-fit LAT position inconsistent with the core of M28. We do not detect significant gamma-ray pulsations at the spin or orbital periods from any other known pulsar in M28, and we place limits on the number of energetic pulsars in the cluster. The derived gamma-ray efficiency, ~2%, is typical of other gamma-ray pulsars with comparable spin-down power, suggesting that the measured spin-down rate ($2.2\times10^{36}$ erg s$^{-1}$) is not appreciably distorted by acceleration in the cluster potential. This confirms PSR B1821-24 as the second very energetic millisecond pulsar in a globular cluster and raises the question of whether these represent a separate class of objects that only form in regions of very high stellar density

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

Einstein@Home Discovery of 24 Pulsars in the Parkes Multi-beam Pulsar Survey

We have conducted a new search for radio pulsars in compact binary systems in the Parkes multi-beam pulsar survey (PMPS) data, employing novel methods to remove the Doppler modulation from binary motion. This has yielded unparalleled sensitivity to pulsars in compact binaries. The required computation time of approximately 17000 CPU core years was provided by the distributed volunteer computing project Einstein@Home, which has a sustained computing power of about 1 PFlop/s. We discovered 24 new pulsars in our search, of which 18 were isolated pulsars, and six were members of binary systems. Despite the wide filterbank channels and relatively slow sampling time of the PMPS data, we found pulsars with very large ratios of dispersion measure (DM) to spin period. Among those is PSR J1748-3009, the millisecond pulsar with the highest known DM (approximately 420 pc/cc). We also discovered PSR J1840-0643, which is in a binary system with an orbital period of 937 days, the fourth largest known. The new pulsar J1750-2536 likely belongs to the rare class of intermediate-mass binary pulsars. Three of the isolated pulsars show long-term nulling or intermittency in their emission, further increasing this growing family. Our discoveries demonstrate the value of distributed volunteer computing for data-driven astronomy and the importance of applying new analysis methods to extensively searched data.

preprint2013arXiv

Einstein@Home discovery of four young gamma-ray pulsars in Fermi LAT data

We report the discovery of four gamma-ray pulsars, detected in computing-intensive blind searches of data from the Fermi Large Area Telescope (LAT). The pulsars were found using a novel search approach, combining volunteer distributed computing via Einstein@Home and methods originally developed in gravitational-wave astronomy. The pulsars PSRs J0554+3107, J1422-6138, J1522-5735, and J1932+1916 are young and energetic, with characteristic ages between 35 and 56 kyr and spin-down powers in the range $6\times10^{34}$ - $10^{36}$ erg s$^{-1}$. They are located in the Galactic plane and have rotation rates of less than 10 Hz, among which the 2.1 Hz spin frequency of PSR J0554+3107 is the slowest of any known gamma-ray pulsar. For two of the new pulsars, we find supernova remnants coincident on the sky and discuss the plausibility of such associations. Deep radio follow-up observations found no pulsations, suggesting that all four pulsars are radio-quiet as viewed from Earth. These discoveries, the first gamma-ray pulsars found by volunteer computing, motivate continued blind pulsar searches of the many other unidentified LAT gamma-ray sources.

preprint2013arXiv

Fermi LAT pulsed detection of PSR J0737-3039A in the double pulsar system

We report the Fermi Large Area Telescope discovery of gamma-ray pulsations from the 22.7 ms pulsar A in the double pulsar system J0737-3039A/B. This is the first mildly recycled millisecond pulsar (MSP) detected in the GeV domain. The 2.7 s companion object PSR J0737-3039B is not detected in gamma rays. PSR J0737-3039A is a faint gamma-ray emitter, so that its spectral properties are only weakly constrained; however, its measured efficiency is typical of other MSPs. The two peaks of the gamma-ray light curve are separated by roughly half a rotation and are well offset from the radio and X-ray emission, suggesting that the GeV radiation originates in a distinct part of the magnetosphere from the other types of emission. From the modeling of the radio and the gamma-ray emission profiles and the analysis of radio polarization data, we constrain the magnetic inclination $α$ and the viewing angle $ζ$ to be close to 90$^\circ$, which is consistent with independent studies of the radio emission from PSR J0737-3039A. A small misalignment angle between the pulsar's spin axis and the system's orbital axis is therefore favored, supporting the hypothesis that pulsar B was formed in a nearly symmetric supernova explosion as has been discussed in the literature already.

preprint2013arXiv

Handbook of LHC Higgs Cross Sections: 3. Higgs Properties

This Report summarizes the results of the activities in 2012 and the first half of 2013 of the LHC Higgs Cross Section Working Group. The main goal of the working group was to present the state of the art of Higgs Physics at the LHC, integrating all new results that have appeared in the last few years. This report follows the first working group report Handbook of LHC Higgs Cross Sections: 1. Inclusive Observables (CERN-2011-002) and the second working group report Handbook of LHC Higgs Cross Sections: 2. Differential Distributions (CERN-2012-002). After the discovery of a Higgs boson at the LHC in mid-2012 this report focuses on refined prediction of Standard Model (SM) Higgs phenomenology around the experimentally observed value of 125-126 GeV, refined predictions for heavy SM-like Higgs bosons as well as predictions in the Minimal Supersymmetric Standard Model and first steps to go beyond these models. The other main focus is on the extraction of the characteristics and properties of the newly discovered particle such as couplings to SM particles, spin and CP-quantum numbers etc.

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: 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

Pathway to the Square Kilometre Array - The German White Paper -

The Square Kilometre Array (SKA) is the most ambitious radio telescope ever planned. With a collecting area of about a square kilometre, the SKA will be far superior in sensitivity and observing speed to all current radio facilities. The scientific capability promised by the SKA and its technological challenges provide an ideal base for interdisciplinary research, technology transfer, and collaboration between universities, research centres and industry. The SKA in the radio regime and the European Extreme Large Telescope (E-ELT) in the optical band are on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI) and have been recognised as the essential facilities for European research in astronomy. This "White Paper" outlines the German science and R&D interests in the SKA project and will provide the basis for future funding applications to secure German involvement in the Square Kilometre Array.

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

PSR J1723-2837: An Eclipsing Binary Radio Millisecond Pulsar

We present a study of PSR J1723-2837, an eclipsing, 1.86 ms millisecond binary radio pulsar discovered in the Parkes Multibeam survey. Radio timing indicates that the pulsar has a circular orbit with a 15 hr orbital period, a low-mass companion, and a measurable orbital period derivative. The eclipse fraction of ~15% during the pulsar's orbit is twice the Roche lobe size inferred for the companion. The timing behavior is significantly affected by unmodeled systematics of astrophysical origin, and higher-order orbital period derivatives are needed in the timing solution to account for these variations. We have identified the pulsar's (non-degenerate) companion using archival ultraviolet, optical, and infrared survey data and new optical photometry. Doppler shifts from optical spectroscopy confirm the star's association with the pulsar and indicate a pulsar-to-companion mass ratio of 3.3 +/- 0.5, corresponding to a companion mass range of 0.4 to 0.7 Msun and an orbital inclination angle range of between 30 and 41 degrees, assuming a pulsar mass range of 1.4-2.0 Msun. Spectroscopy indicates a spectral type of G for the companion and an inferred Roche-lobe-filling distance that is consistent with the distance estimated from radio dispersion. The features of PSR J1723-2837 indicate that it is likely a "redback" system. Unlike the five other Galactic redbacks discovered to date, PSR J1723-2837 has not been detected as a gamma-ray source with Fermi. This may be due to an intrinsic spin-down luminosity that is much smaller than the measured value if the unmeasured contribution from proper motion is large.

preprint2013arXiv

PSR J1840-1419: A very cool neutron star

We present upper limits on the X-ray emission for three neutron stars. For PSR J1840$-$1419, with a characteristic age of 16.5 Myr, we calculate a blackbody temperature upper limit (at 99% confidence) of $kT_{\mathrm{bb}}^{\infty}<24^{+17}_{-10}$ eV, making this one of the coolest neutron stars known. PSRs J1814$-$1744 and J1847$-$0130 are both high magnetic field pulsars, with inferred surface dipole magnetic field strengths of $5.5\times10^{13}$ and $9.4\times10^{13}$ G, respectively. Our temperature upper limits for these stars are $kT_{\mathrm{bb}}^{\infty}<123^{+20}_{-33}$ eV and $kT_{\mathrm{bb}}^{\infty}<115^{+16}_{-33}$ eV, showing that these high magnetic field pulsars are not significantly hotter than those with lower magnetic fields. Finally, we put these limits into context by summarizing all temperature measurements and limits for rotation-driven neutron stars.

preprint2013arXiv

Pulsar searches of Fermi unassociated sources with the Effelsberg telescope

Using the 100-m Effelsberg radio telescope operating at 1.36 GHz, we have performed a targeted radio pulsar survey of 289 unassociated gamma-ray sources discovered by the Large Area Telescope (LAT) aboard the Fermi satellite and published in the 1FGL catalogue (Abdo et al., 2010). This survey resulted in the discovery of millisecond pulsar J1745+1017, which resides in a short-period binary system with a low-mass companion, Mmin ~ 0.0137 Msun, indicative of `Black Widow' type systems. A two-year timing campaign has produced a refined radio ephemeris, accurate enough to allow for phase-folding of the LAT photons, resulting in the detection of a dual-peaked gamma-ray light-curve, proving that PSR J1745+1017 is the source responsible for the gamma-ray emission seen in 1FGL J1745.5 + 1018 (2FGL J1745.6+1015; Nolan et al., 2012). We find the gamma-ray spectrum of PSR J1745+1017 to be well modelled by an exponentially-cut-off power law with cut-off energy 3.2 GeV and photon index 1.6. The observed sources are known to contain a further 10 newly discovered pulsars which were undetected in this survey. Our radio observations of these sources are discussed and in all cases limiting flux densities are calculated. The reasons behind the seemingly low yield of discoveries are also discussed.

preprint2013arXiv

Pulse Broadening Measurements from the Galactic Center Pulsar J1745--2900

We present temporal scattering measurements of single pulses and average profiles of PSR J1745--2900, a magnetar recently discovered only 3 arcsec away from Sagittarius A* (Sgr A*), from 1.2 - 18.95 GHz using the Effelsberg 100-m Radio Telescope, the Nançay Decimetric Radio Telescope, and the Jodrell Bank Lovell Telescope. Single pulse analysis shows that the integrated pulse profile above 2 GHz is dominated by pulse jitter, while below 2 GHz the pulse profile shape is dominated by scattering. The high dispersion measure and rotation measure of the magnetar suggest that it is close to Sgr A* (within ~0.1 pc). This is the first object in the GC with both pulse broadening and angular broadening measurements. We measure a pulse broadening spectral index of alpha = -3.8 +/- 0.2 and a pulse broadening time scale at 1 GHz of tau_GHz = 1.3 +/- 0.2 s, which is several orders of magnitude lower than the scattering predicted by the NE2001 model (Cordes and Lazio 2002). If this scattering timescale is representative of the GC as a whole, then previous surveys should have detected many pulsars. The lack of detections implies either our understanding of scattering in the GC is incomplete or there are fewer pulsars in the GC than previously predicted. Given that magnetars are a rare class of radio pulsar, we believe that there are many canonical and millisecond pulsars in the GC, and not surprisingly, scattering regions in the GC have complex spatial structures.

preprint2013arXiv

Six millisecond pulsars detected by the Fermi Large Area Telescope and the radio/gamma-ray connection of millisecond pulsars

We report on the discovery of gamma-ray pulsations from five millisecond pulsars (MSPs) using the Fermi Large Area Telescope (LAT) and timing ephemerides provided by various radio observatories. We also present confirmation of the gamma-ray pulsations from a sixth source, PSR J2051-0827. Five of these six MSPs are in binary systems: PSRs J1713+0747, J1741+1351, J1600-3053 and the two black widow binary pulsars PSRs J0610-2100 and 2051-0827. The only isolated MSP is the nearby PSR J1024-0719, which is also known to emit X-rays. We present X-ray observations in the direction of PSRs J1600-3053 and J2051-0827. While the latter is firmly detected, we an only give upper limits for the X-ray flux of the former. There are no dedicated X-ray observations available for the other 3 objects. The MSPs mentioned above, together with most of the MSPs detected by Fermi, are used to put together a sample of 30 gamma-ray MSPs which is used to study the morphology and phase connection of radio and gamma-ray pulse profiles. We show that MSPs with pulsed gamma-ray emission which is phase aligned with the radio emission present the steepest radio spectra and the largest magnetic fields at the light cylinder among all MSPs. As well, we also observe a trend towards very low, or undetectable, radio linear polarisation levels. These properties could be attributed to caustic radio emission produced at a range of different altitudes in the magnetosphere. We note that most of these characteristics are also observed in the Crab pulsar, the only other radio pulsar known to exhibit phase-aligned radio and gamma-ray emission.

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

The double pulsar: evidence for neutron star formation without an iron core-collapse supernova

The double pulsar system PSR J0737-3039A/B is a double neutron star binary, with a 2.4-hour orbital period, which has allowed measurement of relativistic orbital perturbations to high precision. The low mass of the second-formed neutron star, as well as the low system eccentricity and proper motion, point to a different evolutionary scenario compared to most other known double neutron star systems. We describe analysis of the pulse profile shape over 6 years of observations, and present the resulting constraints on the system geometry. We find the recycled pulsar in this system, PSR J0737-3039A, to be a near-orthogonal rotator, with an average separation between its spin and magnetic axes of 90 +/- 11 +/- 5 deg. Furthermore, we find a mean 95% upper limit on the misalignment between its spin and orbital angular momentum axes of 3.2 deg, assuming that the observed emission comes from both magnetic poles. This tight constraint lends credence to the idea that the supernova that formed the second pulsar was relatively symmetric, possibly involving electron-capture onto an O-Ne-Mg core.

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 High Time Resolution Universe Pulsar Survey - VII: discovery of five millisecond pulsars and the different luminosity properties of binary and isolated recycled pulsars

This paper presents the discovery and timing parameters for five millisecond pulsars (MSPs), four in binary systems with probable white dwarf companions and one isolated, found in ongoing processing of the High Time Resolution Universe Pulsar Survey (HTRU). We also present high quality polarimetric data on four of them. These further discoveries confirm the high potential of our survey in finding pulsars with very short spin periods. At least two of these five MSPs are excellent candidates to be included in the Pulsar Timing Array projects. Thanks to the wealth of MSP discoveries in the HTRU survey, we revisit the question of whether the luminosity distributions of isolated and binary MSPs are different. Using the Cordes and Lazio distance model and our new and catalogue flux density measurements, we find that 41 of the 42 most luminous MSPs in the Galactic disk are in binaries and a statistical analysis suggests that the luminosity functions differ with 99.9% significance. We conclude that the formation process that leads to solitary MSPs affects their luminosities, despite their period and period derivatives being similar to those of pulsars in binary systems.

preprint2013arXiv

The High Time Resolution Universe Pulsar Survey VIII: The Galactic millisecond pulsar population

We have used millisecond pulsars (MSPs) from the southern High Time Resolution Universe (HTRU) intermediate latitude survey area to simulate the distribution and total population of MSPs in the Galaxy. Our model makes use of the scale factor method, which estimates the ratio of the total number of MSPs in the Galaxy to the known sample. Using our best fit value for the z-height, z=500 pc, we find an underlying population of MSPs of 8.3(\pm 4.2)*10^4 sources down to a limiting luminosity of L_min=0.1 mJy kpc^2 and a luminosity distribution with a steep slope of d\log N/d\log L = -1.45(\pm 0.14). However, at the low end of the luminosity distribution, the uncertainties introduced by small number statistics are large. By omitting very low luminosity pulsars, we find a Galactic population above L_min=0.2 mJy kpc^2 of only 3.0(\pm 0.7)*10^4 MSPs. We have also simulated pulsars with periods shorter than any known MSP, and estimate the maximum number of sub-MSPs in the Galaxy to be 7.8(\pm 5.0)*10^4 pulsars at L=0.1 mJy kpc^2. In addition, we estimate that the high and low latitude parts of the southern HTRU survey will detect 68 and 42 MSPs respectively, including 78 new discoveries. Pulsar luminosity, and hence flux density, is an important input parameter in the model. Some of the published flux densities for the pulsars in our sample do not agree with the observed flux densities from our data set, and we have instead calculated average luminosities from archival data from the Parkes Telescope. We found many luminosities to be very different than their catalogue values, leading to very different population estimates. Large variations in flux density highlight the importance of including scintillation effects in MSP population studies.

preprint2013arXiv

The High Time Resolution Universe Survey - IX: Polarimetry of long-period pulsars

We present a polarimetric analysis of 49 long-period pulsars discovered as part of the High Time Resolution Universe (HTRU) southern survey. The sources exhibit the typical characteristics of "old" pulsars, with low fractional linear and circular polarisation and narrow, multicomponent profiles. Although the position angle swings are generally complex, for two of the analysed pulsars (J1622-3751 and J1710-2616) we obtained an indication of the geometry via the rotating vector model. We were able to determine a value of the rotation measure (RM) for 34 of the sources which, when combined with their dispersion measures (DM), yields an integrated magnetic field strength along the line of sight. With the data presented here, the total number of values of RM associated to pulsars discovered during the HTRU southern survey sums to 51. The RMs are not consistent with the hypothesis of a counter-clockwise direction of the Galactic magnetic field within an annulus included between 4 and 6 kpc from the Galactic centre. A partial agreement with a counter-clockwise sense of the Galactic magnetic field within the spiral arms is however found in the area of the Carina-Sagittarius arm.

preprint2012arXiv

A Shapiro delay detection in the binary system hosting the millisecond pulsar PSR J1910-5959A

PSR J1910-5959A is a binary pulsar with a helium white dwarf companion located about 6 arcmin from the center of the globular cluster NGC6752. Based on 12 years of observations at the Parkes radio telescope, the relativistic Shapiro delay has been detected in this system. We obtain a companion mass Mc = 0.180+/-0.018Msun (1sigma) implying that the pulsar mass lies in the range 1.1Msun <= Mp <= 1.5Msun. We compare our results with previous optical determinations of the companion mass, and examine prospects for using this new measurement for calibrating the mass-radius relation for helium white dwarfs and for investigating their evolution in a pulsar binary system. Finally we examine the set of binary systems hosting a millisecond pulsar and a low mass helium white dwarf for which the mass of both stars has been measured. We confirm that the correlation between the companion mass and the orbital period predicted by Tauris & Savonije reproduces the observed values but find that the predicted Mp - Pb correlation over-estimates the neutron star mass by about 0.5Msun in the orbital period range covered by the observations. Moreover, a few systems do not obey the observed Mp - Pb correlation. We discuss these results in the framework of the mechanisms that inhibit the accretion of matter by a neutron star during its evolution in a low-mass X-ray binary.

preprint2012arXiv

A side-by-side comparison of Daya Bay antineutrino detectors

The Daya Bay Reactor Neutrino Experiment is designed to determine precisely the neutrino mixing angle $θ_{13}$ with a sensitivity better than 0.01 in the parameter sin$^22θ_{13}$ at the 90% confidence level. To achieve this goal, the collaboration will build eight functionally identical antineutrino detectors. The first two detectors have been constructed, installed and commissioned in Experimental Hall 1, with steady data-taking beginning September 23, 2011. A comparison of the data collected over the subsequent three months indicates that the detectors are functionally identical, and that detector-related systematic uncertainties exceed requirements.

preprint2012arXiv

A survey of nulling pulsars using the Giant Meterwave Radio Telescope

Several pulsars show sudden cessation of pulsed emission, which is known as pulsar nulling. In this paper, the nulling behaviour of 15 pulsars is presented. The nulling fractions of these pulsars, along with the degree of reduction in the pulse energy during the null phase, are reported for these pulsars. A quasi-periodic null-burst pattern is reported for PSR J1738-2330. The distributions of lengths of the null and the burst phases as well as the typical nulling time scales are estimated for eight strong pulsars. The nulling pattern of four pulsars with similar nulling fraction are found to be different from each other, suggesting that the fraction of null pulses does not quantify the nulling behaviour of a pulsar in full detail. Analysis of these distributions also indicate that while the null and the burst pulses occur in groups, the underlying distribution of the interval between a transition from the null to the burst phase and vice verse appears to be similar to that of a stochastic Poisson point process.

preprint2012arXiv

Application of the Gaussian mixture model in pulsar astronomy -- pulsar classification and candidates ranking for {\it Fermi} 2FGL catalog

Machine learning, algorithms to extract empirical knowledge from data, can be used to classify data, which is one of the most common tasks in observational astronomy. In this paper, we focus on Bayesian data classification algorithms using the Gaussian mixture model and show two applications in pulsar astronomy. After reviewing the Gaussian mixture model and the related Expectation-Maximization algorithm, we present a data classification method using the Neyman-Pearson test. To demonstrate the method, we apply the algorithm to two classification problems. Firstly, it is applied to the well known period-period derivative diagram, where we find that the pulsar distribution can be modeled with six Gaussian clusters, with two clusters for millisecond pulsars (recycled pulsars) and the rest for normal pulsars. From this distribution, we derive an empirical definition for millisecond pulsars as $\frac{\dot{P}}{10^{-17}} \leq3.23(\frac{P}{100 \textrm{ms}})^{-2.34}$. The two millisecond pulsar clusters may have different evolutionary origins, since the companion stars to these pulsars in the two clusters show different chemical composition. Four clusters are found for normal pulsars. Possible implications for these clusters are also discussed. Our second example is to calculate the likelihood of unidentified \textit{Fermi} point sources being pulsars and rank them accordingly. In the ranked point source list, the top 5% sources contain 50% known pulsars, the top 50% contain 99% known pulsars, and no known active galaxy (the other major population) appears in the top 6%. Such a ranked list can be used to help the future follow-up observations for finding pulsars in unidentified \textit{Fermi} point sources.

preprint2012arXiv

Binary Millisecond Pulsar Discovery via Gamma-Ray Pulsations

Millisecond pulsars, old neutron stars spun-up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.

preprint2012arXiv

Broadband nulling behaviour of PSR B2319+60

Pulse nulling is one of many single pulse phenomena exhibited by radio pulsars. The broadband nature of this phenomena remained unexplained due to lack of coordinated observations. We are reporting broadband nulling behaviour of a well known nulling pulsar PSR B2319+60. The simultaneous observations were carried out at four different frequencies using the Giant Meterwave Radio telescope (325 and 610 MHz), the Westerbork Synthesis Radio Telescope (1420 MHz) and the Effelsberg radio telescope (5100 MHz). The nulling fractions were estimated at all four observed frequencies, which suggest similar degree of nulling across these frequencies. To comment on the broadband behaviour of pulse nulling, we calculated the Pearson cross-correlation coefficients for the occurrence of null and burst pulses between all four frequencies. We conclude that nulling is largely a broadband phenomenon for this pulsar and it points towards a global failure of the magnetospheric currents which produces radio emission at these frequencies.

preprint2012arXiv

Can we see pulsars around Sgr A*? - The latest searches with the Effelsberg telescope

Radio pulsars in relativistic binary systems are unique tools to study the curved space-time around massive compact objects. The discovery of a pulsar closely orbiting the super-massive black hole at the centre of our Galaxy, Sgr A*, would provide a superb test-bed for gravitational physics. To date, the absence of any radio pulsar discoveries within a few arc minutes of Sgr A* has been explained by one principal factor: extreme scattering of radio waves caused by inhomogeneities in the ionized component of the interstellar medium in the central 100 pc around Sgr A*. Scattering, which causes temporal broadening of pulses, can only be mitigated by observing at higher frequencies. Here we describe recent searches of the Galactic centre region performed at a frequency of 18.95 GHz with the Effelsberg radio telescope.

preprint2012arXiv

Changes in Polarization Position Angle across the Eclipse in the Double Pulsar System

We investigate the changes in polarization position angle in radiation from pulsar A around the eclipse in the Double Pulsar system PSR J0737-3039A/B at the 20 cm and 50 cm wavelengths using the Parkes 64-m telescope. The changes are ~2σ during and shortly after the eclipse at 20 cm but less significant at 50 cm. We show that the changes in position angle during the eclipse can be modelled by differential synchrotron absorption in the eclipse regions. Position angle changes after the eclipse are interpreted as Faraday rotation in the magnetotail of pulsar B. Implied charge densities are consistent with the Goldreich-Julian density, suggesting that the particle energies in the magnetotail are mildly relativistic.

preprint2012arXiv

Constraining the optical emission from the double pulsar system J0737-3039

We present the first optical observations of the unique system J0737-3039 (composed of two pulsars, hereafter PSR-A and PSR-B). Ultra-deep optical observations, performed with the High Resolution Camera of the Advanced Camera for Surveys on board the Hubble Space Telescope could not detect any optical emission from the system down to m_F435W=27.0 and m_F606W=28.3. The estimated optical flux limits are used to constrain the three-component (two thermal and one non-thermal) model recently proposed to reproduce the XMM-Newton X-ray spectrum. They suggest the presence of a break at low energies in the non-thermal power law component of PSR-A and are compatible with the expected black-body emission from the PSR-B surface. The corresponding efficiency of the optical emission from PSR-A's magnetosphere would be comparable to that of other Myr-old pulsars, thus suggesting that this parameter may not dramatically evolve over a time-scale of a few Myr.

preprint2012arXiv

Discovery of the millisecond pulsar PSR J2043+1711 in a Fermi source with the Nancay Radio Telescope

We report the discovery of the millisecond pulsar PSR J2043+1711 in a search of a Fermi Large Area Telescope (LAT) source with no known associations, with the Nancay Radio Telescope. The new pulsar, confirmed with the Green Bank Telescope, has a spin period of 2.38 ms, is relatively nearby (d <~ 2 kpc), and is in a 1.48 day orbit around a low mass companion, probably a He-type white dwarf. Pulsed gamma-ray emission was detected in the data recorded by the Fermi LAT. The gamma-ray light curve and spectral properties are typical of other gamma-ray millisecond pulsars seen with Fermi. X-ray observations of the pulsar with Suzaku and the Swift/XRT yielded no detection. At 1.4 GHz we observe strong flux density variations because of interstellar diffractive scintillation, however a sharp peak can be observed at this frequency during bright scintillation states. At 327 MHz the pulsar is detected with a much higher signal-to-noise ratio and its flux density is far more steady. However, at that frequency the Arecibo instrumentation cannot yet fully resolve the pulse profile. Despite that, our pulse time-of-arrival measurements have a post-fit residual rms of 2 \mus. This and the expected stability of this system has made PSR J2043+1711 one of the first new Fermi-selected millisecond pulsars to be added to pulsar gravitational wave timing arrays. It has also allowed a significant measurement of relativistic delays in the times of arrival of the pulses due to the curvature of space-time near the companion, but not yet with enough precision to derive useful masses for the pulsar and the companion. A mass for the pulsar between 1.7 and 2.0 solar masses can be derived if a standard millisecond pulsar formation model is assumed. In this article we also present a comprehensive summary of pulsar searches in Fermi LAT sources with the Nancay Radio Telescope to date.

preprint2012arXiv

Five New Millisecond Pulsars From a Radio Survey of 14 Unidentified Fermi-LAT Gamma-ray Sources

We have discovered five millisecond pulsars (MSPs) in a survey of 14 unidentified Fermi-LAT sources in the southern sky using the Parkes radio telescope. PSRs J0101-6422, J1514-4946, and J1902-5105 reside in binaries, while PSRs J1658-5324 and J1747-4036 are isolated. Using an ephemeris derived from timing observations of PSR J0101-6422 (P =2.57 ms, DM=12 pc cm-3), we have detected γ-ray pulsations and measured its proper motion. Its γ-ray spectrum (a power law of Γ = 0.9 with a cutoff at 1.6 GeV) and efficiency are typical of other MSPs, but its radio and γ-ray light curves challenge simple geometric models of emission. The high success rate of this survey-enabled by selecting γ-ray sources based on their detailed spectral characteristics-and other similarly successful searches indicate that a substantial fraction of the local population of MSPs may soon be known.

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

Handbook of LHC Higgs Cross Sections: 2. Differential Distributions

This Report summarises the results of the second year's activities of the LHC Higgs Cross Section Working Group. The main goal of the working group was to present the state of the art of Higgs Physics at the LHC, integrating all new results that have appeared in the last few years. The first working group report Handbook of LHC Higgs Cross Sections: 1. Inclusive Observables (CERN-2011-002) focuses on predictions (central values and errors) for total Higgs production cross sections and Higgs branching ratios in the Standard Model and its minimal supersymmetric extension, covering also related issues such as Monte Carlo generators, parton distribution functions, and pseudo-observables. This second Report represents the next natural step towards realistic predictions upon providing results on cross sections with benchmark cuts, differential distributions, details of specific decay channels, and further recent developments.

preprint2012arXiv

Improved Measurement of Electron Antineutrino Disappearance at Daya Bay

We report an improved measurement of the neutrino mixing angle $θ_{13}$ from the Daya Bay Reactor Neutrino Experiment. We exclude a zero value for $\sin^22θ_{13}$ with a significance of 7.7 standard deviations. Electron antineutrinos from six reactors of 2.9 GW$_{\rm th}$ were detected in six antineutrino detectors deployed in two near (flux-weighted baselines of 470 m and 576 m) and one far (1648 m) underground experimental halls. Using 139 days of data, 28909 (205308) electron antineutrino candidates were detected at the far hall (near halls). The ratio of the observed to the expected number of antineutrinos assuming no oscillations at the far hall is $0.944\pm 0.007({\rm stat.}) \pm 0.003({\rm syst.})$. An analysis of the relative rates in six detectors finds $\sin^22θ_{13}=0.089\pm 0.010({\rm stat.})\pm0.005({\rm syst.})$ in a three-neutrino framework.

preprint2012arXiv

Observation of electron-antineutrino disappearance at Daya Bay

The Daya Bay Reactor Neutrino Experiment has measured a non-zero value for the neutrino mixing angle $θ_{13}$ with a significance of 5.2 standard deviations. Antineutrinos from six 2.9 GW$_{\rm th}$ reactors were detected in six antineutrino detectors deployed in two near (flux-weighted baseline 470 m and 576 m) and one far (1648 m) underground experimental halls. With a 43,000 ton-GW_{\rm th}-day livetime exposure in 55 days, 10416 (80376) electron antineutrino candidates were detected at the far hall (near halls). The ratio of the observed to expected number of antineutrinos at the far hall is $R=0.940\pm 0.011({\rm stat}) \pm 0.004({\rm syst})$. A rate-only analysis finds $\sin^22θ_{13}=0.092\pm 0.016({\rm stat})\pm0.005({\rm syst})$ in a three-neutrino framework.

preprint2012arXiv

On the origin of a highly-dispersed coherent radio burst

We discuss the possible source of a highly-dispersed radio transient discovered in the Parkes Multi-beam Pulsar Survey (PMPS). The pulse has a dispersion measure of $746\mathrm{cm}^{-3}\mathrm{pc}$, a peak flux density of 400 mJy for the observed pulse width of 7.8 ms, and a flat spectrum across a 288-MHz band centred on 1374 MHz. The flat spectrum suggests that the pulse did not originate from a pulsar, but is consistent with radio-emitting magnetar spectra. The non-detection of subsequent bursts constrains any possible pulsar period to $\gtrsim1$ s, and the pulse energy distribution to being much flatter than typical giant pulse emitting pulsars. The burst is also consistent with the radio signal theorised from an annihilating mini black hole. Extrapolating the PMPS detection rate, provides a limit of $Ω_{BH}\lesssim5\times10^{-14}$ on the density of these objects. We investigate the consistency of these two scenarios, plus several other possible solutions, as potential explanations to the origin of the pulse, as well as for another transient with similar properties: the Lorimer Burst.

preprint2012arXiv

On the origin of LS 5039 and PSR J1825-1446

Context. The gamma-ray binary LS 5039 and the isolated pulsar PSR J1825-1446 were proposed to have been formed in the supernova remnant (SNR) G016.8-01.1. Aims. We aim to obtain the Galactic trajectory of LS 5039 and PSR J1825-1446 to find their origin in the Galaxy, and in particular to check their association with SNR G016.8-01.1 to restrict their age. Methods. By means of radio and optical observations we obtained the proper motion and the space velocity of the sources. Results. The proper motion of PSR J1825-1446 corresponds to a transverse space velocity of 690 km/s at a distance of 5 kpc. Its Galactic velocity at different distances is not compatible with the expected Galactic rotation. The velocity and characteristic age of PSR J1825-1446 make it incompatible with SNR G016.8-01.1. There are no clear OB associations or SNRs crossing the past trajectory of PSR J1825-1446. We estimate the age of the pulsar to be 80-245 kyr, which is compatible with its characteristic age. The proper motion of LS 5039 is 7.09 and -8.82 mas/yr in right ascension and declination, respectively. The association of LS 5039 with SNR G016.8-01.1 is unlikely, although we cannot to discard it. The system would have had to be formed in the association Ser OB2 (at 2.0 kpc) if the age of the system is 1.0-1.2 Myr, or in the association Sct OB3 (distance 1.5-2 kpc) for an age of 0.1-0.2 Myr. If the system were not formed close to Ser OB2, the pseudo-synchronization of the orbit would be unlikely. Conclusions. PSR J1825-1446 is a high-velocity isolated pulsar ejected from the Galaxy. The distance to LS 5039, which needs to be constrained by future astrometric missions such as Gaia, is a key parameter for restricting its origin and age.

preprint2012arXiv

On the Pulse Intensity Modulation of PSR B0823+26

We investigate the radio emission behaviour of PSR B0823+26, a pulsar which is known to undergo pulse nulling, using an 153-d intensive sequence of observations. The pulsar is found to exhibit both short (~min) and unusually long-term (~hours or more) nulls, which not only suggest that the source possesses a distribution of nulling timescales, but that it may also provide a link between conventional nulling pulsars and longer-term intermittent pulsars. Despite seeing evidence for periodicities in the pulsar radio emission, we are uncertain whether they are intrinsic to the source, due to the influence of observation sampling on the periodicity analysis performed. Remarkably, we find evidence to suggest that the pulsar may undergo pre-ignition periods of 'emission flickering', that is rapid changes between radio-on (active) and -off (null) emission states, before transitioning to a steady radio-emitting phase. We find no direct evidence to indicate that the object exhibits any change in spin-down rate between its radio-on and -off emission modes. We do, however, place an upper limit on this variation to be <= 6 % from simulations. This indicates that emission cessation in pulsars does not necessarily lead to large changes in spin-down rate. Moreover, we show that such changes in spin-down rate will not be discernible in the majority of objects which exhibit short-term (<= 1 d) emission cessation. In light of this, we predict that many pulsars could exhibit similar magnetospheric and emission properties to PSR B0823+26, but which have not yet been observed.

preprint2012arXiv

Prospects for probing strong gravity with a pulsar-black hole system

The discovery of a pulsar (PSR) in orbit around a black hole (BH) is expected to provide a superb new probe of relativistic gravity and BH properties. Apart from a precise mass measurement for the BH, one could expect a clean verification of the dragging of space-time caused by the BH spin. In order to measure the quadrupole moment of the BH for testing the no-hair theorem of general relativity (GR), one has to hope for a sufficiently massive BH. In this respect, a PSR orbiting the super-massive BH in the center of our Galaxy would be the ultimate laboratory for gravity tests with PSRs. But even for gravity theories that predict the same properties for BHs as GR, a PSR-BH system would constitute an excellent test system, due to the high grade of asymmetry in the strong field properties of these two components. Here we highlight some of the potential gravity tests that one could expect from different PSR-BH systems, utilizing present and future radio telescopes, like FAST and SKA.

preprint2012arXiv

PSR J1838-0537: Discovery of a young, energetic gamma-ray pulsar

We report the discovery of PSR J1838-0537, a gamma-ray pulsar found through a blind search of data from the Fermi Large Area Telescope (LAT). The pulsar has a spin frequency of 6.9 Hz and a frequency derivative of -2.2e-11 Hz/s, implying a young characteristic age of 4970 years and a large spin-down power of 5.9e36 erg/s. Follow-up observations with radio telescopes detected no pulsations, thus PSR J1838-0537 appears radio-quiet as viewed from Earth. In September 2009 the pulsar suffered the largest glitch so far seen in any gamma-ray-only pulsar, causing a relative increase in spin frequency of about 5.5e-6. After the glitch, during a putative recovery period, the timing analysis is complicated by the sparsity of the LAT photon data, the weakness of the pulsations, and the reduction in average exposure from a coincidental, contemporaneous change in the LAT's sky-survey observing pattern. The pulsar's sky position is coincident with the spatially extended TeV source HESS J1841-055 detected by the High Energy Stereoscopic System (H.E.S.S.). The inferred energetics suggest that HESS J1841-055 contains a pulsar wind nebula powered by the pulsar.

preprint2012arXiv

Radio Detection of the Fermi LAT Blind Search Millisecond Pulsar J1311-3430

We report the detection of radio emission from PSR J1311-3430, the first millisecond pulsar discovered in a blind search of Fermi Large Area Telescope (LAT) gamma-ray data. We detected radio pulsations at 2 GHz, visible for <10% of ~4.5-hrs of observations using the Green Bank Telescope (GBT). Observations at 5 GHz with the GBT and at several lower frequencies with Parkes, Nancay, and the Giant Metrewave Radio Telescope resulted in non-detections. We also report the faint detection of a steep spectrum continuum radio source (0.1 mJy at 5 GHz) in interferometric imaging observations with the Jansky Very Large Array. These detections demonstrate that PSR J1311-3430, is not radio quiet and provides additional evidence that the radio beaming fraction of millisecond pulsars is very large. The radio detection yields a distance estimate of 1.4 kpc for the system, yielding a gamma-ray efficiency of 30%, typical of LAT-detected MSPs. We see apparent excess delay in the radio pulsar as the pulsar appears from eclipse and we speculate on possible mechanisms for the non-detections of the pulse at other orbital phases and observing frequencies.

preprint2012arXiv

Radio emission evolution, polarimetry and multifrequency single pulse analysis of the radio magnetar PSR J1622-4950

Here we report on observations of the radio magnetar PSR J1622-4950 at frequencies from 1.4 to 17 GHz. We show that although its flux density is varying up to a factor of ~10 within a few days, it has on average decreased by a factor of 2 over the last 700 days. At the same time, timing analysis indicates a trend of decreasing spin-down rate over our entire data set, again of about a factor of 2 over 700 days, but also an erratic variability in the spin-down rate within this time span. Integrated pulse profiles are often close to 100 per cent linearly polarized, but large variations in both the profile shape and fractional polarization are regularly observed. Furthermore, the behaviour of the position angle of the linear polarization is very complex - offsets in both the absolute position angle and the phase of the position angle sweep are often seen and the occasional presence of orthogonal mode jumps further complicates the picture. However, model fitting indicates that the magnetic and rotation axes are close to aligned. Finally, a single pulse analysis has been carried out at four observing frequencies, demonstrating that the wide pulse profile is built up of narrow spikes of emission, with widths that scale inversely with observing frequency. All three of the known radio magnetars seem to have similar characteristics, with highly polarized emission, time-variable flux density and pulse profiles, and with spectral indices close to zero.

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

Recycling Pulsars: spins, masses and ages

Although the first millisecond pulsars (MSPs) were discovered 30 years ago we still do not understand all details of their formation process. Here, we present new results from Tauris, Langer & Kramer (2012) on the recycling scenario leading to radio MSPs with helium or carbon-oxygen white dwarf companions via evolution of low- and intermediate mass X-ray binaries (LMXBs, IMXBs). We discuss the location of the spin-up line in the (P,Pdot)-diagram and estimate the amount of accreted mass needed to obtain a given spin period and compare with observations. Finally, we constrain the true ages of observed recycled pulsars via calculated isochrones in the (P,Pdot)-diagram.

preprint2012arXiv

Searches for New Physics: Les Houches Recommendations for the Presentation of LHC Results

We present a set of recommendations for the presentation of LHC results on searches for new physics, which are aimed at providing a more efficient flow of scientific information between the experimental collaborations and the rest of the high energy physics community, and at facilitating the interpretation of the results in a wide class of models. Implementing these recommendations would aid the full exploitation of the physics potential of the LHC.

preprint2012arXiv

The Double Pulsar Eclipses I: Phenomenology and Multi-frequency Analysis

The double pulsar PSR J0737-3039A/B displays short, 30 s eclipses that arise around conjunction when the radio waves emitted by pulsar A are absorbed as they propagate through the magnetosphere of its companion pulsar B. These eclipses offer a unique opportunity to probe directly the magnetospheric structure and the plasma properties of pulsar B. We have performed a comprehensive analysis of the eclipse phenomenology using multi-frequency radio observations obtained with the Green Bank Telescope. We have characterized the periodic flux modulations previously discovered at 820 MHz by McLaughlin et al., and investigated the radio frequency dependence of the duration and depth of the eclipses. Based on their weak radio frequency evolution, we conclude that the plasma in pulsar B's magnetosphere requires a large multiplicity factor (~ 10^5). We also found that, as expected, flux modulations are present at all radio frequencies in which eclipses can be detected. Their complex behavior is consistent with the confinement of the absorbing plasma in the dipolar magnetic field of pulsar B as suggested by Lyutikov & Thompson and such a geometric connection explains that the observed periodicity is harmonically related to pulsar B's spin frequency. We observe that the eclipses require a sharp transition region beyond which the plasma density drops off abruptly. Such a region defines a plasmasphere which would be well inside the magnetospheric boundary of an undisturbed pulsar. It is also two times smaller than the expected standoff radius calculated using the balance of the wind pressure from pulsar A and the nominally estimated magnetic pressure of pulsar B.

preprint2012arXiv

The High Time Resolution Universe Survey - V: Single-pulse energetics and modulation properties of 315 pulsars

We report on the pulse-to-pulse energy distributions and phase-resolved modulation properties for catalogued pulsars in the southern High Time Resolution Universe intermediate-latitude survey. We selected the 315 pulsars detected in a single-pulse search of this survey, allowing a large sample unbiased regarding any rotational parameters of neutron stars. We found that the energy distribution of many pulsars is well-described by a log-normal distribution, with few deviating from a small range in log-normal scale and location parameters. Some pulsars exhibited multiple energy states corresponding to mode changes, and implying that some observed "nulling" may actually be a mode-change effect. PSRJ1900-2600 was found to emit weakly in its previously-identified "null" state. We found evidence for another state-change effect in two pulsars, which show bimodality in their nulling time scales; that is, they switch between a continuous-emission state and a single-pulse-emitting state. Large modulation occurs in many pulsars across the full integrated profile, with increased sporadic bursts at leading and trailing sub-beam edges. Some of these high-energy outbursts may indicate the presence of "giant pulse" phenomena. We found no correlation with modulation and pulsar period, age, or other parameters. Finally, the deviation of integrated pulse energy from its average value was generally quite small, despite the significant phase-resolved modulation in some pulsars; we interpret this as tenuous evidence of energy regulation between distinct pulsar sub-beams.

preprint2012arXiv

The High Time Resolution Universe Survey VI: An Artificial Neural Network and Timing of 75 Pulsars

We present 75 pulsars discovered in the mid-latitude portion of the High Time Resolution Universe survey, 54 of which have full timing solutions. All the pulsars have spin periods greater than 100 ms, and none of those with timing solutions are in binaries. Two display particularly interesting behaviour; PSR J1054-5944 is found to be an intermittent pulsar, and PSR J1809-0119 has glitched twice since its discovery. In the second half of the paper we discuss the development and application of an artificial neural network in the data-processing pipeline for the survey. We discuss the tests that were used to generate scores and find that our neural network was able to reject over 99% of the candidates produced in the data processing, and able to blindly detect 85% of pulsars. We suggest that improvements to the accuracy should be possible if further care is taken when training an artificial neural network; for example ensuring that a representative sample of the pulsar population is used during the training process, or the use of different artificial neural networks for the detection of different types of pulsars.

preprint2012arXiv

The optimal schedule for pulsar timing array observations

In order to maximize the sensitivity of pulsar timing arrays to a stochastic gravitational wave background, we present computational techniques to optimize observing schedules. The techniques are applicable to both single and multi-telescope experiments. The observing schedule is optimized for each telescope by adjusting the observing time allocated to each pulsar while keeping the total amount of observing time constant. The optimized schedule depends on the timing noise characteristics of each individual pulsar as well as the performance of instrumentation. Several examples are given to illustrate the effects of different types of noise. A method to select the most suitable pulsars to be included in a pulsar timing array project is also presented.

preprint2012arXiv

The relativistic pulsar-white dwarf binary PSR J1738+0333 I. Mass determination and evolutionary history

PSR J1738+0333 is one of the four millisecond pulsars known to be orbited by a white dwarf companion bright enough for optical spectroscopy. Of these, it has the shortest orbital period, making it especially interesting for a range of astrophysical and gravity related questions. We present a spectroscopic and photometric study of the white dwarf companion and infer its radial velocity curve, effective temperature, surface gravity and luminosity. We find that the white dwarf has properties consistent with those of low-mass white dwarfs with thick hydrogen envelopes, and use the corresponding mass-radius relation to infer its mass; M_WD = 0.181 +/- +0.007/-0.005 solar masses. Combined with the mass ratio q=8.1 +/- 0.2 inferred from the radial velocities and the precise pulsar timing ephemeris, the neutron star mass is constrained to M_PSR = 1.47 +/- +0.07/-0.06 solar masses. Contrary to expectations, the latter is only slightly above the Chandrasekhar limit. We find that, even if the birth mass of the neutron star was only 1.20 solar masses, more than 60% of the matter that left the surface of the white dwarf progenitor escaped the system. The accurate determination of the component masses transforms this system in a laboratory for fundamental physics by constraining the orbital decay predicted by general relativity. Currently, the agreement is within 1 sigma of the observed decay. Further radio timing observations will allow precise tests of white dwarf models, assuming the validity of general relativity.

preprint2012arXiv

Timing the main-sequence-star binary pulsar J1740-3052

PSR J1740-3052 is a young pulsar in orbit around a companion that is most likely a B-type main-sequence star. Since its discovery more than a decade ago, data have been taken at several frequencies with instruments at the Green Bank, Parkes, Lovell, and Westerbork telescopes. We measure scattering timescales in the pulse profiles and dispersion measure changes as a function of binary orbital phase and present evidence that both of these vary as would be expected due to a wind from the companion star. Using pulse arrival times that have been corrected for the observed periodic dispersion measure changes, we find a timing solution spanning 1997 November to 2011 March. This includes measurements of the advance of periastron and the change in the projected semimajor axis of the orbit and sets constraints on the orbital geometry. From these constraints, we estimate that the pulsar received a kick of at least ~50 km/s at birth. A quasi-periodic signal is present in the timing residuals with a period of 2.2 times the binary orbital period. The origin of this signal is unclear.

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

A search for optical bursts from RRAT J1819-1458: II. Simultaneous ULTRACAM-Lovell Telescope observations

The Rotating RAdio Transient (RRAT) J1819-1458 exhibits ~3 ms bursts in the radio every ~3 min, implying that it is visible for only ~1 per day. Assuming that the optical light behaves in a similar manner, long exposures of the field would be relatively insensitive due to the accumulation of sky photons. A much better way of detecting optical emission from J1819-1458 would then be to observe with a high-speed optical camera simultaneously with radio observations, and co-add only those optical frames coincident with the dispersion-corrected radio bursts. We present the results of such a search, using simultaneous ULTRACAM and Lovell Telescope observations. We find no evidence for optical bursts in J1819-1458 at magnitudes brighter than i'=19.3 (5-sigma limit). This is nearly 3 magnitudes fainter than the previous burst limit, which had no simultaneous radio observations.

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

Detection of the Pulsar Wind Nebula HESS J1825-137 with the Fermi Large Area Telescope

We announce the discovery of 1 - 100 GeV gamma-ray emission from the archetypal TeV pulsar wind nebula HESS J1825-137 using 20 months of survey data from the Fermi Large Area Telescope (LAT). The gamma-ray emission detected by the LAT is significantly spatially extended, with a best-fit rms extension of sigma = 0.56° $\pm$ 0.07° for an assumed Gaussian model. The 1 - 100 GeV LAT spectrum of this source is well described by a power-law with a spectral index of 1.38 $\pm$ 0.12 $\pm$ 0.16 and an integral flux above 1 GeV of (6.50 $\pm$ 0.21 $\pm$ 3.90) x 10^{-9} cm^{-2} s^{-1}. The first errors represent the statistical errors on the fit parameters, while the second ones are the systematic uncertainties. Detailed morphological and spectral analyses bring new constraints on the energetics and magnetic field of the pulsar wind nebula system. The spatial extent and hard spectrum of the GeV emission are consistent with the picture of an inverse Compton origin of the GeV-TeV emission in a cooling-limited nebula powered by the pulsar PSR J1826-1334.

preprint2011arXiv

Discovery of gamma and X-ray pulsations from the young and energetic PSR J1357-6429 with Fermi and XMM-Newton

Since the launch of the Fermi satellite, the number of known gamma-ray pulsars has increased tenfold. Most gamma-ray detected pulsars are young and energetic, and many are associated with TeV sources. PSR J1357-6429 is a high spin-down power pulsar (Edot = 3.1 * 10^36 erg/s), discovered during the Parkes multibeam survey of the Galactic plane, with significant timing noise typical of very young pulsars. In the very-high-energy domain, H.E.S.S. has reported the detection of the extended source HESS J1356-645 (intrinsic Gaussian width of 12') whose centroid lies 7' from PSR J1357-6429. Using a rotational ephemeris obtained with 74 observations made with the Parkes telescope at 1.4 GHz, we phase-fold more than two years of gamma-ray data acquired by Fermi-LAT as well as those collected with XMM-Newton, and perform gamma-ray spectral modeling. Significant gamma and X-ray pulsations are detected from PSR J1357-6429. The light curve in both bands shows one broad peak. Gamma-ray spectral analysis of the pulsed emission suggests that it is well described by a simple power-law of index 1.5 +/- 0.3stat +/- 0.3syst with an exponential cut-off at 0.8 +/- 0.3stat +/- 0.3syst GeV and an integral photon flux above 100 MeV of (6.5 +/- 1.6stat +/- 2.3syst) * 10^-8 cm^-2 s^-1. The X-ray spectra obtained from the new data provide results consistent with those reported by Zavlin (2007). Upper limits on the gamma-ray emission from its potential pulsar wind nebula (PWN) are also reported. Assuming a distance of 2.4 kpc, the Fermi LAT energy flux yields a gamma-ray luminosity for PSR J1357-6429 of L_gamma = (2.13 +/- 0.25stat +/- 0.83syst) * 10^34 erg/s, consistent with an L_gamma \propto sqrt(Edot) relationship. The Fermi non-detection of the pulsar wind nebula associated with HESS J1356-645 provides new constraints on the electron population responsible for the extended TeV emission.

preprint2011arXiv

Discovery of millisecond pulsars in radio searches of southern Fermi LAT sources

Using the Parkes radio telescope we have carried out deep observations of eleven unassociated gamma-ray sources. Periodicity searches of these data have discovered two millisecond pulsars, PSR J1103-5403 (1FGL J1103.9-5355) and PSR J2241-5236 (1FGL J2241.9-5236), and a long period pulsar, PSR J1604-44 (1FGL J1604.7-4443). In addition we searched for but did not detect any radio pulsations from six gammaray pulsars discovered by the Fermi satellite to a level of - 0.04 mJy (for pulsars with a 10% duty cycle). Timing of the millisecond pulsar PSR J1103-5403 has shown that its position is 9' from the centroid of the gamma-ray source. Since these observations were carried out, independent evidence has shown that 1FGL J1103.9-5355 is associated with the flat spectrum radio source PKS 1101-536. It appears certain that the pulsar is not associated with the gamma-ray source, despite the seemingly low probability of a chance detection of a radio millisecond pulsar. We consider that PSR J1604-44 is a chance discovery of a weak, long period pulsar and is unlikely to be associated with 1FGL J1604.7-4443. PSR J2241-5236 has a spin period of 2.2 ms and orbits a very low mass companion with a 3.5 hour orbital period. The relatively high flux density and low dispersion measure of PSR J2241-5236 makes it an excellent candidate for high precision timing experiments. The gamma-rays of 1FGL J2241.9-5236 have a spectrum that is well modelled by a power law with exponential cutoff, and phasebinning with the radio ephemeris results in a multi-peaked gamma-ray pulse profile. Observations with Chandra have identified a coincident X-ray source within 0.1" of the position of the pulsar obtained by radio timing

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

Discovery of two millisecond pulsars in Fermi sources with the Nancay Radio Telescope

We report the discovery of two millisecond pulsars in a search for radio pulsations at the positions of \emph{Fermi Large Area Telescope} sources with no previously known counterparts, using the Nançay radio telescope. The two millisecond pulsars, PSRs J2017+0603 and J2302+4442, have rotational periods of 2.896 and 5.192 ms and are both in binary systems with low-eccentricity orbits and orbital periods of 2.2 and 125.9 days respectively, suggesting long recycling processes. Gamma-ray pulsations were subsequently detected for both objects, indicating that they power the associated \emph{Fermi} sources in which they were found. The gamma-ray light curves and spectral properties are similar to those of previously-detected gamma-ray millisecond pulsars. Detailed modeling of the observed radio and gamma-ray light curves shows that the gamma-ray emission seems to originate at high altitudes in their magnetospheres. Additionally, X-ray observations revealed the presence of an X-ray source at the position of PSR J2302+4442, consistent with thermal emission from a neutron star. These discoveries along with the numerous detections of radio-loud millisecond pulsars in gamma rays suggest that many \emph{Fermi} sources with no known counterpart could be unknown millisecond pulsars.

preprint2011arXiv

Evidence for gravitational quadrupole moment variations in the companion of PSR J2051-0827

We have conducted radio timing observations of the eclipsing millisecond binary pulsar J2051-0827 with the European Pulsar Timing Array network of telescopes and the Parkes radio telescope, spanning over 13 years. The increased data span allows significant measurements of the orbital eccentricity, e = (6.2 {\pm} 1.3) {\times} 10^{-5} and composite proper motion, μ_t = 7.3 {\pm} 0.4 mas/yr. Our timing observations have revealed secular variations of the projected semi-major axis of the pulsar orbit which are much more extreme than those previously published; and of the orbital period of the system. Investigation of the physical mechanisms producing such variations confirm that the variations of the semi-major axis are most probably caused by classical spin-orbit coupling in the binary system, while the variations in orbital period are most likely caused by tidal dissipation leading to changes in the gravitational quadrupole moment of the companion.

preprint2011arXiv

Gravitational wave astronomy of single sources with a pulsar timing array

Abbreviated: We investigate the potential of detecting the gravitational wave from individual binary black hole systems using pulsar timing arrays (PTAs) and calculate the accuracy for determining the GW properties. This is done in a consistent analysis, which at the same time accounts for the measurement of the pulsar distances via the timing parallax. We find that, at low redshift, a PTA is able to detect the nano-Hertz GW from super massive black hole binary systems with masses of $\sim10^8 - 10^{10}\,M_{\sun}$ less than $\sim10^5$\,years before the final merger, and those with less than $\sim10^3 - 10^4$ years before merger may allow us to detect the evolution of binaries. We derive an analytical expression to describe the accuracy of a pulsar distance measurement via timing parallax. We consider five years of bi-weekly observations at a precision of 15\,ns for close-by ($\sim 0.5 - 1$\,kpc) pulsars. Timing twenty pulsars would allow us to detect a GW source with an amplitude larger than $5\times 10^{-17}$. We calculate the corresponding GW and binary orbital parameters and their measurement precision. The accuracy of measuring the binary orbital inclination angle, the sky position, and the GW frequency are calculated as functions of the GW amplitude. We note that the "pulsar term", which is commonly regarded as noise, is essential for obtaining an accurate measurement for the GW source location. We also show that utilizing the information encoded in the GW signal passing the Earth also increases the accuracy of pulsar distance measurements. If the gravitational wave is strong enough, one can achieve sub-parsec distance measurements for nearby pulsars with distance less than $\sim 0.5 - 1$\,kpc.

preprint2011arXiv

High-Precision Timing of 5 Millisecond Pulsars: Space Velocities, Binary Evolution and Equivalence Principles

We present high-precision timing of five millisecond pulsars (MSPs) carried out for more than seven years; four pulsars are in binary systems and one is isolated. We are able to measure the pulsars' proper motions and derive an estimate for their space velocities. The measured two-dimensional velocities are in the range 70-210 km/s, consistent with those measured for other MSPs. We also use all the available proper motion information for isolated and binary MSPs to update the known velocity distribution for these populations. As found by earlier works, we find that the velocity distribution of binary and isolated MSPs are indistinguishable with the current data. Four of the pulsars in our observing program are highly recycled with low-mass white dwarf companions and we are able to derive accurate binary parameters for these systems. For three of these binary systems we are able to place initial constraints on the pulsar masses with best-fit values in the range 1.0-1.6 M_sun. The implications of the results presented here to our understanding of binary pulsar evolution are discussed. The updated parameters for the binary systems studied here, together with recently discovered similar systems, allowed us to update previous limits on the the violation of the strong equivalence principle through the parameter |Delta| to 4.6x10^-3 (95% confidence) and the violation of Lorentz-invariance/momentum-conservation through the parameter |hat{alpha}_3| to 5.5x10^-20 (95% confidence).

preprint2011arXiv

Long term monitoring of mode switching for PSR B0329+54

The mode switching phenomenon of PSR B0329+54 is investigated based on the long-term monitoring from September 2003 to April 2009 made with the Urumqi 25m radio telescope at 1540 MHz. At that frequency, the change of relative intensity between the leading and trailing components is the predominant feature of mode switching. The intensity ratios between the leading and trailing components are measured for the individual profiles averaged over a few minutes. It is found that the ratios follow normal distributions, where the abnormal mode has a wider typical width than the normal mode, indicating that the abnormal mode is less stable than the normal mode. Our data show that 84.9% of the time for PSR B0329+54 was in the normal mode and 15.1% was in the abnormal mode. From the two passages of eight-day quasi-continuous observations in 2004, and supplemented by the daily data observed with 15 m telescope at 610 MHz at Jodrell Bank Observatory, the intrinsic distributions of mode timescales are constrained with the Bayesian inference method. It is found that the gamma distribution with the shape parameter slightly smaller than 1 is favored over the normal, lognormal and Pareto distributions. The optimal scale parameters of the gamma distribution is 31.5 minutes for the abnormal mode and 154 minutes for the normal mode. The shape parameters have very similar values, i.e. 0.75^{+0.22}_{-0.17} for the normal mode and 0.84^{+0.28}_{-0.22} for the abnormal mode, indicating the physical mechanisms in both modes may be the same. No long-term modulation of the relative intensity ratios was found for both the modes, suggesting that the mode switching was stable. The intrinsic timescale distributions, for the first time constrained for this pulsar, provide valuable information to understand the physics of mode switching.

preprint2011arXiv

Observations of Energetic High Magnetic Field Pulsars with the Fermi Large Area Telescope

We report the detection of gamma-ray pulsations from the high-magnetic-field rotation-powered pulsar PSR J1119-6127 using data from the Fermi Large Area Telescope. The gamma-ray light curve of PSR J1119-6127 shows a single, wide peak offset from the radio peak by 0.43 pm 0.02 in phase. Spectral analysis suggests a power law of index 1.0 pm 0.3 with an energy cut-off at 0.8 pm 0.2 GeV. The first uncertainty is statistical and the second is systematic. We discuss the emission models of PSR J1119-6127 and demonstrate that despite the object's high surface magnetic field---near that of magnetars---the field strength and structure in the gamma-ray emitting zone are apparently similar to those of typical young pulsars. Additionally, we present upper limits on the \gam-ray pulsed emission for the magnetically active PSR J1846-0258 in the supernova remnant Kesteven 75 and two other energetic high-B pulsars, PSRs J1718-3718 and J1734-3333. We explore possible explanations for the non-detection of these three objects, including peculiarities in their emission geometry.

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

Precise Gamma-Ray Timing and Radio Observations of 17 Fermi Gamma-Ray Pulsars

We present precise phase-connected pulse timing solutions for 16 gamma-ray-selected pulsars recently discovered using the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope plus one very faint radio pulsar (PSR J1124-5916) that is more effectively timed with the LAT. We describe the analysis techniques including a maximum likelihood method for determining pulse times of arrival from unbinned photon data. A major result of this work is improved position determinations, which are crucial for multi-wavelength follow up. For most of the pulsars, we overlay the timing localizations on X-ray images from Swift and describe the status of X-ray counterpart associations. We report glitches measured in PSRs J0007+7303, J1124-5916, and J1813-1246. We analyze a new 20 ks Chandra ACIS observation of PSR J0633+0632 that reveals an arcminute-scale X-ray nebula extending to the south of the pulsar. We were also able to precisely localize the X-ray point source counterpart to the pulsar and find a spectrum that can be described by an absorbed blackbody or neutron star atmosphere with a hard powerlaw component. Another Chandra ACIS image of PSR J1732-3131 reveals a faint X-ray point source at a location consistent with the timing position of the pulsar. Finally, we present a compilation of new and archival searches for radio pulsations from each of the gamma-ray-selected pulsars as well as a new Parkes radio observation of PSR J1124-5916 to establish the gamma-ray to radio phase offset.

preprint2011arXiv

Present and possible future implications for mSUGRA of the non-discovery of SUSY at the LHC

Both ATLAS and CMS have published results of SUSY searches putting limits on SUSY parameters and masses. A non-discovery of SUSY in the next two years would push these limits further. On the other hand, precision data of low energy measurements and the dark matter relic density favor a light scale of supersymmetry. Therefore we investigate if supersymmetry -- more specifically the highly constraint model mSUGRA -- does at all agree with precision data and LHC exclusions at the same time, and whether the first two years of LHC will be capable of excluding models of supersymmetry. We consider the current non observation of supersymmetry with 35 pb-1 as well as the possible non observation with 1,2 and 7 fb-1 in a global fit using the framework Fittino.

preprint2011arXiv

Profile shape stability and phase jitter analyses of millisecond pulsars

Millisecond pulsars (MSPs) have been studied in detail since their discovery in 1982. The integrated pulse profiles of MSPs appear to be stable, which enables precision monitoring of the pulse times of arrival (TOAs). However, for individual pulses the shape and arrival phase can vary dramatically, which is known as pulse jitter. In this paper, we investigate the stability of integrated pulse profiles for 5 MSPs, and estimate the amount of jitter for PSR J0437-4715. We do not detect intrinsic profile shape variation based on integration times from ~10 to ~100 s with the provided instrumental sensitivity. For PSR J0437-4715 we calculate the jitter parameter to be f_J=0.067+-0.002, and demonstrate that the result is not significantly affected by instrumental TOA uncertainties. Jitter noise is also found to be independent of observing frequency and bandwidth around 1.4 GHz on frequency scales of <100 MHz, which supports the idea that pulses within narrow frequency scale are equally jittered. In addition, we point out that pulse jitter would limit TOA calculation for the timing observations with future telescopes like the Square Kilometre Array and the Five hundred metre Aperture Spherical Telescope. A quantitative understanding of pulse profile stability and the contribution of jitter would enable improved TOA calculations, which are essential for the ongoing endeavours in pulsar timing, such as the detection of the stochastic gravitational wave background.

preprint2011arXiv

Prospects for accurate distance measurements of pulsars with the SKA: enabling fundamental physics

Parallax measurements of pulsars allow for accurate measurements of the interstellar electron density and contribute to accurate tests of general relativity using binary systems. The Square Kilometre Array (SKA) will be an ideal instrument for measuring the parallax of pulsars, because it has a very high sensitivity, as well as baselines extending up to several thousands of kilometres. We performed simulations to estimate the number of pulsars for which the parallax can be measured with the SKA and the distance to which a parallax can be measured. We compare two different methods. The first method measures the parallax directly by utilising the long baselines of the SKA to form high angular resolution images. The second method uses the arrival times of the radio signals of pulsars to fit a transformation between time coordinates in the terrestrial frame and the comoving pulsar frame directly yielding the parallax. We find that with the first method a parallax with an accuracy of 20% or less can be measured up to a maximum distance of 13 kpc, which would include 9,000 pulsars. By timing pulsars with the most stable arrival times for the radio emission, parallaxes can be measured for about 3,600 millisecond pulsars up to a distance of 9 kpc with an accuracy of 20%.

preprint2011arXiv

Prospects for High-Precision Pulsar Timing

Timing pulses of pulsars has proved to be a most powerful technique useful to a host of research areas in astronomy and physics. Importantly, the precision of this timing is not only affected by radiometer noise, but also by intrinsic pulse shape changes, interstellar medium (ISM) evolution, instrumental distortions, etc. In this paper we review the known causes of pulse shape variations and assess their effect on the precision and accuracy of a single measurement of pulse arrival time with current instrumentation. Throughout this analysis we use the brightest and most precisely timed millisecond pulsar (MSP), PSR J0437-4715, as a case study, and develop a set of diagnostic tools to evaluate profile stability in timing observations. We conclude that most causes of distortion can be either corrected by state-of-the-art techniques or taken into account in the estimation of time-of-arrival (TOA) uncertainties. The advent of a new generation of radio telescopes (e.g. the Square Kilometre Array, SKA), and their increase in collecting area has sparked speculation about the timing precision achievable through increases in gain. Based on our analysis of current data, we predict that for normal-brightness MSPs a TOA precision of between 80 and 230 ns can be achieved at 1.4 GHz with 10-minute integrations by the SKA. The actual rms timing residuals for each pulsar will be approximately at the same level only if all the other influences on timing precision (e.g. ISM, spin noise) are either corrected, modelled, or negligible.

preprint2011arXiv

Prospects for Probing the Spacetime of Sgr A* with Pulsars

The discovery of radio pulsars in compact orbits around Sgr A* would allow an unprecedented and detailed investigation of the spacetime of the supermassive black hole. This paper shows that pulsar timing, including that of a single pulsar, has the potential to provide novel tests of general relativity, in particular its cosmic censorship conjecture and no-hair theorem for rotating black holes. These experiments can be performed by timing observations with 100 micro-second precision, achievable with the Square Kilometre Array for a normal pulsar at frequency above 15 GHz. Based on the standard pulsar timing technique, we develop a method that allows the determination of the mass, spin, and quadrupole moment of Sgr A*, and provides a consistent covariance analysis of the measurement errors. Furthermore, we test this method in detailed mock data simulations. It seems likely that only for orbital periods below ~0.3 yr is there the possibility of having negligible external perturbations. For such orbits we expect a ~10^-3 test of the frame dragging and a ~10^-2 test of the no-hair theorem within 5 years, if Sgr A* is spinning rapidly. Our method is also capable of identifying perturbations caused by distributed mass around Sgr A*, thus providing high confidence in these gravity tests. Our analysis is not affected by uncertainties in our knowledge of the distance to the Galactic center, R0. A combination of pulsar timing with the astrometric results of stellar orbits would greatly improve the measurement precision of R0.

preprint2011arXiv

PSR J2030+3641: radio discovery and gamma-ray study of a middle-aged pulsar in the now identified Fermi-LAT source 1FGL J2030.0+3641

In a radio search with the Green Bank Telescope of three unidentified low Galactic latitude Fermi-LAT sources, we have discovered the middle-aged pulsar J2030+3641, associated with 1FGL J2030.0+3641 (2FGL J2030.0+3640). Following the detection of gamma-ray pulsations using a radio ephemeris, we have obtained a phase-coherent timing solution based on gamma-ray and radio pulse arrival times that spans the entire Fermi mission. With a rotation period of 0.2 s, spin-down luminosity of 3e34 erg/s, and characteristic age of 0.5 Myr, PSR J2030+3641 is a middle-aged neutron star with spin parameters similar to those of the exceedingly gamma-ray-bright and radio-undetected Geminga. Its gamma-ray flux is 1% that of Geminga, primarily because of its much larger distance, as suggested by the large integrated column density of free electrons, DM=246 pc/cc. We fit the gamma-ray light curve, along with limited radio polarimetric constraints, to four geometrical models of magnetospheric emission, and while none of the fits have high significance some are encouraging and suggest that further refinements of these models may be worthwhile. We argue that not many more non-millisecond radio pulsars may be detected along the Galactic plane that are responsible for LAT sources, but that modified methods to search for gamma-ray pulsations should be productive -- PSR J2030+3641 would have been found blindly in gamma rays if only >0.8 GeV photons had been considered, owing to its relatively flat spectrum and location in a region of high soft background.

preprint2011arXiv

Pulsed Gamma Rays from the Original Millisecond and Black Widow Pulsars: a case for Caustic Radio Emission?

We report the detection of pulsed gamma-ray emission from the fast millisecond pulsars (MSPs) B1937+21 (also known as J1939+2134) and B1957+20 (J1959+2048) using 18 months of survey data recorded by the \emph{Fermi} Large Area Telescope (LAT) and timing solutions based on radio observations conducted at the Westerbork and Nançay radio telescopes. In addition, we analyzed archival \emph{RXTE} and \emph{XMM-Newton} X-ray data for the two MSPs, confirming the X-ray emission properties of PSR B1937+21 and finding evidence ($\sim 4σ$) for pulsed emission from PSR B1957+20 for the first time. In both cases the gamma-ray emission profile is characterized by two peaks separated by half a rotation and are in close alignment with components observed in radio and X-rays. These two pulsars join PSRs J0034-0534 and J2214+3000 to form an emerging class of gamma-ray MSPs with phase-aligned peaks in different energy bands. The modeling of the radio and gamma-ray emission profiles suggests co-located emission regions in the outer magnetosphere.

preprint2011arXiv

RRATs: New Discoveries, Timing Solutions & Musings

We describe observations of Rotating RAdio Transients (RRATs) that were discovered in a re-analysis of the Parkes Multi-beam Pulsar Survey (PMPS). The sources have now been monitored for sufficiently long to obtain seven new coherent timing solutions, to make a total of 14 now known. Furthermore we announce the discovery of 7 new transient sources, one of which may be extragalactic in origin (with $z\sim0.1$) and would then be a second example of the so-called `Lorimer burst'. The timing solutions allow us to infer neutron star characteristics such as energy-loss rate, magnetic field strength and evolutionary timescales, as well as facilitating multi-wavelength followup by providing accurate astrometry. All of this enables us to consider the question of whether or not RRATs are in any way special, i.e. a distinct and separate population of neutron stars, as has been previously suggested. We see no reason to consider 'RRAT' as anything other than a detection label, the subject of a selection effect in the parameter space searched. However, single-pulse searches can be utilised to great effect to identify pulsars difficult, or impossible, to find by other means, in particular those with long-periods (half of the PMPS RRATs have periods greater than 4 seconds), high-magnetic field strengths ($B\gtrsim 10^{13}$ G) and pulsars approaching the 'death valley'. The detailed nulling properties of such pulsars are unknown but the mounting evidence suggests a broad range of behaviour in the pulsar population. The group of RRATs fit in to the picture where pulsar magnetospheres switch between stable configurations.

preprint2011arXiv

The High Time Resolution Universe Pulsar Survey IV: Discovery and polarimetry of millisecond pulsars

We present the discovery of six millisecond pulsars (MSPs) in the High Time Resolution Universe (HTRU) survey for pulsars and fast transients carried out with the Parkes radio telescope. All six are in binary systems with approximately circular orbits and are likely to have white dwarf companions. PSR J1017-7156 has a high flux density and a narrow pulse width, making it ideal for precision timing experiments. PSRs J1446-4701 and J1125-5825 are coincident with gamma-ray sources, and folding the high-energy photons with the radio timing ephemeris shows evidence of pulsed gamma-ray emission. PSR J1502-6752 has a spin period of 26.7 ms, and its low period derivative implies that it is a recycled pulsar. The orbital parameters indicate it has a very low mass function, and therefore a companion mass much lower than usually expected for such a mildly recycled pulsar. In addition we present polarisation profiles for all 12 MSPs discovered in the HTRU survey to date. Similar to previous observations of MSPs, we find that many have large widths and a wide range of linear and circular polarisation fractions. Their polarisation profiles can be highly complex, and although the observed position angles often do not obey the rotating vector model, we present several examples of those that do. We speculate that the emission heights of MSPs are a substantial fraction of the light cylinder radius in order to explain broad emission profiles, which then naturally leads to a large number of cases where emission from both poles is observed.

preprint2011arXiv

The High Time Resolution Universe Survey - III. Single-pulse searches and preliminary analysis

We present the search methods and initial results for transient radio signals in the High Time Resolution Universe (HTRU) Survey. The HTRU survey's single-pulse search, the software designed to perform the search, and a determination of the HTRU survey's sensitivity to single pulses are described. Initial processing of a small fraction of the survey has produced 11 discoveries, all of which are sparsely-emitting neutron stars, as well as provided confirmation of two previously unconfirmed neutron stars. Most of the newly discovered objects lie in regions surveyed previously, indicating both the improved sensitivity of the HTRU survey observing system and the dynamic nature of the radio sky. The cycles of active and null states in nulling pulsars, rotating radio transients (RRATs), and long-term intermittent pulsars are explored in the context of determining the relationship between these populations, and of the sensitivity of a search to the various radio-intermittent neutron star populations. This analysis supports the case that many RRATs are in fact high-null-fraction pulsars (i. e. with null fraction >~0.95), and indicates that intermittent pulsars appear distinct from nulling pulsars in their activity cycle timescales. We find that in the measured population, there is a deficit of pulsars with typical emission timescales greater than ~300 s that is not readily explained by selection effects. The HTRU low-latitude survey will be capable of addressing whether this deficit is physical. We predict that the HTRU survey will explore pulsars with a broad range of nulling fractions (up to and beyond 0.999), and at its completion is likely to increase the currently known RRATs by a factor of more than two.

preprint2011arXiv

The High Time Resolution Universe Survey II: Discovery of 5 Millisecond Pulsars

We present the discovery of 5 millisecond pulsars found in the mid-Galactic latitude portion of the High Time Resolution Universe (HTRU) Survey. The pulsars have rotational periods from ~2.3 to ~7.5 ms, and all are in binary systems with orbital periods ranging from ~0.3 to ~150 d. In four of these systems, the most likely companion is a white dwarf, with minimum masses of ~0.2 Solar Masses. The other pulsar, J1731-1847, has a very low mass companion and exhibits eclipses, and is thus a member of the "black widow" class of pulsar binaries. These eclipses have been observed in bands centred near frequencies of 700, 1400 and 3000 MHz, from which measurements have been made of the electron density in the eclipse region. These measurements have been used to examine some possible eclipse mechanisms. The eclipse and other properties of this source are used to perform a comparison with the other known eclipsing and "black widow" pulsars. These new discoveries occupy a short-period and high-dispersion measure (DM) region of parameter space, which we demonstrate is a direct consequence of the high time and frequency resolution of the HTRU survey. The large implied distances to our new discoveries makes observation of their companions unlikely with both current optical telescopes and the Fermi Gamma-ray Space Telescope. The extremely circular orbits make any advance of periastron measurements highly unlikely. No relativistic Shapiro delays are obvious in any of the systems, although the low flux densities would make their detection difficult unless the orbits were fortuitously edge-on.

preprint2011arXiv

Transformation of a Star into a Planet in a Millisecond Pulsar Binary

Millisecond pulsars are thought to be neutron stars that have been spun-up by accretion of matter from a binary companion. Although most are in binary systems, some 30% are solitary, and their origin is therefore mysterious. PSR J1719-1438, a 5.7 ms pulsar, was detected in a recent survey with the Parkes 64m radio telescope. We show that it is in a binary system with an orbital period of 2.2 h. Its companion's mass is near that of Jupiter, but its minimum density of 23 g cm$^{-3}$ suggests that it may be an ultra-low mass carbon white dwarf. This system may thus have once been an Ultra Compact Low-Mass X-ray Binary, where the companion narrowly avoided complete destruction.

preprint2011arXiv

Very High Angular Resolution Science with the Square Kilometre Array

Preliminary specifications for the Square Kilometre Array (SKA) call for 25% of the total collecting area of the dish array to be located at distances greater than 180 km from the core, with a maximum baseline of at least 3000 km. The array will provide angular resolution ~ 40 - 2 mas at 0.5 - 10 GHz with image sensitivity reaching < 50 nJy/beam in an 8 hour integration with 500 MHz bandwidth. Given these specifications, the high angular resolution component of the SKA will be capable of detecting brightness temperatures < 200 K with milliarcsecond-scale angular resolution. The aim of this article is to bring together in one place a discussion of the broad range of new and important high angular resolution science that will be enabled by the SKA, and in doing so, address the merits of long baselines as part of the SKA. We highlight the fact that high angular resolution requiring baselines greater than 1000 km provides a rich science case with projects from many areas of astrophysics, including important contributions to key SKA science.

preprint2010arXiv

A 6.5-GHz Multibeam Pulsar Survey

A survey of the Galactic plane in the region $-60\degree \leq l \leq 30\degree$, $|b| \leq 0.25\degree$ was carried out using the seven-beam Parkes Methanol Multibeam (MMB) receiver, which operates at a frequency of 6.5 GHz. Three pulsars were discovered, and 16 previously known pulsars detected. In this paper we present two previously-unpublished discoveries, both with extremely high dispersion measures, one of which is very close, in angular distance, to the Galactic centre. The survey data also contain the first known detection, at radio frequencies, of the radio magnetar PSR J1550-5418. Our survey observation was made 46 days prior to that previously published and places constraints on the beginning of pulsed radio emission from the source. The detection of only three previously undiscovered pulsars argues that there are few pulsars in the direction of the inner Galaxy whose flux density spectrum is governed by a flat power law. However, these pulsars would be likely to remain undetected at lower frequencies due to the large amount of scatter broadening which affects pulsars with high values of dispersion measure. Surveys with future telescopes at high observing frequencies will, therefore, play an important role in the discovery of pulsars at the Galactic centre. By simulating pulsar surveys of the Galaxy with Phase 1 SKA at frequencies of 1.4 GHz and 10 GHz, we find that high-frequency observations are the only way to discover and observe the Galactic-centre pulsar population.

preprint2010arXiv

A white dwarf companion to the relativistic pulsar PSR J1141-6545

Pulsars with compact companions in close eccentric orbits are unique laboratories for testing general relativity and alternative theories of gravity. Moreover, they are excellent targets for future gravitational wave experiments like LISA and they are also highly important for understanding the equation of state of super-dense matter and the evolution of massive binaries. Here we report on optical observations of the 1.02 M_sun companion to the pulsar PSR J1141-6545. We detect an optical counterpart with apparent magnitudes V=25.08(11) and R=24.38(14), consistent with the timing position of the pulsar. We demonstrate that our results are in agreement with a white dwarf companion. However the latter is redder than expected and the inferred values are not consistent with the theoretical cooling tracks, preventing us from deriving the exact age. Our results confirm the importance of the PSR J1141-6545 system for gravitational experiments.

preprint2010arXiv

Discovery of 59ms Pulsations from 1RXS J141256.0+792204 (Calvera)

We report on a multi-wavelength study of the compact object candidate 1RXS J141256.0+792204 (Calvera). Calvera was observed in the X-rays with XMM/EPIC twice for a total exposure time of ~50 ks. The source spectrum is thermal and well reproduced by a two component model composed of either two hydrogen atmosphere models, or two blackbodies (kT_1~ 55/150 eV, kT_2~ 80/250 eV, respectively, as measured at infinity). Evidence was found for an absorption feature at ~0.65 keV; no power-law high-energy tail is statistically required. Using pn and MOS data we discovered pulsations in the X-ray emission at a period P=59.2 ms. The detection is highly significant (> 11 sigma), and unambiguously confirms the neutron star nature of Calvera. The pulse profile is nearly sinusoidal, with a pulsed fraction of ~18%. We looked for the timing signature of Calvera in the Fermi Large Area Telescope (LAT) database and found a significant (~5 sigma) pulsed signal at a period coincident with the X-ray value. The gamma-ray timing analysis yielded a tight upper limit on the period derivative, dP/dt < 5E-18 s/s (dE_rot/dt <1E33 erg/s, B<5E10 G for magneto- dipolar spin-down). Radio searches at 1.36 GHz with the 100-m Effelsberg radio telescope yielded negative results, with a deep upper limit on the pulsed flux of 0.05 mJy. Diffuse, soft (< 1 keV) X-ray emission about 13' west of the Calvera position is present both in our pointed observations and in archive ROSAT all-sky survey images, but is unlikely associated with the X-ray pulsar. Its spectrum is compatible with an old supernova remnant (SNR); no evidence for diffuse emission in the radio and optical bands was found. The most likely interpretations are that Calvera is either a central compact object escaped from a SNR or a mildly recycled pulsar; in both cases the source would be the first ever member of the class detected at gamma-ray energies.

preprint2010arXiv

Long-term timing of four millisecond pulsars

We have timed four millisecond pulses, PSRs J1721-2457, J1745-0952, J1810-2005, and J1918-0642, for up to a total of 10.5 years each using multiple telescopes in the European Pulsar Timing Array network: the Westerbork Synthesis Radio Telescope in The Netherlands, the Nancay Radio Telescope in France and the Lovell telescope at Jodrell Bank in the UK. The long time span has enabled us to measure the proper motions of J1745-0952 and J1918-0642, indicating that they have transverse velocities of 200(50) and 54(7) km/s respectively. We have obtained upper limits on the proper motion of J1721-2457 and J1810-2005, which imply that they have transverse velocities less than 140 and 400 km/s respectively. In all cases, the velocities lie in the range typical of millisecond pulsars. We present pulse profiles for each pulsar taken from observations at multiple frequencies in the range of 350 to 2600 MHz, and show that J1810-2005 shows significant profile evolution in this range. Using our multi-frequency observations, we measured the spectral indices for all four pulsars, and for J1810-2005 it appears to be very flat. The flux density of J1918-0642 shows extensive modulation which we attribute to the combined effects of refractive and diffractive scintillation. We discuss the possible use of including J1721-2457 or J1918-0642 in a pulsar timing array, and find that J1918-0642 will be useful to include when the timing precision of this pulsar is improved over the next few years. We have searched archival optical observations to detect companions of the binary pulsars, but none were detected. However, we provide lower limits on the masses of the white dwarf companions of PSRs J1745-0952 and J1918-0642.

preprint2010arXiv

Measuring the mass of solar system planets using pulsar timing

High-precision pulsar timing relies on a solar-system ephemeris in order to convert times of arrival (TOAs) of pulses measured at an observatory to the solar system barycenter. Any error in the conversion to the barycentric TOAs leads to a systematic variation in the observed timing residuals; specifically, an incorrect planetary mass leads to a predominantly sinusoidal variation having a period and phase associated with the planet's orbital motion about the Sun. By using an array of pulsars (PSRs J0437-4715, J1744-1134, J1857+0943, J1909-3744), the masses of the planetary systems from Mercury to Saturn have been determined. These masses are consistent with the best-known masses determined by spacecraft observations, with the mass of the Jovian system, 9.547921(2)E-4 Msun, being significantly more accurate than the mass determined from the Pioneer and Voyager spacecraft, and consistent with but less accurate than the value from the Galileo spacecraft. While spacecraft are likely to produce the most accurate measurements for individual solar system bodies, the pulsar technique is sensitive to planetary system masses and has the potential to provide the most accurate values of these masses for some planets.

preprint2010arXiv

Observations and Modelling of Relativistic Spin Precession in PSR J1141-6545

Observations of the binary pulsar PSR J1141-6545 using the Parkes radio telescope over 9.3 years show clear time-variations in pulse width, shape and polarization. We interpret these variations in terms of relativistic precession of the pulsar spin axis about the total angular momentum vector of the system. Over the nine years, the pulse width at the 50% level has changed by more than a factor of three. Large variations have also been observed in the 1400-MHz mean flux density. The pulse polarization has been monitored since 2004 April using digital filterbank systems and also shows large and systematic variations in both linear and circular polarization. Position angle variations, both across the pulse profile and over the data span, are complex, with major differences between the central and outer parts of the pulse profile. Modelling of the observed position angle variations by relativistic precession of the pulsar spin axis shows that the spin-orbit misalignment angle is about 110 deg and that the precessional phase has passed through 180 deg during the course of our observations. At the start of our observations, the line-of-sight impact parameter was about 4 deg in magnitude and it reached a minimum very close to 0 deg around early 2007, consistent with the observed pulse width variations. We have therefore mapped approximately one half of the emission beam, showing that it is very asymmetric with respect to the magnetic axis. The derived precessional parameters imply that the pre-supernova star had a mass of about 2 Msun and that the supernova recoil kick velocity was relatively small. With the reversal in the rate of change of the impact parameter, we predict that over the next decade we will see a reversed "replay" of the variations observed in the past decade.

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

Pulsars with the Australian Square Kilometre Array Pathfinder

The Australian Square Kilometre Array Pathfinder (ASKAP) is a 36-element array with a 30-square-degree field of view being built at the proposed SKA site in Western Australia. We are conducting a Design Study for pulsar observations with ASKAP, planning both timing and search observations. We provide an overview of the ASKAP telescope and an update on pulsar-related progress.

preprint2010arXiv

Radio and Gamma-Ray Constraints on the Emission Geometry and Birthplace of PSR J2043+2740

We report on the first year of Fermi gamma-ray observations of pulsed high-energy emission from the old PSR J2043+2740. The study of the gamma-ray efficiency of such old pulsars gives us an insight into the evolution of pulsars' ability to emit in gammma rays as they age. The gamma-ray lightcurve of this pulsar above 0.1 GeV is clearly defined by two sharp peaks, 0.353+/-0.035 periods apart. We have combined the gamma-ray profile characteristics of PSR J2043+2740 with the geometrical properties of the pulsar's radio emission, derived from radio polarization data, and constrained the pulsar-beam geometry in the framework of a Two Pole Caustic and an Outer Gap model. The ranges of magnetic inclination and viewing angle were determined to be {alpha,zeta}~{52-57,61-68} for the Two Pole Caustic model, and {alpha,zeta}~{62-73,74-81} and {alpha,zeta}~{72-83,60-75} for the Outer Gap model. Based on this geometry, we assess possible birth locations for this pulsar and derive a likely proper motion, sufficiently high to be measurable with VLBI. At a characteristic age of 1.2 Myr, PSR J2043+2740 is the third oldest of all discovered, non-recycled, gamma-ray pulsars: it is twice as old as the next oldest, PSR J0357+32, and younger only than the recently discovered PSR J1836+5925 and PSR J2055+25, both of which are at least 5 and 10 times less energetic, respectively.

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

Selection of radio pulsar candidates using artificial neural networks

Radio pulsar surveys are producing many more pulsar candidates than can be inspected by human experts in a practical length of time. Here we present a technique to automatically identify credible pulsar candidates from pulsar surveys using an artificial neural network. The technique has been applied to candidates from a recent re-analysis of the Parkes multi-beam pulsar survey resulting in the discovery of a previously unidentified pulsar.

preprint2010arXiv

The binary companion of PSR J1740-3052

We report on the identification of a near-infrared counterpart to the massive (>11 Msun) binary companion of pulsar J1740-3052. An accurate celestial position of PSR J1740-3052 is determined from interferometric radio observations. Adaptive optics corrected near-infrared imaging observations show a counterpart at the interferometric position of the pulsar. The counterpart has Ks=15.87+-0.10 and J-Ks>0.83. Based on distance and absorption estimates from models of the Galactic electron and dust distributions these observed magnitudes are consistent with those of a main-sequence star as the binary companion. We argue that this counterpart is the binary companion to PSR J1740-3052 and thus rule out a stellar mass black hole as the pulsar companion.

preprint2010arXiv

The High Time Resolution Universe Pulsar Survey I: System configuration and initial discoveries

We have embarked on a survey for pulsars and fast transients using the 13-beam Multibeam receiver on the Parkes radio telescope. Installation of a digital backend allows us to record 400 MHz of bandwidth for each beam, split into 1024 channels and sampled every 64 us. Limits of the receiver package restrict us to a 340 MHz observing band centred at 1352 MHz. The factor of eight improvement in frequency resolution over previous multibeam surveys allows us to probe deeper into the Galactic plane for short duration signals such as the pulses from millisecond pulsars. We plan to survey the entire southern sky in 42641 pointings, split into low, mid and high Galactic latitude regions, with integration times of 4200, 540 and 270 s respectively. Simulations suggest that we will discover 400 pulsars, of which 75 will be millisecond pulsars. With ~30% of the mid-latitude survey complete, we have re-detected 223 previously known pulsars and discovered 27 pulsars, 5 of which are millisecond pulsars. The newly discovered millisecond pulsars tend to have larger dispersion measures than those discovered in previous surveys, as expected from the improved time and frequency resolution of our instrument.

preprint2010arXiv

Three Millisecond Pulsars in FERMI LAT Unassociated Bright Sources

We searched for radio pulsars in 25 of the non-variable, unassociated sources in the Fermi LAT Bright Source List with the Green Bank Telescope at 820 MHz. We report the discovery of three radio and gamma-ray millisecond pulsars (MSPs) from a high Galactic latitude subset of these sources. All of the pulsars are in binary systems, which would have made them virtually impossible to detect in blind gamma-ray pulsation searches. They seem to be relatively normal, nearby (<=2 kpc) millisecond pulsars. These observations, in combination with the Fermi detection of gamma-rays from other known radio MSPs, imply that most, if not all, radio MSPs are efficient gamma-ray producers. The gamma-ray spectra of the pulsars are power-law in nature with exponential cutoffs at a few GeV, as has been found with most other pulsars. The MSPs have all been detected as X-ray point sources. Their soft X-ray luminosities of ~10^{30-31} erg/s are typical of the rare radio MSPs seen in X-rays.

preprint2010arXiv

XMM-Newton X-ray Observation of the High-Magnetic-Field Radio Pulsar PSR J1734--3333

Using observations made with the XMM-Newton Observatory, we report the probable X-ray detection of the high-magnetic-field radio pulsar PSR J1734-3333. This pulsar has an inferred surface dipole magnetic field of B = 5.2e13 G, just below that of one anomalous X-ray pulsar (AXP). We find that the pulsar has an absorbed 0.5-2.0 keV flux of (5-15)e-15 erg/s/cm^2 and that its X-ray luminosity L_X is well below its spin down luminosity E_dot, with L_X < 0.1E_dot. No pulsations were detected in these data although our derived upper limit is unconstraining. Like most of the other high-B pulsars, PSR J1734-3333 is X-ray faint with no sign of magnetar activity. We collect and tabulate the properties of this and all other known high-B radio pulsars with measured X-ray luminosities or luminosity upper limits and plot L_X versus B for them all.

preprint2009arXiv

A new technique for timing the double pulsar system

In 2004, McLaughlin et al. discovered a phenomenon in the radio emission of PSR J0737-3039B (B) that resembles drifting sub-pulses. The repeat rate of the sub-pulses is equal to the spin frequency of PSR J0737-3039A (A); this led to the suggestion that they are caused by incidence upon B's magnetosphere of electromagnetic radiation from A. Here we describe a geometrical model which predicts the delay of B's sub-pulses relative to A's radio pulses. We show that measuring these delays is equivalent to tracking A's rotation from the point of view of an hypothetical observer located near B. This has three main astrophysical applications: (a) to determine the sense of rotation of A relative to its orbital plane; (b) to estimate where in B's magnetosphere the radio sub-pulses are modulated and (c) to provide an independent estimate of the mass ratio of A and B. The latter might improve existing tests of gravitational theories using this system.

preprint2009arXiv

An analysis of the timing irregularities for 366 pulsars

We provide an analysis of timing irregularities observed for 366 pulsars. Observations were obtained using the 76-m Lovell radio telescope at the Jodrell Bank Observatory over the past 36 years. These data sets have allowed us to carry out the first large-scale analysis of pulsar timing noise over time scales of > 10yr, with multiple observing frequencies and for a large sample of pulsars. Our sample includes both normal and recycled pulsars. The timing residuals for the pulsars with the smallest characteristic ages are shown to be dominated by the recovery from glitch events, whereas the timing irregularities seen for older pulsars are quasi-periodic. We emphasise that previous models that explained timing residuals as a low-frequency noise process are not consistent with observation.

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

Discovery of three new pulsars in a 610 MHz pulsar survey with the GMRT

We report on the discovery of three new pulsars in the first blind survey of the north Galactic plane (45 degrees < l < 135 degrees; |b| < 1 degrees with the Giant Meterwave Radio telescope (GMRT) at an intermediate frequency of 610 MHz. The survey covered 106 square degrees with a sensitivity of roughly 1 mJy to long-period pulsars (pulsars with period longer than 1 s). The three new pulsars have periods of 318, 933, and 1056 ms. Their timing parameters and flux densities, obtained in follow up observations with the Lovell Telescope at Jodrell Bank and the GMRT, are presented. We also report on pulse nulling behaviour in one of the newly discovered pulsars, PSR J2208+5500.

preprint2009arXiv

Further Searches for RRATs in the Parkes Multi-Beam Pulsar Survey

We describe the steps involved in performing searches for sources of transient radio emission such as Rotating Radio Transients (RRATs), and present 10 new transient radio sources discovered in a re-analysis of the Parkes Multi-beam Pulsar Survey. Followup observations of each new source as well as one previously known source are also presented. The new sources suggest that the population of transient radio-emitting neutron stars, and hence the neutron star population in general, may be even larger than initially predicted. We highlight the importance of radio frequency interference excision for single-pulse searches. Also, we discuss some interesting properties of individual sources and consider the difficulties involved in precisely defining a RRAT and determining where they fit in with the other known classes of neutron stars.

preprint2009arXiv

Generic tests of the existence of the gravitational dipole radiation and the variation of the gravitational constant

We present results from the high precision timing analysis of the pulsar-white dwarf (WD) binary PSR J1012+5307 using 15 years of multi-telescope data. Observations were performed regularly by the European Pulsar Timing Array (EPTA) network, consisting of Effelsberg, Jodrell Bank, Westerbork and Nançay. All the timing parameters have been improved from the previously published values, most by an order of magnitude. In addition, a parallax measurement of $π= 1.2(3)$ mas is obtained for the first time for PSR J1012+5307, being consistent with the optical estimation from the WD companion. Combining improved 3D velocity information and models for the Galactic potential the complete evolutionary Galactic path of the system is obtained. A new intrinsic eccentricity upper limit of $e<8.4\times 10^{-7}$ is acquired, one of the smallest calculated for a binary system and a measurement of the variation of the projected semi-major axis also constrains the system's orbital orientation for the first time. It is shown that PSR J1012+5307 is an ideal laboratory for testing alternative theories of gravity. The measurement of the change of the orbital period of the system of $\dot{P}_{b} = 5(1)\times 10^{-14}$ is used to set an upper limit on the dipole gravitational wave emission that is valid for a wide class of alternative theories of gravity. Moreover, it is shown that in combination with other binary pulsars PSR J1012+5307 is an ideal system to provide self-consistent, generic limits, based only on millisecond pulsar data, for the dipole radiation and the variation of the gravitational constant $\dot{G}$.

preprint2009arXiv

Phase-resolved Faraday rotation in pulsars

We have detected significant Rotation Measure variations for 9 bright pulsars, as a function of pulse longitude. An additional sample of 10 pulsars showed a rather constant RM with phase, yet a small degree of RM fluctuation is visible in at least 3 of those cases. In all cases, we have found that the rotation of the polarization position angle across our 1.4 GHz observing band is consistent with the wavelength-squared law of interstellar Faraday Rotation. We provide for the first time convincing evidence that RM variations across the pulse are largely due to interstellar scattering, although we cannot exclude that magnetospheric Faraday Rotation may still have a minor contribution; alternative explanations of this phenomenon, like erroneous de-dispersion and the presence of non-orthogonal polarization modes, are excluded. If the observed, phase-resolved RM variations are common amongst pulsars, then many of the previously measured pulsar RMs may be in error by as much as a few tens of rad m-2.

preprint2009arXiv

Polarisation measurements of five pulsars with interpulses

We present polarisation observations of five pulsars whose profiles exhibit two distinct emission regions separated by close to 180 degrees of longitude. We fitted the position angle of the linear polarisation using the rotating vector model and convincingly show that all the pulsars have the angle between their magnetic and rotation axes close to 90 degrees. The simplest interpretation of the results is that we see `main pulse' emission from one pole and `interpulse' emission from the opposite pole. We have attempted to produce emission maps of the magnetosphere above the polar caps for each pulsar and find that the maps support the view that the emission region in pulsars is complex, even when the profile appears simple. For three pulsars, we can derive emission heights and polar maps which are consistent with emission regions located symmetrically about the magnetic axis and confined to the open field lines. For two pulsars, we find that either the emission arises from `closed' field lines or that the profiles are highly asymmetric with respect to the magnetic axis.

preprint2009arXiv

Precision Astrometry with the Very Long Baseline Array: Parallaxes and Proper Motions for 14 Pulsars

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced observations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 0.13+-0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+1.3-1.1 kpc. We detail our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. Based on our new model-independent estimates for distance and transverse velocity, we investigate the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distribution. We find that young pulsars are moving away from the Galactic plane, as expected, and that age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable exceptions. Given its present trajectory, the pulsar B2045-16 was plausibly born in the open cluster NGC 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the parallax distances by a factor of two, while in others the estimates agree with or are larger than the parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The VLBA astrometric results for the recycled pulsar J1713+0747 are consistent with two independent estimates from pulse timing, enabling a consistency check between the different reference frames.

preprint2009arXiv

Pulsar science with the Five hundred metre Aperture Spherical Telescope

With a collecting area of 70 000 m^2, the Five hundred metre Aperture Spherical Telescope (FAST) will allow for great advances in pulsar astronomy. We have performed simulations to estimate the number of previously unknown pulsars FAST will find with its 19-beam or possibly 100-beam receivers for different survey strategies. With the 19-beam receiver, a total of 5200 previously unknown pulsars could be discovered in the Galactic plane, including about 460 millisecond pulsars (MSPs). Such a survey would take just over 200 days with eight hours survey time per day. We also estimate that, with about 80 six-hour days, a survey of M31 and M33 could yield 50--100 extra-Galactic pulsars. A 19-beam receiver would produce just under 500 MB of data per second and requires about 9 tera-ops to perform the major part of a real time analysis. We also simulate the logistics of high-precision timing of MSPs with FAST. Timing of the 50 brightest MSPs to a signal-to-noise of 500 would take about 24 hours per epoch.

preprint2009arXiv

Radio detection of LAT PSRs J1741-2054 and J2032+4127: no longer just gamma-ray pulsars

Sixteen pulsars have been discovered so far in blind searches of photons collected with the Large Area Telescope on the Fermi Gamma-ray Space Telescope. We here report the discovery of radio pulsations from two of them. PSR J1741-2054, with period P=413ms, was detected in archival Parkes telescope data and subsequently has been detected at the Green Bank Telescope (GBT). Its received flux varies greatly due to interstellar scintillation and it has a very small dispersion measure of DM=4.7pc/cc, implying a distance of ~0.4kpc and possibly the smallest luminosity of any known radio pulsar. At this distance, for isotropic emission, its gamma-ray luminosity above 0.1GeV corresponds to 25% of the spin-down luminosity of dE/dt=9.4e33erg/s. The gamma-ray profile occupies 1/3 of pulse phase and has three closely-spaced peaks with the first peak lagging the radio pulse by delta=0.29P. We have also identified a soft Swift source that is the likely X-ray counterpart. In many respects PSR J1741-2054 resembles the Geminga pulsar. The second source, PSR J2032+4127, was detected at the GBT. It has P=143ms, and its DM=115pc/cc suggests a distance of ~3.6kpc, but we consider it likely that it is located within the Cyg OB2 stellar association at half that distance. The radio emission is nearly 100% linearly polarized, and the main radio peak precedes by delta=0.15P the first of two narrow gamma-ray peaks that are separated by Delta=0.50P. Faint, diffuse X-ray emission in a Chandra image is possibly its pulsar wind nebula. PSR J2032+4127 likely accounts for the EGRET source 3EG J2033+4118, while its pulsar wind is responsible for the formerly unidentified HEGRA source TeV J2032+4130.

preprint2009arXiv

The international pulsar timing array project: using pulsars as a gravitational wave detector

The International Pulsar Timing Array project combines observations of pulsars from both Northern and Southern hemisphere observatories with the main aim of detecting ultra-low frequency (~10^-9 to 10^-8 Hz) gravitational waves. Here we introduce the project, review the methods used to search for gravitational waves emitted from coalescing supermassive binary black-hole systems in the centres of merging galaxies and discuss the status of the project.

preprint2009arXiv

Timing Observations of Rotating Radio Transients

We present radio timing measurements of six rotating radio transient (RRAT) sources discovered in the Parkes Multibeam Pulsar Survey. These provide four new phase-connected timing solutions and two updated ones, making a total of seven of the original 11 reported RRATs now with high-precision rotational and astrometric parameters. Three of these seven RRATs have magnetic fields greater than 10^13 G, with spin-down properties similar to those of the magnetars and X-ray detected isolated neutron stars. Another two of these RRATs have long periods and large characteristic ages, and lie near the `death-line' for radio pulsar emission. The remaining two RRATs with timing solutions have properties typical of the bulk of the pulsar population. The new solutions offer insights into what might be responsible for the unusual emission properties. We demonstrate that the RRATs have significantly longer periods and higher magnetic fields than normal radio pulsars, and find no correlation with other spin-down parameters. These solutions also provide precise positions, which will facilitate follow-up studies at high energies, crucial for relating these sources with other neutron star populations.

preprint2009arXiv

Unusual glitch activity in the RRAT J1819-1458: an exhausted magnetar?

We present an analysis of regular timing observations of the high-magnetic-field Rotating Radio Transient (RRAT) J1819$-$1458 obtained using the 64-m Parkes and 76-m Lovell radio telescopes over the past five years. During this time, the RRAT has suffered two significant glitches with fractional frequency changes of $0.6\times10^{-6}$ and $0.1\times10^{-6}$. Glitches of this magnitude are a phenomenon displayed by both radio pulsars and magnetars. However, the behaviour of J1819$-$1458 following these glitches is quite different to that which follows glitches in other neutron stars, since the glitch activity resulted in a significant long-term net decrease in the slow-down rate. If such glitches occur every 30 years, the spin-down rate, and by inference the magnetic dipole moment, will drop to zero on a timescale of a few thousand years. There are also significant increases in the rate of pulse detection and in the radio pulse energy immediately following the glitches.

preprint2008arXiv

High-Resolution Timing Observations of Spin-Powered Pulsars with the AGILE Gamma-Ray Telescope

AGILE is a small gamma-ray astronomy satellite mission of the Italian Space Agency dedicated to high-energy astrophysics launched in 2007 April. Its 1 microsecond absolute time tagging capability coupled with a good sensitivity in the 30 MeV-30 GeV range, with simultaneous X-ray monitoring in the 18-60 keV band, makes it perfectly suited for the study of gamma-ray pulsars following up on the CGRO/EGRET heritage. In this paper we present the first AGILE timing results on the known gamma-ray pulsars Vela, Crab, Geminga and B 1706-44. The data were collected from 2007 July to 2008 April, exploiting the mission Science Verification Phase, the Instrument Timing Calibration and the early Observing Pointing Program. Thanks to its large field of view, AGILE collected a large number of gamma-ray photons from these pulsars (about 10,000 pulsed counts for Vela) in only few months of observations. The coupling of AGILE timing capabilities, simultaneous radio/X-ray monitoring and new tools aimed at precise photon phasing, exploiting also timing noise correction, unveiled new interesting features at sub-millisecond level in the pulsars' high-energy light-curves.

preprint2008arXiv

On the nature of the intermittent pulsar PSR B1931+24

PSR B1931+24 is the first intermittent radio pulsar discovered to date, characterized by a 0.8 s pulsation which turns on and off quasi-periodically every ~35 days, with a duty cycle of ~10%. We present here X-ray and optical observations of PSR B1931+24 performed with the Chandra X-ray Observatory and Isaac Newton Telescope, respectively. Simultaneous monitoring from the Jodrell Bank Observatory showed that this intermittent pulsar was in the radio-on phase during our observations. We do not find any X-ray or optical counterpart to PSR B1931+24 translating into an upper limit of ~2x10^{31} erg/s on the X-ray luminosity, and of g' > 22.6 on the optical magnitude. If the pulsar is isolated, these limits cannot constrain the dim X-ray and optical emission expected for a pulsar of that age (~1.6 Myr). We discuss the possibility that the quasi-periodic intermittent behavior of PSR B1931+24 is due to the presence of a low mass companion star or gaseous planet, tight with the pulsar in an eccentric orbit. In order to constrain the parameters of this putative binary system we re-analysed the pulsar radio timing residuals and we found that (if indeed hosted in a binary system), PSR B1931+24 should have a very low mass companion and an orbit of low inclination.

preprint2008arXiv

Pulsar searches and timing with the SKA

The Square Kilometre Array (SKA) is a planned multi purpose radio telescope with a collecting area approaching 1 million square metres. One of the key science objectives of the SKA is to provide exquisite strong-field tests of gravitational physics by finding and timing pulsars in extreme binary systems such as a pulsar-black hole binary. To find out how three preliminary SKA configurations will affect a pulsar survey, we have simulated SKA pulsar surveys for each configuration. We estimate that the total number of normal pulsars the SKA will detect, using only the 1-km core and 30 minutes integration time, is around 14000 normal pulsar and 6000 millisecond pulsars. We describe a simple strategy for follow-up timing observations and find that, depending on the configuration, it would take 1-6 days to obtain a single timing point for 14000 pulsars. Obtaining a single timing point for the high-precision timing projects of the SKA, will take less than 14 hours, 2 days, or 3 days, depending on the configuration. The presence of aperture arrays will be of great benefit here. We also study the computational requirements for beam forming and data analysis for a pulsar survey. Beam forming of the full field of view of the single-pixel feed 15-m dishes using the 1-km core of the SKA requires about 2.2*10^15 operations per second. The corresponding data rate from such a pulsar survey is about 4.7*10^11 bytes per second. The required computational power for a deep real time analysis is estimated to be 1.2*10^16 operations per second. For an aperture array or dishes equipped with phased array feeds, the survey can be performed faster, but the computational requirements and data rates will go up.

preprint2007arXiv

Tevatron-for-LHC Report: Top and Electroweak Physics

The top quark and electroweak bosons (W and Z) represent the most massive fundamental particles yet discovered, and as such refer directly to the Standard Model's greatest remaining mystery: the mechanism by which all particles gained mass. This report summarizes the work done within the top-ew group of the Tevatron-for-LHC workshop. It represents a collection of both Tevatron results, and LHC predictions. The hope is that by considering and comparing both machines, the LHC program can be improved and aided by knowledge from the Tevatron, and that particle physics as a whole can be enriched. The report includes measurements of the top quark mass, searches for single top quark production, and physics of the electroweak bosons at hadron colliders.

preprint2006arXiv

Strange baryon resonance production in $\sqrt{s_{NN}} = 200$ GeV $p+p$ and $Au+Au$ collisions

We report the measurements of $Σ(1385)$ and $Λ(1520)$ production in $p+p$ and $Au+Au$ collisions at $\sqrt{s_{NN}} = 200$ GeV from the STAR collaboration. The yields and the $p_{T}$ spectra are presented and discussed in terms of chemical and thermal freeze-out conditions and compared to model predictions. Thermal and microscopic models do not adequately describe the yields of all the resonances produced in central $Au+Au$ collisions. Our results indicate that there may be a time-span between chemical and thermal freeze-out during which elastic hadronic interactions occur.

preprint2005arXiv

A statistical study of 233 pulsar proper motions

We present and analyse a catalogue of 233 pulsars with proper motion measurements. The sample contains a wide variety of pulsars including recycled objects and those associated with globular clusters or supernova remnants. After taking the most precise proper motions for those pulsars for which multiple measurements are available, the majority of the proper motions (58%) are derived from pulsar timing methods, 41% using interferometers and the remaining 1% using optical telescopes. Many of the 1-D and 2-D speeds (referring to speeds measured in one coordinate only and the magnitudes of the transverse velocities respectively) derived from these measurements are somewhat lower than earlier estimates because of the use of the most recent electron density model in determining pulsar distances. The mean 1-D speeds for the normal and recycled pulsars are 152(10) and 54(6) km/s respectively. The corresponding mean 2-D speeds are 246(22) and 87(13) km/s. PSRs B2011+38 and B2224+64 have the highest inferred 2-D speeds of ~1600 km/s. We study the mean speeds for different subsamples and find that, in general, they agree with previous results. Applying a novel deconvolution technique to the sample of 73 pulsars with characteristic ages less than 3 Myr, we find the mean 3-D pulsar birth velocity to be 400(40) km/s. The distribution of velocities is well described by a Maxwellian distribution with 1-D rms sigma=265 km/s. There is no evidence for a bimodal velocity distribution. The proper motions for PSRs B1830-08 and B2334+61 are consistent with their proposed associations with the supernova remnants W41 and G114.3+0.3 respectively.

preprint2005arXiv

Pulsars as Tools for Fundamental Physics and Astrophysics

The sheer number of pulsars discovered by the SKA, in combination with the exceptional timing precision it can provide, will revolutionize the field of pulsar astrophysics. The SKA will provide a complete census of pulsars in both the Galaxy and in Galactic globular clusters that can be used to provide a detailed map of the electron density and magnetic fields, the dynamics of the systems, and their evolutionary histories. This complete census will provide examples of nearly every possible outcome of the evolution of massive stars, including the discovery of very exotic systems such as pulsar black-hole systems and sub-millisecond pulsars, if they exist. These exotic systems will allow unique tests of the strong field limit of relativistic gravity and the equation of state at extreme densities. Masses of pulsars and their binary companions -- planets, white dwarfs, other neutron stars, and black holes -- will be determined to $\sim 1$% for hundreds of objects. With the SKA we can discover and time highly-stable millisecond pulsars that comprise a pulsar-timing array for the detection of low-frequency gravitational waves. The SKA will also provide partial censuses of nearby galaxies through periodicity and single-pulse detections, yielding important information on the intergalactic medium.

preprint2004arXiv

A Double-Pulsar System - A Rare Laboratory for Relativistic Gravity and Plasma Physics

The clock-like properties of pulsars moving in the gravitational fields of their unseen neutron-star companions have allowed unique tests of general relativity and provided evidence for gravitational radiation. We report here the detection of the 2.8-sec pulsar J0737-3039B as the companion to the 23-ms pulsar J0737-3039A in a highly-relativistic double-neutron-star system, allowing unprecedented tests of fundamental gravitational physics. We observe a short eclipse of J0737-3039A by J0737-3039B and orbital modulation of the flux density and pulse shape of J0737-3039B, probably due to the influence of J0737-3039A's energy flux upon its magnetosphere. These effects will allow us to probe magneto-ionic properties of a pulsar magnetosphere.

preprint2004arXiv

The Parkes multibeam pulsar survey: IV. Discovery of 180 pulsars and parameters for 281 previously known pulsars

The Parkes multibeam pulsar survey has led to the discovery of more than 700 pulsars. In this paper, we provide timing solutions, flux densities and pulse profiles for 180 of these new discoveries. Two pulsars, PSRs J1736-2843 and J1847-0130 have rotational periods P > 6s and are therefore among the slowest rotating radio pulsars known. Conversely, with P = 1.8ms, PSR J1843-1113 has the third shortest period of pulsars currently known. This pulsar and PSR J1905+0400 (P = 3.8ms) are both solitary. We also provide orbital parameters for a new binary system, PSR J1420-5625, which has P = 34ms, an orbital period of 40 days and a minimum companion mass of 0.4 solar masses. The 10 degree-wide strip along the Galactic plane that was surveyed is known to contain 264 radio pulsars that were discovered prior to the multibeam pulsar survey. We have redetected almost all of these pulsars and provide new dispersion measure values and flux densities at 20cm for the redetected pulsars.