Search for: All records

Abstract In 2007 we were part of a team that discovered the so-called “Lorimer Burst”, the first example of a new class of objects now known as fast radio bursts (FRBs). These enigmatic events are only a few ms in duration and occur at random locations on the sky at a rate of a few thousand per day. Several thousand FRBs are currently known. While it is now well established that they have a cosmological origin, and about 10% of all currently known sources have been seen to exhibit multiple bursts, the origins of these enigmatic sources are currently poorly understood. In this article, we review the discovery of FRBs and present some of the highlights from the vast body of work by an international community. Following a brief overview of the scale of the visible Universe in §1, we describe the key moments in radio astronomy (§2) that led up to the discovery of the Lorimer burst (§3). Early efforts to find more FRBs are described in §4 which led to the discovery of the first repeating source (§5). In §6, as we close out on the second decade of FRBs, we outline some of the many open questions in the field and look ahead to the coming years where many surprises are surely in store. 

Abstract A pulsar’s scintillation bandwidth is inversely proportional to the scattering delay, making accurate measurements of scintillation bandwidth critical to characterize unmitigated delays in efforts to measure low-frequency gravitational waves with pulsar timing arrays. In this pilot work, we searched for a subset of known pulsars within ∼97% of the data taken with the Puerto Rico Ultimate Pulsar Processing Instrument for the AO327 survey with the Arecibo telescope, attempting to measure the scintillation bandwidths in the data set by fitting to the 2D autocorrelation function of their dynamic spectra. We successfully measured 38 bandwidths from 23 pulsars (six without prior literature values), finding that: almost all of the measurements are larger than the predictions from NE2001 and YMW16 (two popular galactic models); NE2001 is more consistent with our measurements than YMW16; Gaussian fits to the bandwidth are more consistent with both electron density models than Lorentzian ones; and for the 17 pulsars with prior literature values, the measurements between various sources often vary by a few factors. The success of Gaussian fits may be due to the use of Gaussian fits to train models in previous work. The variance of literature values over time could relate to the scaling factor used to compare measurements, but also seems consistent with time-varying interstellar medium parameters. This work can be extended to the rest of AO327 to further investigate these trends, highlighting the continuing importance of large archival data sets for projects beyond their initial conception. 

Abstract 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. 

Abstract Although neutron star–black hole binaries have been identified through mergers detected in gravitational waves, a pulsar–black hole binary has yet to be detected. While short-period binaries are detectable due to a clear signal in the pulsar’s timing residuals, effects from a long-period binary could be masked by other timing effects, allowing them to go undetected. In particular, a long-period binary measured over a small subset of its orbital period could manifest via time derivatives of the spin frequency incompatible with isolated pulsar properties. We assess the possibility of pulsars having unknown companions in long-period binaries and put constraints on the range of binary properties that may remain undetected in current data, but that may be detectable with further observations. We find that for 35% of canonical pulsars with published higher-order derivatives, the precision of measurements is not enough to confidently reject binarity (period ≳2 kyr), and that a black hole binary companion could not be ruled out for a sample of pulsars without published constraints if the period is >1 kyr. While we find no convincing cases in the literature, we put more stringent limits on orbital period and longitude of periastron for the few pulsars with published higher-order frequency derivatives (n≥ 3). We discuss the detectability of candidates and find that a sample pulsar in a 100 yr orbit could be detectable within 5–10 yr. 

Abstract 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. 

Abstract We simulate scattering delays from the interstellar medium to examine the effectiveness of three estimators in recovering these delays in pulsar timing data. Two of these estimators use the more traditional process of fitting autocorrelation functions to pulsar dynamic spectra to extract scintillation bandwidths, while the third estimator uses the newer technique of cyclic spectroscopy on baseband pulsar data to recover the interstellar medium’s impulse response function. We find that either fitting a Lorentzian or Gaussian distribution to an autocorrelation function or recovering the impulse response function from the cyclic spectrum are, on average, accurate in recovering scattering delays, although autocorrelation function estimators have a large variance, even at high signal-to-noise ratio (S/N). We find that, given sufficient S/N, cyclic spectroscopy is more accurate than both Gaussian and Lorentzian fitting for recovering scattering delays at specific epochs, suggesting that cyclic spectroscopy is a superior method for scattering estimation in high-quality data. 

