Simultaneous Multi-band Radio & X-ray Observations of the Galactic Center Magnetar SGR 1745−-2900

We report on multi-frequency, wideband radio observations of the Galactic Center magnetar (SGR 1745$-$2900) with the Green Bank Telescope for $\sim$100 days immediately following its initial X-ray outburst in April 2013. We made multiple simultaneous observations at 1.5, 2.0, and 8.9 GHz, allowing us to examine the magnetar's flux evolution, radio spectrum, and interstellar medium parameters (such as the dispersion measure (DM), the scattering timescale and its index). During two epochs, we have simultaneous observations from the Chandra X-ray Observatory, which permitted the absolute alignment of the radio and X-ray profiles. As with the two other radio magnetars with published alignments, the radio profile lies within the broad peak of the X-ray profile, preceding the X-ray profile maximum by $\sim$0.2 rotations. We also find that the radio spectral index $\gamma$ is significantly negative between $\sim$2 and 9 GHz; during the final $\sim$30 days of our observations $\gamma \sim -1.4$, which is typical of canonical pulsars. The radio flux has not decreased during this outburst, whereas the long-term trends in the other radio magnetars show concomitant fading of the radio and X-ray fluxes. Finally, our wideband measurements of the DMs taken in adjacent frequency bands in tandem are stochastically inconsistent with one another. Based on recent theoretical predictions, we consider the possibility that the dispersion measure is frequency-dependent. Despite having several properties in common with the other radio magnetars, such as $L_{\textrm{X,qui}}/L_{\textrm{rot}} \lesssim 1$, an increase in the radio flux during the X-ray flux decay has not been observed thus far in other systems.Comment: 15 pages, 9 figures, 3 tables; accepted to Ap

The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.Comment: 10 pages, 11 figures. Accepted by Astrophysical Journal. Please send any comments/questions to S. J. Vigeland ([email protected]

Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.Comment: 14 pages, 6 figures. Accepted by Ap

Testing Theories of Gravitation Using 21-Year Timing of Pulsar Binary J1713+0747

We report 21-year timing of one of the most precise pulsars: PSR J1713+0747. Its pulse times of arrival are well modeled by a comprehensive pulsar binary model including its three-dimensional orbit and a noise model that incorporates short-and long-timescale correlated noise such as jitter and red noise. Its timing residuals have weighted root mean square similar to 92 ns. The new data set allows us to update and improve previous measurements of the system properties, including the masses of the neutron star (1.31 +/- 0.11 M-circle dot) and the companion white dwarf (0.286 +/- 0.012 M-circle dot) as well as their parallax distance 1.15 +/- 0.03 kpc. We measured the intrinsic change in orbital period, (P) over dot(b)(Int), is -0.20 +/- 0.17 ps s(-1), which is not distinguishable from zero. This result, combined with the measured (P) over dot(b)(Int) of other pulsars, can place a generic limit on potential changes in the gravitational constant G. We found that (G) over dot/G is consistent with zero [(-0.6 +/- 1.1) x 10(-12) yr(-1), 95% confidence] and changes at least a factor of 31 (99.7% confidence) more slowly than the average expansion rate of the universe. This is the best (G) over dot/G limit from pulsar binary systems. The (P) over dot(b)(Int) of pulsar binaries can also place limits on the putative coupling constant for dipole gravitational radiation kappa(D) = (-0.9 +/- 3.3) 10(-4) (95% confidence). Finally, the nearly circular orbit of this pulsar binary allows us to constrain statistically the strong-field post-Newtonian parameters Delta, which describes the violation of strong equivalence principle, and (alpha) over cap (3), which describes a breaking of both Lorentz invariance in gravitation and conservation of momentum. We found, at 95% confidence, Delta <0.01 and (3) <2 x 10(-20) based on PSR J1713+0747

The NANOGrav Nine-year Data Set:Noise Budget for Pulsar Arrival Times on Intraday Timescales

The use of pulsars as astrophysical clocks for gravitational wave (GW) 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 GWs 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 GW detection pipelines

The NANOGrav Nine-year Data Set: Measurement and Analysis of Variations in Dispersion Measures

We analyze dispersion measure (DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. 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 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. 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 interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs 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 we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause

The NANOGrav Nine-year Data Set: Excess Noise in Millisecond Pulsar Arrival Times

Gravitational wave (GW) astronomy using a pulsar timing array requires high-quality millisecond pulsars (MSPs), 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 MSPs 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. We briefly discuss our results in terms of detection of GWs at nanohertz frequencies

The NANOGrav Nine-year Data Set:Observations, Arrival Time Measurements, and Analysis of 37 Millisecond Pulsars

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

The NANOGrav Nine-year Data Set:Mass and Geometric Measurements of Binary Millisecond Pulsars

We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make 14 significant measurements of the Shapiro delay, including new detections in four pulsar-binary systems (PSRs J0613−0200, J2017+0603, J2302+4442, and J2317+1439), and derive estimates of the binary-component masses and orbital inclination for these MSP-binary systems. We find a wide range of binary pulsar masses, with values as low as ${m}_{{\rm{p}}}={1.18}_{-0.09}^{+0.10}\,{M}_{\odot }$ for PSR J1918−0642 and as high as ${m}_{{\rm{p}}}={1.928}_{-0.017}^{+0.017}\,{M}_{\odot }$ for PSR J1614−2230 (both 68.3% credibility). We make an improved measurement of the Shapiro timing delay in the PSR J1918−0642 and J2043+1711 systems, measuring the pulsar mass in the latter system to be ${m}_{{\rm{p}}}={1.41}_{-0.18}^{+0.21}\,{M}_{\odot }$ (68.3% credibility) for the first time. We measure secular variations of one or more orbital elements in many systems, and use these measurements to further constrain our estimates of the pulsar and companion masses whenever possible. In particular, we used the observed Shapiro delay and periastron advance due to relativistic gravity in the PSR J1903+0327 system to derive a pulsar mass of ${m}_{{\rm{p}}}={1.65}_{-0.02}^{+0.02}\,{M}_{\odot }$ (68.3% credibility). We discuss the implications that our mass measurements have on the overall neutron-star mass distribution, and on the "mass/orbital-period" correlation due to extended mass transfer

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 pulsar timing array (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 ionized 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 similar to 60 per cent compared to a model that includes DM variations and spin-noise only