Many physically motivated extensions to general relativity (GR) predict substantial deviations in the properties of spacetime surrounding massive neutron stars. We report the measurement of a 2.01 ± 0.04 solar mass (M⊙) pulsar in a 2.46-hour orbit with a 0.172 ± 0.003 M⊙ white dwarf. The high pulsar mass and the compact orbit make this system a sensitive laboratory of a previously untested strong-field gravity regime. Thus far, the observed orbital decay agrees with GR, supporting its validity even for the extreme conditions present in the system. The resulting constraints on deviations support the use of GR-based templates for ground-based gravitational wave detectors. Additionally, the system strengthens recent constraints on the properties of dense matter and provides insight to binary stellar astrophysics and pulsar recycling.

Radio pulsars with millisecond spin periods are thought to have been spun up by the transfer of matter and angular momentum from a low-mass companion star during an x-ray-emitting phase. The spin periods of the neutron stars in several such low-mass x-ray binary (LMXB) systems have been shown to be in the millisecond regime, but no radio pulsations have been detected. Here we report on detection and follow-up observations of a nearby radio millisecond pulsar (MSP) in a circular binary orbit with an optically identified companion star. Optical observations indicate that an accretion disk was present in this system within the past decade. Our optical data show no evidence that one exists today, suggesting that the radio MSP has turned on after a recent LMXB phase.

Gravitationally bound three-body systems have been studied for hundreds of years 1, 2 and are common in our Galaxy 3, 4 . They show complex orbital interactions, which can constrain the compositions, masses, and interior structures of the bodies 5 and test theories of gravity 6 , if sufficiently precise measurements are available. A triple system containing a radio pulsar could provide such measurements, but the only previously known such system, B1620−26 7, 8 (with a millisecond pulsar, a white dwarf, and a planetary-mass object in an orbit of several decades), shows only weak interactions. Here we report precision timing and multiwavelength observations of PSR J0337+1715, a millisecond pulsar in a hierarchical triple system with two other stars. Strong gravitational interactions are apparent and provide the masses of the pulsar (1.4378(13) M , where M is the solar mass and the parentheses contain the uncertainty in the final decimal places) and the two white dwarf companions (0.19751(15) M and 0.4101(3) M ), as well as the inclinations of the orbits (both ∼39.2• ). The unexpectedly coplanar and nearly circular orbits indicate a complex and exotic evolutionary past that differs from those of known stellar systems. The gravitational field of the outer white dwarf strongly accelerates the inner binary containing the neutron star, and the system will thus provide an ideal laboratory in which to test the strong equivalence principle of general relativity. Millisecond pulsars (MSPs) are neutron stars that rotate hundreds of times per second and emit radio waves in a lighthouse-like fashion. They are thought to form in binary systems 9 and their rotation rates and orbital properties can be measured with exquisite precision via the unambiguous pulse-counting methodology known as pulsar timing. As part of a large-scale pulsar survey 10, 11 with the Green Bank Telescope (GBT), we have discovered the only known MSP in a stellar triple system. The pulsar has a spin period of 2.73 ms, is relatively bright (∼2 mJy at 1.4 GHz), and has a complex radio pulse profile with multiple narrow components.Though initial timing observations showed a seemingly typical binary MSP system with a 1.6-day circular orbit and a 0.1−0.2 M white dwarf (WD) companion, large timing systematics quickly appeared, strongly suggesting the presence of a third body. There are two other MSPs known to have multiple companions: the famous pulsar B1257+12 which hosts at least 3 low-mass planets 12,13 , and the MSP triple system B1620−26 in globular cluster M4 with a WD inner companion and a roughly Jupiter-mass outer companion 7, 8 . The timing perturbations from J0337+1715 were much too large to be caused by a planetary mass companion.We began an intensive multi-frequency radio timing campaign (Methods) using the GBT, the Arecibo telescope, and the Westerbork Synthesis Radio Telescope (WSRT) to constrain the system's position, orbital parameters, and the nature of the third body. At Arecibo, we achieve median arrival time uncertainties of 0.8 µs...

We present first results from a LOFAR census of non-recycled pulsars. The census includes almost all such pulsars known (194 sources) at declinations Dec > 8• and Galactic latitudes |Gb| > 3• , regardless of their expected flux densities and scattering times. Each pulsar was observed for ≥20 min 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 × 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 lowfrequency 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.

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 through a search for isolated dispersed pulses. All of these objects are Galactic and have measured periods between 0.4 and 4.7 s. One of the new discoveries has a duty cycle of 0.01%, smaller than that of any other radio pulsar. We discuss the impact of selection effects on the detectability and classification of intermittent sources, and compare the efficiencies of periodicity and single-pulse (SP) searches for various pulsar classes. For some cases we find that the apparent intermittency is likely to be caused by off-axis detection or a short time window that selects only a few bright pulses and favors detection with our SP algorithm. In other cases, the intermittency appears to be intrinsic to the source. No transients were found with DMs large enough to require that they originate from sources outside our Galaxy. Accounting for the on-axis gain of the ALFA system, as well as the low gain but large solid-angle coverage of far-out sidelobes, we use the results of the survey so far to place limits on the amplitudes and event rates of transients of arbitrary origin.

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 ∼2 years of observing, 405 cosmic-ray events in the energy range of 10 16 −10 18 eV have been detected in the band from 30−80 MHz. Each of these air showers is registered with up to ∼1000 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.

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-240 MHz, 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.

4Cosmic 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 1 comes from accelerators capable of producing cosmic rays of these energies 2 . 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, X max 3 , or the composition of shower particles reaching the ground 4 .Current measurements 5 suffer from either low precision, and/or a low duty cycle. Radio detection of cosmic rays 6-8 is a rapidly developing technique 9 , suitable for determination of X max 10, 11 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 6, 12 . Here we report radio measurements of X max with a mean precision of 16 g/cm 2 between 10 17 − 10 17.5 eV. Because of the high resolution in X max 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.Observations were made with the Low Frequency Array (LOFAR 13 ), a radio telescope consisting of thousands of crossed dipoles, with built-in air shower detection capability 14 . LOFAR records the radio signals from air showers continuously while running astronomical observations simultaneously. It comprises a scintillator array (LORA), that triggers the readout of buffers, stor-5 ing the full waveforms received by all antennas.We have selected air showers from the period June 2011 -January 2015 with radio pulses in at least 192 antennas. The total uptime was ∼150 days, limited by construction and commissioning of the telescope. Showers that occurred within an hour from lightning activity, or have a polarisation pattern that is indicative of influences from atmospheric electric fields are excluded from the sample 15 .Radio intensity patterns from air showers are asymmetric due to the interference between geomagnetic and charge excess radiation. They can be reproduced from first principles by summing the radio contributions of all electrons and positrons in the shower. We use the radio simulation code CoREAS 16 , a plug-in of CORSIKA 17 , which follows this approach.It has been shown that X max can be accurately reconstructed from densely sampled radio measurements 18 . We use a hybrid approach, simultaneously fitting the radio and particle data. The radio component is very sensitive to X max , while the particle component is used for the energy measurement.The fit contains four free parameters: the shower core position (x, y), and scaling factors for the partic...