Search for: All records

Broadband radio waves emitted from pulsars are distorted and delayed as they propagate toward the Earth due to interactions with the free electrons that compose the interstellar medium (ISM), with lower radio frequencies being more impacted than higher frequencies. Multipath propagation in the ISM results in both later times of arrival for the lower frequencies and causes the observed pulse to arrive with a broadened tail described via the pulse broadening function. We employ the CLEAN deconvolution technique to recover the pulse broadening timescale and by proxy the intrinsic pulse shape. This work expands upon previous descriptions of CLEAN deconvolution used in pulse broadening analyses by parameterizing the efficacy on simulated data and developing a suite of tests to establish which of a set of figures of merit leads to an automatic and consistent determination of the scattering timescale and its uncertainty. We compare our algorithm to the cyclic spectroscopy method of estimating the scattering timescale, specifically to the simulations performed in Dolch et al. (2021). We test our improved algorithm on the highly scattered millisecond pulsar J1903+0327, showing the scattering timescale to change over years, consistent with estimates of the refractive timescale of the pulsar.

The Galactic electron density model NE2001 describes the multicomponent ionized structure of the Milky Way interstellar medium. NE2001 forward models the dispersion and scattering of compact radio sources, including pulsars, fast radio bursts, active galactic nuclei, and masers, and the model is routinely used to predict the distances of radio sources lacking independent distance measures. Here we present the open-source package NE2001p, a fully Python implementation of NE2001. The model parameters are identical to NE2001 but the computational architecture is optimized for Python, yielding small (<1%) numerical differences between NE2001p and the Fortran code. NE2001p can be used on the command-line and through Python scripts available on PyPI. Future package releases will include modular extensions aimed at providing short-term improvements to model accuracy, including a modified thick disk scale height and additional clumps and voids. This implementation of NE2001 is a springboard to a next-generation Galactic electron density model now in development.

Four supernova remnants and four anomalous X-ray pulsars were previously observed with the Parkes telescope in a campaign to detect pulsed radio emission from associated neutron stars. No signals were detected in the original searches of these data. I have reprocessed the data with the more recently developed HEIMDALL and FETCH software packages, which are optimized for single-pulse detection and classification. In this new analysis, no astrophysical pulses were detected having a signal-to-noise ratio above 7 from any of the targets at dispersion measures ranging from 0 to 10⁴pc cm⁻³. I include calculated fluence limits on single radio pulses from these targets.

We use the upgraded Giant Metrewave Radio Telescope (uGMRT) to measure scintillation arc properties in six bright canonical pulsars with simultaneous dual-frequency coverage. These observations, at frequencies from 300 to 750 MHz, allowed for detailed analysis of arc evolution across frequency and epoch. We perform more robust determinations of frequency dependence for arc curvature, scintillation bandwidth, and scintillation timescale, and comparison between arc curvature and pseudo-curvature than allowed by single-frequency-band-per-epoch measurements, which we find to agree with theory and previous literature. We find a strong correlation between arc asymmetry and arc curvature, which we have replicated using simulations, and attribute to a bias in the Hough transform approach to scintillation arc analysis. Possible evidence for an approximately week-long timescale over which a given scattering screen dominates signal propagation was found by tracking visible scintillation arcs in each epoch in PSR J1136+1551. The inclusion of a 155-minute observation allowed us to resolve the scale of scintillation variations on short timescales, which we find to be directly tied to the amount of interstellar medium sampled over the observation. Some of our pulsars showed either consistent or emerging asymmetries in arc curvature, indicating instances of refraction across their lines of sight. Significant features in various pulsars, such as multiple scintillation arcs in PSR J1136+1551 and flat arclets in PSR J1509+5531, that have been found in previous works, were also detected. The simultaneous multiple-band observing capability of the upgraded GMRT shows excellent promise for future pulsar scintillation work.

Observations of pulsar scintillation are among the few astrophysical probes of very small-scale (≲ au) phenomena in the interstellar medium (ISM). In particular, characterization of scintillation arcs, including their curvature and intensity distributions, can be related to interstellar turbulence and potentially overpressurized plasma in local ISM inhomogeneities, such as supernova remnants, H ii regions, and bow shocks. Here we present a survey of eight pulsars conducted at the Five-hundred-metre Aperture Spherical Telescope (FAST), revealing a diverse range of scintillation arc characteristics at high sensitivity. These observations reveal more arcs than measured previously for our sample. At least nine arcs are observed toward B1929+10 at screen distances spanning $\sim 90~{{\ \rm per\ cent}}$ of the pulsar’s 361 pc path length to the observer. Four arcs are observed toward B0355+54, with one arc yielding a screen distance as close as ∼105 au (<1 pc) from either the pulsar or the observer. Several pulsars show highly truncated, low-curvature arcs that may be attributable to scattering near the pulsar. The scattering screen constraints are synthesized with continuum maps of the local ISM and other well-characterized pulsar scintillation arcs, yielding a three-dimensional view of the scattering media in context.

