Engineering Papers⌕ Search

Engineering topics

Stairs, I. H.

Publications and source records attributed to Stairs, I. H..

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

Observations indicate that nearly all galaxies contain supermassive black holes 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. We have searched the North American Nanohertz Observatory for Gravitational Waves 11 yr 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. At f(gw) = 8 nHz, we placed a sky-averaged upper limit of h(0) < 7.3(3) × 10(exp −15). 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. From these upper limits, we ruled out SMBHBs emitting GWs f(gw) = 8 nHz within 120 Mpc for M = 10(exp 9) Solar Mass, and within 5.5 Gpc for M= 10(exp 10) Solar Mass at our most sensitive sky location. We also determined that there are no SMBHBs with M > 1.6 x 10(exp 9) Solar Mass emitting GWs with f(gw) = 2.8–317.8 nHz in the Virgo Cluster. Finally, we compared our strain upper limits to simulated populations of SMBHBs, based on galaxies in the Two Micron All-Sky Survey and merger rates from the Illustris cosmological simulation project, and found that only 34 out of 75,000 realizations of the local universe contained a detectable source.

Aggarwal, K.↗

The NANOGrav 12.5 yr Data Set: The Frequency Dependence of Pulse Jitter in Precision Millisecond Pulsars

Low-frequency gravitational-wave experiments require the highest timing precision from an array of the moststable millisecond pulsars. Several known sources of noise on short timescales in single radio pulsar observationsare well described by a simple model of three components: template fitting from a finite signal-to-noise ratio, pulsephase/amplitude jitter from single-pulse stochasticity, and scintillation errors from short-timescale interstellarscattering variations. Currently template-fitting errors dominate, but as radio telescopes push toward higher signalto-noise ratios, jitter becomes the next dominant term for most millisecond pulsars. Understanding the statistics ofjitter becomes crucial for properly characterizing arrival time uncertainties. We characterize the radio frequencydependence of jitter using data on 48 pulsars in the North American Nanohertz Observatory for GravitationalWaves timing program. We detect significant jitter in 43 of the pulsars and test several functional forms for itsfrequency dependence; we find significant frequency dependence for 30 pulsars. We find moderate correlations ofrms jitter with pulse width (R = 0.62) and number of profile components (R = 0.40); the single-pulse rms jitter istypically 1% of pulse phase. The average frequency dependence for all pulsars using a power-law model hasindex 0.42. We investigate the jitter variations for the interpulse of PSR B1937+21 and find no significantdeviations from the main pulse rms jitter. We also test the time variation of jitter in two pulsars and find thatsystematics likely bias the results for high-precision pulsars. Pulsar timing array analyses must properly modeljitter as a significant component of the noise within the detector.

Lam, M. T.↗

The NANOGrav 11 yr Data Set: Solar Wind Sounding through Pulsar Timing

The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has observed dozens ofmillisecond 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 ofsight cross through the solar wind (SW), which produces correlated DM fluctuations in all pulsars. We develop andapply techniques for extracting the imprint of the SW from the full collection of DM measurements in the recentlyreleased NANOGrav 11 yr data set. We filter out long-timescale DM fluctuations attributable to structure in theinterstellar medium and carry out a simultaneous analysis of all pulsars in our sample that can differentiate thecorrelated signature of the wind from signals unique to individual lines of sight. When treating the SW asspherically symmetric and constant in time, we find the electron number density at 1au to be 7.9 ± 0.2 cm(exp -3). Wefind our data to be insensitive to long-term variation in the density of the wind. We argue that our techniques pairedwith a high-cadence, low-radio-frequency observing campaign of near-ecliptic pulsars would be capable ofmapping out large-scale latitudinal structure in the wind.

Madison, D.R.↗

The NANOGrav 11 Year Data Set: Pulsar-Timing Constraints on the Stochastic Gravitational-Wave Background

We search for an isotropic stochastic gravitational-wave background (GWB) in the newly released 11 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). While we find no evidence for a GWB, we place constraints on a population of inspiraling supermassive black hole (SMBH) binaries, a network of decaying cosmic strings, and a primordial GWB. For the first time, we find that the GWB constraints 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 pulsar-timing array (PTA) constraints that are robust against SSE errors. We thus place a 95% upper limit on the GW-strain amplitude of A (sub GWB) < 1.45 × 10 (exp -15) at a frequency of f=1 yr(exp -1) for a fiducial f (exp -2/3) power-law spectrum and with interpulsar correlations modeled. This is a factor of approximately 2 improvement over the NANOGrav nine-year limit calculated using the same procedure. Previous PTA upper limits on the GWB (as well as their astrophysical and cosmological interpretations) will need revision in light of SSE systematic errors. 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 the cosmic-string tension using recent simulations, yielding an SSE-marginalized 95% upper limit of G (sub mu) < 5.3 × 10(exp -11) - a factor of approximately 2 better than the published NANOGrav nine-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 (sub GWB)(f) h (exp 2) < 3.4 × 10 (exp -10).

