sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Arzoumanian, Z.

Normalized to: Arzoumanian, Z.

131 article(s) in total. 1697 co-authors, from 1 to 53 common article(s). Median position in authors list is 5,0.

[1]  oai:arXiv.org:2006.03074  [pdf] - 2128368
A NICER look at the state transitions of the black hole candidate MAXI J1535-571 during its reflares
Comments: Accepted for publication in MNRAS main journal
Submitted: 2020-06-04, last modified: 2020-07-06
The black hole candidate and X-ray binary MAXI J1535-571 was discovered in September 2017. During the decay of its discovery outburst, and before returning to quiescence, the source underwent at least four reflaring events, with peak luminosities of $\sim$10$^{35-36}$ erg s$^{-1}$ (d/4.1 kpc)$^2$. To investigate the nature of these flares, we analysed a sample of NICER observations taken with almost daily cadence. In this work we present the detailed spectral and timing analysis of the evolution of the four reflares. The higher sensitivity of NICER at lower energies, in comparison with other X-ray detectors, allowed us to constrain the disc component of the spectrum at $\sim$0.5 keV. We found that during each reflare the source appears to trace out a q-shaped track in the hardness-intensity diagram similar to those observed in black hole binaries during full outbursts. MAXI J1535-571 transits between the hard state (valleys) and softer states (peaks) during these flares. Moreover, the Comptonised component is undetected at the peak of the first reflare, while the disc component is undetected during the valleys. Assuming the most likely distance of 4.1 kpc, we find that the hard-to-soft transitions take place at the lowest luminosities ever observed in a black hole transient, while the soft-to-hard transitions occur at some of the lowest luminosities ever reported for such systems.
[2]  oai:arXiv.org:2006.11358  [pdf] - 2119143
A possible polar origin for the FRB associated with a Galactic magnetar
Comments: Submitted to Nature
Submitted: 2020-06-19
Fast Radio Bursts (FRBs) are millisecond-long radio pulses of extragalactic origin with peak luminosities far exceeding any Milky Way sources. The prevalent invocation for the FRB origin involves magnetars: young, magnetically powered neutron stars with the strongest magnetic fields in the Universe. A magnetar-defining signature is the emission of bright, hard X-ray bursts of sub-second duration. These occur in isolation or during a burst storm, when several hundred are observed within minutes to hours. On April 27th 2020, the Galactic magnetar SGR J1935+2154 entered an active period, emitting hundreds of X-ray bursts in a few hours. Remarkably, only one of these temporally coincided with an FRB-like radio burst. Here we report on the spectral and temporal analyses of 24 X-ray bursts emitted 13 hours prior to the FRB and seen simultaneously with NASA NICER and Fermi/GBM missions in their combined energy range. We demonstrate that the FRB-associated X-ray burst is very similar temporally, albeit strikingly different spectrally, from the 24 NICER/GBM bursts. If the FRB-associated burst were drawn from this magnetar burst population, its occurrence rate would be at most around 1 in 7000. This rarity combined with the unusual X-ray burst spectrum is perhaps indicative of an uncommon locale for the origin of the FRB-associated burst. We suggest that this unique event originated in quasi-polar open or closed magnetic field lines extending to high altitudes where radio emission can be generated, possibly from a collimated plasma flow.
[3]  oai:arXiv.org:2005.06495  [pdf] - 2107873
The NANOGrav 12.5-year Data Set: Wideband Timing of 47 Millisecond Pulsars
Comments: 64 pages, 55 figures, 6 tables, 3 appendices. \n Submitted to ApJS. Data available at http://nanograv.org/data/. Replacement: v1 authorship list caused the paper to be incorrectly referenced as Faisal Alam et al.; it should now show up as Alam et al
Submitted: 2020-05-13, last modified: 2020-06-04
We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.
[4]  oai:arXiv.org:2005.06490  [pdf] - 2107872
The NANOGrav 12.5-year Data Set: Observations and Narrowband Timing of 47 Millisecond Pulsars
Comments: 53 pages, 52 figures, 7 tables, 1 appendix. Submitted to The Astrophysical Journal Supplement. Data available at http://nanograv.org/data/. Replacement: v1 authorship list caused the paper to be incorrectly referenced as Faisal Alam et al.; it should now be corrected referenced as Alam et al
Submitted: 2020-05-13, last modified: 2020-06-04
We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.
[5]  oai:arXiv.org:2005.08410  [pdf] - 2096528
Bright X-ray and Radio Pulses from a Recently Reactivated Magnetar
Comments: 54 pages, 8 figures, Submitted for publication in Nature Astronomy
Submitted: 2020-05-17
Magnetars are young, rotating neutron stars that possess larger magnetic fields ($B$ $\approx$ $10^{13}$-$10^{15}$ G) and longer rotational periods ($P$ $\approx$ 1-12 s) than ordinary pulsars. In contrast to rotation-powered pulsars, magnetar emission is thought to be fueled by the evolution and decay of their powerful magnetic fields. They display highly variable radio and X-ray emission, but the processes responsible for this behavior remain a mystery. We report the discovery of bright, persistent individual X-ray pulses from XTE J1810-197, a transient radio magnetar, using the Neutron star Interior Composition Explorer (NICER) following its recent radio reactivation. Similar behavior has only been previously observed from a magnetar during short time periods following a giant flare. However, the X-ray pulses presented here were detected outside of a flaring state. They are less energetic and display temporal structure that differs from the impulsive X-ray events previously observed from the magnetar class, such as giant flares and short X-ray bursts. Our high frequency radio observations of the magnetar, carried out simultaneously with the X-ray observations, demonstrate that the relative alignment between the X-ray and radio pulses varies on rotational timescales. No correlation was found between the amplitudes or temporal structure of the X-ray and radio pulses. The magnetar's 8.3 GHz radio pulses displayed frequency structure, which was not observed in the pulses detected simultaneously at 31.9 GHz. Many of the radio pulses were also not detected simultaneously at both frequencies, which indicates that the underlying emission mechanism producing these pulses is not broadband. We find that the radio pulses from XTE J1810-197 share similar characteristics to radio bursts detected from fast radio burst (FRB) sources, some of which are now thought to be produced by active magnetars.
[6]  oai:arXiv.org:2005.07123  [pdf] - 2095173
Multi-Messenger Gravitational Wave Searches with Pulsar Timing Arrays: Application to 3C66B Using the NANOGrav 11-year Data Set
Comments: 14 pages, 6 figures. Submitted to ApJ
Submitted: 2020-05-14
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.
[7]  oai:arXiv.org:1912.00482  [pdf] - 2076756
The NANOGrav 11-year Data Set: Constraints on Planetary Masses Around 45 Millisecond Pulsars
Comments: Revised and accepted by ApJ Letters
Submitted: 2019-12-01, last modified: 2020-03-24
We search for extrasolar planets around millisecond pulsars using pulsar timing data and seek to determine the minimum detectable planetary masses as a function of orbital period. Using the 11-year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), we look for variations from our models of pulse arrival times due to the presence of exoplanets. No planets are detected around the millisecond pulsars in the NANOGrav 11-year data set, but taking into consideration the noise levels of each pulsar and the sampling rate of our observations, we develop limits that show we are sensitive to planetary masses as low as that of the moon. We analyzed potential planet periods, P, in the range 7 days < P < 2000 days, with somewhat smaller ranges for some binary pulsars. The planetary mass limit for our median-sensitivity pulsar within this period range is 1 M_moon (P / 100 days)^(-2/3).
[8]  oai:arXiv.org:2003.01012  [pdf] - 2065459
A rapid change in X-ray variability and a jet ejection in the black hole transient MAXI J1820+070
Comments: Accepted for publication in ApJ Letters
Submitted: 2020-03-02
We present Neutron Star Interior Composition Explorer X-ray and Arcminute Microkelvin Imager Large Array radio observations of a rapid hard-to-soft state transition in the black hole X-ray transient MAXI J1820+070. During the transition from the hard state to the soft state a switch between two particular types of quasiperiodic oscillations (QPOs) was seen in the X-ray power density spectra, from type-C to type-B, along with a drop in the strength of the broadband X-ray variability and a brief flare in the 7-12 keV band. Soon after this switch (~1.5-2.5 hr) a strong radio flare was observed that corresponded to the launch of superluminal ejecta. Although hints of a connection between QPO transitions and radio flares have been seen in other black hole X-ray transients, our observations constitute the strongest observational evidence to date for a link between the appearance of type-B QPOs and the launch of discrete jet ejections.
[9]  oai:arXiv.org:2002.09691  [pdf] - 2057945
The Soft State of the Black Hole Transient Source MAXI J1820+070: Emission from the Edge of the Plunge Region?
Comments: 8 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2020-02-22
The Galactic black hole X-ray binary MAXI J1820+070 had a bright outburst in 2018 when it became the second brightest X-ray source in the Sky. It was too bright for X-ray CCD instruments such as XMM-Newton and Chandra, but was well observed by photon-counting instruments such as NICER and NuSTAR. We report here on the discovery of an excess emission component during the soft state. It is best modelled with a blackbody spectrum in addition to the regular disk emission, modelled either as diskbb or kerrbb. Its temperature varies from about 0.9 to 1.1 keV which is about 30 to 80 per cent higher than the inner disc temperature of diskbb. Its flux varies between 4 and 12 percent of the disc flux. Simulations of magnetised accretion discs have predicted the possibility of excess emission associated with a non-zero torque at the Innermost Stable Circular Orbit (ISCO) about the black hole, which from other NuSTAR studies lies at about 5 gravitational radii or about 60 km (for a black hole mass is 8 M). In this case the emitting region at the ISCO has a width varying between 1.3 and 4.6 km and would encompass the start of the plunge region where matter begins to fall freely into the black hole.
[10]  oai:arXiv.org:2002.04003  [pdf] - 2050444
A strongly changing accretion morphology during the outburst decay of the neutron star X-ray binary 4U 1608-52
Comments: Accepted for publication in MNRAS
Submitted: 2020-02-10
It is commonly assumed that the properties and geometry of the accretion flow in transient low-mass X-ray binaries (LMXBs) significantly change when the X-ray luminosity decays below $\sim 10^{-2}$ of the Eddington limit ($L_{\rm Edd}$). However, there are few observational cases where the evolution of the accretion flow is tracked in a single X-ray binary over a wide dynamic range. In this work, we use NuSTAR and NICER observations obtained during the 2018 accretion outburst of the neutron star LMXB 4U 1608-52, to study changes in the reflection spectrum. We find that the broad Fe-K$\alpha$ line and Compton hump, clearly seen during the peak of the outburst when the X-ray luminosity is $\sim 10^{37}$ erg/s ($\sim 0.05$ $L_{\rm Edd}$), disappear during the decay of the outburst when the source luminosity drops to $\sim 4.5 \times 10^{35}$ erg/s ($\sim 0.002$ $L_{\rm Edd}$). We show that this non-detection of the reflection features cannot be explained by the lower signal-to-noise at lower flux, but is instead caused by physical changes in the accretion flow. Simulating synthetic NuSTAR observations on a grid of inner disk radius, disk ionisation, and reflection fraction, we find that the disappearance of the reflection features can be explained by either increased disk ionisation ($\log \xi \geq 4.1$) or a much decreased reflection fraction. A changing disk truncation alone, however, cannot account for the lack of reprocessed Fe-K$\alpha$ emission. The required increase in ionisation parameter could occur if the inner accretion flow evaporates from a thin disk into a geometrically thicker flow, such as the commonly assumed formation of an radiatively inefficient accretion flow at lower mass accretion rates.
[11]  oai:arXiv.org:1910.04124  [pdf] - 2046278
The SUrvey for Pulsars and Extragalactic Radio Bursts IV: Discovery and polarimetry of a 12.1-second radio pulsar
Comments: Accepted for publication in MNRAS, 13 pages, 8 figures
Submitted: 2019-10-09, last modified: 2020-01-31
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.
[12]  oai:arXiv.org:2001.11513  [pdf] - 2076872
A NICER View of Spectral and Profile Evolution for Three X-ray Emitting Millisecond Pulsars
Comments: 16 pages, 13 figures, accepted for publication in ApJ
Submitted: 2020-01-30
We present two years of Neutron star Interior Composition Explorer (NICER) X-ray observations of three energetic rotation-powered millisecond pulsars (MSPs): PSRs B1937+21, B1821-24, and J0218+4232. We fit Gaussians and Lorentzians to the pulse profiles for different energy sub-bands of the soft X-ray regime to measure the energy dependence of pulse separation and width. We find that the separation between pulse components of PSR J0218+4232 decreases with increasing energy at $\gt 3\sigma$ confidence. The 95% upper limit on pulse separation evolution for PSRs B1937+21 and B1821-24 is less than 2 milliperiods per keV. Our phase-resolved spectral results provide updated constraints on the non-thermal X-ray emission of these three pulsars. The photon indices of the modeled X-ray emission spectra for each pulse component of PSR B1937+21 are inconsistent with each other at the 90% confidence level, suggesting different emission origins for each pulse. We find that the PSR B1821-24 and PSR J0218+4232 emission spectra are invariant with phase at the 90% confidence level. We describe the implications of our profile and spectral results in the context of equatorial current sheet emission models for these three MSPs with non-thermal, magnetospheric X-ray emission.
[13]  oai:arXiv.org:2001.00595  [pdf] - 2085100
Modeling the uncertainties of solar-system ephemerides for robust gravitational-wave searches with pulsar timing arrays
Comments: Fixed typo in author list. Code that supports all calculations and figures is available at github.com/nanograv/11yr_stochastic_analysis/tree/master/bayesephem
Submitted: 2020-01-02, last modified: 2020-01-06
The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.
[14]  oai:arXiv.org:1912.05703  [pdf] - 2013922
A NICER view of PSR J0030+0451: Implications for the dense matter equation of state
Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State, 17 pages, 5 figures
Submitted: 2019-12-11
Both the mass and radius of the millisecond pulsar PSR J0030+0451 have been inferred via pulse-profile modeling of X-ray data obtained by NASA's NICER mission. In this Letter we study the implications of the mass-radius inference reported for this source by Riley et al. (2019) for the dense matter equation of state (EOS), in the context of prior information from nuclear physics at low densities. Using a Bayesian framework we infer central densities and EOS properties for two choices of high-density extensions: a piecewise-polytropic model and a model based on assumptions of the speed of sound in dense matter. Around nuclear saturation density these extensions are matched to an EOS uncertainty band obtained from calculations based on chiral effective field theory interactions, which provide a realistic description of atomic nuclei as well as empirical nuclear matter properties within uncertainties. We further constrain EOS expectations with input from the current highest measured pulsar mass; together, these constraints offer a narrow Bayesian prior informed by theory as well as laboratory and astrophysical measurements. The NICER mass-radius likelihood function derived by Riley et al. (2019) using pulse-profile modeling is consistent with the highest-density region of this prior. The present relatively large uncertainties on mass and radius for PSR J0030+0451 offer, however, only a weak posterior information gain over the prior. We explore the sensitivity to the inferred geometry of the heated regions that give rise to the pulsed emission, and find a small increase in posterior gain for an alternative (but less preferred) model. Lastly, we investigate the hypothetical scenario of increasing the NICER exposure time for PSR J0030+0451.
[15]  oai:arXiv.org:1912.05708  [pdf] - 2015849
NICER X-ray Observations of Seven Nearby Rotation-Powered Millisecond Pulsars
Comments: 25 pages, 11 tables, 4 figures. In press in The Astrophysical Journal Letters
Submitted: 2019-12-11
NICER observed several rotation-powered millisecond pulsars to search for or confirm the presence of X-ray pulsations. When broad and sine-like, these pulsations may indicate thermal emission from hot polar caps at the magnetic poles on the neutron star surface. We report confident detections ($\ge4.7\sigma$ after background filtering) of X-ray pulsations for five of the seven pulsars in our target sample: PSR J0614-3329, PSR J0636+5129, PSR J0751+1807, PSR J1012+5307, and PSR J2241-5236, while PSR J1552+5437 and PSR J1744-1134 remain undetected. Of those, only PSR J0751+1807 and PSR J1012+5307 had pulsations previously detected at the 1.7$\sigma$ and almost 3$\sigma$ confidence levels, respectively, in XMM-Newton data. All detected sources exhibit broad sine-like pulses, which are indicative of surface thermal radiation. As such, these MSPs are promising targets for future X-ray observations aimed at constraining the neutron star mass-radius relation and the dense matter equation of state using detailed pulse profile modeling. Furthermore, we find that three of the detected millisecond pulsars exhibit a significant phase offset between their X-ray and radio pulses.
[16]  oai:arXiv.org:1912.05705  [pdf] - 2015848
PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter
Comments: 49 pages, 16 figures, part of The Astrophysical Journal Letters focus issue on the Neutron Star Interior Composition Explorer
Submitted: 2019-12-11
Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.
[17]  oai:arXiv.org:1912.05707  [pdf] - 2013926
Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. II. Emission from Hot Spots on a Rapidly Rotating Neutron Star
Comments: 31 pages, 14 figures; accepted for publication in the Astrophysical Journal Letters
Submitted: 2019-12-11
We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.
[18]  oai:arXiv.org:1912.05706  [pdf] - 2013925
Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. I. The Millisecond Pulsar X-ray Data Set
Comments: 28 pages, 9 figures; accepted for publication in the Astrophysical Journal Letters
Submitted: 2019-12-11
We present the set of deep Neutron Star Interior Composition Explorer (NICER) X-ray timing observations of the nearby rotation-powered millisecond pulsars PSRs J0437-4715, J0030+0451, J1231-1411, and J2124-3358, selected as targets for constraining the mass-radius relation of neutron stars and the dense matter equation of state via modeling of their pulsed thermal X-ray emission. We describe the instrument, observations, and data processing/reduction procedures, as well as the series of investigations conducted to ensure that the properties of the data sets are suitable for parameter estimation analyses to produce reliable constraints on the neutron star mass-radius relation and the dense matter equation of state. We find that the long-term timing and flux behavior and the Fourier-domain properties of the event data do not exhibit any anomalies that could adversely affect the intended measurements. From phase-selected spectroscopy, we find that emission from the individual pulse peaks is well described by a single-temperature hydrogen atmosphere spectrum, with the exception of PSR J0437-4715, for which multiple temperatures are required.
[19]  oai:arXiv.org:1912.05704  [pdf] - 2013923
A NICER view of PSR J0030+0451: evidence for a global-scale multipolar magnetic field
Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State
Submitted: 2019-12-11
Recent modeling of NICER observations of thermal X-ray pulsations from the surface of the isolated millisecond pulsar PSR J0030+0451 suggests that the hot emitting regions on the pulsar's surface are far from antipodal, which is at odds with the classical assumption that the magnetic field in the pulsar magnetosphere is predominantly that of a centered dipole. Here, we review these results and examine previous attempts to constrain the magnetospheric configuration of PSR J0030+0451. To the best of our knowledge, there is in fact no direct observational evidence that PSR J0030+0451's magnetic field is a centered dipole. Developing models of physically motivated, non-canonical magnetic field configurations and the currents that they can support poses a challenging task. However, such models may have profound implications for many aspects of pulsar research, including pulsar braking, estimates of birth velocities, and interpretations of multi-wavelength magnetospheric emission.
[20]  oai:arXiv.org:1912.05702  [pdf] - 2013921
A NICER View of PSR J0030+0451: Millisecond Pulsar Parameter Estimation
Comments: Appears in ApJ Letters Focus Issue on NICER Constraints on the Dense Matter Equation of State; 76 pages, 24 figures, 7 tables, 8 figure sets (available in the online journal or from the authors)
Submitted: 2019-12-11
We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSR J0030$+$0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer (NICER) X-ray spectral-timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar's surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply-connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass $M$ and equatorial radius $R_\mathrm{eq}$ are, respectively, $1.34_{-0.16}^{+0.15}$ M$_{\odot}$ and $12.71_{-1.19}^{+1.14}$ km, whilst the compactness $GM/R_\mathrm{eq}c^2 = 0.156_{-0.010}^{+0.008}$ is more tightly constrained; the credible interval bounds reported here are approximately the $16\%$ and $84\%$ quantiles in marginal posterior mass.
[21]  oai:arXiv.org:1911.08488  [pdf] - 2042263
The NANOGrav 11-Year Data Set: Limits on Gravitational Wave Memory
Comments: 10 pages, 6 figures, submitted to ApJ
Submitted: 2019-11-19, last modified: 2019-12-06
The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.
[22]  oai:arXiv.org:1911.11829  [pdf] - 2026447
Deciphering the Nature of the Pulsar Wind Nebula CTB 87 with XMM-Newton
Comments: 9 pages, 10 figures. Accepted for publication in MNRAS
Submitted: 2019-11-26
CTB 87 (G74.9+1.2) is an evolved supernova remnant (SNR) which hosts a peculiar pulsar wind nebula (PWN). The X-ray peak is offset from that observed in radio and lies towards the edge of the radio nebula. The putative pulsar, CXOU~J201609.2+371110, was first resolved with \textit{Chandra} and is surrounded by a compact and a more extended X-ray nebula. Here we use a deep {\textit{XMM-Newton}} observation to examine the morphology and evolutionary stage of the PWN and to search for thermal emission expected from a supernova shell or reverse shock interaction with supernova ejecta. We do not find evidence of thermal X-ray emission from the SNR and place an upper limit on the electron density of 0.05~cm$^{-3}$ for a plasma temperature $kT\sim 0.8$ keV. The morphology and spectral properties are consistent with a $\sim$20~kyr-old relic PWN expanding into a stellar wind-blown bubble. We also present the first X-ray spectral index map from the PWN and show that we can reproduce its morphology by means of 2D axisymmetric relativistic hydrodynamical simulations.
[23]  oai:arXiv.org:1910.13628  [pdf] - 2026196
A pulsar-based timescale from the International Pulsar Timing Array
Comments: accepted by MNRAS
Submitted: 2019-10-29
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.
[24]  oai:arXiv.org:1910.04174  [pdf] - 1977454
A Black Hole X-ray Binary at $\sim$100 Hz: Multiwavelength Timing of MAXI J1820+070 with HiPERCAM and NICER
Comments: MNRAS Letters, In Press
Submitted: 2019-10-09
We report on simultaneous sub-second optical and X-ray timing observations of the low mass X-ray binary black hole candidate MAXI J1820+070. The bright 2018 outburst rise allowed simultaneous photometry in five optical bands ($ugriz_s$) with HiPERCAM/GTC (Optical) at frame rates over 100 Hz, together with NICER/ISS observations (X-rays). Intense (factor of two) red flaring activity in the optical is seen over a broad range of timescales down to $\sim$10 ms. Cross-correlating the bands reveals a prominent anti-correlation on timescales of $\sim$seconds, and a narrow sub-second correlation at a lag of $\approx$+165 ms (optical lagging X-rays). This lag increases with optical wavelength, and is approximately constant over Fourier frequencies of $\sim$0.3-10 Hz. These features are consistent with an origin in the inner accretion flow and jet base within $\sim$5000 Gravitational radii. An additional $\sim$+5 s lag feature may be ascribable to disc reprocessing. MAXI J1820+070 is the third black hole transient to display a clear $\sim$0.1s optical lag, which may be common feature in such objects. The sub-second lag $variation$ with wavelength is novel, and may allow constraints on internal shock jet stratification models.
[25]  oai:arXiv.org:1910.03062  [pdf] - 1975523
Timing the pulsations of the accreting millisecond pulsar SAX J1808.4-3658 during its 2019 outburst
Comments: 12 pages, 7 figures, 2 tables. Submitted to ApJ
Submitted: 2019-10-07
In this paper we present a coherent timing analysis of the 401 Hz pulsations of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its 2019 outburst. Using observations collected with the Neutron Star Interior Composition Explorer (NICER), we establish the pulsar spin frequency and orbital phase during its latest epoch. We find that the 2019 outburst shows a pronounced evolution in pulse phase over the course of the outburst. These phase shifts are found to correlate with the source flux, and are interpreted in terms of hot-spot drift on the stellar surface, driven by changes in the mass accretion rate. Additionally, we find that the long-term evolution of the pulsar spin frequency shows evidence for a modulation at the Earth's orbital period, enabling pulsar timing based astrometry of this accreting millisecond pulsar.
[26]  oai:arXiv.org:1909.08644  [pdf] - 2054030
The NANOGrav 11-Year Data Set: Evolution of Gravitational Wave Background Statistics
Comments: 14 pages, 13 figures, fixed typo in abstract of earlier version
Submitted: 2019-09-18, last modified: 2019-09-20
An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.
[27]  oai:arXiv.org:1909.07338  [pdf] - 1986841
An evolving broad iron line from the first Galactic ultraluminous X-ray pulsar Swift J0243.6+6124
Comments: 11 pages, 8 Figures, Accepted for publication in ApJ
Submitted: 2019-09-16
We present a spectral study of the ultraluminous Be/X-ray transient pulsar Swift J0243.6+6124 using Neutron Star Interior Composition Explorer (NICER) observations during the system's 2017--2018 giant outburst. The 1.2--10~keV energy spectrum of the source can be approximated with an absorbed cut-off power law model. We detect strong, luminosity-dependent emission lines in the 6--7 keV energy range. A narrow 6.42 keV line, observed in the sub-Eddington regime, is seen to evolve into a broad Fe-line profile in the super-Eddington regime. Other features are found at 6.67 and 6.97 keV in the Fe-line complex. An asymmetric broad line profile, peaking at 6.67 keV, is possibly due to Doppler effects and gravitational redshift. The 1.2--79 keV broadband spectrum from NuSTAR and NICER observations at the outburst peak is well described by an absorbed cut-off power law plus multiple Gaussian lines and a blackbody component. Physical reflection models are also tested to probe the broad iron line feature. Depending on the mass accretion rate, we found emission sites that are evolving from ~5000 km to a range closer to the surface of the neutron star. Our findings are discussed in the framework of the accretion disk and its implication on the magnetic field, the presence of optically thick accretion curtain in the magnetosphere, jet emission, and the massive, ultra-fast outflow expected at super-Eddington accretion rates. We do not detect any signatures of a cyclotron absorption line in the NICER or NuSTAR data.
[28]  oai:arXiv.org:1904.06759  [pdf] - 1961559
Relativistic Shapiro delay measurements of an extremely massive millisecond pulsar
Comments: 12 pages, 3 figures, 1 table, published in Nature Astronomy (2019)
Submitted: 2019-04-14, last modified: 2019-09-13
Despite its importance to our understanding of physics at supranuclear densities, the equation of state (EoS) of matter deep within neutron stars remains poorly understood. Millisecond pulsars (MSPs) are among the most useful astrophysical objects in the Universe for testing fundamental physics, and place some of the most stringent constraints on this high-density EoS. Pulsar timing - the process of accounting for every rotation of a pulsar over long time periods - can precisely measure a wide variety of physical phenomena, including those that allow the measurement of the masses of the components of a pulsar binary system (Lorimer & Kramer 2005). One of these, called relativistic Shapiro delay (Shapiro 1964), can yield precise masses for both an MSP and its companion; however, it is only easily observed in a small subset of high-precision, highly inclined (nearly edge-on) binary pulsar systems. By combining data from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 12.5-year data set with recent orbital-phase-specific observations using the Green Bank Telescope, we have measured the mass of the MSP J0740+6620 to be $\mathbf{2.14^{+0.10}_{-0.09}}$ solar masses (68.3% credibility interval; 95.4% credibility interval is $\mathbf{2.14^{+0.20}_{-0.18}}$ solar masses). It is highly likely to be the most massive neutron star yet observed, and serves as a strong constraint on the neutron star interior EoS.
[29]  oai:arXiv.org:1909.04534  [pdf] - 1983956
The International Pulsar Timing Array: Second data release
Comments: Submitted to MNRAS and in review, 23 pages, 5 figures
Submitted: 2019-09-10
In this paper, we describe the International Pulsar Timing Array second data release, which includes recent pulsar timing data obtained by three regional consortia: the European Pulsar Timing Array, the North American Nanohertz Observatory for Gravitational Waves, and the Parkes Pulsar Timing Array. We analyse and where possible combine high-precision timing data for 65 millisecond pulsars which are regularly observed by these groups. A basic noise analysis, including the processes which are both correlated and uncorrelated in time, provides noise models and timing ephemerides for the pulsars. We find that the timing precisions of pulsars are generally improved compared to the previous data release, mainly due to the addition of new data in the combination. The main purpose of this work is to create the most up-to-date IPTA data release. These data are publicly available for searches for low-frequency gravitational waves and other pulsar science.
[30]  oai:arXiv.org:1909.03595  [pdf] - 1994193
A NICER thermonuclear burst from the millisecond X-ray pulsar SAX J1808.4-3658
Comments: 10 pages, 5 figures, submitted to ApJ Letters
Submitted: 2019-09-08
The Neutron Star Interior Composition Explorer (NICER) has extensively monitored the August 2019 outburst of the 401 Hz millisecond X-ray pulsar SAX J1808.4-3658. In this Letter, we report on the detection of a bright helium-fueled Type I X-ray burst. With a bolometric peak flux of $(2.3\pm0.1)\times 10^{-7}$ erg/cm^2/s, this was the brightest X-ray burst among all bursting sources observed with NICER to date. The burst shows a remarkable two-stage evolution in flux, emission lines at $1.0$ keV and $6.7$ keV, and burst oscillations at the known pulsar spin frequency, with $\approx4$\% fractional sinusoidal amplitude. We interpret the burst flux evolution as the detection of the local Eddington limits associated with the hydrogen and helium layers of the neutron star envelope. The emission lines are likely associated with Fe, due to reprocessing of the burst emission in the accretion disk.
[31]  oai:arXiv.org:1908.08023  [pdf] - 1994140
A NICER Look at Strong X-ray Obscuration in the Seyfert-2 Galaxy NGC 4388
Comments: Accepted for publication in ApJ
Submitted: 2019-08-21
We present an analysis of the time-averaged spectrum of the Seyfert-2 active galaxy NGC 4388, obtained by NICER. The intrinsic strength of the reflection spectrum in NGC 4388, the large collecting area and favorable pass band of NICER, and a net exposure of 105.6 ks yielded an exceptionally sensitive spectrum. Using two independent families of models, the intrinsic spectrum from the central engine is found to be highly obscured but not Compton-thick. Enforcing physical self-consistency within each model, the independent treatments give formally consistent results: N_H = 2.67 (-0.03,+0.02) E+23 cm^-2 or N_H = 2.64 (-0.03, +0.03) E+23 cm^-2. Past measurements made with Suzaku and XMM-Newton are in broad agreement with these column density values. A more recent measurement with NuSTAR (in late 2013) recorded a column density about twice as large; the robustness of this variability is reinforced by the use of consistent models and procedures. The neutral Fe K-alpha line in the NICER spectrum is nominally resolved and consistent with an origin in the optical broad line region (BLR). The data also require ionized absorption in the Fe K band, similar to the "warm absorbers" detected in Seyfert-1 active galactic nuclei (AGN). The low-energy spectrum is consistent with a set of ionized plasma components. We discuss these findings and note that the geometric inferences that derive from this analysis can be tested with XRISM and Athena.
[32]  oai:arXiv.org:1908.03373  [pdf] - 1964698
NICER observes a secondary peak in the decay of a thermonuclear burst from 4U 1608-52
Comments: 8 pages, 6 figures, Accepted for publication in ApJ
Submitted: 2019-08-09
We report for the first time below 1.5 keV, the detection of a secondary peak in an Eddington-limited thermonuclear X-ray burst observed by the Neutron Star Interior Composition Explorer (NICER) from the low-mass X-ray binary 4U 1608-52. Our time-resolved spectroscopy of the burst is consistent with a model consisting of a varying-temperature blackbody, and an evolving persistent flux contribution, likely attributed to the accretion process. The dip in the burst intensity before the secondary peak is also visible in the bolometric flux. Prior to the dip, the blackbody temperature reached a maximum of $\approx3$ keV. Our analysis suggests that the dip and secondary peak are not related to photospheric expansion, varying circumstellar absorption, or scattering. Instead, we discuss the observation in the context of hydrodynamical instabilities, thermonuclear flame spreading models, and re-burning in the cooling tail of the burst.
[33]  oai:arXiv.org:1908.01206  [pdf] - 1928259
NICER Observation of Unusual Burst Oscillations in 4U 1728-34
Comments: 14 pages, 6 figures, published in ApJ 878, 2, 145
Submitted: 2019-08-03
The Neutron Star Interior Composition Explorer (NICER) has observed seven thermonuclear X-ray bursts from the Low Mass X-ray Binary (LMXB) neutron star 4U 1728-34 from the start of the mission's operations until February of 2019. Three of these bursts show oscillations in their decaying tail with frequencies that are within 1 Hz of the previously detected burst oscillations from this source. Two of these burst oscillations have unusual properties: They have large fractional rms amplitudes of $ 48 \pm 9 \%$ and $ 46 \pm 9 \%$, and they are detected only at photon energies above 6 keV. By contrast, the third detected burst oscillation is compatible with previous observations of this source, with a fractional rms amplitude of $7.7 \pm 1.5\%$ rms in the 0.3 to 6.2 keV energy band. We discuss the implications of these large-amplitude burst oscillations, finding they are difficult to explain with the current theoretical models for X-ray burst tail oscillations.
[34]  oai:arXiv.org:1908.00539  [pdf] - 1966922
Observations of the Ultra-compact X-ray Binary 4U 1543-624 in Outburst with NICER, INTEGRAL, Swift, and ATCA
Comments: 13 pages, 9 figures, 2 tables, accepted for publication in ApJ
Submitted: 2019-08-01
We report on X-ray and radio observations of the ultra-compact X-ray binary 4U 1543-624 taken in August 2017 during an enhanced accretion episode. We obtained NICER monitoring of the source over a $\sim10$ day period during which target-of-opportunity observations were also conducted with Swift, INTEGRAL, and ATCA. Emission lines were measured in the NICER X-ray spectrum at $\sim0.64$ keV and $\sim6.4$ keV that correspond to O and Fe, respectively. By modeling these line components, we are able to track changes in the accretion disk throughout this period. The innermost accretion flow appears to move inwards from hundreds of gravitational radii ($R_{g}=GM/c^{2}$) at the beginning of the outburst to $<8.7$ $R_{g}$ at peak intensity. We do not detect the source in radio, but are able to place a $3\sigma$ upper limit on the flux density at $27$ $\mu$Jy beam$^{-1}$. Comparing the radio and X-ray luminosities, we find that the source lies significantly away from the range typical of black holes in the ${L}_{{r}}$-${L}_{{x}}$ plane, suggesting a neutron star (NS) primary. This adds to the evidence that NSs do not follow a single track in the ${L}_{{r}}$-${L}_{{x}}$ plane, limiting its use in distinguishing between different classes of NSs based on radio and X-ray observations alone.
[35]  oai:arXiv.org:1902.08507  [pdf] - 1907525
Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015-2017 LIGO Data
The LIGO Scientific Collaboration; the Virgo Collaboration; Abbott, B. P.; Abbott, R.; Abbott, T. D.; Abraham, S.; Acernese, F.; Ackley, K.; Adams, C.; Adhikari, R. X.; Adya, V. B.; Affeldt, C.; Agathos, M.; Agatsuma, K.; Aggarwal, N.; Aguiar, O. D.; Aiello, L.; Ain, A.; Ajith, P.; Allen, G.; Allocca, A.; Aloy, M. A.; Altin, P. A.; Amato, A.; Ananyeva, A.; Anderson, S. B.; Anderson, W. G.; Angelova, S. V.; Antier, S.; Appert, S.; Arai, K.; Araya, M. C.; Areeda, J. S.; Arène, M.; Arnaud, N.; Ascenzi, S.; Ashton, G.; Aston, S. M.; Astone, P.; Aubin, F.; Aufmuth, P.; AultONeal, K.; Austin, C.; Avendano, V.; Avila-Alvarez, A.; Babak, S.; Bacon, P.; Badaracco, F.; Bader, M. K. M.; Bae, S.; Baker, P. T.; Baldaccini, F.; Ballardin, G.; Ballmer, S. W.; Banagiri, S.; Barayoga, J. C.; Barclay, S. E.; Barish, B. C.; Barker, D.; Barkett, K.; Barnum, S.; Barone, F.; Barr, B.; Barsotti, L.; Barsuglia, M.; Barta, D.; Bartlett, J.; Bartos, I.; Bassiri, R.; Basti, A.; Bawaj, M.; Bayley, J. C.; Bazzan, M.; Bécsy, B.; Bejger, M.; Belahcene, I.; Bell, A. S.; Beniwal, D.; Berger, B. K.; Bergmann, G.; Bernuzzi, S.; Bero, J. J.; Berry, C. P. L.; Bersanetti, D.; Bertolini, A.; Betzwieser, J.; Bhandare, R.; Bidler, J.; Bilenko, I. A.; Bilgili, S. A.; Billingsley, G.; Birch, J.; Birney, R.; Birnholtz, O.; Biscans, S.; Biscoveanu, S.; Bisht, A.; Bitossi, M.; Bizouard, M. A.; Blackburn, J. K.; Blair, C. D.; Blair, D. G.; Blair, R. M.; Bloemen, S.; Bode, N.; Boer, M.; Boetzel, Y.; Bogaert, G.; Bondu, F.; Bonilla, E.; Bonnand, R.; Booker, P.; Boom, B. A.; Booth, C. D.; Bork, R.; Boschi, V.; Bose, S.; Bossie, K.; Bossilkov, V.; Bosveld, J.; Bouffanais, Y.; Bozzi, A.; Bradaschia, C.; Brady, P. R.; Bramley, A.; Branchesi, M.; Brau, J. E.; Briant, T.; Briggs, J. H.; Brighenti, F.; Brillet, A.; Brinkmann, M.; Brisson, V.; Brockill, P.; Brooks, A. F.; Brown, D. D.; Brunett, S.; Buikema, A.; Bulik, T.; Bulten, H. J.; Buonanno, A.; Buskulic, D.; Buy, C.; Byer, R. L.; Cabero, M.; Cadonati, L.; Cagnoli, G.; Cahillane, C.; Bustillo, J. Calderón; Callister, T. A.; Calloni, E.; Camp, J. B.; Campbell, W. A.; Canepa, M.; Cannon, K. C.; Cao, H.; Cao, J.; Capocasa, E.; Carbognani, F.; Caride, S.; Carney, M. F.; Carullo, G.; Diaz, J. Casanueva; Casentini, C.; Caudill, S.; Cavaglià, M.; Cavalier, F.; Cavalieri, R.; Cella, G.; Cerdá-Durán, P.; Cerretani, G.; Cesarini, E.; Chaibi, O.; Chakravarti, K.; Chamberlin, S. J.; Chan, M.; Chao, S.; Charlton, P.; Chase, E. A.; Chassande-Mottin, E.; Chatterjee, D.; Chaturvedi, M.; Cheeseboro, B. D.; Chen, H. Y.; Chen, X.; Chen, Y.; Cheng, H. -P.; Cheong, C. K.; Chia, H. Y.; Chincarini, A.; Chiummo, A.; Cho, G.; Cho, H. S.; Cho, M.; Christensen, N.; Chu, Q.; Chua, S.; Chung, K. W.; Chung, S.; Ciani, G.; Ciobanu, A. A.; Ciolfi, R.; Cipriano, F.; Cirone, A.; Clara, F.; Clark, J. A.; Clearwater, P.; Cleva, F.; Cocchieri, C.; Coccia, E.; Cohadon, P. -F.; Cohen, D.; Colgan, R.; Colleoni, M.; Collette, C. G.; Collins, C.; Cominsky, L. R.; Constancio, M.; Conti, L.; Cooper, S. J.; Corban, P.; Corbitt, T. R.; Cordero-Carrión, I.; Corley, K. R.; Cornish, N.; Corsi, A.; Cortese, S.; Costa, C. A.; Cotesta, R.; Coughlin, M. W.; Coughlin, S. B.; Coulon, J. -P.; Countryman, S. T.; Couvares, P.; Covas, P. B.; Cowan, E. E.; Coward, D. M.; Cowart, M. J.; Coyne, D. C.; Coyne, R.; Creighton, J. D. E.; Creighton, T. D.; Cripe, J.; Croquette, M.; Crowder, S. G.; Cullen, T. J.; Cumming, A.; Cunningham, L.; Cuoco, E.; Canton, T. Dal; Dálya, G.; Danilishin, S. L.; D'Antonio, S.; Danzmann, K.; Dasgupta, A.; Costa, C. F. Da Silva; Datrier, L. E. H.; Dattilo, V.; Dave, I.; Davier, M.; Davis, D.; Daw, E. J.; DeBra, D.; Deenadayalan, M.; Degallaix, J.; De Laurentis, M.; Deléglise, S.; Del Pozzo, W.; DeMarchi, L. M.; Demos, N.; Dent, T.; De Pietri, R.; Derby, J.; De Rosa, R.; De Rossi, C.; DeSalvo, R.; de Varona, O.; Dhurandhar, S.; Díaz, M. C.; Dietrich, T.; Di Fiore, L.; Di Giovanni, M.; Di Girolamo, T.; Di Lieto, A.; Ding, B.; Di Pace, S.; Di Palma, I.; Di Renzo, F.; Dmitriev, A.; Doctor, Z.; Donovan, F.; Dooley, K. L.; Doravari, S.; Dorrington, I.; Downes, T. P.; Drago, M.; Driggers, J. C.; Du, Z.; Ducoin, J. -G.; Dupej, P.; Dwyer, S. E.; Easter, P. J.; Edo, T. B.; Edwards, M. C.; Effler, A.; Ehrens, P.; Eichholz, J.; Eikenberry, S. S.; Eisenmann, M.; Eisenstein, R. A.; Essick, R. C.; Estelles, H.; Estevez, D.; Etienne, Z. B.; Etzel, T.; Evans, M.; Evans, T. M.; Fafone, V.; Fair, H.; Fairhurst, S.; Fan, X.; Farinon, S.; Farr, B.; Farr, W. M.; Fauchon-Jones, E. J.; Favata, M.; Fays, M.; Fazio, M.; Fee, C.; Feicht, J.; Fejer, M. M.; Feng, F.; Fernandez-Galiana, A.; Ferrante, I.; Ferreira, E. C.; Ferreira, T. A.; Ferrini, F.; Fidecaro, F.; Fiori, I.; Fiorucci, D.; Fishbach, M.; Fisher, R. P.; Fishner, J. M.; Fitz-Axen, M.; Flaminio, R.; Fletcher, M.; Flynn, E.; Fong, H.; Font, J. A.; Forsyth, P. W. F.; Fournier, J. -D.; Frasca, S.; Frasconi, F.; Frei, Z.; Freise, A.; Frey, R.; Frey, V.; Fritschel, P.; Frolov, V. V.; Fulda, P.; Fyffe, M.; Gabbard, H. A.; Gadre, B. U.; Gaebel, S. M.; Gair, J. R.; Gammaitoni, L.; Ganija, M. R.; Gaonkar, S. G.; Garcia, A.; García-Quirós, C.; Garufi, F.; Gateley, B.; Gaudio, S.; Gaur, G.; Gayathri, V.; Gemme, G.; Genin, E.; Gennai, A.; George, D.; George, J.; Gergely, L.; Germain, V.; Ghonge, S.; Ghosh, Abhirup; Ghosh, Archisman; Ghosh, S.; Giacomazzo, B.; Giaime, J. A.; Giardina, K. D.; Giazotto, A.; Gill, K.; Giordano, G.; Glover, L.; Godwin, P.; Goetz, E.; Goetz, R.; Goncharov, B.; González, G.; Castro, J. M. Gonzalez; Gopakumar, A.; Gorodetsky, M. L.; Gossan, S. E.; Gosselin, M.; Gouaty, R.; Grado, A.; Graef, C.; Granata, M.; Grant, A.; Gras, S.; Grassia, P.; Gray, C.; Gray, R.; Greco, G.; Green, A. C.; Green, R.; Gretarsson, E. M.; Groot, P.; Grote, H.; Grunewald, S.; Gruning, P.; Guidi, G. M.; Gulati, H. K.; Guo, Y.; Gupta, A.; Gupta, M. K.; Gustafson, E. K.; Gustafson, R.; Haegel, L.; Halim, O.; Hall, B. R.; Hall, E. D.; Hamilton, E. Z.; Hammond, G.; Haney, M.; Hanke, M. M.; Hanks, J.; Hanna, C.; Hannam, M. D.; Hannuksela, O. A.; Hanson, J.; Hardwick, T.; Haris, K.; Harms, J.; Harry, G. M.; Harry, I. W.; Haster, C. -J.; Haughian, K.; Hayes, F. J.; Healy, J.; Heidmann, A.; Heintze, M. C.; Heitmann, H.; Hello, P.; Hemming, G.; Hendry, M.; Heng, I. S.; Hennig, J.; Heptonstall, A. W.; Vivanco, Francisco Hernandez; Heurs, M.; Hild, S.; Hinderer, T.; Hoak, D.; Hochheim, S.; Hofman, D.; Holgado, A. M.; Holland, N. A.; Holt, K.; Holz, D. E.; Hopkins, P.; Horst, C.; Hough, J.; Howell, E. J.; Hoy, C. G.; Hreibi, A.; Huerta, E. A.; Huet, D.; Hughey, B.; Hulko, M.; Husa, S.; Huttner, S. H.; Huynh-Dinh, T.; Idzkowski, B.; Iess, A.; Ingram, C.; Inta, R.; Intini, G.; Irwin, B.; Isa, H. N.; Isac, J. -M.; Isi, M.; Iyer, B. R.; Izumi, K.; Jacqmin, T.; Jadhav, S. J.; Jani, K.; Janthalur, N. N.; Jaranowski, P.; Jenkins, A. C.; Jiang, J.; Johnson, D. S.; Jones, A. W.; Jones, D. I.; Jones, R.; Jonker, R. J. G.; Ju, L.; Junker, J.; Kalaghatgi, C. V.; Kalogera, V.; Kamai, B.; Kandhasamy, S.; Kang, G.; Kanner, J. B.; Kapadia, S. J.; Karki, S.; Karvinen, K. S.; Kashyap, R.; Kasprzack, M.; Katsanevas, S.; Katsavounidis, E.; Katzman, W.; Kaufer, S.; Kawabe, K.; Keerthana, N. V.; Kéfélian, F.; Keitel, D.; Kennedy, R.; Key, J. S.; Khalili, F. Y.; Khan, H.; Khan, I.; Khan, S.; Khan, Z.; Khazanov, E. A.; Khursheed, M.; Kijbunchoo, N.; Kim, Chunglee; Kim, J. C.; Kim, K.; Kim, W.; Kim, W. S.; Kim, Y. -M.; Kimball, C.; King, E. J.; King, P. J.; Kinley-Hanlon, M.; Kirchhoff, R.; Kissel, J. S.; Kleybolte, L.; Klika, J. H.; Klimenko, S.; Knowles, T. D.; Koch, P.; Koehlenbeck, S. M.; Koekoek, G.; Koley, S.; Kondrashov, V.; Kontos, A.; Koper, N.; Korobko, M.; Korth, W. Z.; Kowalska, I.; Kozak, D. B.; Kringel, V.; Krishnendu, N.; Królak, A.; Kuehn, G.; Kumar, A.; Kumar, P.; Kumar, R.; Kumar, S.; Kuo, L.; Kutynia, A.; Kwang, S.; Lackey, B. D.; Lai, K. H.; Lam, T. L.; Landry, M.; Lane, B. B.; Lang, R. N.; Lange, J.; Lantz, B.; Lanza, R. K.; Lartaux-Vollard, A.; Lasky, P. D.; Laxen, M.; Lazzarini, A.; Lazzaro, C.; Leaci, P.; Leavey, S.; Lecoeuche, Y. K.; Lee, C. H.; Lee, H. K.; Lee, H. M.; Lee, H. W.; Lee, J.; Lee, K.; Lehmann, J.; Lenon, A.; Leroy, N.; Letendre, N.; Levin, Y.; Li, J.; Li, K. J. L.; Li, T. G. F.; Li, X.; Lin, F.; Linde, F.; Linker, S. D.; Littenberg, T. B.; Liu, J.; Liu, X.; Lo, R. K. L.; Lockerbie, N. A.; London, L. T.; Longo, A.; Lorenzini, M.; Loriette, V.; Lormand, M.; Losurdo, G.; Lough, J. D.; Lousto, C. O.; Lovelace, G.; Lower, M. E.; Lück, H.; Lumaca, D.; Lundgren, A. P.; Lynch, R.; Ma, Y.; Macas, R.; Macfoy, S.; MacInnis, M.; Macleod, D. M.; Macquet, A.; Magaña-Sandoval, F.; Zertuche, L. Magaña; Magee, R. M.; Majorana, E.; Maksimovic, I.; Malik, A.; Man, N.; Mandic, V.; Mangano, V.; Mansell, G. L.; Manske, M.; Mantovani, M.; Marchesoni, F.; Marion, F.; Márka, S.; Márka, Z.; Markakis, C.; Markosyan, A. S.; Markowitz, A.; Maros, E.; Marquina, A.; Marsat, S.; Martelli, F.; Martin, I. W.; Martin, R. M.; Martynov, D. V.; Mason, K.; Massera, E.; Masserot, A.; Massinger, T. J.; Masso-Reid, M.; Mastrogiovanni, S.; Matas, A.; Matichard, F.; Matone, L.; Mavalvala, N.; Mazumder, N.; McCann, J. J.; McCarthy, R.; McClelland, D. E.; McCormick, S.; McCuller, L.; McGuire, S. C.; McIver, J.; McManus, D. J.; McRae, T.; McWilliams, S. T.; Meacher, D.; Meadors, G. D.; Mehmet, M.; Mehta, A. K.; Meidam, J.; Melatos, A.; Mendell, G.; Mercer, R. A.; Mereni, L.; Merilh, E. L.; Merzougui, M.; Meshkov, S.; Messenger, C.; Messick, C.; Metzdorff, R.; Meyers, P. M.; Miao, H.; Michel, C.; Middleton, H.; Mikhailov, E. E.; Milano, L.; Miller, A. L.; Miller, A.; Millhouse, M.; Mills, J. C.; Milovich-Goff, M. C.; Minazzoli, O.; Minenkov, Y.; Mishkin, A.; Mishra, C.; Mistry, T.; Mitra, S.; Mitrofanov, V. P.; Mitselmakher, G.; Mittleman, R.; Mo, G.; Moffa, D.; Mogushi, K.; Mohapatra, S. R. P.; Montani, M.; Moore, C. J.; Moraru, D.; Moreno, G.; Morisaki, S.; Mours, B.; Mow-Lowry, C. M.; Mukherjee, Arunava; Mukherjee, D.; Mukherjee, S.; Mukund, N.; Mullavey, A.; Munch, J.; Muñiz, E. A.; Muratore, M.; Murray, P. G.; Nagar, A.; Nardecchia, I.; Naticchioni, L.; Nayak, R. K.; Neilson, J.; Nelemans, G.; Nelson, T. J. N.; Nery, M.; Neunzert, A.; Ng, K. Y.; Ng, S.; Nguyen, P.; Nichols, D.; Nissanke, S.; Nocera, F.; North, C.; Nuttall, L. K.; Obergaulinger, M.; Oberling, J.; O'Brien, B. D.; O'Dea, G. D.; Ogin, G. H.; Oh, J. J.; Oh, S. H.; Ohme, F.; Ohta, H.; Okada, M. A.; Oliver, M.; Oppermann, P.; Oram, Richard J.; O'Reilly, B.; Ormiston, R. G.; Ortega, L. F.; O'Shaughnessy, R.; Ossokine, S.; Ottaway, D. J.; Overmier, H.; Owen, B. J.; Pace, A. E.; Pagano, G.; Page, M. A.; Pai, A.; Pai, S. A.; Palamos, J. R.; Palashov, O.; Palomba, C.; Pal-Singh, A.; Pan, Huang-Wei; Pang, B.; Pang, P. T. H.; Pankow, C.; Pannarale, F.; Pant, B. C.; Paoletti, F.; Paoli, A.; Parida, A.; Parker, W.; Pascucci, D.; Pasqualetti, A.; Passaquieti, R.; Passuello, D.; Patil, M.; Patricelli, B.; Pearlstone, B. L.; Pedersen, C.; Pedraza, M.; Pedurand, R.; Pele, A.; Penn, S.; Perez, C. J.; Perreca, A.; Pfeiffer, H. P.; Phelps, M.; Phukon, K. S.; Piccinni, O. J.; Pichot, M.; Piergiovanni, F.; Pillant, G.; Pinard, L.; Pirello, M.; Pitkin, M.; Poggiani, R.; Pong, D. Y. T.; Ponrathnam, S.; Popolizio, P.; Porter, E. K.; Powell, J.; Prajapati, A. K.; Prasad, J.; Prasai, K.; Prasanna, R.; Pratten, G.; Prestegard, T.; Privitera, S.; Prodi, G. A.; Prokhorov, L. G.; Puncken, O.; Punturo, M.; Puppo, P.; Pürrer, M.; Qi, H.; Quetschke, V.; Quinonez, P. J.; Quintero, E. A.; Quitzow-James, R.; Raab, F. J.; Radkins, H.; Radulescu, N.; Raffai, P.; Raja, S.; Rajan, C.; Rajbhandari, B.; Rakhmanov, M.; Ramirez, K. E.; Ramos-Buades, A.; Rana, Javed; Rao, K.; Rapagnani, P.; Raymond, V.; Razzano, M.; Read, J.; Regimbau, T.; Rei, L.; Reid, S.; Reitze, D. H.; Ren, W.; Ricci, F.; Richardson, C. J.; Richardson, J. W.; Ricker, P. M.; Riles, K.; Rizzo, M.; Robertson, N. A.; Robie, R.; Robinet, F.; Rocchi, A.; Rolland, L.; Rollins, J. G.; Roma, V. J.; Romanelli, M.; Romano, R.; Romel, C. L.; Romie, J. H.; Rose, K.; Rosińska, D.; Rosofsky, S. G.; Ross, M. P.; Rowan, S.; Rüdiger, A.; Ruggi, P.; Rutins, G.; Ryan, K.; Sachdev, S.; Sadecki, T.; Sakellariadou, M.; Salconi, L.; Saleem, M.; Samajdar, A.; Sammut, L.; Sanchez, E. J.; Sanchez, L. E.; Sanchis-Gual, N.; Sandberg, V.; Sanders, J. R.; Santiago, K. A.; Sarin, N.; Sassolas, B.; Saulson, P. R.; Sauter, O.; Savage, R. L.; Schale, P.; Scheel, M.; Scheuer, J.; Schmidt, P.; Schnabel, R.; Schofield, R. M. S.; Schönbeck, A.; Schreiber, E.; Schulte, B. W.; Schutz, B. F.; Schwalbe, S. G.; Scott, J.; Scott, S. M.; Seidel, E.; Sellers, D.; Sengupta, A. S.; Sennett, N.; Sentenac, D.; Sequino, V.; Sergeev, A.; Setyawati, Y.; Shaddock, D. A.; Shaffer, T.; Shahriar, M. S.; Shaner, M. B.; Shao, L.; Sharma, P.; Shawhan, P.; Shen, H.; Shink, R.; Shoemaker, D. H.; Shoemaker, D. M.; ShyamSundar, S.; Siellez, K.; Sieniawska, M.; Sigg, D.; Silva, A. D.; Singer, L. P.; Singh, N.; Singhal, A.; Sintes, A. M.; Sitmukhambetov, S.; Skliris, V.; Slagmolen, B. J. J.; Slaven-Blair, T. J.; Smith, J. R.; Smith, R. J. E.; Somala, S.; Son, E. J.; Sorazu, B.; Sorrentino, F.; Souradeep, T.; Sowell, E.; Spencer, A. P.; Srivastava, A. K.; Srivastava, V.; Staats, K.; Stachie, C.; Standke, M.; Steer, D. A.; Steinke, M.; Steinlechner, J.; Steinlechner, S.; Steinmeyer, D.; Stevenson, S. P.; Stocks, D.; Stone, R.; Stops, D. J.; Strain, K. A.; Stratta, G.; Strigin, S. E.; Strunk, A.; Sturani, R.; Stuver, A. L.; Sudhir, V.; Summerscales, T. Z.; Sun, L.; Sunil, S.; Suresh, J.; Sutton, P. J.; Swinkels, B. L.; Szczepańczyk, M. J.; Tacca, M.; Tait, S. C.; Talbot, C.; Talukder, D.; Tanner, D. B.; Tápai, M.; Taracchini, A.; Tasson, J. D.; Taylor, R.; Thies, F.; Thomas, M.; Thomas, P.; Thondapu, S. R.; Thorne, K. A.; Thrane, E.; Tiwari, Shubhanshu; Tiwari, Srishti; Tiwari, V.; Toland, K.; Tonelli, M.; Tornasi, Z.; Torres-Forné, A.; Torrie, C. I.; Töyrä, D.; Travasso, F.; Traylor, G.; Tringali, M. C.; Trovato, A.; Trozzo, L.; Trudeau, R.; Tsang, K. W.; Tse, M.; Tso, R.; Tsukada, L.; Tsuna, D.; Tuyenbayev, D.; Ueno, K.; Ugolini, D.; Unnikrishnan, C. S.; Urban, A. L.; Usman, S. A.; Vahlbruch, H.; Vajente, G.; Valdes, G.; van Bakel, N.; van Beuzekom, M.; Brand, J. F. J. van den; Broeck, C. Van Den; Vander-Hyde, D. C.; van Heijningen, J. V.; van der Schaaf, L.; van Veggel, A. A.; Vardaro, M.; Varma, V.; Vass, S.; Vasúth, M.; Vecchio, A.; Vedovato, G.; Veitch, J.; Veitch, P. J.; Venkateswara, K.; Venugopalan, G.; Verkindt, D.; Vetrano, F.; Viceré, A.; Viets, A. D.; Vine, D. J.; Vinet, J. -Y.; Vitale, S.; Vo, T.; Vocca, H.; Vorvick, C.; Vyatchanin, S. P.; Wade, A. R.; Wade, L. E.; Wade, M.; Walet, R.; Walker, M.; Wallace, L.; Walsh, S.; Wang, G.; Wang, H.; Wang, J. Z.; Wang, W. H.; Wang, Y. F.; Ward, R. L.; Warden, Z. A.; Warner, J.; Was, M.; Watchi, J.; Weaver, B.; Wei, L. -W.; Weinert, M.; Weinstein, A. J.; Weiss, R.; Wellmann, F.; Wen, L.; Wessel, E. K.; Weßels, P.; Westhouse, J. W.; Wette, K.; Whelan, J. T.; Whiting, B. F.; Whittle, C.; Wilken, D. M.; Williams, D.; Williamson, A. R.; Willis, J. L.; Willke, B.; Wimmer, M. H.; Winkler, W.; Wipf, C. C.; Wittel, H.; Woan, G.; Woehler, J.; Wofford, J. K.; Worden, J.; Wright, J. L.; Wu, D. S.; Wysocki, D. M.; Xiao, L.; Yamamoto, H.; Yancey, C. C.; Yang, L.; Yap, M. J.; Yazback, M.; Yeeles, D. W.; Yu, Hang; Yu, Haocun; Yuen, S. H. R.; Yvert, M.; Zadrożny, A. K.; Zanolin, M.; Zelenova, T.; Zendri, J. -P.; Zevin, M.; Zhang, J.; Zhang, L.; Zhang, T.; Zhao, C.; Zhou, M.; Zhou, Z.; Zhu, X. J.; Zucker, M. E.; Zweizig, J.; Arzoumanian, Z.; Bogdanov, S.; Cognard, I.; Corongiu, A.; Enoto, T.; Freire, P.; Gendreau, K. C.; Guillemot, L.; Harding, A. K.; Jankowski, F.; Keith, M. J.; Kerr, M.; Lyne, A.; Palfreyman, J.; Possenti, A.; Ridolfi, A.; Stappers, B.; Theureau, G.; Weltervrede, P.
Comments: 39 pages, 9 figures, 2 tables. Accepted for publication in The Astrophysical Journal. Results from the tables are contained in the ancillary file
Submitted: 2019-02-22, last modified: 2019-06-27
We present a search for gravitational waves from 221 pulsars with rotation frequencies $\gtrsim 10$ Hz. We use advanced LIGO data from its first and second observing runs spanning 2015-2017, which provides the highest-sensitivity gravitational-wave data so far obtained. In this search we target emission from both the $l = m = 2$ mass quadrupole mode, with a frequency at twice that of the pulsar's rotation, and from the $l = 2$, $m = 1$ mode, with a frequency at the pulsar rotation frequency. The search finds no evidence for gravitational-wave emission from any pulsar at either frequency. For the $l = m = 2$ mode search, we provide updated upper limits on the gravitational-wave amplitude, mass quadrupole moment, and fiducial ellipticity for 167 pulsars, and the first such limits for a further 55. For 20 young pulsars these results give limits that are below those inferred from the pulsars' spin-down. For the Crab and Vela pulsars our results constrain gravitational-wave emission to account for less than 0.017% and 0.18% of the spin-down luminosity, respectively. For the recycled millisecond pulsar J0711-6830 our limits are only a factor of 1.3 above the spin-down limit, assuming the canonical value of $10^{38}$ kg m$^2$ for the star's moment of inertia, and imply a gravitational-wave-derived upper limit on the star's ellipticity of $1.2\!\times\!10^{-8}$. We also place new limits on the emission amplitude at the rotation frequency of the pulsars.
[36]  oai:arXiv.org:1906.02806  [pdf] - 1925083
A joint NICER and XMM-Newton view of the "Magnificent" thermally emitting X-ray Isolated Neutron Star RX J1605.3+3249
Comments: Accepted for publication in ApJ
Submitted: 2019-06-06
Thermally emitting X-ray isolated neutron stars represent excellent targets for testing cooling surface emission and atmosphere models, which are used to infer physical parameters of the neutron star. Among the seven known members of this class, RX J1605.3+3249 is the only one that still lacks confirmation of its spin period. Here we analyze NICER and XMM-Newton observations of RX J1605.3+3249, in order to address its timing and spectral behavior. Contrary to a previous tentative detection, but in agreement with the recent work by Pires et al. (2019), we find no significant pulsation with pulsed fraction higher than 1.3% (3{\sigma}) for periods above 150 ms. We also find a limit of 2.6% for periods above 2 ms, despite searches in different energy bands. The X-ray spectrum can be fit by either a double-blackbody model or by a single-temperature magnetized atmosphere model, both modified by a Gaussian absorption line at ~0.44 keV. The origin of the absorption feature as a proton cyclotron line or as an atomic transition in the neutron star atmosphere is discussed. The predictions of the best-fit X-ray models extended to IR, optical and UV bands are compared with archival data. Our results are interpreted in the framework of a fallback disk scenario.
[37]  oai:arXiv.org:1906.00974  [pdf] - 1901987
NICER Discovers Spectral Lines During Photospheric Radius Expansion Bursts from 4U 1820-30: Evidence for Burst-driven Winds
Comments: 24 pages, 9 figures, accepted for publication in The Astrophysical Journal Letters
Submitted: 2019-06-03
We report the discovery with the Neutron Star Interior Composition Explorer (NICER) of narrow emission and absorption lines during photospheric radius expansion (PRE) X-ray bursts from the ultracompact binary 4U 1820-30. NICER observed the source in 2017 August accumulating about 60 ks of exposure. Five thermonuclear X-ray bursts were detected of which four showed clear signs of PRE. We extracted spectra during the PRE phases and fit each to a model that includes a comptonized component to describe the accretion-driven emission, and a black body for the burst thermal radiation. The temperature and spherical emitting radius of the fitted black body are used to assess the strength of PRE in each burst. The two strongest PRE bursts (burst pair 1) had black body temperatures of approximately 0.6 keV and emitting radii of 100 km (at a distance of 8.4 kpc). The other two bursts (burst pair 2) had higher temperatures (~0.67 keV) and smaller radii (75 km). All of the PRE bursts show evidence of narrow line emission near 1 keV. By co-adding the PRE phase spectra of burst pairs 1 and, separately, 2 we find, in both co-added spectra, significant, narrow, spectral features near 1.0 (emission), 1.7 and 3.0 keV (both in absorption). Remarkably, all the fitted line centroids in the co-added spectrum of burst pair 1 appear systematically blue-shifted by a factor of $1.046 \pm 0.006$ compared to the centroids of pair 2, strongly indicative of a gravitational red-shift, a wind-induced blue-shift, or more likely some combination of both effects. The observed shifts are consistent with this scenario in that the stronger PRE bursts in pair 1 reach larger photospheric radii, and thus have weaker gravitational red-shifts, and they generate faster outflows, yielding higher blue-shifts. We discuss possible elemental identifications for the observed features in the context of recent burst-driven wind models.
[38]  oai:arXiv.org:1811.09218  [pdf] - 1920810
Anti-glitches in the Ultraluminous Accreting Pulsar NGC 300 ULX-1 Observed with NICER
Comments: ApJ, accepted
Submitted: 2018-11-22, last modified: 2019-05-28
We present evidence for three spin-down glitches (or `anti-glitches') in the ultraluminous accreting X-ray pulsar NGC 300 ULX-1 in timing observations made with the Neutron Star Interior Composition Explorer (NICER). Our timing analysis reveals three sudden spin-down events of magnitudes $\Delta\nu = -23, -30,$ and $-43 \,\mu$Hz (fractional amplitudes $\Delta\nu/\nu = -4.4, -5.5,$ and $-7.7 \times 10^{-4}$). We determined fully phase-coherent timing solutions through the first two glitches, giving us high confidence in their detection, while the third candidate glitch is somewhat less secure. These are larger in magnitude (and opposite in sign) than any known radio pulsar glitch. This may be caused by the prolonged rapid spin-up of the pulsar causing a sudden transfer of angular momentum between the superfluid and non-superfluid components of the star. We find no evidence for profile or spectral changes at the epochs of the glitches, supporting the conclusion that these are due to the same process as in normal pulsar glitches, but in reverse.
[39]  oai:arXiv.org:1905.12060  [pdf] - 1899410
Discovery of Soft X-ray Pulsations from PSR J1231-1411 using NICER
Comments: ApJ Letters, in press
Submitted: 2019-05-28
We report the discovery of soft X-ray pulsations from the nearby millisecond pulsar PSR J1231$-$1411 using NICER. The pulsed emission is characterized by a broad and asymmetric main pulse and a much fainter secondary interpulse, with a total pulsed count rate of 0.055 c s$^{-1}$ in the 0.35-1.5 keV band. We analyzed Fermi LAT data to update the pulse timing model covering 10 years of data and used that model to coherently combine NICER data over a year of observations. Spectral modeling suggests that the flux is dominated by thermal emission from a hot spot (or spots) on the neutron star surface. The phase relationship between the X-ray pulse and the radio and $\gamma$ rays provides insight into the geometry of the system.
[40]  oai:arXiv.org:1812.11585  [pdf] - 1969541
The NANOGrav 11-Year Data Set: Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries
Comments: 10 pages, 11 figures. Accepted by Astrophysical Journal. Please send any comments/questions to S. J. Vigeland (vigeland@uwm.edu)
Submitted: 2018-12-30, last modified: 2019-05-21
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.
[41]  oai:arXiv.org:1905.04440  [pdf] - 1894361
NICER Observations of the 2018 Outburst of XTE J1810-197
Comments: Accepted for publication in the Astrophysical Journal Letters
Submitted: 2019-05-10
We present the earliest available soft X-ray observations of XTE J1810-187, the prototypical transient magnetar, obtained 75--84 days after its 2018 outburst with the Neutron Star Interior Composition Explorer (NICER). Using a series of observations covering eight days we find that its decreasing X-ray flux is well-described by either a blackbody plus power-law or a two-blackbody spectral model. The 2-10 keV flux of the source varied from (1.206+/-0.007)x10^{-10} to (1.125+/-0.004)x10^{-10} erg s^{-1} cm^{-2}, a decrease of about 7% within our observations and 44% from that measured 7-14 days after the outburst with NuSTAR. We confirm that the pulsed fraction and spin pulse phase of the neutron star are energy dependent up to at least 8 keV. Phase resolved spectroscopy of the pulsar suggests magnetospheric variations relative to the line of sight.
[42]  oai:arXiv.org:1904.11311  [pdf] - 1894334
On the curious pulsation properties of the accreting millisecond pulsar IGR J17379-3747
Comments: 16 pages, 8 figures, 3 tables. Accepted for publication in ApJ
Submitted: 2019-04-25
We report on the Neutron Star Interior Composition Explorer (NICER) monitoring campaign of the 468 Hz accreting millisecond X-ray pulsar IGR J17379-3747. From a detailed spectral and timing analysis of the coherent pulsations we find that they show a strong energy dependence, with soft thermal emission lagging about 640 microseconds behind the hard, Comptonized emission. Additionally, we observe uncommonly large pulse fractions, with measured amplitudes in excess of 20% sinusoidal fractional amplitude across the NICER passband and fluctuations of up to ~70%. Based on a phase-resolved spectral analysis, we suggest that these extreme properties might be explained if the source has an unusually favorable viewing geometry with a large magnetic misalignment angle. Due to these large pulse fractions, we were able to detect pulsations down to quiescent luminosities (~5 x 10^33 erg s^-1). We discuss these low-luminosity pulsations in the context of transitional millisecond pulsars.
[43]  oai:arXiv.org:1903.04648  [pdf] - 1847053
Determining the Equation of State of Cold, Dense Matter with X-ray Observations of Neutron Stars
Comments: Submitted in response to the Astronomy and Astrophysics Decadal Survey (Astro2020) call for Science White Papers; 8 pages, 2 figures
Submitted: 2019-03-11
The unknown state of matter at ultra-high density, large proton/neutron number asymmetry, and low temperature is a major long-standing problem in modern physics. Neutron stars provide the only known setting in the Universe where matter in this regime can stably exist. Valuable information about the interior structure of neutron stars can be extracted via sensitive observations of their exteriors. There are several complementary techniques that require different combinations of high time resolution, superb spectral resolution, and high spatial resolution. In the upcoming decade and beyond, measurements of the masses and radii of an ensemble of neutron stars using these techniques, based on data from multiple proposed next-generation X-ray telescopes, can produce definitive empirical constraints on the allowed dense matter equation of state.
[44]  oai:arXiv.org:1903.03035  [pdf] - 1845076
STROBE-X: X-ray Timing and Spectroscopy on Dynamical Timescales from Microseconds to Years
Ray, Paul S.; Arzoumanian, Zaven; Ballantyne, David; Bozzo, Enrico; Brandt, Soren; Brenneman, Laura; Chakrabarty, Deepto; Christophersen, Marc; DeRosa, Alessandra; Feroci, Marco; Gendreau, Keith; Goldstein, Adam; Hartmann, Dieter; Hernanz, Margarita; Jenke, Peter; Kara, Erin; Maccarone, Tom; McDonald, Michael; Nowak, Michael; Phlips, Bernard; Remillard, Ron; Stevens, Abigail; Tomsick, John; Watts, Anna; Wilson-Hodge, Colleen; Wood, Kent; Zane, Silvia; Ajello, Marco; Alston, Will; Altamirano, Diego; Antoniou, Vallia; Arur, Kavitha; Ashton, Dominic; Auchettl, Katie; Ayres, Tom; Bachetti, Matteo; Balokovic, Mislav; Baring, Matthew; Baykal, Altan; Begelman, Mitch; Bhat, Narayana; Bogdanov, Slavko; Briggs, Michael; Bulbul, Esra; Bult, Petrus; Burns, Eric; Cackett, Ed; Campana, Riccardo; Caspi, Amir; Cavecchi, Yuri; Chenevez, Jerome; Cherry, Mike; Corbet, Robin; Corcoran, Michael; Corsi, Alessandra; Degenaar, Nathalie; Drake, Jeremy; Eikenberry, Steve; Enoto, Teruaki; Fragile, Chris; Fuerst, Felix; Gandhi, Poshak; Garcia, Javier; Goldstein, Adam; Gonzalez, Anthony; Grefenstette, Brian; Grinberg, Victoria; Grossan, Bruce; Guillot, Sebastien; Guver, Tolga; Haggard, Daryl; Heinke, Craig; Heinz, Sebastian; Hemphill, Paul; Homan, Jeroen; Hui, Michelle; Huppenkothen, Daniela; Ingram, Adam; Irwin, Jimmy; Jaisawal, Gaurava; Jaodand, Amruta; Kalemci, Emrah; Kaplan, David; Keek, Laurens; Kennea, Jamie; Kerr, Matthew; van der Klis, Michiel; Kocevski, Daniel; Koss, Mike; Kowalski, Adam; Lai, Dong; Lamb, Fred; Laycock, Silas; Lazio, Joseph; Lazzati, Davide; Longcope, Dana; Loewenstein, Michael; Maitra, Dipankair; Majid, Walid; Maksym, W. Peter; Malacaria, Christian; Margutti, Raffaella; Martindale, Adrian; McHardy, Ian; Meyer, Manuel; Middleton, Matt; Miller, Jon; Miller, Cole; Motta, Sara; Neilsen, Joey; Nelson, Tommy; Noble, Scott; O'Brien, Paul; Osborne, Julian; Osten, Rachel; Ozel, Feryal; Palliyaguru, Nipuni; Pasham, Dheeraj; Patruno, Alessandro; Pelassa, Vero; Petropoulou, Maria; Pilia, Maura; Pohl, Martin; Pooley, David; Prescod-Weinstein, Chanda; Psaltis, Dimitrios; Raaijmakers, Geert; Reynolds, Chris; Riley, Thomas E.; Salvesen, Greg; Santangelo, Andrea; Scaringi, Simone; Schanne, Stephane; Schnittman, Jeremy; Smith, David; Smith, Krista Lynne; Snios, Bradford; Steiner, Andrew; Steiner, Jack; Stella, Luigi; Strohmayer, Tod; Sun, Ming; Tauris, Thomas; Taylor, Corbin; Tohuvavohu, Aaron; Vacchi, Andrea; Vasilopoulos, Georgios; Veledina, Alexandra; Walsh, Jonelle; Weinberg, Nevin; Wilkins, Dan; Willingale, Richard; Wilms, Joern; Winter, Lisa; Wolff, Michael; Zand, Jean in 't; Zezas, Andreas; Zhang, Bing; Zoghbi, Abdu
Comments: 50 pages, Probe class mission concept study report submitted to NASA for Astro2020 Decadal Survey
Submitted: 2019-03-07, last modified: 2019-03-08
We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.
[45]  oai:arXiv.org:1902.07130  [pdf] - 1863991
High-Precision X-ray Timing of Three Millisecond Pulsars with NICER: Stability Estimates and Comparison with Radio
Comments: 30 pages, 15 figures; accepted to ApJ
Submitted: 2019-02-19
The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years.
[46]  oai:arXiv.org:1902.00144  [pdf] - 1894219
Neutron Star Interior Composition Explorer X-ray Timing of the Radio and Gamma-ray Quiet Pulsars PSR J1412+7922 AND PSR J1849-0001
Comments: 12 pages, 4 figures; submitted to the Astrophysical Journal
Submitted: 2019-01-31, last modified: 2019-02-04
We present new timing and spectral analyses of PSR J1412+7922 (Calvera) and PSR J1849-0001, which are only seen as pulsars in X-rays, based on observations conducted with the Neutron Star Interior Composition Explorer (NICER). We obtain updated and substantially improved pulse ephemerides compared to previous X-ray studies, as well as spectra that can be well-fit by simple blackbodies and/or a power law. Our refined timing measurements enable deeper searches for pulsations at other wavelengths and sensitive targeted searches by LIGO/Virgo for continuous gravitational waves from these neutron stars. Using the sensitivity of LIGO's first observing run, we estimate constraints that a gravitational wave search of these pulsars would obtain on the size of their mass deformation and r-mode fluid oscillation.
[47]  oai:arXiv.org:1901.03877  [pdf] - 1818876
The corona contracts in a black-hole transient
Comments: Published in Nature January 9, 2019
Submitted: 2019-01-12
The geometry of the accretion flow around stellar-mass black holes can change on timescales of days to months. When a black hole emerges from quiescence (that is, it "turns on" after accreting material from its companion) it has a very hard (high-energy) X-ray spectrum produced by a hot corona positioned above its accretion disk, and then transitions to a soft (lower-energy) spectrum dominated by emission from the geometrically thin accretion disk, which extends to the innermost stable circular orbit. Much debate persists over how this transition occurs and whether it is driven largely by a reduction in the truncation radius of the disk or by a reduction in the spatial extent of the corona. Observations of X-ray reverberation lags in supermassive black-hole systems suggest that the corona is compact and that the disk extends nearly to the central black hole. Observations of stellar-mass black holes, however, reveal equivalent (mass-scaled) reverberation lags that are much larger, leading to the suggestion that the accretion disk in the hard X-ray state of stellar-mass black holes is truncated at a few hundreds of gravitational radii from the black hole. Here we report X-ray observations of the black-hole transient MAXI J1820+070. We find that the reverberation time lags between the continuum-emitting corona and the irradiated accretion disk are 6 to 20 times shorter than previously seen. The timescale of the reverberation lags shortens by an order of magnitude over a period of weeks, whereas the shape of the broadened iron K emission line remains remarkably constant. This suggests a reduction in the spatial extent of the corona, rather than a change in the inner edge of the accretion disk.
[48]  oai:arXiv.org:1901.03731  [pdf] - 1814298
A new class of flares from accreting supermassive black holes
Comments: Published in Nature Astronomy
Submitted: 2019-01-11
Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed, including tidal disruption events, 'changing look' AGNs and other extreme variability objects. The physics behind the 're-ignition', enhancement and 'shut-down' of accretion onto SMBHs is not entirely understood. Here we present a rapid increase in ultraviolet-optical emission in the centre of a nearby galaxy, marking the onset of sudden increased accretion onto a SMBH. The optical spectrum of this flare, dubbed AT 2017bgt, exhibits a mix of emission features. Some are typical of luminous, unobscured AGNs, but others are likely driven by Bowen fluorescence - robustly linked here with high-velocity gas in the vicinity of the accreting SMBH. The spectral features and increased ultraviolet flux show little evolution over a period of at least 14 months. This disfavours the tidal disruption of a star as their origin, and instead suggests a longer-term event of intensified accretion. Together with two other recently reported events with similar properties, we define a new class of SMBH-related flares. This has important implications for the classification of different types of enhanced accretion onto SMBHs.
[49]  oai:arXiv.org:1809.05064  [pdf] - 1814977
PSR J2234+0611: A new laboratory for stellar evolution
Comments: 16 pages, 5 figures, accepted by ApJ
Submitted: 2018-09-13, last modified: 2018-11-26
We report timing results for PSR J2234+0611, a 3.6-ms pulsar in a 32-day, eccentric (e = 0.13) orbit with a helium white dwarf companion discovered as part of the Arecibo Observatory 327 MHz drift scan survey. The precise timing and the eccentric nature of the orbit allow precise measurements of an unusual number of parameters: a) a precise proper motion of 27.10(3) mas/yr and a parallax of 1.05(4) mas resulting in a pulsar distance of 0.95(4) kpc; this allows a precise estimate of the transverse velocity, 123(5) km/s. Together with previously published spectroscopic measurements of the systemic radial velocity, this allows a full 3-D determination of the system's velocity; b) precise measurements of the rate of advance of periastron, which after subtraction of the contribution of the proper motion yields a total system mass of $1.6518^{+0.0033}_{-0.0035}$ solar masses; c) a Shapiro delay measurement, h_3 = $82 \pm 14$ ns despite the orbital inclination not being near 90 deg; combined with the measurement of the total mass, this yields a pulsar mass of $1.353^{+0.014}_{-0.017}$ solar masses and a companion mass of $0.298^{+0.015}_{-0.012}$ solar masses; d) we measure precisely the secular variation of the projected semi-major axis and detect significant annual orbital parallax; together these allow a determination of the full 3-D orbital geometry, including an unambiguous orbital inclination (i = $138.7^{+2.5}_{-2.2}$ deg) and a position angle for the line of nodes (Omega = $44^{+5}_{-4}$ deg). We discuss the component masses to investigate hypotheses previously advanced to explain the origin of eccentric MSPs. The unprecedented determination of the full 3-D position, motion and orbital orientation of the system, plus the precisely measured pulsar and WD mass and the latter's optical detection make this system an unique test of our understanding of white dwarfs and their atmospheres.
[50]  oai:arXiv.org:1810.08269  [pdf] - 1791213
The NANOGrav 11-year Data Set: Pulse Profile Variability
Comments: 38 pages, 25 figures
Submitted: 2018-10-18, last modified: 2018-11-14
Access to 50 years of data has led to the discovery of pulsar emission and rotation variability on timescales of months and years. Most of this long-term variability has been seen in long-period pulsars, with relatively little focus on recycled millisecond pulsars. We have analyzed a 38-pulsar sub-set of the 45 millisecond pulsars in the NANOGrav 11-year data set, in order to review their pulse profile stability. The most variability, on any timescale, is seen in PSRs J1713+0747, B1937+21 and J2145-0750. The strongest evidence for long-timescale pulse profile changes is seen in PSRs B1937+21 and J1643-1224. We have focused our analyses on these four pulsars in an attempt to elucidate the causes of their profile variability. Effects of scintillation seem to be responsible for the profile modifications of PSR J2145-0750. We see evidence that imperfect polarization calibration contributes to the profile variability of PSRs J1713+0747 and B1937+21, along with radio frequency interference around 2 GHz, but find that propagation effects also have an influence. The changes seen in PSR J1643-1224 have been reported previously, yet elude explanation beyond their astrophysical nature. Regardless of cause, unmodeled pulse profile changes are detrimental to the accuracy of pulsar timing and must be incorporated into the timing models where possible.
[51]  oai:arXiv.org:1809.10744  [pdf] - 1767703
Studying the solar system with the International Pulsar Timing Array
Comments: Accepted for publication by the Monthly Notices of the Royal Astronomical Society
Submitted: 2018-09-27
Pulsar-timing analyses are sensitive to errors in the solar-system ephemerides (SSEs) that timing models utilise to estimate the location of the solar-system barycentre, the quasi-inertial reference frame to which all recorded pulse times-of-arrival are referred. Any error in the SSE will affect all pulsars, therefore pulsar timing arrays (PTAs) are a suitable tool to search for such errors and impose independent constraints on relevant physical parameters. We employ the first data release of the International Pulsar Timing Array to constrain the masses of the planet-moons systems and to search for possible unmodelled objects (UMOs) in the solar system. We employ ten SSEs from two independent research groups, derive and compare mass constraints of planetary systems, and derive the first PTA mass constraints on asteroid-belt objects. Constraints on planetary-system masses have been improved by factors of up to 20 from the previous relevant study using the same assumptions, with the mass of the Jovian system measured at 9.5479189(3)$\times10^{-4}$ $M_{\odot}$. The mass of the dwarf planet Ceres is measured at 4.7(4)$\times10^{-10}$ $M_{\odot}$. We also present the first sensitivity curves using real data that place generic limits on the masses of UMOs, which can also be used as upper limits on the mass of putative exotic objects. For example, upper limits on dark-matter clumps are comparable to published limits using independent methods. While the constraints on planetary masses derived with all employed SSEs are consistent, we note and discuss differences in the associated timing residuals and UMO sensitivity curves.
[52]  oai:arXiv.org:1809.07556  [pdf] - 1763115
A NICER Discovery of a Low-Frequency Quasi-Periodic Oscillation in the Soft-Intermediate State of MAXI J1535-571
Comments: Accepted for publication in ApJ Letters
Submitted: 2018-09-20, last modified: 2018-09-24
We present the discovery of a low-frequency $\approx 5.7$ Hz quasi-periodic oscillation (QPO) feature in observations of the black hole X-ray binary MAXI J1535-571 in its soft-intermediate state, obtained in September-October 2017 by the Neutron Star Interior Composition Explorer (NICER). The feature is relatively broad (compared to other low-frequency QPOs; quality factor $Q\approx 2$) and weak (1.9% rms in 3-10 keV), and is accompanied by a weak harmonic and low-amplitude broadband noise. These characteristics identify it as a weak Type A/B QPO, similar to ones previously identified in the soft-intermediate state of the transient black hole X-ray binary XTE J1550-564. The lag-energy spectrum of the QPO shows increasing soft lags towards lower energies, approaching 50 ms at 1 keV (with respect to a 3-10 keV continuum). This large phase shift has similar amplitude but opposite sign to that seen in Rossi X-ray Timing Explorer data for a Type B QPO from the transient black hole X-ray binary GX 339-4. Previous phase-resolved spectroscopy analysis of the Type B QPO in GX 339-4 pointed towards a precessing jet-like corona illuminating the accretion disk as the origin of the QPO signal. We suggest that this QPO in MAXI J1535-571 may have the same origin, with the different lag sign depending on the scale height of the emitting region and the observer inclination angle.
[53]  oai:arXiv.org:1809.03058  [pdf] - 1842358
The NANOGrav 12.5-Year Data Set: The Frequency Dependence of Pulse Jitter in Precision Millisecond Pulsars
Comments: 18 pages, 10 figures, submitted to ApJ
Submitted: 2018-09-09
Low-frequency gravitational-wave experiments require the highest timing precision from an array of the most stable millisecond pulsars. Several known sources of noise on short timescales in single radio-pulsar observations are well described by a simple model of three components: template-fitting from a finite signal-to-noise ratio, pulse phase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellar scattering variations. Currently template-fitting errors dominate, but as radio telescopes push towards higher signal-to-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics of jitter becomes crucial for properly characterizing arrival-time uncertainties. We characterize the radio-frequency dependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for its frequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations of rms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter is typically ~1% of pulse phase. The average frequency dependence for all pulsars using a power-law model has index -0.42. We investigate the jitter variations for the interpulse of PSR~B1937+21 and find no significant deviations from the main pulse rms jitter. We also test the time-variation of jitter in two pulsars and find that systematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly model jitter as a significant component of the noise within the detector.
[54]  oai:arXiv.org:1808.10195  [pdf] - 1759371
NuSTAR and NICER reveal IGR J17591-2342 as a new accreting millisecond X-ray pulsar
Comments: 6 pages, 5 figures, 2 tables. Submitted to A&A
Submitted: 2018-08-30
We report on the discovery by the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Neutron Star Interior Composition Explorer (NICER) of the accreting millisecond X-ray pulsar IGR J17591-2342, detecting coherent X-ray pulsations around 527.4 Hz (1.9 ms) with a clear Doppler modulation. This implies an orbital period of ~8.8 hours and a projected semi-major axis of ~1.23 lt-s. From the binary mass function, we estimate a minimum companion mass of 0.42 solar masses, obtained assuming a neutron star mass of 1.4 solar masses and an inclination angle lower than 60 degrees, as suggested by the absence of eclipses or dips in the light-curve of the source. The broad-band energy spectrum is dominated by Comptonisation of soft thermal seed photons with a temperature of ~0.7 keV by electrons heated to 21 keV. We also detect black-body-like thermal direct emission compatible with an emission region of a few kilometers and temperature compatible with the seed source of Comptonisation. A weak Gaussian line centered on the iron K-alpha; complex can be interpreted as a signature of disc reflection. A similar spectrum characterises the NICER spectra, measured during the outburst fading.
[55]  oai:arXiv.org:1808.07078  [pdf] - 1838174
The NANOGrav 11-year Data Set: Solar Wind Sounding Through Pulsar Timing
Comments: 15 pages, submitted to ApJ
Submitted: 2018-08-21
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens of millisecond pulsars for over a decade. We have accrued a large collection of dispersion measure (DM) measurements sensitive to the total electron content between Earth and the pulsars at each observation. All lines of sight cross through the solar wind which produces correlated DM fluctuations in all pulsars. We develop and apply techniques for extracting the imprint of the solar wind from the full collection of DM measurements in the recently released NANOGrav 11-yr data set. We filter out long time scale DM fluctuations attributable to structure in the interstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate the correlated signature of the wind from signals unique to individual lines of sight. When treating the solar wind as spherically symmetric and constant in time, we find the electron number density at 1 A.U. to be $7.9\pm0.2$ cm$^{-3}$. Our data shows little evidence of long-term variation in the density of the wind. We argue that our techniques paired with a high cadence, low radio frequency observing campaign of near-ecliptic pulsars would be capable of mapping out large-scale latitudinal structure in the wind.
[56]  oai:arXiv.org:1808.06442  [pdf] - 1735589
NICER Detection of Strong Photospheric Expansion during a Thermonuclear X-Ray Burst from 4U 1820-30
Comments: 6 pages, 4 figures, 1 table, published by ApJ Letters
Submitted: 2018-08-20
The Neutron Star Interior Composition Explorer (NICER) on the International Space Station (ISS) observed strong photospheric expansion of the neutron star in 4U 1820-30 during a Type I X-ray burst. A thermonuclear helium flash in the star's envelope powered a burst that reached the Eddington limit. Radiation pressure pushed the photosphere out to ~200 km, while the blackbody temperature dropped to 0.45 keV. Previous observations of similar bursts were performed with instruments that are sensitive only above 3 keV, and the burst signal was weak at low temperatures. NICER's 0.2-12 keV passband enables the first complete detailed observation of strong expansion bursts. The strong expansion lasted only 0.6 s, and was followed by moderate expansion with a 20 km apparent radius, before the photosphere finally settled back down at 3 s after the burst onset. In addition to thermal emission from the neutron star, the NICER spectra reveal a second component that is well fit by optically thick Comptonization. During the strong expansion, this component is six times brighter than prior to the burst, and it accounts for 71% of the flux. In the moderate expansion phase, the Comptonization flux drops, while the thermal component brightens, and the total flux remains constant at the Eddington limit. We speculate that the thermal emission is reprocessed in the accretion environment to form the Comptonization component, and that changes in the covering fraction of the star explain the evolution of the relative contributions to the total flux.
[57]  oai:arXiv.org:1808.06441  [pdf] - 1735588
NICER Observes the Effects of an X-Ray Burst on the Accretion Environment in Aql X-1
Comments: 6 pages, 5 figures, published by ApJ Letters
Submitted: 2018-08-20
Accretion disks around neutron stars regularly undergo sudden strong irradiation by Type I X-ray bursts powered by unstable thermonuclear burning on the stellar surface. We investigate the impact on the disk during one of the first X-ray burst observations with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. The burst is seen from Aql X-1 during the hard spectral state. In addition to thermal emission from the neutron star, the burst spectrum exhibits an excess of soft X-ray photons below 1 keV, where NICER's sensitivity peaks. We interpret the excess as a combination of reprocessing by the strongly photoionized disk and enhancement of the pre-burst persistent flux, possibly due to Poynting Robertson drag or coronal reprocessing. This is the first such detection for a short sub-Eddington burst. As these bursts are observed frequently, NICER will be able to study how X-ray bursts affect the disk and corona for a range of accreting neutron star systems and disk states.
[58]  oai:arXiv.org:1808.04392  [pdf] - 1732879
NICER Discovers the Ultracompact Orbit of the Accreting Millisecond Pulsar IGR J17062-6143
Comments: 19 pages, 7 figures, published in the Astrophysical Journal Letters
Submitted: 2018-08-13
We present results of recent Neutron Star Interior Composition Explorer observations of the accreting millisecond X-ray pulsar IGR J17062-6143 that show that it resides in a circular, ultracompact binary with a 38 minute orbital period. NICER observed the source for approximately 26 ksec over a 5.3 day span in 2017 August, and again for 14 and 11 ksec in 2017 October and November, respectively. A power spectral analysis of the August exposure confirms the previous detection of pulsations at 163.656 Hz in Rossi X-ray Timing Explorer data, and reveals phase modulation due to orbital motion of the neutron star. A coherent search for the orbital solution using the Z^2 method finds a best-fitting circular orbit with a period of 2278.21 s (37.97 min), a projected semi-major axis of 0.00390 lt-sec, and a barycentric pulsar frequency of 163.6561105 Hz. This is currently the shortest known orbital period for an AMXP. The mass function is 9.12 e-8} solar masses, presently the smallest known for a stellar binary. The minimum donor mass ranges from about 0.005 - 0.007 solar masses, for a neutron star mass from 1.2 - 2 solar masses. Assuming mass transfer is driven by gravitational radiation, we find donor mass and binary inclination bounds of 0.0175 - 0.0155 solar masses and 19 deg < i < 27.5 deg, where the lower and upper bounds correspond to 1.4 and 2 solar mass neutron stars, respectively. Folding the data accounting for the orbital modulation reveals a sinusoidal profile with fractional amplitude 2.04 +- 0.11 % (0.3 - 3.2 keV).
[59]  oai:arXiv.org:1808.04294  [pdf] - 1763058
NICER Discovers mHz Oscillations in the "Clocked" Burster GS 1826-238
Comments: 26 pages, 10 figures, accepted for publication in the Astrophysical Journal
Submitted: 2018-08-13
We report the discovery with the Neutron Star Interior Composition Explorer (NICER) of mHz X-ray brightness oscillations from the "clocked burster." NICER observed the source in the periods 2017 June 20 - 29, July 11 - 13, and September 9 - 15, for a total useful exposure of 34 ks. Two consecutive dwells obtained on 2017 September 9 revealed highly significant oscillations at a frequency of 8 mHz. The fractional, sinusoidal modulation amplitude increases from 0.7 % at 1 keV to approximately 2 % at 6 keV. Similar oscillations were also detected at lower significance in three additional dwells. The oscillation frequency and amplitude are consistent with those of mHz QPOs reported in other accreting neutron star systems. A thermonuclear X-ray burst was also observed on 2017 June 22. The burst properties and X-ray colors are both consistent with the source being in a soft spectral state during these observations, findings that are confirmed by ongoing monitoring with MAXI and SWIFT-BAT. Assuming that the mHz oscillations are associated with black body emission from the neutron star surface, modeling of the phase-resolved spectra shows that the oscillation is consistent with being produced by modulation of the temperature component of this emission. In this interpretation, the black body normalization, proportional to the emitting surface area, is consistent with being constant through the oscillation cycle. We place the observations in the context of the current theory of marginally stable burning and briefly discuss the potential for constraining neutron star properties using mHz oscillations.
[60]  oai:arXiv.org:1805.11906  [pdf] - 1721149
NICER Detects a Soft X-ray Kilohertz Quasi-periodic Oscillation in 4U 0614+09
Comments: proof correction applied. 7 pages, 3 figures, 1 table, accepted for publication in ApJ Letters
Submitted: 2018-05-30, last modified: 2018-07-24
We report on the detection of a kilohertz quasi-periodic oscillation (QPO) with the Neutron Star Interior Composition Explorer (NICER). Analyzing approximately 165 ks of NICER exposure on the X-ray burster 4U 0614+09, we detect multiple instances of a single-peak upper kHz QPO, with centroid frequencies that range from 400 Hz to 750 Hz. We resolve the kHz QPO as a function of energy, and measure, for the first time, the QPO amplitude below 2 keV. We find the fractional amplitude at 1 keV is on the order of 2% rms, and discuss the implications for the QPO emission process in the context of Comptonization models.
[61]  oai:arXiv.org:1807.09100  [pdf] - 1743826
On the 2018 outburst of the accreting millisecond X-ray pulsar Swift J1756.9-2508 as seen with NICER
Comments: 9 pages, 4 figures, 2 tables. Accepted for publication in ApJ
Submitted: 2018-07-24
We report on the coherent timing analysis of the 182 Hz accreting millisecond X-ray pulsar Swift J1756.9$-$2508 during its 2018 outburst as observed with the Neutron Star Interior Composition Explorer (NICER). Combining our NICER observations with Rossi X-ray Timing Explorer observations of the 2007 and 2009 outbursts, we also studied the long-term spin and orbital evolution of this source. We find that the binary system is well described by a constant orbital period model, with an upper limit on the orbital period derivative of $|\dot P_b| < 7.4\times 10^{-13}$ s s$^{-1}$. Additionally, we improve upon the source coordinates through astrometric analysis of the pulse arrival times, finding R.A. = $17^{\rm h}56^{\rm m}57.18^{\rm s}\pm0.08^{\rm s}$ and Decl. = $-25^{\circ}06'27.8''\pm3.5''$, while simultaneously measuring the long-term spin frequency derivative as $\dot\nu = -7.3\times 10^{-16}$ Hz s$^{-1}$. We briefly discuss the implications of these measurements in the context of the wider population of accreting millisecond pulsars.
[62]  oai:arXiv.org:1712.03651  [pdf] - 1714763
A Second Chromatic Timing Event of Interstellar Origin toward PSR J1713+0747
Comments: 9 pages, 3 figures, published in ApJ
Submitted: 2017-12-11, last modified: 2018-07-12
The frequency dependence of radio pulse arrival times provides a probe of structures in the intervening media. Demorest et al. 2013 was the first to show a short-term (~100-200 days) reduction in the electron content along the line of sight to PSR J1713+0747 in data from 2008 (approximately MJD 54750) based on an apparent dip in the dispersion measure of the pulsar. We report on a similar event in 2016 (approximately MJD 57510), with average residual pulse-arrival times of approximately 3.0,-1.3, and -0.7 microseconds at 820, 1400, and 2300 MHz, respectively. Timing analyses indicate possible departures from the standard nu^-2 dispersive-delay dependence. We discuss and rule out a wide variety of potential interpretations. We find the likeliest scenario to be lensing of the radio emission by some structure in the interstellar medium, which causes multiple frequency-dependent pulse arrival-time delays.
[63]  oai:arXiv.org:1807.01179  [pdf] - 1709616
STROBE-X: A probe-class mission for X-ray spectroscopy and timing on timescales from microseconds to years
Comments: 20 pages, 10 figures, to appear in Proc. SPIE Astronomical Telescopes and Instrumentation 2018
Submitted: 2018-07-03
We describe the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept that will provide an unprecedented view of the X-ray sky, performing timing and spectroscopy over both a broad energy band (0.2-30 keV) and a wide range of timescales from microseconds to years. STROBE-X comprises two narrow-field instruments and a wide field monitor. The soft or low-energy band (0.2-12 keV) is covered by an array of lightweight optics (3-m focal length) that concentrate incident photons onto small solid-state detectors with CCD-level (85-175 eV) energy resolution, 100 ns time resolution, and low background rates. This technology has been fully developed for NICER and will be scaled up to take advantage of the longer focal length of STROBE-X. The higher-energy band (2-30 keV) is covered by large-area, collimated silicon drift detectors that were developed for the European LOFT mission concept. Each instrument will provide an order of magnitude improvement in effective area over its predecessor (NICER in the soft band and RXTE in the hard band). Finally, STROBE-X offers a sensitive wide-field monitor (WFM), both to act as a trigger for pointed observations of X-ray transients and also to provide high duty-cycle, high time-resolution, and high spectral-resolution monitoring of the variable X-ray sky. The WFM will boast approximately 20 times the sensitivity of the RXTE All-Sky Monitor, enabling multi-wavelength and multi-messenger investigations with a large instantaneous field of view. This mission concept will be presented to the 2020 Decadal Survey for consideration.
[64]  oai:arXiv.org:1807.00708  [pdf] - 1729100
The NANOGrav 11-Year Data Set: Arecibo Observatory Polarimetry And Pulse Microcomponents
Comments: 41 pages, 20 figures
Submitted: 2018-07-02
We present the polarization pulse profiles for 28 pulsars observed with the Arecibo Observatory by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) timing project at 2.1 GHz, 1.4 GHz, and 430 MHz. These profiles represent some of the most sensitive polarimetric millisecond pulsar profiles to date, revealing the existence of microcomponents (that is, pulse components with peak intensities much lower than the total pulse peak intensity). Although microcomponents have been detected in some pulsars previously, we present microcomponents for PSRs B1937+21, J1713+0747, and J2234+0944 for the first time. These microcomponents can have an impact on pulsar timing, geometry, and flux density determination. We present rotation measures for all 28 pulsars, determined independently at different observation frequencies and epochs, and find the Galactic magnetic fields derived from these rotation measures to be consistent with current models. These polarization profiles were made using measurement equation template matching, which allows us to generate the polarimetric response of the Arecibo Observatory on an epoch-by-epoch basis. We use this method to describe its time variability, and find that the polarimetric responses of the Arecibo Observatory's 1.4 and 2.1 GHz receivers vary significantly with time.
[65]  oai:arXiv.org:1806.10094  [pdf] - 1732756
NICER and Fermi GBM Observations of the First Galactic Ultraluminous X-ray Pulsar Swift J0243.6+6124
Comments: 35 pages, 7 figures (main text), 5 figures (appendix); Accepted for publication in ApJ on June 19, 2018
Submitted: 2018-06-26
Swift J0243.6+6124 is a newly discovered Galactic Be/X-ray binary, revealed in late September 2017 in a giant outburst with a peak luminosity of 2E+39 (d/7 kpc)^2 erg/s (0.1-10 keV), with no formerly reported activity. At this luminosity, Swift J0243.6+6124 is the first known galactic ultraluminous X-ray pulsar. We describe Neutron star Interior Composition Explorer (NICER)} and Fermi Gamma-ray Burst Monitor (GBM) timing and spectral analyses for this source. A new orbital ephemeris is obtained for the binary system using spin-frequencies measured with GBM and 15-50 keV fluxes measured with the Neil Gehrels Swift Observatory Burst Alert Telescope to model the system's intrinsic spin-up. Power spectra measured with NICER show considerable evolution with luminosity, including a quasi-periodic oscillation (QPO) near 50 mHz that is omnipresent at low luminosity and has an evolving central frequency. Pulse profiles measured over the combined 0.2-100 keV range show complex evolution that is both luminosity and energy dependent. Near the critical luminosity of L~1E+38 erg/s, the pulse profiles transition from single-peaked to double peaked, the pulsed fraction reaches a minimum in all energy bands, and the hardness ratios in both NICER and GBM show a turn-over to softening as the intensity increases. This behavior repeats as the outburst rises and fades, indicating two distinct accretion regimes. These two regimes are suggestive of the accretion structure on the neutron star surface transitioning from a Coulomb collisional stopping mechanism at lower luminosities to a radiation-dominated stopping mechanism at higher luminosities. This is the highest observed (to date) value of the critical luminosity, suggesting a magnetic field of B ~1E+13 G.
[66]  oai:arXiv.org:1806.04115  [pdf] - 1713009
A NICER Spectrum of MAXI J1535-571: Near-Maximal Black Hole Spin and Potential Disk Warping
Comments: Accepted for publication in ApJ Letters
Submitted: 2018-06-11
We report on a NICER observation of the Galactic X-ray binary and stellar-mass black hole candidate, MAXI J1535-571. The source was likely observed in an "intermediate" or "very high" state, with important contributions from both an accretion disk and hard X-ray corona. The 2.3-10 keV spectrum shows clear hallmarks of relativistic disk reflection. Fits with a suitable model strongly indicate a near-maximal spin parameter of a = cJ/GM^2 = 0.994(2) and a disk that extends close to the innermost stable circular orbit, r/r_ISCO = 1.08(8) (1-sigma statistical errors). In addition to the relativistic spectrum from the innermost disk, a relatively narrow Fe K emission line is also required. The resolution of NICER reveals that the narrow line may be asymmetric, indicating a specific range of emission radii. Fits with a relativistic line model suggest an inner radius of r = 144 (+140,-60) GM/c^2 for the putative second reflection geometry; full reflection models suggest that radii a few times larger are possible. The origin of the narrow line is uncertain but a warp likely provides the most physically plausible explanation. We discuss our results in terms of the potential for NICER to reveal new features of the inner and intermediate accretion disk around black holes.
[67]  oai:arXiv.org:1801.02617  [pdf] - 1696252
The NANOGrav 11-year Data Set: Pulsar-timing Constraints On The Stochastic Gravitational-wave Background
Comments: 21 pages, 12 figures, 9 tables. Published in The Astrophysical Journal. Please send any comments/questions to S. R. Taylor (srtaylor@caltech.edu)
Submitted: 2018-01-08, last modified: 2018-06-07
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, cosmic strings, and 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. We developed an approach that bridges systematic SSE differences, producing the first PTA constraints that are robust against SSE uncertainties. We thus 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 use our constraints to characterize the combined influence on the GWB of the stellar mass-density in galactic cores, the eccentricity of SMBH binaries, and SMBH--galactic-bulge scaling relationships. We constrain cosmic-string tension using 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 published NANOGrav $9$-year constraints. Our SSE-marginalized $95\%$ upper limit on the energy density of a primordial GWB (for a radiation-dominated post-inflation Universe) is $\Omega_\mathrm{GWB}(f)h^2<3.4\times10^{-10}$.
[68]  oai:arXiv.org:1806.02342  [pdf] - 1709519
A Persistent Disk Wind in GRS 1915+105 with NICER
Comments: Accepted for publication in ApJL. Comments welcome
Submitted: 2018-06-06
The bright, erratic black hole X-ray binary GRS 1915+105 has long been a target for studies of disk instabilities, radio/infrared jets, and accretion disk winds, with implications that often apply to sources that do not exhibit its exotic X-ray variability. With the launch of NICER, we have a new opportunity to study the disk wind in GRS 1915+105 and its variability on short and long timescales. Here we present our analysis of 39 NICER observations of GRS 1915+105 collected during five months of the mission data validation and verification phase, focusing on Fe XXV and Fe XXVI absorption. We report the detection of strong Fe XXVI in 32 (>80%) of these observations, with another four marginal detections; Fe XXV is less common, but both likely arise in the well-known disk wind. We explore how the properties of this wind depends on broad characteristics of the X-ray lightcurve: mean count rate, hardness ratio, and fractional RMS variability. The trends with count rate and RMS are consistent with an average wind column density that is fairly steady between observations but varies rapidly with the source on timescales of seconds. The line dependence on spectral hardness echoes known behavior of disk winds in outbursts of Galactic black holes; these results clearly indicate that NICER is a powerful tool for studying black hole winds.
[69]  oai:arXiv.org:1805.06850  [pdf] - 1684686
A NICER look at the Aql X-1 hard state
Comments: 7 pages, 4 figures, 2 tables, accepted for publication in ApJL
Submitted: 2018-05-17
We report on a spectral-timing analysis of the neutron star low-mass X-ray binary Aql~X-1 with the Neutron Star Interior Composition Explorer (NICER) on the International Space Station. Aql~X-1 was observed with NICER during a dim outburst in 2017 July, collecting approximately $50$ ks of good exposure. The spectral and timing properties of the source correspond to that of an (hard) extreme island state in the atoll classification. We find that the fractional amplitude of the low frequency ($<0.3$ Hz) band-limited noise shows a dramatic turnover as a function of energy: it peaks at 0.5 keV with nearly 5\% rms, drops to $12\%$ rms at 2 keV, and rises to $15\%$ rms at 10 keV. Through the analysis of covariance spectra, we demonstrate that band-limited noise exists in both the soft thermal emission and the power-law emission. Additionally, we measure hard time lags, indicating the thermal emission at $0.5$ keV leads the power-law emission at 10 keV on a timescale of $\sim100$ ms at $0.3$ Hz to $\sim10$ ms at $3$ Hz. Our results demonstrate that the thermal emission in the hard state is intrinsically variable, and driving the modulation of the higher energy power-law. Interpreting the thermal spectrum as disk emission, we find our results are consistent with the disk propagation model proposed for accretion onto black holes.
[70]  oai:arXiv.org:1804.10214  [pdf] - 1682696
Detection of Reflection Features in the Neutron Star Low-Mass X-ray Binary Serpens X-1 with NICER
Comments: 6 pages, 4 figures, 1 table, accepted for publication in ApJL
Submitted: 2018-04-26, last modified: 2018-05-02
We present Neutron Star Interior Composition Explorer (NICER) observations of the neutron star low-mass X-ray binary Serpens X-1 during the early mission phase in 2017. With the high spectral sensitivity and low-energy X-ray passband of NICER, we are able to detect the Fe L line complex in addition to the signature broad, asymmetric Fe K line. We confirm the presence of these lines by comparing the NICER data to archival observations with XMM-Newton/RGS and NuSTAR. Both features originate close to the innermost stable circular orbit (ISCO). When modeling the lines with the relativistic line model RELLINE, we find the Fe L blend requires an inner disk radius of $1.4_{-0.1}^{+0.2}$ $R_{\mathrm{ISCO}}$ and Fe K is at $1.03_{-0.03}^{+0.13}$ $R_{\mathrm{ISCO}}$ (errors quoted at 90%). This corresponds to a position of $17.3_{-1.2}^{+2.5}$ km and $12.7_{-0.4}^{+1.6}$ km for a canonical neutron star mass ($M_{\mathrm{NS}}=1.4\ M_{\odot}$) and dimensionless spin value of $a=0$. Additionally, we employ a new version of the RELXILL model tailored for neutron stars and determine that these features arise from a dense disk and supersolar Fe abundance.
[71]  oai:arXiv.org:1801.01837  [pdf] - 1664520
The NANOGrav 11-year Data Set: High-precision timing of 45 Millisecond Pulsars
Comments: Published in the Astrophysical Journal Supplement Series
Submitted: 2017-12-23, last modified: 2018-04-10
We present high-precision timing data over time spans of up to 11 years for 45 millisecond pulsars observed as part of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project, aimed at detecting and characterizing low-frequency gravitational waves. The pulsars were observed with the Arecibo Observatory and/or the Green Bank Telescope at frequencies ranging from 327 MHz to 2.3 GHz. Most pulsars were observed with approximately monthly cadence, with six high--timing-precision pulsars observed weekly, and all were observed at widely separated frequencies at each observing epoch in order to fit for time-variable dispersion delays. We describe our methods for data processing, time-of-arrival (TOA) calculation, and the implementation of a new, automated method for removing outlier TOAs. We fit a timing model for each pulsar that includes spin, astrometric, and, if necessary, binary parameters, in addition to time-variable dispersion delays and parameters that quantify pulse-profile evolution with frequency. The new timing solutions provide three new parallax measurements, two new Shapiro delay measurements, and two new measurements of large orbital-period variations. We fit models that characterize sources of noise for each pulsar. We find that 11 pulsars show significant red noise, with generally smaller spectral indices than typically measured for non-recycled pulsars, possibly suggesting a different origin. Future papers will use these data to constrain or detect the signatures of gravitational-wave signals.
[72]  oai:arXiv.org:1802.09206  [pdf] - 1779499
Tests of Gravitational Symmetries with Pulsar Binary J1713+0747
Comments: 11 pages, 6 figures, and 1 table; submitted to MNRAS
Submitted: 2018-02-26
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.
[73]  oai:arXiv.org:1711.08371  [pdf] - 1736232
Searching for X-ray Pulsations from Neutron Stars Using NICER
Comments: 4 pages, to appear in Proceedings of IAU Symposium 337: Pulsar Astrophysics - The Next 50 Years
Submitted: 2017-11-22
The Neutron Star Interior Composition Explorer (NICER) presents an exciting new capability for exploring the modulation properties of X-ray emitting neutron stars, including large area, low background, extremely precise absolute event time stamps, superb low-energy response and flexible scheduling. The Pulsation Searches and Multiwavelength Coordination working group has designed a 2.5 Ms observing program to search for emission and characterize the modulation properties of about 30 known or suspected neutron star sources across a number of source categories. A key early goal will be to search for pulsations from millisecond pulsars that might exhibit thermal pulsations from the surface suitable for pulse profile modeling to constrain the neutron star equation of state. In addition, we will search for pulsations from transitional millisecond pulsars, isolated neutron stars, low-mass X-ray binaries (LMXBs), accretion-powered millisecond pulsars, central compact objects and other sources. We present our science plan and initial results from the first months of the NICER mission, including the discovery of pulsations from the millisecond pulsar J1231-1411.
[74]  oai:arXiv.org:1709.03494  [pdf] - 1938266
STROBE-X: X-ray Timing and Spectroscopy on Dynamical Timescales from Microseconds to Years
Comments: Accepted for Publication in Results in Physics
Submitted: 2017-09-08
The Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X) probes strong gravity for stellar mass to supermassive black holes and ultradense matter with unprecedented effective area, high time-resolution, and good spectral resolution, while providing a powerful time-domain X-ray observatory.
[75]  oai:arXiv.org:1612.03187  [pdf] - 1580860
The NANOGrav Nine-Year Data Set: Measurement and Interpretation of Variations in Dispersion Measures
Comments: 20 pages, 13 figures, submitted to ApJ
Submitted: 2016-12-09
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.
[76]  oai:arXiv.org:1610.01731  [pdf] - 1531913
The NANOGrav Nine-Year Data Set: Excess Noise in Millisecond Pulsar Arrival Times
Comments: 20 pages, 16 figures, to be published in ApJ
Submitted: 2016-10-06, last modified: 2016-12-02
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.
[77]  oai:arXiv.org:1603.00545  [pdf] - 1530508
The NANOGrav Nine-year Data Set: Mass and Geometric Measurements of Binary Millisecond Pulsars
Comments: Fixed typos, added Table 5 and analysis of PSR J2145-0750. Accepted by ApJ on 24 September 2016
Submitted: 2016-03-01, last modified: 2016-10-04
We analyze 24 binary radio pulsars in the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) nine-year data set. We make fourteen significant measurements of 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.10}_{-0.09}\text{M}_{\odot}$ for PSR J1918$-$0642 and as high as $m_{\rm p} = 1.928^{+0.017}_{-0.017}\text{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.21}_{-0.18}\text{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}\text{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.
[78]  oai:arXiv.org:1512.08326  [pdf] - 1488368
The NANOGrav Nine-Year Data Set: Noise Budget for Pulsar Arrival Times on Intraday Timescales
Comments: 22 pages, 11 figures, published in ApJ
Submitted: 2015-12-28, last modified: 2016-09-28
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.
[79]  oai:arXiv.org:1604.00131  [pdf] - 1447711
PSR J1024-0719: A Millisecond Pulsar in an Unusual Long-Period Orbit
Comments: Accepted by ApJ. 14 pages, 8 figures
Submitted: 2016-04-01, last modified: 2016-07-15
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.
[80]  oai:arXiv.org:1602.05570  [pdf] - 1370407
From Spin Noise to Systematics: Stochastic Processes in the First International Pulsar Timing Array Data Release
Comments: 29 pages. 16 figures. Accepted for publication in MNRAS
Submitted: 2016-02-16
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.
[81]  oai:arXiv.org:1602.03640  [pdf] - 1362729
The International Pulsar Timing Array: First Data Release
Comments: 25 pages, 6 tables, 5 figures. Accepted for publication in MNRAS
Submitted: 2016-02-11
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.
[82]  oai:arXiv.org:1601.04490  [pdf] - 1370349
The NANOGrav Nine-year Data Set: Monitoring Interstellar Scattering Delays
Comments: 14 pages, 5 figures, 3 tables. Accepted for publication in ApJ
Submitted: 2016-01-18
We report on an effort to extract and monitor interstellar scintillation parameters in regular timing observations collected for the NANOGrav pulsar timing array. Scattering delays are measured by creating dynamic spectra for each pulsar and observing epoch of wide-band observations centered near 1500 MHz and carried out at the Green Bank Telescope and the Arecibo Observatory. The ~800-MHz wide frequency bands imply dramatic changes in scintillation bandwidth across the bandpass, and a stretching routine has been included to account for this scaling. For most of the 10 pulsars for which the scaling has been measured, the bandwidths scale with frequency less steeply than expected for a Kolmogorov medium. We find estimated scattering delay values that vary with time by up to an order of magnitude. The mean measured scattering delays are similar to previously published values and slightly higher than predicted by interstellar medium models. We investigate the possibility of increasing the timing precision by mitigating timing errors introduced by the scattering delays. For most of the pulsars, the uncertainty in the time of arrival of a single timing point is much larger than the maximum variation of the scattering delay, suggesting that diffractive scintillation remains only a negligible part of their noise budget.
[83]  oai:arXiv.org:1505.07540  [pdf] - 1335431
The NANOGrav Nine-year Data Set: Observations, Arrival Time Measurements, and Analysis of 37 Millisecond Pulsars
Comments: 37 pages, 43 figures, this revision matches the final ApJ accepted version
Submitted: 2015-05-27, last modified: 2016-01-05
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.
[84]  oai:arXiv.org:1509.08982  [pdf] - 1358884
The NANOGrav Nine-year Data Set: Astrometric Measurements of 37 Millisecond Pulsars
Comments: Accepted by ApJ
Submitted: 2015-09-29, last modified: 2016-01-02
Using the nine-year radio-pulsar timing data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), collected at Arecibo Observatory and the Green Bank Telescope, we have measured the positions, proper motions, and parallaxes for 37 millisecond pulsars. We report twelve significant parallax measurements and distance measurements, and eighteen lower limits on distance. We compare these measurements to distances predicted by the NE2001 interstellar electron density model and find them to be in general agreement. We use measured orbital-decay rates and spin-down rates to confirm two of the parallax distances and to place distance upper limits on other sources; these distance limits agree with the parallax distances with one exception, PSR J1024-0719, which we discuss at length. Using the proper motions of the 37 NANOGrav pulsars in combination with other published measurements, we calculate the velocity dispersion of the millisecond pulsar population in Galactocentric coordinates. We find the radial, azimuthal, and perpendicular dispersions to be 46, 40, and 24 km s-1, respectively, in a model that allows for high-velocity outliers; or 81, 58, and 62 km s-1 for the full population. These velocity dispersions are far smaller than those of the canonical pulsar population, and are similar to older Galactic disk populations. This suggests that millisecond pulsar velocities are largely attributable to their being an old population rather than being artifacts of their birth and evolution as neutron star binary systems. The components of these velocity dispersions follow similar proportions to other Galactic populations, suggesting that our results are not biased by selection effects.
[85]  oai:arXiv.org:1512.03067  [pdf] - 1528011
Measuring Neutron Star Radii via Pulse Profile Modeling with NICER
Comments: Submitted to ApJ
Submitted: 2015-12-09
The Neutron-star Interior Composition Explorer (NICER) is an X-ray astrophysics payload that will be placed on the International Space Station. Its primary science goal is to measure with high accuracy the pulse profiles that arise from the non-uniform thermal surface emission of rotation-powered pulsars. Modeling general relativistic effects on the profiles will lead to measuring the radii of these neutron stars and to constraining their equation of state. Achieving this goal will depend, among other things, on accurate knowledge of the source, sky, and instrument backgrounds. We use here simple analytic estimates to quantify the level at which these backgrounds need to be known in order for the upcoming measurements to provide significant constraints on the properties of neutron stars. We show that, even in the minimal-information scenario, knowledge of the background at a few percent level for a background-to-source countrate ratio of 0.2 allows for a measurement of the neutron star compactness to better than 10% uncertainty for most of the parameter space. These constraints improve further when more realistic assumptions are made about the neutron star emission and spin, and when additional information about the source itself, such as its mass or distance, are incorporated.
[86]  oai:arXiv.org:1504.00662  [pdf] - 1264695
Testing Theories of Gravitation Using 21-Year Timing of Pulsar Binary J1713+0747
Comments: 14 pages, 6 figures. Published on ApJ, 809, 41 (2015)
Submitted: 2015-04-02, last modified: 2015-08-19
We report 21-yr 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 correlated noise such as jitter and red noise. Its timing residuals have weighted root mean square $\sim 92$ ns. The new dataset allows us to update and improve previous measurements of the system properties, including the masses of the neutron star ($1.31\pm0.11$ $M_{\odot}$) and the companion white dwarf ($0.286\pm0.012$ $M_{\odot}$) and the parallax distance $1.15\pm0.03$ kpc. We measured the intrinsic change in orbital period, $\dot{P}^{\rm Int}_{\rm b}$, is $-0.20\pm0.17$ ps s$^{-1}$, which is not distinguishable from zero. This result, combined with the measured $\dot{P}^{\rm Int}_{\rm b}$ of other pulsars, can place a generic limit on potential changes in the gravitational constant $G$. We found that $\dot{G}/G$ is consistent with zero [$(-0.6\pm1.1)\times10^{-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 $\dot{G}/G$ limit from pulsar binary systems. The $\dot{P}^{\rm Int}_{\rm b}$ of pulsar binaries can also place limits on the putative coupling constant for dipole gravitational radiation $\kappa_D=(-0.9\pm3.3)\times10^{-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 $\hat{\alpha}_3$, which describes a breaking of both Lorentz invariance in gravitation and conservation of momentum. We found, at 95\% confidence, $\Delta<0.01$ and $\hat{\alpha}_3<2\times10^{-20}$ based on PSR J1713+0747.
[87]  oai:arXiv.org:1508.03024  [pdf] - 1387105
The NANOGrav Nine-year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background
Comments: 21 pages, 12 figures. Please send any comments/questions to justin.ellis18@gmail.com
Submitted: 2015-08-12
We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. We set upper limits for a GWB from supermassive black hole binaries under power law, broken power law, and free spectral coefficient GW spectrum models. We place a 95\% upper limit on the strain amplitude (at a frequency of yr$^{-1}$) in the power law model of $A_{\rm gw} < 1.5\times 10^{-15}$. For a broken power law model, we place priors on the strain amplitude derived from simulations of Sesana (2013) and McWilliams et al. (2014). We find that the data favor a broken power law to a pure power law with odds ratios of 22 and 2.2 to one for the McWilliams and Sesana prior models, respectively. The McWilliams model is essentially ruled out by the data, and the Sesana model is in tension with the data under the assumption of a pure power law. Using the broken power-law analysis we construct posterior distributions on environmental factors that drive the binary to the GW-driven regime including the stellar mass density for stellar-scattering, mass accretion rate for circumbinary disk interaction, and orbital eccentricity for eccentric binaries, marking the first time that the shape of the GWB spectrum has been used to make astrophysical inferences. We then place the most stringent limits so far on the energy density of relic GWs, $\Omega_\mathrm{gw}(f)\,h^2 < 4.2 \times 10^{-10}$, yielding a limit on the Hubble parameter during inflation of $H_*=1.6\times10^{-2}~m_{Pl}$, where $m_{Pl}$ is the Planck mass. Our limit on the cosmic string GWB, $\Omega_\mathrm{gw}(f)\, h^2 < 2.2 \times 10^{-10}$, translates to a conservative limit of $G\mu<3.3\times 10^{-8}$ - a factor of 4 better than the joint Planck and high-$l$ CMB data from other experiments.
[88]  oai:arXiv.org:1501.05343  [pdf] - 1280782
NANOGrav Constraints on Gravitational Wave Bursts with Memory
Comments: 12 pages, 7 figures. We will submit this article to ApJ two weeks from now. Questions or comments should be directed to Dustin Madison (drm252@cornell.edu)
Submitted: 2015-01-21
Among efforts to detect gravitational radiation, pulsar timing arrays are uniquely poised to detect "memory" signatures, permanent perturbations in spacetime from highly energetic astrophysical events such as mergers of supermassive black hole binaries. The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) observes dozens of the most stable millisecond pulsars using the Arecibo and Green Bank radio telescopes in an effort to study, among other things, gravitational wave memory. We herein present the results of a search for gravitational wave bursts with memory (BWMs) using the first five years of NANOGrav observations. We develop original methods for dramatically speeding up searches for BWM signals. In the directions of the sky where our sensitivity to BWMs is best, we would detect mergers of binaries with reduced masses of $10^9$ $M_\odot$ out to distances of 30 Mpc; such massive mergers in the Virgo cluster would be marginally detectable. We find no evidence for BWMs. However, with our non-detection, we set upper limits on the rate at which BWMs of various amplitudes could have occurred during the time spanned by our data--e.g., BWMs with amplitudes greater than $10^{-13}$ must occur at a rate less than 1.5 yr$^{-1}$.
[89]  oai:arXiv.org:1404.1267  [pdf] - 1208811
NANOGrav Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries in Circular Orbits
Comments: 21 pages, 9 figures. Matches version submitted to ApJ
Submitted: 2014-04-04, last modified: 2014-05-19
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) project currently observes 43 pulsars using the Green Bank and Arecibo radio telescopes. In this work we use a subset of 17 pulsars timed for a span of roughly five years (2005--2010). We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Within the timing data, we perform a search for continuous gravitational waves from individual supermassive black hole binaries in circular orbits using robust frequentist and Bayesian techniques. We find that there is no evidence for the presence of a detectable continuous gravitational wave; however, we can use these data to place the most constraining upper limits to date on the strength of such gravitational waves. Using the full 17 pulsar dataset we place a 95% upper limit on the sky-averaged strain amplitude of $h_0\lesssim 3.8\times 10^{-14}$ at a frequency of 10 nHz. Furthermore, we place 95% \emph{all sky} lower limits on the luminosity distance to such gravitational wave sources finding that the $d_L \gtrsim 425$ Mpc for sources at a frequency of 10 nHz and chirp mass $10^{10}{\rm M}_{\odot}$. We find that for gravitational wave sources near our best timed pulsars in the sky, the sensitivity of the pulsar timing array is increased by a factor of $\sim$4 over the sky-averaged sensitivity. Finally we place limits on the coalescence rate of the most massive supermassive black hole binaries.
[90]  oai:arXiv.org:1206.3324  [pdf] - 525400
Fermi-LAT constraints on the Pulsar Wind Nebula nature of HESS J1857+026
Comments: A&A - Year: 2012 Accepted: 05/23/12
Submitted: 2012-06-14
Since its launch, the Fermi satellite has firmly identified 5 pulsar wind nebulae plus a large number of candidates, all powered by young and energetic pulsars. HESS J1857+026 is a spatially extended gamma-ray source detected by H.E.S.S. and classified as a possible pulsar wind nebula candidate powered by PSR J1856+0245. We search for gamma-ray pulsations from PSR J1856+0245 and explore the characteristics of its associated pulsar wind nebula. Using a rotational ephemeris obtained from the Lovell telescope at Jodrell Bank Observatory at 1.5 GHz, we phase-fold 36 months of gamma-ray data acquired by the Large Area Telescope (LAT) aboard Fermi. We also perform a complete gamma-ray spectral and morphological analysis. No gamma-ray pulsations were detected from PSR J1856+0245. However, significant emission is detected at a position coincident with the TeV source HESS J1857+026. The gamma-ray spectrum is well described by a simple power-law with a spectral index of 1.53 \pm 0.11_{\rm stat} \pm 0.55_{\rm syst}$ and an energy flux of $G(0.1$--100 GeV$)=(2.71 \pm 0.52_{\rm stat} \pm 1.51_{\rm syst}) \times 10^{-11}$ ergs cm$^{-2}$ s$^{-1}$. The $\gamma$-ray luminosity is $L_{PWN}^{\gamma} (0.1$--100 GeV$)=(2.5 \pm 0.5_{stat} \pm 1.5_{syst}) \times 10^{35} (\frac{d}{9 kpc})^2$ ergs s$^{-1}$, assuming a distance of 9 kpc. This implies a $\gamma-$ray efficiency of $\sim$ 5% for $\dot{E}=4.6 \times 10^{36}$ erg $s^{-1}$, in the range expected for pulsar wind nebulae. Detailed multi-wavelength modeling provides new constraints on its pulsar wind nebula nature.
[91]  oai:arXiv.org:1201.6641  [pdf] - 1116285
Limits on the Stochastic Gravitational Wave Background from the North American Nanohertz Observatory for Gravitational Waves
Comments: To be submitted to ApJ
Submitted: 2012-01-31
We present an analysis of high-precision pulsar timing data taken as part of the North American Nanohertz Observatory for Gravitational waves (NANOGrav) project. We have observed 17 pulsars for a span of roughly five years using the Green Bank and Arecibo radio telescopes. We analyze these data using standard pulsar timing models, with the addition of time-variable dispersion measure and frequency-variable pulse shape terms. Sub-microsecond timing residuals are obtained in nearly all cases, and the best root-mean-square timing residuals in this set are ~30-50 ns. We present methods for analyzing post-fit timing residuals for the presence of a gravitational wave signal with a specified spectral shape. These optimally take into account the timing fluctuation power removed by the model fit, and can be applied to either data from a single pulsar, or to a set of pulsars to detect a correlated signal. We apply these methods to our dataset to set an upper limit on the strength of the nHz-frequency stochastic supermassive black hole gravitational wave background of h_c (1 yr^-1) < 7x10^-15 (95%). This result is dominated by the timing of the two best pulsars in the set, PSRs J1713+0747 and J1909-3744.
[92]  oai:arXiv.org:1105.3185  [pdf] - 1076662
Discovery of an Energetic Pulsar Associated with SNR G76.9+1.0
Comments: 8 pages in emulateapj format. Minor changes (including a shortened abstract) to reflect the version accepted for publication
Submitted: 2011-05-16, last modified: 2011-06-10
We report the discovery of PSR J2022+3842, a 24 ms radio and X-ray pulsar in the supernova remnant G76.9+1.0, in observations with the Chandra X-ray telescope, the Robert C. Byrd Green Bank Radio Telescope, and the Rossi X-ray Timing Explorer (RXTE). The pulsar's spin-down rate implies a rotation-powered luminosity Edot = 1.2 x 10^{38} erg/s, a surface dipole magnetic field strength B_s = 1.0 x 10^{12} G, and a characteristic age of 8.9 kyr. PSR J2022+3842 is thus the second-most energetic Galactic pulsar known, after the Crab, as well as the most rapidly-rotating young, radio-bright pulsar known. The radio pulsations are highly dispersed and broadened by interstellar scattering, and we find that a large (delta-f / f ~= 1.9 x 10^{-6}) spin glitch must have occurred between our discovery and confirmation observations. The X-ray pulses are narrow (0.06 cycles FWHM) and visible up to 20 keV, consistent with magnetospheric emission from a rotation-powered pulsar. The Chandra X-ray image identifies the pulsar with a hard, unresolved source at the midpoint of the double-lobed radio morphology of SNR G76.9+1.0 and embedded within faint, compact X-ray nebulosity. The spatial relationship of the X-ray and radio emissions is remarkably similar to extended structure seen around the Vela pulsar. The combined Chandra and RXTE pulsar spectrum is well-fitted by an absorbed power-law model with column density N_H = (1.7\pm0.3) x 10^{22} cm^{-2} and photon index Gamma = 1.0\pm0.2; it implies that the Chandra point-source flux is virtually 100% pulsed. For a distance of 10 kpc, the X-ray luminosity of PSR J2022+3842 is L_X(2-10 keV) = 7.0 x 10^{33} erg s^{-1}. Despite being extraordinarily energetic, PSR J2022+3842 lacks a bright X-ray wind nebula and has an unusually low conversion efficiency of spin-down power to X-ray luminosity, L_X/Edot = 5.9 x 10^{-5}.
[93]  oai:arXiv.org:1012.3507  [pdf] - 1042608
Pulsars with the Australian Square Kilometre Array Pathfinder
Comments: To appear in proceedings of "Radio Pulsars: An astrophysical key to unlock the secrets of the Universe"
Submitted: 2010-12-15
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.
[94]  oai:arXiv.org:0911.5206  [pdf] - 172800
The international pulsar timing array project: using pulsars as a gravitational wave detector
Comments: accepted by Classical and Quantum Gravity. Review talk for the Amaldi8 conference series
Submitted: 2009-11-27
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.
[95]  oai:arXiv.org:0909.1058  [pdf] - 28042
The North American Nanohertz Observatory for Gravitational Waves
Comments: This document is the NANOGrav consortium's submission to Astro2010's Program Prioritization Panel on Particle Astrophysics and Gravitation
Submitted: 2009-09-05
The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) is a consortium of astronomers whose goal is the creation of a galactic scale gravitational wave observatory sensitive to gravitational waves in the nHz-microHz band. It is just one component of an international collaboration involving similar organizations of European and Australian astronomers who share the same goal. Gravitational waves, a prediction of Einstein's general theory of relativity, are a phenomenon of dynamical space-time generated by the bulk motion of matter, and the dynamics of space-time itself. They are detectable by the small disturbance they cause in the light travel time between some light source and an observer. NANOGrav exploits radio pulsars as both the light (radio) source and the clock against which the light travel time is measured. In an array of radio pulsars gravitational waves manifest themselves as correlated disturbances in the pulse arrival times. The timing precision of today's best measured pulsars is less than 100 ns. With improved instrumentation and signal-to-noise it is widely believed that the next decade could see a pulsar timing network of 100 pulsars each with better than 100 ns timing precision. Such a pulsar timing array (PTA), observed with a regular cadence of days to weeks, would be capable of observing supermassive black hole binaries following galactic mergers, relic radiation from early universe phenomena such as cosmic strings, cosmic superstrings, or inflation, and more generally providing a vantage on the universe whose revolutionary potential has not been seen in the 400 years since Galileo first turned a telescope to the heavens.
[96]  oai:arXiv.org:0902.3264  [pdf] - 21536
X-ray Timing of Neutron Stars, Astrophysical Probes of Extreme Physics
Comments: 8 pages. Science white paper submitted to Astro2010: The Astronomy and Astrophysics Decadal Survey
Submitted: 2009-02-18
The characteristic physical timescales near stellar-mass compact objects are measured in milliseconds. These timescales -- the free-fall time, the fastest stable orbital period, and stellar spin periods -- encode the fundamental physical properties of compact objects: mass, radius, and angular momentum. The characteristic temperature of matter in the vicinity of neutron stars is such that the principal electromagnetic window into their realms is the X-ray band. Because of these connections to the fundamental properties of neutron stars, X-ray timing studies remain today the most direct means of probing their structure and dynamics. While current X-ray observatories have revealed many relevant and fascinating phenomena, they lack the sensitivity to fully exploit them to uncover the fundamental properties of compact objects and their extreme physics. With this white paper, we summarize and highlight the science opportunities that will accompany an order-of-magnitude improvement in X-ray timing sensitivity, a goal attainable in the coming decade.
[97]  oai:arXiv.org:0902.2891  [pdf] - 21455
Constraining the Bulk Properties of Dense Matter by Measuring Millisecond Pulsar Masses - A White Paper for the Astronomy and Astrophysics Decadal Survey, CFP Panel
Comments: & pages + cover page, 3 figures
Submitted: 2009-02-17
More than four decades after the discovery of pulsars, the composition of matter at their cores is still a mystery. This white paper summarizes how recent high-precision measurements of millisecond pulsar masses have introduced new experimental constraints on the properties of super-dense matter, and how continued timing of intriguing new objects, coupled with radio telescope surveys to discover more pulsars, might introduce significantly more stringent constraints.
[98]  oai:arXiv.org:0807.0811  [pdf] - 14227
DA495 - an aging pulsar wind nebula
Comments: 25 pages, accepted by ApJ
Submitted: 2008-07-04, last modified: 2008-07-08
We present a radio continuum study of the pulsar wind nebula (PWN) DA 495 (G65.7+1.2), including images of total intensity and linear polarization from 408 to 10550 MHz based on the Canadian Galactic Plane Survey and observations with the Effelsberg 100-m Radio Telescope. Removal of flux density contributions from a superimposed \ion{H}{2} region and from compact extragalactic sources reveals a break in the spectrum of DA 495 at 1.3 GHz, with a spectral index ${\alpha}={-0.45 \pm 0.20}$ below the break and ${\alpha}={-0.87 \pm 0.10}$ above it (${S}_\nu \propto{\nu^{\alpha}}$). The spectral break is more than three times lower in frequency than the lowest break detected in any other PWN. The break in the spectrum is likely the result of synchrotron cooling, and DA 495, at an age of $\sim$20,000 yr, may have evolved from an object similar to the Vela X nebula, with a similarly energetic pulsar. We find a magnetic field of $\sim$1.3 mG inside the nebula. After correcting for the resulting high internal rotation measure, the magnetic field structure is quite simple, resembling the inner part of a dipole field projected onto the plane of the sky, although a toroidal component is likely also present. The dipole field axis, which should be parallel to the spin axis of the putative pulsar, lies at an angle of ${\sim}50\degr$ east of the North Celestial Pole and is pointing away from us towards the south-west. The upper limit for the radio surface brightness of any shell-type supernova remnant emission around DA 495 is $\Sigma_{1 GHz} \sim 5.4 \times 10^{-23}$ OAWatt m$^{-2}$ Hz$^{-1}$ sr$^{-1}$ (assuming a radio spectral index of $\alpha = -0.5$), lower than the faintest shell-type remnant known to date.
[99]  oai:arXiv.org:0806.3766  [pdf] - 13820
Chandra Confirmation of a Pulsar Wind Nebula in DA 495
Comments: 26 pages, 5 figures, AASTeX preprint style. Accepted for publication in The Astrophysical Journal
Submitted: 2008-06-23
As part of a multiwavelength study of the unusual radio supernova remnant DA 495, we present observations made with the Chandra X-ray Observatory. Imaging and spectroscopic analysis confirms the previously detected X-ray source at the heart of the annular radio nebula, establishing the radiative properties of two key emission components: a soft unresolved source with a blackbody temperature of 1 MK consistent with a neutron star, surrounded by a nonthermal nebula 40'' in diameter exhibiting a power-law spectrum with photon index Gamma = 1.6+/-0.3, typical of a pulsar wind nebula. The implied spin-down luminosity of the neutron star, assuming a conversion efficiency to nebular flux appropriate to Vela-like pulsars, is ~10^{35} ergs/s, again typical of objects a few tens of kyr old. Morphologically, the nebular flux is slightly enhanced along a direction, in projection on the sky, independently demonstrated to be of significance in radio polarization observations; we argue that this represents the orientation of the pulsar spin axis. At smaller scales, a narrow X-ray feature is seen extending out 5'' from the point source, a distance consistent with the sizes of resolved wind termination shocks around many Vela-like pulsars. Finally, we argue based on synchrotron lifetimes in the estimated nebular magnetic field that DA 495 represents a rare pulsar wind nebula in which electromagnetic flux makes up a significant part, together with particle flux, of the neutron star's wind, and that this high magnetization factor may account for the nebula's low luminosity.
[100]  oai:arXiv.org:0805.2396  [pdf] - 12688
An Eccentric Binary Millisecond Pulsar in the Galactic Plane
Comments: 28 pages, 4 figures inc Supplementary On-Line Material. Accepted for publication in Science, published on Science Express: 10.1126/science.1157580
Submitted: 2008-05-15
Binary pulsar systems are superb probes of stellar and binary evolution and the physics of extreme environments. In a survey with the Arecibo telescope, we have found PSR J1903+0327, a radio pulsar with a rotational period of 2.15 ms in a highly eccentric (e = 0.44) 95-day orbit around a solar mass companion. Infrared observations identify a possible main-sequence companion star. Conventional binary stellar evolution models predict neither large orbital eccentricities nor main-sequence companions around millisecond pulsars. Alternative formation scenarios involve recycling a neutron star in a globular cluster then ejecting it into the Galactic disk or membership in a hierarchical triple system. A relativistic analysis of timing observations of the pulsar finds its mass to be 1.74+/-0.04 Msun, an unusually high value.
[101]  oai:arXiv.org:astro-ph/0701103  [pdf] - 88232
Correlations between pulsed X-ray flux and radio arrival time in the Vela pulsar
Comments: Proceedings of the 363. WE-Heraeus Seminar on: Neutron Stars and Pulsars (Posters and contributed talks) Physikzentrum Bad Honnef, Germany, May.14-19, 2006, eds. W.Becker, H.H.Huang, MPE Report 291, pp.16-19
Submitted: 2007-01-04
We report the results of simultaneous observations of the Vela pulsar in X-rays and radio from the RXTE satellite and the Mount Pleasant Radio Observatory in Tasmania. We sought correlations between the Vela's X-ray and radio flux densities and radio arrival times on a pulse by pulse basis. We found significantly higher flux density in Vela's main X-ray peak during radio pulses that arrived early. This excess flux shifts to the 'trough' following the 2nd X-ray peak during radio pulses that arrive later. We suggest that the mechanism producing the radio pulses is intimately connected to the mechanism producing X-rays. Current models using resonant absorption in the outer magnetosphere as a cause of the radio emission, and less directly of the X-ray emission, are explored as a possible explanation for the correlation.
[102]  oai:arXiv.org:astro-ph/0611450  [pdf] - 86836
Correlation between X-ray Lightcurve Shape and Radio Arrival Time in the Vela Pulsar
Comments: 6 pages, 5 figures, accepted by ApJ
Submitted: 2006-11-14
We report the results of simultaneous observations of the Vela pulsar in X-rays and radio from the RXTE satellite and the Mount Pleasant Radio Observatory in Tasmania. We sought correlations between the Vela's X-ray emission and radio arrival times on a pulse by pulse basis. At a confidence level of 99.8% we have found significantly higher flux density in Vela's main X-ray peak during radio pulses that arrived early. This excess flux shifts to the 'trough' following the 2nd X-ray peak during radio pulses that arrive later. Our results suggest that the mechanism producing the radio pulses is intimately connected to the mechanism producing X-rays. Current models using resonant absorption of radio emission in the outer magnetosphere as a cause of the X-ray emission are explored as a possible explanation for the correlation.
[103]  oai:arXiv.org:astro-ph/0604392  [pdf] - 81461
Arecibo and the ALFA Pulsar Survey
Comments: 8 pages, invited review for the Hanas Pulsar Symposium, ChJAA in press
Submitted: 2006-04-19
The recently started Arecibo L-band Feed Array (ALFA) pulsar survey aims to find ~1000 new pulsars. Due to its high time and frequency resolution the survey is especially sensitive to millisecond pulsars, which have the potential to test gravitational theories, detect gravitational waves and probe the neutron-star equation of state. Here we report the results of our preliminary analysis: in the first months we have discovered 21 new pulsars. One of these, PSR J1906+0746, is a young 144-ms pulsar in a highly relativistic 3.98-hr low-eccentricity orbit. The 2.61 +- 0.02 solar-mass system is expected to coalesce in ~300 Myr and contributes significantly to the computed cosmic inspiral rate of compact binary systems.
[104]  oai:arXiv.org:physics/0601012  [pdf] - 119140
Development of Ground-testable Phase Fresnel Lenses in Silicon
Comments: Presented at GammaWave05: "Focusing Telescopes in Nuclear Astrophysics", Bonifacio, Corsica, September 2005, to be published in Experimental Astronomy, 8 pages, 3 figures
Submitted: 2006-01-03
Diffractive/refractive optics, such as Phase Fresnel Lenses (PFL's), offer the potential to achieve excellent imaging performance in the x-ray and gamma-ray photon regimes. In principle, the angular resolution obtained with these devices can be diffraction limited. Furthermore, improvements in signal sensitivity can be achieved as virtually the entire flux incident on a lens can be concentrated onto a small detector area. In order to verify experimentally the imaging performance, we have fabricated PFL's in silicon using gray-scale lithography to produce the required Fresnel profile. These devices are to be evaluated in the recently constructed 600-meter x-ray interferometry testbed at NASA/GSFC. Profile measurements of the Fresnel structures in fabricated PFL's have been performed and have been used to obtain initial characterization of the expected PFL imaging efficiencies.
[105]  oai:arXiv.org:astro-ph/0511523  [pdf] - 77881
The young, highly relativistic binary pulsar J1906+0746
Comments: 19 pages, 3 figures, accepted for publication in ApJ
Submitted: 2005-11-17
We report the discovery of PSR J1906+0746, a young 144-ms pulsar in a highly relativistic 3.98-hr orbit with an eccentricity of 0.085 and expected gravitational wave coalescence time of 300 Myr. The new pulsar was found during precursor survey observations with the Arecibo 1.4-GHz feed array system and retrospectively detected in the Parkes Multibeam plane pulsar survey data. From radio follow-up observations with Arecibo, Jodrell Bank, Green Bank, and Parkes, we have measured the spin-down and binary parameters of the pulsar and its basic spectral and polarization properties. We also present evidence for pulse profile evolution, which is likely due to geodetic precession, a relativistic effect caused by the misalignment of the pulsar spin and total angular momentum vectors. Our measurements show that PSR J1906+0746 is a young object with a characteristic age of 112 kyr. From the measured rate of orbital periastron advance 7.57+/-0.03 deg/yr, we infer a total system mass of 2.61+/-0.02 Msun. While these parameters suggest that the PSR J1906+0746 binary system might be a younger version of the double pulsar system, intensive searches for radio pulses from the companion have so far been unsuccessful. It is therefore not known whether the companion is another neutron star or a massive white dwarf. Regardless of the nature of the companion, a simple calculation suggests that the Galactic birth rate of binaries similar to PSR J1906+0746 is 60/Myr. This implies that PSR J1906+0746 will make a significant contribution to the computed cosmic inspiral rate of compact binary systems.
[106]  oai:arXiv.org:astro-ph/0509732  [pdf] - 76313
Arecibo Pulsar Survey Using ALFA. I. Survey Strategy and First Discoveries
Comments: 10 pp, 9 figures, accepted by the Astrophysical Journal
Submitted: 2005-09-23
We report results from the initial stage of a long-term pulsar survey of the Galactic plane using the Arecibo L-band Feed Array (ALFA), a seven-beam receiver operating at 1.4 GHz with 0.3 GHz bandwidth. The search targets Galactic latitudes |b| < 5 deg in the longitude ranges 32 deg < l < 77 deg and 168 deg < l < 77 deg. Data discussed here were collected over a 100 MHz passband centered on 1.42 GHz using a spectrometer that recorded 256 channels every 64 microsec. In a preliminary, standard period-DM analysis, we have detected 29 previously known pulsars and discovered 11 new ones. One of these, with a period of 69 ms and a low characteristic age of 82 kyr, is a plausible candidate for association with the unidentified EGRET source 3EG J1928+1733. Another is a non-recycled pulsar in a relativistic binary with orbital period of 3.98 hr. We also search the data for isolated dispersed pulses, a technique that yielded discovery of an extremely sporadic radio emitter with a spin period of 1.2 s. Simulations we have carried out indicate that about 1000 new pulsars will be found in the ALFA survey. In addition to providing a large sample for use in population analyses and for probing the magnetoionic interstellar medium, the survey maximizes the chances of finding rapidly spinning millisecond pulsars and pulsars in compact binary systems. Our search algorithms will exploit the multiple data streams from ALFA to discriminate between radio frequency interference and celestial signals, including pulsars and possibly new classes of transient radio sources.
[107]  oai:arXiv.org:astro-ph/0508320  [pdf] - 75209
Arecibo timing and single-pulse observations of 17 pulsars
Comments: 9 pages, 6 figures. Accepted for publication in MNRAS
Submitted: 2005-08-15, last modified: 2005-08-16
We report on timing and single-pulse observations of 17 pulsars discovered at the Arecibo observatory. The highlights of our sample are the recycled pulsars J1829+2456, J1944+0907 and the drifting subpulses observed in PSR J0815+0939. For the double neutron star binary J1829+2456, in addition to improving upon our existing measurement of relativistic periastron advance, we have now measured the pulsar's spin period derivative. This new result sets an upper limit on the transverse speed of 120 km/s and a lower limit on the characteristic age of 12.4 Gyr. From our measurement of proper motion of the isolated 5.2-ms pulsar J1944+0907, we infer a transverse speed of 188 +/- 65 km/s. This is higher than that of any other isolated millisecond pulsar. An estimate of the speed, using interstellar scintillation, of 235 +/- 45 km/s indicates that the scattering medium along the line of sight is non-uniform. We discuss the drifting subpulses detected from three pulsars in the sample, in particular the remarkable drifting subpulse properties of the 645-ms pulsar J0815+0939. Drifting is observed in all four components of the pulse profile, with the sense of drift varying among the different components. This unusual `bi-drifting'' behaviour challenges standard explanations of the drifting subpulse phenomenon.
[108]  oai:arXiv.org:astro-ph/0411208  [pdf] - 68807
GBT Observations of Very Low Mass Binary Millisecond Pulsars: A Search for Eclipses
Comments: 6 pages, to appear in "Binary Radio Pulsars", eds F. A. Rasio & I. H. Stairs, PASP conf series
Submitted: 2004-11-08
We searched for eclipses of two millisecond pulsars, PSR J1807-2459 and PSR B1908+00. These pulsars are in very low mass binary systems with orbital parameters similar to those of eclipsing binaries. Observations were made with the GBT at frequencies as low as 575 MHz. No eclipses were detected in either system. Observations of well-established eclipsing binary J2051-0827 found eclipses to be substantially weaker than previously seen, with the pulsar detected throughout the eclipse region at 575 MHz and with an electron column density an order of magnitude smaller than previously measured.
[109]  oai:arXiv.org:astro-ph/0408457  [pdf] - 66948
Measurement of gravitational spin-orbit coupling in a binary pulsar system
Comments: 4 pages, accepted to PRL. Version with high-resolution figure 2 available at http://www.astro.ubc.ca/people/stairs/papers/sta04b.ps.gz
Submitted: 2004-08-24
In relativistic gravity, a spinning pulsar will precess as it orbits a compact companion star. We have measured the effect of such precession on the average shape and polarization of the radiation from PSR B1534+12. We have also detected, with limited precision, special-relativistic aberration of the revolving pulsar beam due to orbital motion. Our observations fix the system geometry, including the misalignment between the spin and orbital angular momenta, and yield a measurement of the precession timescale consistent with the predictions of General Relativity.
[110]  oai:arXiv.org:astro-ph/0405146  [pdf] - 64666
A multi-wavelength search for a counterpart of the unidentified gamma-ray source 3EG J2020+4017 (2CG078+2)
Comments: Accepted for publication in ApJ. A PDF version of the accepted paper with higher resolution images is available at ftp://www.xray.mpe.mpg.de/people/web/gamma-Cygni
Submitted: 2004-05-07, last modified: 2004-07-08
In search of the counterpart to the brightest unidentified gamma-ray source 3EG J2020+4017 we report on new X-ray and radio observations of the gamma-Cygni field with the Chandra X-ray Observatory and with the Green Bank Telescope. We also report on reanalysis of archival ROSAT data. With Chandra it became possible for the first time to measure the position of the putative gamma-ray counterpart RX J2020.2+4026 with sub-arcsec accuracy and to deduce its X-ray spectral characteristics. These observations demonstrate that RX J2020.2+4026 is associated with a K field star and therefore is unlikely to be the counterpart of the bright gamma-ray source 2CG078+2 in the SNR G78.2+2.1 as had been previously suggested. The Chandra observation detected 37 additional X-ray sources which were correlated with catalogs of optical and infrared data. Subsequent GBT radio observations covered the complete 99% EGRET likelihood contour of 3EG J2020+4017 with a sensitivity limit of L_820 ~ 0.1 mJy kpc^2 which is lower than most of the recent deep radio search limits. If there is a pulsar operating in 3EG J2020+4017 this sensitivity limit suggests that the pulsar either does not produce significant amounts of radio emission or that its geometry is such that the radio beam does not intersect with the line of sight. Finally, reanalysis of archival ROSAT data leads to a flux upper limit of f_x < 1.8 x 10^-13 erg s^-1 cm^-2 for a putative point-like X-ray source located within the 68% confidence contour of 3EG J2020+4017. Adopting the SNR age of 5400 yrs and assuming a spin-down to X-ray energy conversion factor of 10^-4 this upper limit constrains the parameters of a putative neutron star as a counterpart for 3EG J2020+4017 to be P >= 160 (d/1.5 kpc)^-1 ms, P_dot >= 5 x 10^-13 (d/1.5 kpc)^-1 s s^-1 and B_perp >= 9 x 10^12 (d/1.5 kpc)^-1 G.
[111]  oai:arXiv.org:astro-ph/0404181  [pdf] - 64089
New Binary and Millisecond Pulsars from Arecibo Drift-Scan Searches
Comments: 6 pages, 3 figures, to appear in proceedings of Aspen Center for Physics Conference on ``Binary Radio Pulsars'' Eds. F. Rasio and I. Stairs
Submitted: 2004-04-08
We discuss four recycled pulsars found in Arecibo drift-scan searches. PSR J1944+0907 has a spin period of 5.2 ms and is isolated. The 5.8-ms pulsar J1453+19 may have a low-mass companion. We discuss these pulsars in the context of isolated millisecond pulsar formation and the minimum spin period of neutron stars. The isolated 56-ms pulsar J0609+2130 is possibly the remnant of a disrupted double neutron star binary. The 41-ms pulsar J1829+2456 is in a relativistic orbit. Its companion is most likely another neutron star, making this the eighth known double neutron star binary system.
[112]  oai:arXiv.org:astro-ph/0403553  [pdf] - 63744
PSR J1829+2456: a relativistic binary pulsar
Comments: Accepted MNRAS, 5 pages, 3 figures
Submitted: 2004-03-23
We report the discovery of a new binary pulsar, PSR J1829+2456, found during a mid-latitude drift-scan survey with the Arecibo telescope. Our initial timing observations show the 41-ms pulsar to be in a 28-hr, slightly eccentric, binary orbit. The advance of periastron, omegadot = 0.28 +/- 0.01 deg/yr is derived from our timing observations spanning 200 days. Assuming that the advance of periastron is purely relativistic and a reasonable range of neutron star masses for PSR J1829+2456 we constrain the companion mass to be between 1.22 Msun and 1.38 Msun, making it likely to be another neutron star. We also place a firm upper limit on the pulsar mass of 1.38 Msun. The expected coalescence time due to gravitational-wave emission is long (~60 Gyr) and this system will not significantly impact upon calculations of merger rates that are relevant to upcoming instruments such as LIGO.
[113]  oai:arXiv.org:astro-ph/0312044  [pdf] - 61257
Pulsar Parallaxes at 5 GHz with the Very Long Baseline Array
Comments: 8 pages, including 3 figures and 3 tables; emulateapj; ApJ submitted
Submitted: 2003-12-01
We present the first pulsar parallaxes measured with phase-referenced pulsar VLBI observations at 5 GHz. Due to the steep spectra of pulsars, previous astrometric measurements have been at lower frequencies. However, the strongest pulsars can be observed at 5 GHz, offering the benefit of lower combined ionospheric and tropospheric phase errors, which usually limit VLBI astrometric accuracy. The pulsars B0329+54, B0355+54 and B1929+10 were observed for 7 epochs spread evenly over 2 years. For B0329+54, large systematic errors lead to only an upper limit on the parallax (pi < 1.5 mas). A new proper motion and parallax were measured for B0355+54 (pi = 0.91 +- 0.16 mas), implying a distance of 1.04+0.21-0.16 kpc and a transverse velocity of 61+12-9 km/s. The parallax and proper motion for B1929+10 were significantly improved (pi = 2.77 +- 0.07 mas), yielding a distance of 361+10-8 pc and a transverse velocity of 177+4-5 km/s. We demonstrate that the astrometric errors are correlated with the angular separation between the phase reference calibrator and the target source, with significantly lower errors at 5 GHz compared to 1.6 GHz. Finally, based on our new distance determinations for B1929+10 and B0355+54, we derive or constrain the luminosities of each pulsar at high energies. We show that, for thermal emission models, the emitting area for X-rays from PSR B1929+10 is roughly consistent with the canonical size for a heated polar cap, and that the conversion of spin-down power to gamma-ray luminosity in B0355+54 must be low. The new proper motion for B1929+10 also implies that its progenitor is unlikely to have been the binary companion of the runaway O-star zeta-Ophiuchi.
[114]  oai:arXiv.org:astro-ph/0308299  [pdf] - 58578
PSR J0609+2130: A disrupted binary pulsar?
Comments: 5 pages, 2 figures, accepted for publication in MNRAS (letters)
Submitted: 2003-08-18, last modified: 2003-11-12
We report the discovery and initial timing observations of a 55.7-ms pulsar, J0609+2130, found during a 430-MHz drift-scan survey with the Arecibo radio telescope. With a spin-down rate of $3.1 \times 10^{-19}$ s s$^{-1}$ and an inferred surface dipole magnetic field of only $4.2 \times 10^{9}$ G, J0609+2130 has very similar spin parameters to the isolated pulsar J2235+1506 found by Camilo, Nice & Taylor (1993). While the origin of these weakly magnetized isolated neutron stars is not fully understood, one intriguing possibility is that they are the remains of high-mass X-ray binary systems which were disrupted by the supernova explosion of the secondary star.
[115]  oai:arXiv.org:astro-ph/0310454  [pdf] - 60117
New Pulsars Discovered in Arecibo Drift-Scan Searches
Comments: 2 pages, 1 figure. To appear in "Young Neutron Stars and Their Environments" (IAU Symposium 218, ASP Conference Proceedings), eds F. Camilo and B. M. Gaensler
Submitted: 2003-10-16
We report on new pulsars discovered in Arecibo drift-scan data. Processing of 2200 square degrees of data has resulted in the detection of 41 known and 12 new pulsars. New pulsars include two millisecond pulsars, one solitary and one binary recycled pulsar, and one pulsar with very unusual pulse profile morphology and complex drifting subpulse behavior.
[116]  oai:arXiv.org:astro-ph/0212180  [pdf] - 53570
Improved Bounds on Violation of the Strong Equivalence Principle
Comments: 6 pages, 2 figures, to appear in proceedings: "Radio Pulsars," Chania, Crete, August 2002 (ASP conference series)
Submitted: 2002-12-06
I describe a unique, 20-year-long timing program for the binary pulsar B0655+64, the stalwart control experiment for measurements of gravitational radiation damping in relativistic neutron-star binaries. Observed limits on evolution of the B0655+64 orbit provide new bounds on the existence of dipolar gravitational radiation, and hence on violation of the Strong Equivalence Principle.
[117]  oai:arXiv.org:astro-ph/0211261  [pdf] - 52996
New Pulsars from Arecibo Drift-Scan Searches
Comments: For Proceedings of the "Radio Pulsars" conference in Crete in August, 2002
Submitted: 2002-11-12
We report new pulsars discovered in drift-scan data taken by two collaborations (Berkeley/Cornell and STScI/NAIC) during the latter stages of the Arecibo upgrade period. The data were taken with the Penn State Pulsar Machine and are being processed on the COBRA cluster at Jodrell Bank. Processing is roughly 70% complete and has resulted in the detection of 10 new and 31 known pulsars, in addition to a number of pulsar candidates. The 10 new pulsars include one pulsar with a spin-period of 55 ms and another with a spin period of 5.8 ms. At the completion of the processing, we expect to have discovered roughly 20 new pulsars. All new pulsars are being subjected to a program of followup observations at Arecibo to determine spin and astrometric parameters.
[118]  oai:arXiv.org:astro-ph/0208281  [pdf] - 51097
Probing the Masses of the PSR J0621+1002 Binary System Through Relativistic Apsidal Motion
Comments: Accepted by ApJ, 10 pages, 5 figures
Submitted: 2002-08-14
Orbital, spin and astrometric parameters of the millisecond pulsar PSR J0621+1002 have been determined through six years of timing observations at three radio telescopes. The chief result is a measurement of the rate of periastron advance, omega_dot = 0.0116 +/- 0.0008 deg/yr. Interpreted as a general relativistic effect, this implies the sum of the pulsar mass, m_1, and the companion mass, m_2, to be M = m_1 + m_2 = 2.81 +/- 0.30 msun. The Keplerian parameters rule out certain combinations of m_1 and m_2, as does the non-detection of Shapiro delay in the pulse arrival times. These constraints, together with the assumption that the companion is a white dwarf, lead to the 68% confidence maximum likelihood values of m_1 = 1.70(+0.32 -0.29) msun and m_2 =0.97(+0.27 - 0.15) msun and to the 95% confidence maximum likelihood values of m_1 = 1.70(+0.59 -0.63) msun and m_2 = 0.97(+0.43 -0.24) msun. The other major finding is that the pulsar experiences dramatic variability in its dispersion measure (DM), with gradients as steep as 0.013 pc cm^{-3} / yr. A structure function analysis of the DM variations uncovers spatial fluctuations in the interstellar electron density that cannot be fit to a single power law, unlike the Kolmogorov turbulent spectrum that has been seen in the direction of other pulsars. Other results from the timing analysis include the first measurements of the pulsar's proper motion, mu = 3.5 +/- 0.3 mas / yr, and of its spin-down rate, dP/dt = 4.7 x 10^{-20}, which, when corrected for kinematic biases and combined with the pulse period, P = 28.8 ms, gives a characteristic age of 1.1 x 10^{10} yr and a surface magnetic field strength of 1.2 x 10^{9} G.
[119]  oai:arXiv.org:astro-ph/0204219  [pdf] - 48764
Discovery of radio pulsations from the X-ray pulsar J0205+6449 in supernova remnant 3C58 with the Green Bank Telescope
Comments: 5 pages, 2 included figures, accepted for publication by ApJ Letters
Submitted: 2002-04-12
We report the discovery with the 100 m Green Bank Telescope of 65 ms radio pulsations from the X-ray pulsar J0205+6449 at the center of supernova remnant 3C58, making this possibly the youngest radio pulsar known. From our observations at frequencies of 820 and 1375 MHz, the free electron column density to PSR J0205+6449 is found to be 140.7 +- 0.3 pc/cc. The barycentric pulsar period P and period derivative determined from a phase-coherent timing solution are consistent with the values previously measured from X-ray observations. The averaged radio profile of PSR J0205+6449 consists of one sharp pulse of width ~ 3 ms ~ 0.05 P. The pulsar is an exceedingly weak radio source, with pulse-averaged flux density in the 1400 MHz band of 0.045 mJy and a spectral index of ~ -2.1. Its radio luminosity of ~ 0.5 mJy kpc^2 at 1400 MHz is lower than that of ~ 99% of known pulsars and is the lowest among known young pulsars.
[120]  oai:arXiv.org:astro-ph/0106159  [pdf] - 42955
The Velocity Distribution of Isolated Radio Pulsars
Comments: 16 pages, 4 figures, ApJ (in press). Revisions to reflect refereed version include a new sub(sub)section on effects of the 'standard candle" assumption, and a revised Figure 2. Essential results unchanged
Submitted: 2001-06-08, last modified: 2001-11-29
(Abridged) We infer the velocity distribution of radio pulsars by modelling their birth, evolution, and detection in large-scale 0.4 GHz pulsar surveys, and by comparing model distributions of measurable pulsar properties with survey data using a likelihood function. We test models that characterize a population's birth rate, luminosity, shutoff of radio emission, birth locations, and birth velocities. We infer that the radio beam luminosity (i) is comparable to the energy flux of relativistic particles in models for spin-driven magnetospheres, signifying that radio emission losses reach nearly 100% for the oldest pulsars; and (ii) scales approximately as sqrt(Edot) which, in magnetosphere models, is proportional to the voltage drop available for acceleration of particles. We find that a two-component velocity distribution with characteristic velocities of 90 km/s and 500 km/s is greatly preferred to any one-component distribution. We explore some consequences of the preferred birth velocity distribution: (i) roughly 50% of pulsars in the solar neighborhood will escape the Galaxy, while ~15% have velocities greater than 1000 km/s; (ii) an important low-velocity population exists that increases the fraction of neutron stars retained by globular clusters and is consistent with the number of old objects that accrete from the interstellar medium; (iii) under standard assumptions for supernova remnant expansion and pulsar spin-down, ~10% of pulsars younger than 20 kyr will appear to lie outside of their host remnants. Finally, we comment on the ramifications of our birth velocity distribution for binary survival and the population of inspiraling binary neutron stars relevant to some GRB models and potential sources for LIGO.
[121]  oai:arXiv.org:astro-ph/0107145  [pdf] - 43516
A Search for Supernova-Remnant Masers Toward Unidentified EGRET Sources
Comments: 5 pages, 1 figure. To appear in Proceedings, GAMMA2001 (Baltimore, MD, April 4-6, 2001), eds. N. Gehrels, C. Shrader, and S. Ritz
Submitted: 2001-07-08
Supernova remnants expanding into adjacent molecular clouds are believed to be sites of cosmic ray acceleration and sources of energetic gamma-rays. Under certain environmental conditions, such interactions also give rise to unusual OH masers in which the 1720 MHz satellite line dominates over the more common 1665/7 MHz emission. Motivated by the apparent coincidence of a handful of EGRET sources with OH(1720 MHz) maser-producing supernova remnants, we have carried out a search using the Very Large Array for new OH(1720 MHz) masers within the error regions of 11 unidentified EGRET sources at low Galactic latitude. While a previously known maser associated with an HII region was serendipitously detected, initial results indicate that no new masers were found down to a limiting flux of, typically, 50 mJy. We discuss the implications of this result on the nature of the unidentified Galactic EGRET sources.
[122]  oai:arXiv.org:astro-ph/0106371  [pdf] - 43167
PSR J1740+1000: A Young Pulsar Well Out Of The Galactic Plane
Comments: 27 pages, 9 figures, submitted to ApJ
Submitted: 2001-06-20
We discuss PSR J1740+1000, one of five pulsars recently discovered in a search of 470 square degrees at 430 MHz during the upgrade of the 305-m Arecibo telescope. The period of 154 ms and period derivative of 2.1 x 10^-14 s/s imply a spin-down age of 114 kyr that is smaller than 95% of all known pulsars. The youth and proximity of this pulsar make it a good candidate for detection at X-ray and gamma-ray energies. Its high Galactic latitude of 20.4 degrees suggests a very high velocity if the pulsar was born in the midplane of the Galaxy and if its kinematic age equals its spindown age. Interstellar scintillations, however, suggest a much lower velocity. We discuss possible explanations for this discrepancy, taking into account (a) possible birth sites away from the midplane; (b) contributions from the unmeasured radial velocity; (c) a kinematic age different from the spin-down age; and (d) biasing of the scintillation velocity by enhanced scattering from the North Polar Spur.
[123]  oai:arXiv.org:astro-ph/0008054  [pdf] - 37370
New Pulsars from an Arecibo Drift Scan Search
Comments: 21 pages, 4 figures, Accepted for Publication in ApJ
Submitted: 2000-08-02
We report the discovery of pulsars J0030+0451, J0711+0931, and J1313+0931 that were found in a search of 470 square degrees at 430 MHz using the 305m Arecibo telescope. The search has an estimated sensitivity for long period, low dispersion measure, low zenith angle, and high Galactic latitude pulsars of ~1 mJy, comparable to previous Arecibo searches. Spin and astrometric parameters for the three pulsars are presented along with polarimetry at 430 MHz. PSR J0030+0451, a nearby pulsar with a period of 4.8 ms, belongs to the less common category of isolated millisecond pulsars. We have measured significant polarization in PSR J0030+0451 over more than 50% of the period, and use these data for a detailed discussion of its magnetospheric geometry. Scintillation observations of PSR J0030+0451 provide an estimate of the plasma turbulence level along the line of sight through the local interstellar medium.
[124]  oai:arXiv.org:astro-ph/9911198  [pdf] - 109326
Geodetic Precession in PSR B1534+12
Comments: 4 pages, 2 figures. Corrected a reference. To appear in proceedings of IAU Colloquium 177: Pulsar Astronomy - 2000 and Beyond
Submitted: 1999-11-11, last modified: 1999-11-24
We present Arecibo observations of PSR B1534+12 which confirm previous suggestions that the pulse profile is evolving secularly. This effect is similar to that seen in PSR B1913+16, and is almost certainly due to general relativistic precession of the pulsar's spin axis.
[125]  oai:arXiv.org:astro-ph/9911211  [pdf] - 109339
Binary Eclipsing Millisecond Pulsars: A Decade of Timing
Comments: Six pages, for the proceedings of Pulsar Astronomy--2000 and Beyond
Submitted: 1999-11-11
We present results of long-term timing of eclipsing binaries PSR B1744-24A and PSR B1957+20 at Arecibo, the VLA, and Green Bank. Both pulsars exhibit irregularities in pulsar rotation and orbital motion. Increases and decreases of the orbital period of PSR B1957+20 are of order Delta-P_b/P_b ~ 10^-7, varying on a time scale of a few years. Over a decade of observations, the orbital period of PSR B1744-24A has only decreased, with time scale |P_b/P_b-dot| ~ 200 Myr. When the effects of orbital motion are removed from the timing data, long-term trends remain in the pulse phase residuals, with amplitudes of order 30 and 500 us, respectively, for B1957+20 and B1744-24A. Such large "timing noise" is not seen in other spun-up pulsars (isolated or binary), leading us to conclude that it is a consequence of mass flow in the system. Possible causes include variations in the rotation of the pulsars and movement of the binary systems along the line of sight (perhaps due to gravitational interactions with outflowing matter).
[126]  oai:arXiv.org:astro-ph/9903227  [pdf] - 105631
The Triple Pulsar System PSR B1620-26 in M4
Comments: 12 pages, 6 figures, 3 tables. Very minor clarifications and rewording. Accepted for publication in the Astrophys. J
Submitted: 1999-03-15, last modified: 1999-05-12
The millisecond pulsar PSR B1620-26, in the globular cluster M4, has a white dwarf companion in a half-year orbit. Anomalously large variations in the pulsar's apparent spin-down rate have suggested the presence of a second companion in a much wider orbit. Using timing observations made on more than seven hundred days spanning eleven years, we confirm this anomalous timing behavior. We explicitly demonstrate, for the first time, that a timing model consisting of the sum of two non-interacting Keplerian orbits can account for the observed signal. Both circular and elliptical orbits are allowed, although highly eccentric orbits require improbable orbital geometries. The motion of the pulsar in the inner orbit is very nearly a Keplerian ellipse, but the tidal effects of the outer companion cause variations in the orbital elements. We have measured the change in the projected semi-major axis of the orbit, which is dominated by precession-driven changes in the orbital inclination. This measurement, along with limits on the rate of change of other orbital elements, can be used to significantly restrict the properties of the outer orbit. We find that the second companion most likely has a mass m~0.01 Msun --- it is almost certainly below the hydrogen burning limit (m<0.036 Msun, 95% confidence) --- and has a current distance from the binary of ~35 AU and orbital period of order one hundred years. Circular (and near-circular) orbits are allowed only if the pulsar magnetic field is ~3x10^9 G, an order of magnitude higher than a typical millisecond pulsar field strength. In this case, the companion has mass m~1.2x10^-3 Msun and orbital period ~62 years.
[127]  oai:arXiv.org:astro-ph/9811323  [pdf] - 103949
Pulsar Spin-up, Kinematics, and the Birthrate of Neutron Star Binaries
Comments: 9 pages (LaTeX2e with emulateapj.sty and apjfonts.sty), 3 PostScript figures. Submitted to ApJ
Submitted: 1998-11-20
From considerations of spin evolution and kinematics in the galactic potential, we argue that the pulsars B1913+16, B1534+12, and B2127+11C may be younger than previously assumed, and we find that a lower bound on the formation and merger rate of close double neutron star binaries is $10^{-6.7} f_b^{-1} f_d^{-1}$ yr$^{-1}$, where $f_b$ is the beaming fraction and $f_d$ accounts for the possibility that the known NS-NS binaries are atypical of the underlying population (e.g., if most such binaries are born with shorter orbital periods). If we assume no prior knowledge of the detectable lifetimes of such systems, the rate could be as large as $\simeq 10^{-5.0} f_b^{-1} f_d^{-1}$ yr$^{-1}$. From both plausible bounds on $f_b$ and $f_d$, and a revision of the independently derived limit proposed by Bailes (1996), we argue that a firm upper bound is $10^{-4}$ yr$^{-1}$. We also present a unifying empirical overview of the spin-up of massive binary pulsars based on their distribution in spin period $P$ and spin-down rate $\dot P$, finding evidence for two distinct spin-up processes, one dominated by disk accretion, the other by wind accretion. We argue that the positions of binary pulsars in the $P$-$\dot P$ diagram can be understood if (1) there exists a Galactic population of pulsars in double neutron star systems with combinations of spin and orbital periods that have prevented their detection in surveys to date; (2) recycled pulsars in wide-orbit binaries are not born near the canonical spin-up line in the $P$-$\dot P$ diagram because they were predominantly spun up through wind accretion; and (3) there exists a disfavored evolutionary endpoint for radio pulsars with spin periods 5--30 ms and $\dot P > 10^{-19}$ s-s$^{-1}$.
[128]  oai:arXiv.org:astro-ph/9712296  [pdf] - 99785
Measurement of Relativistic Orbital Decay in the PSR B1534+12 Binary System
Comments: 17 pages, 4 figures, submitted to the ApJ
Submitted: 1997-12-22
We have made timing observations of binary pulsar PSR B1534+12 with radio telescopes at Arecibo, Green Bank, and Jodrell Bank. By combining our new observations with data collected up to seven years earlier, we obtain a significantly improved solution for the astrometric, spin, and orbital parameters of the system. For the first time in any binary pulsar system, no fewer than five relativistic or "post-Keplerian" orbital parameters are measurable with useful accuracies in a theory-independent way. We find the orbital period of the system to be decreasing at a rate close to that expected from gravitational radiation damping, according to general relativity, although the precision of this test is limited to about 15% by the otherwise poorly known distance to the pulsar. The remaining post-Keplerian parameters are all consistent with one another and all but one of them have fractional accuracies better than 1%. By assuming that general relativity is the correct theory of gravity, at least to the accuracy demanded by this experiment, we find the masses of the pulsar and companion star each to be 1.339+-0.003 Msun and the system's distance to be d = 1.1+-0.2 kpc, marginally larger than the d ~ 0.7 kpc estimated from the dispersion measure. The increased distance reduces estimates of the projected rate of coalescence of double neutron-star systems in the universe, a quantity of considerable interest for experiments with terrestrial gravitational wave detectors such as LIGO.
[129]  oai:arXiv.org:astro-ph/9605141  [pdf] - 94697
Orbital Parameters of the PSR B1620-26 Triple System
Comments: uses PASP LaTeX macros, 6 pages, to appear in IAU Colloquium 160, Pulsars: Problems and Progress, also available from http://ensor.mit.edu/~rasio/
Submitted: 1996-05-22
Previous timing data for PSR~B1620$-$26 were consistent with a second companion mass $m_2$ anywhere in the range $\sim10^{-3}-1\,M_\odot$, i.e., from a Jupiter-type planet to a star. We present the latest timing parameters for the system, including a significant change in the projected semi-major axis of the inner binary, a marginal detection of the fourth time derivative of the pulse frequency, and the pulsar proper motion (which is in agreement with published values for the proper motion of M4), and use them to further constrain the mass $m_2$ and the orbital parameters. Using the observed value of $\stackrel{\ldots.}{f}$, we obtain a one-parameter family of solutions, all with $m_2 \lo 10^{-2}\,M_\odot$, i.e., excluding stellar masses. Varying $\stackrel{\ldots.}{f}$ within its formal $1\sigma$ error bar does not affect the mass range significantly. However, if we vary $\stackrel{\ldots.}{f}$ within a $4\sigma$ error bar, we find that stellar-mass solutions are still possible. We also calculate the predicted rate of change of the projected semi-major axis of the inner binary and show that it agrees with the measured value.
[130]  oai:arXiv.org:astro-ph/9312032  [pdf] - 91133
Orbital Variability in the Eclipsing Pulsar Binary PSR B1957+20
Comments: 9 pages, 3 figures, LaTeX, submitted ApJL 13 Dec. 1993, arz-001
Submitted: 1993-12-15
We have conducted timing observations of the eclipsing millisecond binary pulsar PSR~B1957+20, extending the span of data on this pulsar to more than five years. During this time the orbital period of the system has varied by roughly $\Delta P_b/P_b = 1.6 \times 10^{-7}$, changing quadratically with time and displaying an orbital period second derivative $\ddot P_b = (1.43 \pm 0.08) \times 10^{-18}\,$s$^{-1}$. The previous measurement of a large negative orbital period derivative reflected only the short-term behavior of the system during the early observations; the orbital period derivative is now positive and increasing rapidly. If, as we suspect, the PSR~B1957+20 system is undergoing quasi-cyclic orbital period variations similar to those found in other close binaries such as Algol and RS CVn, then the $0.025\,M{_\odot}$ companion to PSR~B1957+20 is most likely non-degenerate, convective, and magnetically active.
[131]  oai:arXiv.org:astro-ph/9303002  [pdf] - 90710
The Masses of Two Binary Neutron Star Systems
Comments: 10 pages plus 3 figures, latex, accepted ApJ Letters, sets-001
Submitted: 1993-03-01, last modified: 1993-03-03
The measurement or constraint of the masses of neutron stars and their binary companions tests theories of neutron star structure and of pulsar formation and evolution. We have measured the rate of the general relativistic advance of the longitude of periastron for the pulsar PSR B1802$-$07: $\dot\omega=0\fdg060\pm0\fdg009\,\mbox{yr}^{-1}$, which implies a total system mass, pulsar plus companion star, of $M=1.7\pm0.4\,\Msun$. We also present a much improved measurement of the rate of periastron advance for PSR B2303+46: $\dot\omega = 0\fdg0099\pm0\fdg0002\,\mbox{yr}^{-1}$, implying $M=2.53\pm0.08\,\Msun$ for this system. We discuss the available constraints on distribution of mass between the pulsars and their companions, and we compare the pulsar masses with other determinations of neutron star masses.