We present the discovery and timing solutions of four millisecond pulsars (MSPs) discovered in the Arecibo 327 MHz Drift-Scan Pulsar Survey. Three of these pulsars are in binary systems, consisting of a redback (PSR J2055+1545), a black widow (PSR J1630+3550), and a neutron star–white dwarf binary (PSR J2116+1345). The fourth MSP, PSR J2212+2450, is isolated. We present the multiyear timing solutions as well as polarization properties across a range of radio frequencies for each pulsar. We perform a multiwavelength search for emission from these systems and find an optical counterpart for PSR J2055+1545 in Gaia DR3, as well as a gamma-ray counterpart for PSR J2116+1345 with the Fermi-LAT telescope. Despite the close colocation of PSR J2055+1545 with a Fermi source, we are unable to detect gamma-ray pulsations, likely due to the large orbital variability of the system. This work presents the first two binaries found by this survey with orbital periods shorter than a day; we expect to find more in the 40% of the survey data that have yet to be searched.

Recently we found compelling evidence for a gravitational-wave background with Hellings and Downs (HD) correlations in our 15 yr data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes, which produce different interpulsar correlations. In this work, we search the NANOGrav 15 yr data set for evidence of a gravitational-wave background with quadrupolar HD and scalar-transverse (ST) correlations. We find that HD correlations are the best fit to the data and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors ∼2 when comparing HD to ST correlations, and ∼1 for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis.

Pulsar distances are notoriously difficult to measure, and play an important role in many fundamental physics experiments, such as pulsar timing arrays. Here, we perform a cross-match between International PTA pulsars (IPTA) and Gaia's Data Release 2 (DR2) and Data Release 3 (DR3). We then combine the IPTA pulsar’s parallax with its binary companion’s parallax, found in Gaia, to improve the distance measurement to the binary. We find seven cross-matched IPTA pulsars in Gaia DR2, and when using Gaia DR3 we find six IPTA pulsar cross-matches but with seven Gaia objects. Moving from Gaia DR2 to Gaia DR3, we find that the Gaia parallaxes for the successfully cross-matched pulsars improved by 53%, and pulsar distances improved by 29%. Finally, we find that binary companions with a <3.0σdetection are unreliable associations, setting a high bar for successful cross-matches.

Dual and off-nucleus active supermassive black holes are expected to be common in the hierarchical structure formation paradigm, but their identification at parsec scales remains a challenge due to strict angular resolution requirements. We conducted a systematic study using the Very Long Baseline Array (VLBA) to examine 23 radio-bright candidate dual and off-nucleus quasars. The targets are selected by a novel astrometric technique (varstrometry) from Gaia, aiming to identify dual or off-nucleus quasars at (sub)kiloparsec scales. Among these quasars, eight exhibit either multiple radio components or significant (>3σ) positional offsets between the VLBA and Gaia positions. The radio emission from the three candidates, which exhibit multiple radio components, is likely to originate from small-scale jets based on their morphology. Among the remaining five candidates with significant VLBA-Gaia offsets, three are identified as potential dual quasars at parsec scales, one is likely attributed to small-scale jets, and the origin of the last candidate remains unclear. We explore alternative explanations for the observed VLBA-Gaia offsets. We find no evidence for optical jets at kiloparsec scales, nor any contamination to Gaia astrometric noise from the host galaxy; misaligned coordinate systems are unlikely to account for our offsets. Our study highlights the promise of the varstrometry technique in discovering candidate dual or off-nucleus quasars and emphasizes the need for further confirmation and investigation to validate and understand these intriguing candidates.

The standard Bayesian technique for searching pulsar timing data for gravitational-wave bursts with memory (BWMs) using Markov Chain Monte Carlo (MCMC) sampling is very computationally expensive to perform. In this paper, we explain the implementation of an efficient Bayesian technique for searching for BWMs. This technique makes use of the fact that the signal model for Earth-term BWMs (BWMs passing over the Earth) is fully factorizable. We estimate that this implementation reduces the computational complexity by a factor of 100. We also demonstrate that this technique gives upper limits consistent with published results using the standard Bayesian technique, and may be used to perform all of the same analyses of BWMs that standard MCMC techniques can perform.