-
We present timing solutions for ten pulsars discovered in 350 MHz searches
with the Green Bank Telescope. Nine of these were discovered in the Green Bank
Northern Celestial Cap survey and one was discovered by students in the Pulsar
Search Collaboratory program in analysis of drift-scan data. Following
discovery and confirmation with the Green Bank Telescope, timing has yielded
phase-connected solutions with high precision measurements of rotational and
astrometric parameters. Eight of the pulsars are slow and isolated, including
PSR J0930$-$2301, a pulsar with nulling fraction lower limit of $\sim$30\% and
nulling timescale of seconds to minutes. This pulsar also shows evidence of
mode changing. The remaining two pulsars have undergone recycling, accreting
material from binary companions, resulting in higher spin frequencies. PSR
J0557$-$2948 is an isolated, 44 \rm{ms} pulsar that has been partially recycled
and is likely a former member of a binary system which was disrupted by a
second supernova. The paucity of such so-called `disrupted binary pulsars'
(DRPs) compared to double neutron star (DNS) binaries can be used to test
current evolutionary scenarios, especially the kicks imparted on the neutron
stars in the second supernova. There is some evidence that DRPs have larger
space velocities, which could explain their small numbers. PSR J1806+2819 is a
15 \rm{ms} pulsar in a 44 day orbit with a low mass white dwarf companion. We
did not detect the companion in archival optical data, indicating that it must
be older than 1200 Myr.
-
Symmetries play an important role in modern theories of gravity. The strong
equivalence principle (SEP) constitutes a collection of gravitational
symmetries which are all implemented by general relativity. Alternative
theories, however, are generally expected to violate some aspects of SEP. We
test three aspects of SEP using observed change rates in the orbital period and
eccentricity of binary pulsar J1713+0747: 1. the gravitational constant's
constancy as part of locational invariance of gravitation; 2. the
post-Newtonian parameter $\hat{\alpha}_3$ in gravitational Lorentz invariance;
3. the universality of free fall (UFF) for strongly self-gravitating bodies.
Based on the pulsar timing result of the combined dataset from the North
American Nanohertz Gravitational Observatory (NANOGrav) and the European Pulsar
Timing Array (EPTA), we find $\dot{G}/G = (-0.1 \pm 0.9) \times 10^{-12}\,{\rm
yr}^{-1}$, which is weaker than Solar system limits, but applies for strongly
self-gravitating objects. Furthermore, we obtain the constraints $|\Delta|<
0.002$ for the UFF test and $-3\times10^{-20} < \hat{\alpha}_3 <
4\times10^{-20}$ at 95% confidence. These are the first direct UFF and
$\hat{\alpha}_3$ tests based on pulsar binaries, and they overcome various
limitations of previous tests.
-
We report the discovery and initial follow-up of a double neutron star (DNS)
system, PSR J1946$+$2052, with the Arecibo L-Band Feed Array pulsar (PALFA)
survey. PSR J1946$+$2052 is a 17-ms pulsar in a 1.88-hour, eccentric ($e \, =\,
0.06$) orbit with a $\gtrsim 1.2 \, M_\odot$ companion. We have used the Jansky
Very Large Array to localize PSR J1946$+$2052 to a precision of 0.09 arcseconds
using a new phase binning mode. We have searched multiwavelength catalogs for
coincident sources but did not find any counterparts. The improved position
enabled a measurement of the spin period derivative of the pulsar ($\dot{P} \,
= \, 9\,\pm \, 2 \,\times 10^{-19}$); the small inferred magnetic field
strength at the surface ($B_S \, = \, 4 \, \times \, 10^9 \, \rm G$) indicates
that this pulsar has been recycled. This and the orbital eccentricity lead to
the conclusion that PSR J1946$+$2052 is in a DNS system. Among all known radio
pulsars in DNS systems, PSR J1946$+$2052 has the shortest orbital period and
the shortest estimated merger timescale, 46 Myr; at that time it will display
the largest spin effects on gravitational wave waveforms of any such system
discovered to date. We have measured the advance of periastron passage for this
system, $\dot{\omega} \, = \, 25.6 \, \pm \, 0.3\, \deg \rm yr^{-1}$, implying
a total system mass of only 2.50 $\pm$ 0.04 $M_\odot$, so it is among the
lowest mass DNS systems. This total mass measurement combined with the minimum
companion mass constrains the pulsar mass to $\lesssim 1.3 \, M_\odot$.
-
We search for an isotropic stochastic gravitational-wave background (GWB) in
the newly released $11$-year dataset from the North American Nanohertz
Observatory for Gravitational Waves (NANOGrav). While we find no significant
evidence for a GWB, we place constraints on a GWB from a population of
supermassive black-hole binaries, from cosmic strings, and from a primordial
GWB. For the first time, we find that the GWB upper limits and detection
statistics are sensitive to the Solar System ephemeris (SSE) model used, and
that SSE errors can mimic a GWB signal. To mitigate this effect, we developed
and implemented a novel approach that bridges systematic SSE differences,
producing the first PTA constraints that are robust against SSE uncertainties.
We place a $95\%$ upper limit on the GW strain amplitude of
$A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a
fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations
modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year
limit, calculated using the same procedure. Previous PTA upper limits on the
GWB will need revision in light of SSE systematic uncertainties. We also
characterize the combined influence of the mass-density of stars in galactic
cores, the eccentricity of binaries at formation, and the relation between the
mass of the central supermassive black hole and the galactic bulge (the
$M_\mathrm{BH}-M_\mathrm{bulge}$ relation). We constrain cosmic-string tension
on the basis of recent simulations, yielding an SSE-marginalized 95\% upper
limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of
$\sim 2$ better than the NANOGrav $9$-year constraints. We then use our new
Bayesian SSE model to limit the energy density of primordial GWBs,
corresponding to $\Omega_\mathrm{GWB}(f)h^2<3.4 \times 10^{-10}$ for a
radiation-dominated inflationary era. [ABRIDGED]
-
We describe an image-based method that uses two radio criteria, compactness
and spectral index, to identify promising pulsar candidates among Fermi Large
Area Telescope (LAT) unassociated sources. These criteria are applied to those
radio sources from the Giant Metrewave Radio Telescope all-sky survey at 150
MHz (TGSS ADR1) found within the error ellipses of unassociated sources from
the 3FGL catalog and a preliminary source list based on 7 years of LAT data.
After follow-up interferometric observations to identify extended or variable
sources, a list of 16 compact, steep-spectrum candidates is generated. An
ongoing search for pulsations in these candidates, in gamma rays and radio, has
found six millisecond pulsars and one normal pulsar. A comparison of this
method with existing selection criteria based on gamma-ray spectral and
variability properties suggests that the pulsar discovery space using Fermi may
be larger than previously thought. Radio imaging is a hitherto underutilized
source selection method that can be used, as with other multi-wavelength
techniques, in the search for Fermi pulsars.
-
In this work, we report the discovery and characterization of PSR J1411+2551,
a new binary pulsar discovered in the Arecibo 327 MHz Drift Pulsar Survey. Our
timing observations of the radio pulsar in the system span a period of about
2.5 years. This timing campaign allowed a precise measurement of its spin
period (62.4 ms) and its derivative (9.6 $\pm$ 0.7) $\times 10^{-20}\, \rm s\,
s^{-1}$; from these, we derive a characteristic age of $\sim 10\,$Gyr and a
surface magnetic field strength of 2.5 $\times 10^{9}$ G. These numbers
indicate that this pulsar was mildly recycled by accretion of matter from the
progenitor of the companion star. The system has an eccentric ($e\, = \, 0.17$)
2.61 day orbit. This eccentricity allows a highly significant measurement of
the rate of advance of periastron, $\dot{\omega} = 0.07686 \pm 0.00046
^{\circ}~{\rm yr}^{-1}$. Assuming general relativity accurately models the
orbital motion, this implies a total system mass M = $2.538 \pm 0.022
M_{\odot}$. The minimum companion mass is $0.92\, M_{\odot}$ and the maximum
pulsar mass is $1.62\, M_{\odot}$. The large companion mass and the orbital
eccentricity suggest that PSR J1411+2551 is a double neutron star system; the
lightest known to date including the DNS merger GW 170817. Furthermore, the
relatively low orbital eccentricity and small proper motion limits suggest that
the second supernova had a relatively small associated kick; this and the low
system mass suggest that it was an ultra-stripped supernova.
-
We report on a search for Fast Radio Bursts (FRBs) with the Green Bank
Northern Celestial Cap (GBNCC) Pulsar Survey at 350 MHz. Pointings amounting to
a total on-sky time of 61 days were searched to a DM of 3000 pc cm$^{-3}$ while
the rest (23 days; 29% of the total time) were searched to a DM of 500 pc
cm$^{-3}$. No FRBs were detected in the pointings observed through May 2016. We
estimate a 95% confidence upper limit on the FRB rate of $3.6\times 10^3$ FRBs
sky$^{-1}$ day$^{-1}$ above a peak flux density of 0.63 Jy at 350 MHz for an
intrinsic pulse width of 5 ms. We place constraints on the spectral index
$\alpha$ by running simulations for different astrophysical scenarios and
cumulative flux density distributions. The non-detection with GBNCC is
consistent with the 1.4-GHz rate reported for the Parkes surveys for $\alpha >
+0.35 $ in the absence of scattering and free-free absorption and $\alpha >
-0.3$ in the presence of scattering, for a Euclidean flux distribution. The
constraints imply that FRBs exhibit either a flat spectrum or a spectral
turnover at frequencies above 400 MHz. These constraints also allow estimation
of the number of bursts that can be detected with current and upcoming surveys.
We predict that CHIME may detect anywhere from several to $\sim$50 FRBs a day
(depending on model assumptions), making it well suited for interesting
constraints on spectral index, the log $N$-log $S$ slope and pulse profile
evolution across its bandwidth (400-800 MHz).
-
We report on the discovery and timing observations of 29 distant long-period
pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Following
discovery with the Arecibo Telescope, confirmation and timing observations of
these pulsars over several years at Jodrell Bank Observatory have yielded
high-precision positions and measurements of rotation and radiation properties.
We have used multi-frequency data to measure the interstellar scattering
properties of some of these pulsars. Most of the pulsars have properties that
mirror those of the previously known pulsar population, although four show some
notable characteristics. PSRs J1907+0631 and J1925+1720 are young and are
associated with supernova remnants or plerionic nebulae: J1907+0631 lies close
to the center of SNR G40.5-0.5, while J1925+1720 is coincident with a
high-energy Fermi gamma-ray source. One pulsar, J1932+1500, is in a
surprisingly eccentric, 199-day binary orbit with a companion having a minimum
mass of 0.33 solar masses. Several of the sources exhibit timing noise, and
two, PSRs J0611+1436 and J1907+0631, have both suffered large glitches, but
with very different post-glitch rotation properties. In particular, the
rotational period of PSR J0611+1436 will not recover to its pre-glitch value
for about 12 years, a far greater recovery timescale than seen following any
other large glitches.
-
We analyze dispersion measure (DM) variations of 37 millisecond pulsars in
the 9-year NANOGrav data release and constrain the sources of these variations.
Variations are significant for nearly all pulsars, with characteristic
timescales comparable to or even shorter than the average spacing between
observations. Five pulsars have periodic annual variations, 14 pulsars have
monotonically increasing or decreasing trends, and 13 pulsars show both
effects. Several pulsars show correlations between DM excesses and lines of
sight that pass close to the Sun. Mapping of the DM variations as a function of
the pulsar trajectory can identify localized ISM features and, in one case, an
upper limit to the size of the dispersing region of 13.2 AU. Finally, five
pulsars show very nearly quadratic structure functions, which could be
indicative of an underlying Kolmogorov medium. Four pulsars show roughly
Kolmogorov structure functions and another four show structure functions less
steep than Kolmogorov. One pulsar has too large an uncertainty to allow
comparisons. We discuss explanations for apparent departures from a
Kolmogorov-like spectrum, and show that the presence of other trends in the
data is the most likely cause.
-
We report the discovery of two long-term intermittent radio pulsars in the
ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo
Telescope, extended observations of these pulsars over several years at Jodrell
Bank Observatory have revealed the details of their rotation and radiation
properties. PSRs J1910+0517 and J1929+1357 show long-term extreme bi-modal
intermittency, switching between active (ON) and inactive (OFF) emission states
and indicating the presence of a large, hitherto unrecognised, underlying
population of such objects. For PSR J1929+1357, the initial duty cycle was
fON=0.008, but two years later this changed, quite abruptly, to fON=0.16. This
is the first time that a significant evolution in the activity of an
intermittent pulsar has been seen and we show that the spin-down rate of the
pulsar is proportional to the activity. The spin-down rate of PSR J1929+1357 is
increased by a factor of 1.8 when it is in active mode, similar to the increase
seen in the other three known long-term intermittent pulsars.
-
Gravitational wave astronomy using a pulsar timing array requires
high-quality millisecond pulsars, correctable interstellar propagation delays,
and high-precision measurements of pulse times of arrival. Here we identify
noise in timing residuals that exceeds that predicted for arrival time
estimation for millisecond pulsars observed by the North American Nanohertz
Observatory for Gravitational Waves. We characterize the excess noise using
variance and structure function analyses. We find that 26 out of 37 pulsars
show inconsistencies with a white-noise-only model based on the short timescale
analysis of each pulsar and we demonstrate that the excess noise has a red
power spectrum for 15 pulsars. We also decompose the excess noise into
chromatic (radio-frequency-dependent) and achromatic components. Associating
the achromatic red-noise component with spin noise and including additional
power-spectrum-based estimates from the literature, we estimate a scaling law
in terms of spin parameters (frequency and frequency derivative) and data-span
length and compare it to the scaling law of Shannon \& Cordes (2010). We
briefly discuss our results in terms of detection of gravitational waves at
nanohertz frequencies.
-
The use of pulsars as astrophysical clocks for gravitational wave experiments
demands the highest possible timing precision. Pulse times of arrival (TOAs)
are limited by stochastic processes that occur in the pulsar itself, along the
line of sight through the interstellar medium, and in the measurement process.
On timescales of seconds to hours, the TOA variance exceeds that from
template-fitting errors due to additive noise. We assess contributions to the
total variance from two additional effects: amplitude and phase jitter
intrinsic to single pulses and changes in the interstellar impulse response
from scattering. The three effects have different dependencies on time,
frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the
North American Nanohertz Observatory for Gravitational Waves to assess the
individual contributions to the overall intraday noise budget for each pulsar.
We detect jitter in 22 pulsars and estimate the average value of RMS jitter in
our pulsars to be $\sim 1\%$ of pulse phase. We examine how jitter evolves as a
function of frequency and find evidence for evolution. Finally, we compare our
measurements with previous noise parameter estimates and discuss methods to
improve gravitational wave detection pipelines.
-
Utilizing the all-sky imaging capabilities of the LWA1 radio telescope along
with a host of all-sky optical cameras, we have now observed 44 optical meteor
counterparts to radio afterglows. Combining these observations we have
determined the geographic positions of all 44 afterglows. Comparing the number
of radio detections as a function of altitude above sea level to the number of
expected bright meteors we find a strong altitudinal dependence characterized
by a cutoff below $\sim$ 90 km, below which no radio emission occurs, despite
the fact that many of the observed optical meteors penetrated well below this
altitude. This cutoff suggests that wave damping from electron collisions is an
important factor for the evolution of radio afterglows, which agrees with the
hypothesis that the emission is the result of electron plasma wave emission.
-
We report the discovery and timing results for five millisecond pulsars
(MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726,
J1938+2516, and J1957+2516. Timing observations of the 5 pulsars were conducted
with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3
yr. All of the MSPs except one (PSR J1914+0659) are in binary systems with low
eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ~0.1
$M_\odot$ companion and possible eclipses that last ~10% of the orbit. The
position of PSR J1957+2516 is also coincident with a NIR source. All 5 MSPs are
distant (>3.1 kpc) as determined from their dispersion measures, and none of
them show evidence of $\gamma$-ray pulsations in a search of Fermi Gamma-Ray
Space Telescope data. These 5 MSPs bring the total number of MSPs discovered by
the PALFA survey to 26 and further demonstrate the power of this survey in
finding distant, highly dispersed MSPs deep in the Galactic plane.
-
We report here the Einstein@Home discovery of PSR J1913+1102, a 27.3-ms
pulsar found in data from the ongoing Arecibo PALFA pulsar survey. The pulsar
is in a 4.95-hr double neutron star (DNS) system with an eccentricity of 0.089.
From radio timing with the Arecibo 305-m telescope, we measure the rate of
advance of periastron to be 5.632(18) deg/yr. Assuming general relativity
accurately models the orbital motion, this corresponds to a total system mass
of 2.875(14) solar masses, similar to the mass of the most massive DNS known to
date, B1913+16, but with a much smaller eccentricity. The small eccentricity
indicates that the second-formed neutron star (the companion of PSR J1913+1102)
was born in a supernova with a very small associated kick and mass loss. In
that case this companion is likely, by analogy with other systems, to be a
light (1.2 solar mass) neutron star; the system would then be highly
asymmetric. A search for radio pulsations from the companion yielded no
plausible detections, so we can't yet confirm this mass asymmetry. By the end
of 2016, timing observations should permit the detection of two additional
post-Keplerian parameters: the Einstein delay, which will enable precise mass
measurements and a verification of the possible mass asymmetry of the system,
and the orbital decay due to the emission of gravitational waves, which will
allow another test of the radiative properties of gravity. The latter effect
will cause the system to coalesce in ~0.5 Gyr.
-
Dense, continuous pulsar timing observations over a 24-hr period provide a
method for probing intermediate gravitational wave (GW) frequencies from 10
microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North
American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes
Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array
(IPTA) all use millisecond pulsar observations to detect or constrain GWs
typically at nanohertz frequencies. In the case of the IPTA's nine-telescope
24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the
intermediate frequency regime can be produced. The negligible change in
dispersion measure during the observation minimizes red noise in the timing
residuals, constraining any contributions from GWs due to individual sources.
At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in
the pulsar's direction.
-
We present a large survey of giant pulses from the Crab Pulsar as observed
with the first station of the Long Wavelength Array. Automated methods for
detecting giant pulses at low frequencies where scattering becomes prevalent
are also explored. More than 1400 pulses were detected across four frequency
bands between 20 - 84 MHz over a seven month period beginning in 2013, with
additional followup observations in late 2014 and early 2015. A handful of
these pulses were detected simultaneously across all four frequency bands. We
examine pulse characteristics, including pulse broadening and power law indices
for amplitude distributions. We find that the flux density increases toward
shorter wavelengths, consistent with a spectral turnover at 100 MHz. Our
observations uniquely span multiple scattering epochs, manifesting as a notable
trend in the number of detections per observation. These results are
characteristic of the variable interface between the synchrotron nebula and the
surrounding interstellar medium.
-
We report on simultaneous X-ray and radio observations of the mode-switching
pulsar PSR B0943+10 obtained with the XMM-Newton satellite and the LOFAR, LWA
and Arecibo radio telescopes in November 2014. We confirm the synchronous
X-ray/radio switching between a radio-bright (B) and a radio-quiet (Q) mode, in
which the X-ray flux is a factor ~2.4 higher than in the B-mode. We discovered
X-ray pulsations, with pulsed fraction of 38+/-5% (0.5-2 keV), during the
B-mode, and confirm their presence in Q-mode, where the pulsed fraction
increases with energy from ~20% up to ~65% at 2 keV. We found marginal evidence
for an increase in the X-ray pulsed fraction during B-mode on a timescale of
hours. The Q-mode X-ray spectrum requires a fit with a two-component model
(either a power-law plus blackbody or the sum of two blackbodies), while the
B-mode spectrum is well fit by a single blackbody (a single power-law is
rejected). With a maximum likelihood analysis, we found that in Q-mode the
pulsed emission has a thermal blackbody spectrum with temperature ~3.4x10^6 K
and the unpulsed emission is a power-law with photon index ~2.5, while during
B-mode both the pulsed and unpulsed emission can be fit by either a blackbody
or a power law with similar values of temperature and photon index. A Chandra
image shows no evidence for diffuse X-ray emission. These results support a
scenario in which both unpulsed non-thermal emission, likely of magnetospheric
origin, and pulsed thermal emission from a small polar cap (~1500 m^2) with a
strong non-dipolar magnetic field (~10^{14} G), are present during both radio
modes and vary in intensity in a correlated way. This is broadly consistent
with the predictions of the partially screened gap model and does not
necessarily imply global magnetospheric rearrangements to explain the mode
switching.
-
PSR J1024$-$0719 is a millisecond pulsar that was long thought to be
isolated. However, puzzling results concerning its velocity, distance, and low
rotational period derivative have led to reexamination of its properties. We
present updated radio timing observations along with new and archival optical
data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr)
binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z
\approx -0.9\,$ dex) main sequence star. Such a system can explain most of the
anomalous properties of this pulsar. We suggest that this system formed through
a dynamical exchange in a globular cluster that ejected it into a halo orbit,
consistent with the low observed metallicity for the stellar companion. Further
astrometric and radio timing observations such as measurement of the third
period derivative could strongly constrain the range of orbital parameters.
-
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.
-
We present Clusterrank, a new algorithm for identifying dispersed
astrophysical pulses. Such pulses are commonly detected from Galactic pulsars
and rotating radio transients (RRATs), which are neutron stars with sporadic
radio emission. More recently, isolated, highly dispersed pulses dubbed fast
radio bursts (FRBs) have been identified as the potential signature of an
extragalactic cataclysmic radio source distinct from pulsars and RRATs.
Clusterrank helped us discover 14 pulsars and 8 RRATs in data from the Arecibo
327 MHz Drift Pulsar Survey (AO327). The new RRATs have DMs in the range $23.5
- 86.6$ pc cm$^{-3}$ and periods in the range $0.172 - 3.901$ s. The new
pulsars have DMs in the range $23.6 - 133.3$ pc cm$^{-3}$ and periods in the
range $1.249 - 5.012$ s, and include two nullers and a mode-switching object.
We estimate an upper limit on the all-sky FRB rate of $10^5$ day$^{-1}$ for
bursts with a width of 10 ms and flux density $\gtrsim 83$ mJy. The DMs of all
new discoveries are consistent with a Galactic origin. In comparing statistics
of the new RRATs with sources from the RRATalog, we find that both sets are
drawn from the same period distribution. In contrast, we find that the period
distribution of the new pulsars is different from the period distributions of
canonical pulsars in the ATNF catalog or pulsars found in AO327 data by a
periodicity search. This indicates that Clusterrank is a powerful complement to
periodicity searches and uncovers a subset of the pulsar population that has so
far been underrepresented in survey results and therefore in Galactic pulsar
population models.
-
We analyse the stochastic properties of the 49 pulsars that comprise the
first International Pulsar Timing Array (IPTA) data release. We use Bayesian
methodology, performing model selection to determine the optimal description of
the stochastic signals present in each pulsar. In addition to spin-noise and
dispersion-measure (DM) variations, these models can include timing noise
unique to a single observing system, or frequency band. We show the improved
radio-frequency coverage and presence of overlapping data from different
observing systems in the IPTA data set enables us to separate both system and
band-dependent effects with much greater efficacy than in the individual PTA
data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM
variations, significant band-dependent noise that is coherent between PTAs
which we interpret as coming from time-variable scattering or refraction in the
ionised interstellar medium. Failing to model these different contributions
appropriately can dramatically alter the astrophysical interpretation of the
stochastic signals observed in the residuals. In some cases, the spectral
exponent of the spin noise signal can vary from 1.6 to 4 depending upon the
model, which has direct implications for the long-term sensitivity of the
pulsar to a stochastic gravitational-wave (GW) background. By using a more
appropriate model, however, we can greatly improve a pulsar's sensitivity to
GWs. For example, including system and band-dependent signals in the PSR
J0437$-$4715 data set improves the upper limit on a fiducial GW background by
$\sim 60\%$ compared to a model that includes DM variations and spin-noise
only.
-
The 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.
-
We present high-precision timing observations spanning up to nine years for
37 millisecond pulsars monitored with the Green Bank and Arecibo radio
telescopes as part of the North American Nanohertz Observatory for
Gravitational Waves (NANOGrav) project. We describe the observational and
instrumental setups used to collect the data, and methodology applied for
calculating pulse times of arrival; these include novel methods for measuring
instrumental offsets and characterizing low signal-to-noise ratio timing
results. The time of arrival data are fit to a physical timing model for each
source, including terms that characterize time-variable dispersion measure and
frequency-dependent pulse shape evolution. In conjunction with the timing model
fit, we have performed a Bayesian analysis of a parameterized timing noise
model for each source, and detect evidence for excess low-frequency, or "red,"
timing noise in 10 of the pulsars. For 5 of these cases this is likely due to
interstellar medium propagation effects rather than intrisic spin variations.
Subsequent papers in this series will present further analysis of this data set
aimed at detecting or limiting the presence of nanohertz-frequency
gravitational wave signals.
-
We present dynamic spectra from the LWA1 telescope of two large meteors
(fireballs) observed to emit between 37 and 54 MHz. These spectra show the
first ever recorded broadband measurements of this newly discovered VHF
emission. The spectra show that the emission is smooth and steep, getting very
bright at lower frequencies. We suggest that this signal is possibly emission
of Langmuir waves and that these waves could be excited by a weak electron beam
within the trail. The spectra of one fireball displays broadband temporal
frequency sweeps. We suggest that these sweeps are evidence of individual
expanding clumps of emitting plasma. While some of these proposed clumps may
have formed at the very beginning of the fireball event, others must have
formed seconds after the initial event.
We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released $11$-year dataset from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no significant evidence for a GWB, we place constraints on a GWB from a population of supermassive black-hole binaries, from cosmic strings, and from a primordial GWB. For the first time, we find that the GWB upper limits and detection statistics are sensitive to the Solar System ephemeris (SSE) model used, and that SSE errors can mimic a GWB signal. To mitigate this effect, we developed and implemented a novel approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We place a $95\%$ upper limit on the GW strain amplitude of $A_\mathrm{GWB}<1.45\times 10^{-15}$ at a frequency of $f=1$ yr$^{-1}$ for a fiducial $f^{-2/3}$ power-law spectrum, and with inter-pulsar correlations modeled. This is a factor of $\sim 2$ improvement over the NANOGrav $9$-year limit, calculated using the same procedure. Previous PTA upper limits on the GWB will need revision in light of SSE systematic uncertainties. We also characterize the combined influence of the mass-density of stars in galactic cores, the eccentricity of binaries at formation, and the relation between the mass of the central supermassive black hole and the galactic bulge (the $M_\mathrm{BH}-M_\mathrm{bulge}$ relation). We constrain cosmic-string tension on the basis of recent simulations, yielding an SSE-marginalized 95\% upper limit on the cosmic string tension of $G\mu < 5.3\times 10^{-11}$---a factor of $\sim 2$ better than the NANOGrav $9$-year constraints. We then use our new Bayesian SSE model to limit the energy density of primordial GWBs, corresponding to $\Omega_\mathrm{GWB}(f)h^2<3.4 \times 10^{-10}$ for a radiation-dominated inflationary era. [ABRIDGED]
We report on a search for Fast Radio Bursts (FRBs) with the Green Bank Northern Celestial Cap (GBNCC) Pulsar Survey at 350 MHz. Pointings amounting to a total on-sky time of 61 days were searched to a DM of 3000 pc cm$^{-3}$ while the rest (23 days; 29% of the total time) were searched to a DM of 500 pc cm$^{-3}$. No FRBs were detected in the pointings observed through May 2016. We estimate a 95% confidence upper limit on the FRB rate of $3.6\times 10^3$ FRBs sky$^{-1}$ day$^{-1}$ above a peak flux density of 0.63 Jy at 350 MHz for an intrinsic pulse width of 5 ms. We place constraints on the spectral index $\alpha$ by running simulations for different astrophysical scenarios and cumulative flux density distributions. The non-detection with GBNCC is consistent with the 1.4-GHz rate reported for the Parkes surveys for $\alpha > +0.35 $ in the absence of scattering and free-free absorption and $\alpha > -0.3$ in the presence of scattering, for a Euclidean flux distribution. The constraints imply that FRBs exhibit either a flat spectrum or a spectral turnover at frequencies above 400 MHz. These constraints also allow estimation of the number of bursts that can be detected with current and upcoming surveys. We predict that CHIME may detect anywhere from several to $\sim$50 FRBs a day (depending on model assumptions), making it well suited for interesting constraints on spectral index, the log $N$-log $S$ slope and pulse profile evolution across its bandwidth (400-800 MHz).
We report on the discovery and timing observations of 29 distant long-period pulsars discovered in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, confirmation and timing observations of these pulsars over several years at Jodrell Bank Observatory have yielded high-precision positions and measurements of rotation and radiation properties. We have used multi-frequency data to measure the interstellar scattering properties of some of these pulsars. Most of the pulsars have properties that mirror those of the previously known pulsar population, although four show some notable characteristics. PSRs J1907+0631 and J1925+1720 are young and are associated with supernova remnants or plerionic nebulae: J1907+0631 lies close to the center of SNR G40.5-0.5, while J1925+1720 is coincident with a high-energy Fermi gamma-ray source. One pulsar, J1932+1500, is in a surprisingly eccentric, 199-day binary orbit with a companion having a minimum mass of 0.33 solar masses. Several of the sources exhibit timing noise, and two, PSRs J0611+1436 and J1907+0631, have both suffered large glitches, but with very different post-glitch rotation properties. In particular, the rotational period of PSR J0611+1436 will not recover to its pre-glitch value for about 12 years, a far greater recovery timescale than seen following any other large glitches.
We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the 9-year NANOGrav data release and constrain the sources of these variations. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 13 pulsars show both effects. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized ISM features and, in one case, an upper limit to the size of the dispersing region of 13.2 AU. Finally, five pulsars show very nearly quadratic structure functions, which could be indicative of an underlying Kolmogorov medium. Four pulsars show roughly Kolmogorov structure functions and another four show structure functions less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and show that the presence of other trends in the data is the most likely cause.
We report the discovery of two long-term intermittent radio pulsars in the ongoing Arecibo PALFA pulsar survey. Following discovery with the Arecibo Telescope, extended observations of these pulsars over several years at Jodrell Bank Observatory have revealed the details of their rotation and radiation properties. PSRs J1910+0517 and J1929+1357 show long-term extreme bi-modal intermittency, switching between active (ON) and inactive (OFF) emission states and indicating the presence of a large, hitherto unrecognised, underlying population of such objects. For PSR J1929+1357, the initial duty cycle was fON=0.008, but two years later this changed, quite abruptly, to fON=0.16. This is the first time that a significant evolution in the activity of an intermittent pulsar has been seen and we show that the spin-down rate of the pulsar is proportional to the activity. The spin-down rate of PSR J1929+1357 is increased by a factor of 1.8 when it is in active mode, similar to the increase seen in the other three known long-term intermittent pulsars.
Gravitational wave astronomy using a pulsar timing array requires high-quality millisecond pulsars, correctable interstellar propagation delays, and high-precision measurements of pulse times of arrival. Here we identify noise in timing residuals that exceeds that predicted for arrival time estimation for millisecond pulsars observed by the North American Nanohertz Observatory for Gravitational Waves. We characterize the excess noise using variance and structure function analyses. We find that 26 out of 37 pulsars show inconsistencies with a white-noise-only model based on the short timescale analysis of each pulsar and we demonstrate that the excess noise has a red power spectrum for 15 pulsars. We also decompose the excess noise into chromatic (radio-frequency-dependent) and achromatic components. Associating the achromatic red-noise component with spin noise and including additional power-spectrum-based estimates from the literature, we estimate a scaling law in terms of spin parameters (frequency and frequency derivative) and data-span length and compare it to the scaling law of Shannon \& Cordes (2010). We briefly discuss our results in terms of detection of gravitational waves at nanohertz frequencies.
The use of pulsars as astrophysical clocks for gravitational wave experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting errors due to additive noise. We assess contributions to the total variance from two additional effects: amplitude and phase jitter intrinsic to single pulses and changes in the interstellar impulse response from scattering. The three effects have different dependencies on time, frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the North American Nanohertz Observatory for Gravitational Waves to assess the individual contributions to the overall intraday noise budget for each pulsar. We detect jitter in 22 pulsars and estimate the average value of RMS jitter in our pulsars to be $\sim 1\%$ of pulse phase. We examine how jitter evolves as a function of frequency and find evidence for evolution. Finally, we compare our measurements with previous noise parameter estimates and discuss methods to improve gravitational wave detection pipelines.
Utilizing the all-sky imaging capabilities of the LWA1 radio telescope along with a host of all-sky optical cameras, we have now observed 44 optical meteor counterparts to radio afterglows. Combining these observations we have determined the geographic positions of all 44 afterglows. Comparing the number of radio detections as a function of altitude above sea level to the number of expected bright meteors we find a strong altitudinal dependence characterized by a cutoff below $\sim$ 90 km, below which no radio emission occurs, despite the fact that many of the observed optical meteors penetrated well below this altitude. This cutoff suggests that wave damping from electron collisions is an important factor for the evolution of radio afterglows, which agrees with the hypothesis that the emission is the result of electron plasma wave emission.
We report the discovery and timing results for five millisecond pulsars (MSPs) from the Arecibo PALFA survey: PSRs J1906+0055, J1914+0659, J1933+1726, J1938+2516, and J1957+2516. Timing observations of the 5 pulsars were conducted with the Arecibo and Lovell telescopes for time spans ranging from 1.5 to 3.3 yr. All of the MSPs except one (PSR J1914+0659) are in binary systems with low eccentricities. PSR J1957+2516 is likely a redback pulsar, with a ~0.1 $M_\odot$ companion and possible eclipses that last ~10% of the orbit. The position of PSR J1957+2516 is also coincident with a NIR source. All 5 MSPs are distant (>3.1 kpc) as determined from their dispersion measures, and none of them show evidence of $\gamma$-ray pulsations in a search of Fermi Gamma-Ray Space Telescope data. These 5 MSPs bring the total number of MSPs discovered by the PALFA survey to 26 and further demonstrate the power of this survey in finding distant, highly dispersed MSPs deep in the Galactic plane.
We report here the Einstein@Home discovery of PSR J1913+1102, a 27.3-ms pulsar found in data from the ongoing Arecibo PALFA pulsar survey. The pulsar is in a 4.95-hr double neutron star (DNS) system with an eccentricity of 0.089. From radio timing with the Arecibo 305-m telescope, we measure the rate of advance of periastron to be 5.632(18) deg/yr. Assuming general relativity accurately models the orbital motion, this corresponds to a total system mass of 2.875(14) solar masses, similar to the mass of the most massive DNS known to date, B1913+16, but with a much smaller eccentricity. The small eccentricity indicates that the second-formed neutron star (the companion of PSR J1913+1102) was born in a supernova with a very small associated kick and mass loss. In that case this companion is likely, by analogy with other systems, to be a light (1.2 solar mass) neutron star; the system would then be highly asymmetric. A search for radio pulsations from the companion yielded no plausible detections, so we can't yet confirm this mass asymmetry. By the end of 2016, timing observations should permit the detection of two additional post-Keplerian parameters: the Einstein delay, which will enable precise mass measurements and a verification of the possible mass asymmetry of the system, and the orbital decay due to the emission of gravitational waves, which will allow another test of the radiative properties of gravity. The latter effect will cause the system to coalesce in ~0.5 Gyr.
Dense, continuous pulsar timing observations over a 24-hr period provide a method for probing intermediate gravitational wave (GW) frequencies from 10 microhertz to 20 millihertz. The European Pulsar Timing Array (EPTA), the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), the Parkes Pulsar Timing Array (PPTA), and the combined International Pulsar Timing Array (IPTA) all use millisecond pulsar observations to detect or constrain GWs typically at nanohertz frequencies. In the case of the IPTA's nine-telescope 24-Hour Global Campaign on millisecond pulsar J1713+0747, GW limits in the intermediate frequency regime can be produced. The negligible change in dispersion measure during the observation minimizes red noise in the timing residuals, constraining any contributions from GWs due to individual sources. At 10$^{-5}$Hz, the 95% upper limit on strain is 10$^{-11}$ for GW sources in the pulsar's direction.
We present a large survey of giant pulses from the Crab Pulsar as observed with the first station of the Long Wavelength Array. Automated methods for detecting giant pulses at low frequencies where scattering becomes prevalent are also explored. More than 1400 pulses were detected across four frequency bands between 20 - 84 MHz over a seven month period beginning in 2013, with additional followup observations in late 2014 and early 2015. A handful of these pulses were detected simultaneously across all four frequency bands. We examine pulse characteristics, including pulse broadening and power law indices for amplitude distributions. We find that the flux density increases toward shorter wavelengths, consistent with a spectral turnover at 100 MHz. Our observations uniquely span multiple scattering epochs, manifesting as a notable trend in the number of detections per observation. These results are characteristic of the variable interface between the synchrotron nebula and the surrounding interstellar medium.
We report on simultaneous X-ray and radio observations of the mode-switching pulsar PSR B0943+10 obtained with the XMM-Newton satellite and the LOFAR, LWA and Arecibo radio telescopes in November 2014. We confirm the synchronous X-ray/radio switching between a radio-bright (B) and a radio-quiet (Q) mode, in which the X-ray flux is a factor ~2.4 higher than in the B-mode. We discovered X-ray pulsations, with pulsed fraction of 38+/-5% (0.5-2 keV), during the B-mode, and confirm their presence in Q-mode, where the pulsed fraction increases with energy from ~20% up to ~65% at 2 keV. We found marginal evidence for an increase in the X-ray pulsed fraction during B-mode on a timescale of hours. The Q-mode X-ray spectrum requires a fit with a two-component model (either a power-law plus blackbody or the sum of two blackbodies), while the B-mode spectrum is well fit by a single blackbody (a single power-law is rejected). With a maximum likelihood analysis, we found that in Q-mode the pulsed emission has a thermal blackbody spectrum with temperature ~3.4x10^6 K and the unpulsed emission is a power-law with photon index ~2.5, while during B-mode both the pulsed and unpulsed emission can be fit by either a blackbody or a power law with similar values of temperature and photon index. A Chandra image shows no evidence for diffuse X-ray emission. These results support a scenario in which both unpulsed non-thermal emission, likely of magnetospheric origin, and pulsed thermal emission from a small polar cap (~1500 m^2) with a strong non-dipolar magnetic field (~10^{14} G), are present during both radio modes and vary in intensity in a correlated way. This is broadly consistent with the predictions of the partially screened gap model and does not necessarily imply global magnetospheric rearrangements to explain the mode switching.
PSR J1024$-$0719 is a millisecond pulsar that was long thought to be isolated. However, puzzling results concerning its velocity, distance, and low rotational period derivative have led to reexamination of its properties. We present updated radio timing observations along with new and archival optical data that show PSR J1024$-$0719 is most likely in a long period (2$-$20 kyr) binary system with a low-mass ($\approx 0.4\,M_\odot$) low-metallicity ($Z \approx -0.9\,$ dex) main sequence star. Such a system can explain most of the anomalous properties of this pulsar. We suggest that this system formed through a dynamical exchange in a globular cluster that ejected it into a halo orbit, consistent with the low observed metallicity for the stellar companion. Further astrometric and radio timing observations such as measurement of the third period derivative could strongly constrain the range of orbital parameters.
We present Clusterrank, a new algorithm for identifying dispersed astrophysical pulses. Such pulses are commonly detected from Galactic pulsars and rotating radio transients (RRATs), which are neutron stars with sporadic radio emission. More recently, isolated, highly dispersed pulses dubbed fast radio bursts (FRBs) have been identified as the potential signature of an extragalactic cataclysmic radio source distinct from pulsars and RRATs. Clusterrank helped us discover 14 pulsars and 8 RRATs in data from the Arecibo 327 MHz Drift Pulsar Survey (AO327). The new RRATs have DMs in the range $23.5 - 86.6$ pc cm$^{-3}$ and periods in the range $0.172 - 3.901$ s. The new pulsars have DMs in the range $23.6 - 133.3$ pc cm$^{-3}$ and periods in the range $1.249 - 5.012$ s, and include two nullers and a mode-switching object. We estimate an upper limit on the all-sky FRB rate of $10^5$ day$^{-1}$ for bursts with a width of 10 ms and flux density $\gtrsim 83$ mJy. The DMs of all new discoveries are consistent with a Galactic origin. In comparing statistics of the new RRATs with sources from the RRATalog, we find that both sets are drawn from the same period distribution. In contrast, we find that the period distribution of the new pulsars is different from the period distributions of canonical pulsars in the ATNF catalog or pulsars found in AO327 data by a periodicity search. This indicates that Clusterrank is a powerful complement to periodicity searches and uncovers a subset of the pulsar population that has so far been underrepresented in survey results and therefore in Galactic pulsar population models.
We analyse the stochastic properties of the 49 pulsars that comprise the first International Pulsar Timing Array (IPTA) data release. We use Bayesian methodology, performing model selection to determine the optimal description of the stochastic signals present in each pulsar. In addition to spin-noise and dispersion-measure (DM) variations, these models can include timing noise unique to a single observing system, or frequency band. We show the improved radio-frequency coverage and presence of overlapping data from different observing systems in the IPTA data set enables us to separate both system and band-dependent effects with much greater efficacy than in the individual PTA data sets. For example, we show that PSR J1643$-$1224 has, in addition to DM variations, significant band-dependent noise that is coherent between PTAs which we interpret as coming from time-variable scattering or refraction in the ionised interstellar medium. Failing to model these different contributions appropriately can dramatically alter the astrophysical interpretation of the stochastic signals observed in the residuals. In some cases, the spectral exponent of the spin noise signal can vary from 1.6 to 4 depending upon the model, which has direct implications for the long-term sensitivity of the pulsar to a stochastic gravitational-wave (GW) background. By using a more appropriate model, however, we can greatly improve a pulsar's sensitivity to GWs. For example, including system and band-dependent signals in the PSR J0437$-$4715 data set improves the upper limit on a fiducial GW background by $\sim 60\%$ compared to a model that includes DM variations and spin-noise only.
The 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.
We present high-precision timing observations spanning up to nine years for 37 millisecond pulsars monitored with the Green Bank and Arecibo radio telescopes as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project. We describe the observational and instrumental setups used to collect the data, and methodology applied for calculating pulse times of arrival; these include novel methods for measuring instrumental offsets and characterizing low signal-to-noise ratio timing results. The time of arrival data are fit to a physical timing model for each source, including terms that characterize time-variable dispersion measure and frequency-dependent pulse shape evolution. In conjunction with the timing model fit, we have performed a Bayesian analysis of a parameterized timing noise model for each source, and detect evidence for excess low-frequency, or "red," timing noise in 10 of the pulsars. For 5 of these cases this is likely due to interstellar medium propagation effects rather than intrisic spin variations. Subsequent papers in this series will present further analysis of this data set aimed at detecting or limiting the presence of nanohertz-frequency gravitational wave signals.
We present dynamic spectra from the LWA1 telescope of two large meteors (fireballs) observed to emit between 37 and 54 MHz. These spectra show the first ever recorded broadband measurements of this newly discovered VHF emission. The spectra show that the emission is smooth and steep, getting very bright at lower frequencies. We suggest that this signal is possibly emission of Langmuir waves and that these waves could be excited by a weak electron beam within the trail. The spectra of one fireball displays broadband temporal frequency sweeps. We suggest that these sweeps are evidence of individual expanding clumps of emitting plasma. While some of these proposed clumps may have formed at the very beginning of the fireball event, others must have formed seconds after the initial event.