Supermassive Black Holes↗

A Search for Fast Radio Bursts with the GBNCC Pulsar Survey

We report on a search for fast radio bursts (FRBs) with the Green Bank Northern Celestial Cap (GBNCC) Pulsar Survey at 350 MHz. Pointings amounting to a total on-sky time of 61 days were searched to a dispersion measure (DM) of 3000 pc cm{sup −3}, while the rest (23 days; 29% of the total time) were searched to a DM of 500 pc cm{sup −3}. No FRBs were detected in the pointings observed through 2016 May. We estimate a 95% confidence upper limit on the FRB rate of 3.6×10{sup 3} FRBs sky{sup −1} day{sup −1} above a peak flux density of 0.63 Jy at 350 MHz for an intrinsic pulse width of 5 ms. We place constraints on the spectral index α by running simulations for different astrophysical scenarios and cumulative flux density distributions. The nondetection with GBNCC is consistent with the 1.4 GHz rate reported for the Parkes surveys for α > +0.35 in the absence of scattering and free–free absorption and α > −0.3 in the presence of scattering, for a Euclidean flux distribution. The constraints imply that FRBs exhibit either a flat spectrum or a spectral turnover at frequencies above 400 MHz. These constraints also allow estimation of the number of bursts that can be detected with current and upcoming surveys. We predict that CHIME may detect anywhere from several to ∼50 FRBs per day (depending on model assumptions), making it well suited for interesting constraints on spectral index, the log N–log S slope, and pulse profile evolution across its bandwidth (400–800 MHz).

79 ASTRONOMY AND ASTROPHYSICS↗

The Nanograv Nine-Year Data Set: Measurement and Analysis of Variations in Dispersion Measures