Abstract Noise characterization for pulsar-timing applications accounts for interstellar dispersion by assuming a known frequency dependence of the delay it introduces in the times of arrival (TOAs). However, calculations of this delay suffer from misestimations due to other chromatic effects in the observations. The precision in modeling dispersion is dependent on the observed bandwidth. In this work, we calculate the offsets in infinite-frequency TOAs due to misestimations in the modeling of dispersion when using varying bandwidths at the Green Bank Telescope. We use a set of broadband observations of PSR J1643−1224, a pulsar with unusual chromatic timing behavior. We artificially restricted these observations to a narrowband frequency range, then used both the broad- and narrowband data sets to calculate residuals with a timing model that does not account for time variations in the dispersion. By fitting the resulting residuals to a dispersion model and comparing the fits, we quantify the error introduced in the timing parameters due to using a reduced frequency range. Moreover, by calculating the autocovariance function of the parameters, we obtained a characteristic timescale over which the dispersion misestimates are correlated. For PSR J1643−1224, which has one of the highest dispersion measures (DM) in the NANOGrav pulsar timing array, we find that the infinite-frequency TOAs suffer from a systematic offset of ∼22μs due to incomplete frequency sampling, with correlations over about one month. For lower-DM pulsars, the offset is ∼7μs. This error quantification can be used to provide more robust noise modeling in the NANOGrav data, thereby increasing the sensitivity and improving the parameter estimation in gravitational wave searches. 

Abstract We present an analysis of a densely repeating sample of bursts from the first repeating fast radio burst, FRB 121102. We reanalyzed the data used by Gourdji et al. and detected 93 additional bursts using our single-pulse search pipeline. In total, we detected 133 bursts in three hours of data at a center frequency of 1.4 GHz using the Arecibo telescope, and develop robust modeling strategies to constrain the spectro-temporal properties of all of the bursts in the sample. Most of the burst profiles show a scattering tail, and burst spectra are well modeled by a Gaussian with a median width of 230 MHz. We find a lack of emission below 1300 MHz, consistent with previous studies of FRB 121102. We also find that the peak of the log-normal distribution of wait times decreases from 207 to 75 s using our larger sample of bursts, as compared to that of Gourdji et al. Our observations do not favor either Poissonian or Weibull distributions for the burst rate distribution. We searched for periodicity in the bursts using multiple techniques, but did not detect any significant period. The cumulative burst energy distribution exhibits a broken power-law shape, with the lower- and higher-energy slopes of −0.4 ± 0.1 and −1.8 ± 0.2, with the break at (2.3 ± 0.2) × 10 37 erg. We provide our burst fitting routines as a Python package burstfit 4 4 https://github.com/thepetabyteproject/burstfit that can be used to model the spectrogram of any complex fast radio burst or pulsar pulse using robust fitting techniques. All of the other analysis scripts and results are publicly available. 5 5 https://github.com/thepetabyteproject/FRB121102 

Abstract With strong evidence of a common-spectrum stochastic process in the most recent data sets from the NANOGrav Collaboration, the European Pulsar Timing Array (PTA), Parkes PTA, and the International PTA, it is crucial to assess the effects of the several astrophysical and cosmological sources that could contribute to the stochastic gravitational wave background (GWB). Using the same data set creation and injection techniques as in Pol et al., we assess the separability of multiple GWBs by creating single and multiple GWB source data sets. We search for these injected sources using Bayesian PTA analysis techniques to assess recovery and separability of multiple astrophysical and cosmological backgrounds. For a GWB due to supermassive black hole binaries and an underlying weaker background due to primordial gravitational waves with a GW energy-density ratio of Ω_PGW/Ω_SMBHB= 0.5, the Bayes’ factor for a second process exceeds unity at 17 yr, and increases with additional data. At 20 yr of data, we are able to constrain the spectral index and amplitude of the weaker GWB at this density ratio to a fractional uncertainty of 64% and 110%, respectively, using current PTA methods and techniques. Using these methods and findings, we outline a basic protocol to search for multiple backgrounds in future PTA data sets.