We analyze dispersion measure(DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause.

pulsars: general↗

The NANOGrav Nine-Year Data Set: Limits on the Isotropic Stochastic Gravitational Wave Background

We compute upper limits on the nanohertz-frequency isotropic stochastic gravitational wave background (GWB) using the 9 year data set from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) collaboration. Well-tested Bayesian techniques are used to set upper limits on the dimensionless strain amplitude (at a frequency of 1 yr(exp -1) for a GWB from supermassive black hole binaries of A(sub gw) less than 1.5 x 10(exp -15). We also parameterize the GWB spectrum with a broken power-law model by placing priors on the strain amplitude derived from simulations of Sesana and McWilliams et al. Using Bayesian model selection we find that the data favor a broken power law to a pure power law with odds ratios of 2.2 and 22 to one for the Sesana and McWilliams prior models, respectively. 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. Returning to a power-law model, we place stringent limits on the energy density of relic GWs, OMEGA(sub gw) (f) h squared less than 4.2 x 10(exp -10). Our limit on the cosmic string GWB, OMEGA(sub gw) (f) h squared less than 2.2 x 10(exp -10), translates to a conservative limit on the cosmic string tension with G mu less than 3.3 x 10(exp -8), a factor of four better than the joint Planck and high-l‚ cosmic microwave background data from other experiments.

Arzoumanian, Z.↗

Testing theories of gravitation using 21-year timing of pulsar binary J1713+0747

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

79 ASTRONOMY AND ASTROPHYSICS↗

TIMING OF FIVE MILLISECOND PULSARS DISCOVERED IN THE PALFA SURVEY

We present the discovery of five millisecond pulsars (MSPs) from the PALFA Galactic plane survey using Arecibo. Four of these (PSRs J0557+1551, J1850+0244, J1902+0300, and J1943+2210) are binary pulsars whose companions are likely white dwarfs, and one (PSR J1905+0453) is isolated. Phase-coherent timing solutions, ranging from ∼1 to ∼3 yr in length, and based on observations from the Jodrell Bank and Arecibo telescopes, provide precise determinations of spin, orbital, and astrometric parameters. All five pulsars have large dispersion measures (>100 pc cm{sup –3}, within the top 20% of all known Galactic field MSPs) and are faint (1.4 GHz flux density ≲0.1 mJy, within the faintest 5% of all known Galactic field MSPs), illustrating PALFA's ability to find increasingly faint, distant MSPs in the Galactic plane. In particular, PSR J1850+0244 has a dispersion measure of 540 pc cm{sup –3}, the highest of all known MSPs. Such distant, faint MSPs are important input for accurately modeling the total Galactic MSP population.

79 ASTRONOMY AND ASTROPHYSICS↗

THE BINARY COMPANION OF YOUNG, RELATIVISTIC PULSAR J1906+0746

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

79 ASTRONOMY AND ASTROPHYSICS↗

Realistic modeling of the pulse profile of PSR J0737-3039A

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

79 ASTRONOMY AND ASTROPHYSICS↗

X-ray observations of black widow pulsars

We describe the first X-ray observations of five short orbital period (P{sub B} < 1 day), γ-ray emitting, binary millisecond pulsars (MSPs). Four of these—PSRs J0023+0923, J1124–3653, J1810+1744, and J2256–1024—are 'black-widow' pulsars, with degenerate companions of mass <<0.1 M {sub ☉}, three of which exhibit radio eclipses. The fifth source, PSR J2215+5135, is an eclipsing 'redback' with a near Roche-lobe filling ∼0.2 solar mass non-degenerate companion. Data were taken using the Chandra X-Ray Observatory and covered a full binary orbit for each pulsar. Two pulsars, PSRs J2215+5135 and J2256–1024, show significant orbital variability while PSR J1124–3653 shows marginal orbital variability. The lightcurves for these three pulsars have X-ray flux minima coinciding with the phases of the radio eclipses. This phenomenon is consistent with an intrabinary shock emission interpretation for the X-rays. The other two pulsars, PSRs J0023+0923 and J1810+1744, are fainter and do not demonstrate variability at a level we can detect in these data. All five spectra are fit with three separate models: a power-law model, a blackbody model, and a combined model with both power-law and blackbody components. The preferred spectral fits yield power-law indices that range from 1.3 to 3.2 and blackbody temperatures in the hundreds of eV. The spectrum for PSR J2215+5135 shows a significant hard X-ray component, with a large number of counts above 2 keV, which is additional evidence for the presence of intrabinary shock emission. This is similar to what has been detected in the low-mass X-ray binary to MSP transition object PSR J1023+0038.

79 ASTRONOMY AND ASTROPHYSICS↗

Searching for pulsars using image pattern recognition

In the modern era of big data, many fields of astronomy are generating huge volumes of data, the analysis of which can sometimes be the limiting factor in research. Fortunately, computer scientists have developed powerful data-mining techniques that can be applied to various fields. In this paper, we present a novel artificial intelligence (AI) program that identifies pulsars from recent surveys by using image pattern recognition with deep neural nets—the PICS (Pulsar Image-based Classification System) AI. The AI mimics human experts and distinguishes pulsars from noise and interference by looking for patterns from candidate plots. Different from other pulsar selection programs that search for expected patterns, the PICS AI is taught the salient features of different pulsars from a set of human-labeled candidates through machine learning. The training candidates are collected from the Pulsar Arecibo L-band Feed Array (PALFA) survey. The information from each pulsar candidate is synthesized in four diagnostic plots, which consist of image data with up to thousands of pixels. The AI takes these data from each candidate as its input and uses thousands of such candidates to train its ∼9000 neurons. The deep neural networks in this AI system grant it superior ability to recognize various types of pulsars as well as their harmonic signals. The trained AI's performance has been validated with a large set of candidates from a different pulsar survey, the Green Bank North Celestial Cap survey. In this completely independent test, the PICS ranked 264 out of 277 pulsar-related candidates, including all 56 previously known pulsars and 208 of their harmonics, in the top 961 (1%) of 90,008 test candidates, missing only 13 harmonics. The first non-pulsar candidate appears at rank 187, following 45 pulsars and 141 harmonics. In other words, 100% of the pulsars were ranked in the top 1% of all candidates, while 80% were ranked higher than any noise or interference. The performance of this system can be improved over time as more training data are accumulated. This AI system has been integrated into the PALFA survey pipeline and has discovered six new pulsars to date.

79 ASTRONOMY AND ASTROPHYSICS↗

Fermi LAT Pulsed Detection of PSR J0737-3039A in the Double Pulsar System

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

PULSED DETECTION↗

Measurement of Gravitational Spin-Orbit Coupling in a Binary Pulsar System

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.

Stairs, I. H.↗

Discovery of Radio Pulsations from the X-ray Pulsar JO205+6449 in Supernova Remnant 3C58 with the Green Bank Telescope

We report the discovery with the 100m 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 USSR J0205+6449 is found to be 140.7 +/- 0.3/cc pc. The barycentric pulsar period P and P(dot) determined from a phase-coherent timing solution are consistent with the values previously measured from X-ray observations. The averaged radio profile of USSR 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 approximately 45 micro-Jy and a spectral index of approximately -2.1. Its radio luminosity of approximately 0.5 may kpc(exp 2) at 1400 MHz is lower than that of approximately 99% of known pulsar and is the lowest among known young pulsars.

Camilo, F.↗