Search for: All records

Abstract We explore possible advantages of cyclic spectroscopy for observations of pulsars in instances where full cyclic deconvolution is not possible. We compute cyclic merits and full-deconvolution regime boundaries for pulsars observed by NANOGrav and discuss which sources stand to benefit the most from using cyclic spectroscopy when observed with the Green Bank Telescope and DSA-2000 in a given frequency range. We compare data products, namely the wavefield, in both full-deconvolution and partial-deconvolution regimes to demonstrate what can be accomplished with incomplete phase retrieval. Additionally, we show how some phase retrieval can still be achieved in the partial-deconvolution regime and how this allows for additional information in scintillation-based data products, like the dynamic wavefield power, compared to what can be found in traditional dynamic spectra. An examination of dynamic wavefield phase as a function of observing frequency reveals more complete phase retrieval is achieved the closer one gets to the full-deconvolution regime, agreeing with the expectations of cyclic merit. While we demonstrate that fragmentary recovery of the secondary wavefield can be accomplished in the partial-deconvolution regime, we advocate for a synergistic approach with phase retrieval methods like theθ−θtransform, although we also provide discussion about shortcomings of this strategy. Finally, we use the combination of modest cyclic merit and lack of discernible results for PSR J1903+0327 to motivate the creation of an updated “cyclic merit 2.0,” which relies on scintillation bandwidth instead of observing bandwidth. 

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 We use cyclic spectroscopy to perform high-frequency resolution analyses of multihour baseband Arecibo observations of the millisecond pulsar PSR B1937+21. This technique allows for the examination of scintillation features in far greater detail than is otherwise possible under most pulsar timing array observing setups. We measure scintillation bandwidths and timescales in each of eight subbands across a 200 MHz observing band in each observation. Through these measurements we obtain intra-epoch estimates of the frequency scalings for scintillation bandwidth and timescale. Thanks to our high-frequency resolution and the narrow scintles of this pulsar, we resolve scintillation arcs in the secondary spectra due to the increased Nyquist limit, which would not have been resolved at the same observing frequency with a traditional filterbank spectrum using NANOGrav’s current time and frequency resolutions, and the frequency-dependent evolution of scintillation arc features within individual observations. We observe the dimming of prominent arc features at higher frequencies, possibly due to a combination of decreasing flux density and the frequency dependence of the plasma refractive index of the interstellar medium. We also find agreement with arc curvature frequency dependence predicted by Stinebring et al. in some epochs. Thanks to the frequency-resolution improvement provided by cyclic spectroscopy, these results show strong promise for future such analyses with millisecond pulsars, particularly for pulsar timing arrays, where such techniques can allow for detailed studies of the interstellar medium in highly scattered pulsars without sacrificing the timing resolution that is crucial to their gravitational-wave detection efforts. 

Abstract We report on findings from scintillation analyses using high-cadence observations of eight canonical pulsars with observing baselines ranging from 1–3 yr. We obtain scintillation bandwidth and timescale measurements for all pulsars in our survey and scintillation arc curvature measurements for four, and we detect multiple arcs for two. We find evidence of a previously undocumented scattering screen along the line of sight (LOS) to PSR J1645−0317, as well as evidence that a scattering screen along the LOS to PSR J2313+4253 may reside somewhere within the Milky Way’s Orion–Cygnus arm. We report evidence of a significant change in the scintillation pattern in PSR J2022+5154 from the previous two decades of literature, wherein both the scintillation bandwidth and timescale decreased by an order of magnitude relative to earlier observations at the same frequencies, potentially as a result of a different screen dominating the observed scattering. By augmenting the results of previous studies, we find general agreement with estimations of scattering delays from pulsar observations and predictions by the NE2001 electron density model but not for the newest data we have collected, providing some evidence of changes in the ISM along various LOSs over the timespans considered. In a similar manner, we find additional evidence of a correlation between a pulsar’s dispersion measure and the overall variability of its scattering delays over time. The plethora of interesting science obtained through these observations demonstrates the capabilities of the Green Bank Observatory’s 20 m telescope to contribute to pulsar-based studies of the interstellar medium. 

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 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 Pulsar timing arrays (PTAs) have recently entered the detection era, quickly moving beyond the goal of simply improving sensitivity at the lowest frequencies for the sake of observing the stochastic gravitational wave background (GWB), and focusing on its accurate spectral characterization. While all PTA collaborations around the world use Fourier-domain Gaussian processes to model the GWB and intrinsic long time-correlated (red) noise, techniques to model the time-correlated radio-frequency-dependent (chromatic) processes have varied from collaboration to collaboration. Here we test a new class of models for PTA data, Gaussian processes based on time-domain kernels that model the statistics of the chromatic processes starting from the covariance matrix. As we will show, these models can be effectively equivalent to Fourier-domain models in mitigating chromatic noise. This work presents a method for Bayesian model selection across the various choices of kernel as well as deterministic chromatic models for nonstationary chromatic events and the solar wind. As PTAs turn toward high frequency (>1 yr⁻¹) sensitivity, the size of the basis used to model these processes will need to increase, and these time-domain models present some computational efficiencies compared to Fourier-domain models. 

Abstract The observed gravitational-wave background (GWB) spectrum is higher in amplitude than model predictions by a factor of 2–3. Using a semi-analytic model, we evaluate the effect of a high-scatter supermassive black hole (SMBH) scaling relation (M_BH–M_bulge) on models of the nanohertz GWB. By implementing an intrinsic scatter of theM_BH–M_bulgerelation, which is larger at higher redshift, but matches local observations, we find that the amplitude of GWB models increases to be consistent with the low-frequency end of the GWB spectrum. This amplitude increase is not uniform across frequencies, a strongly evolving scatter preferentially increases the number density of the most massive SMBHs which, in the GWB spectrum, minimizes the strength of the low-frequency turnover. Our models with positively evolving intrinsic scatter can reproduce the electromagnetically observed overmassive SMBHs at 4 < z < 6 without changing theM_BH–M_bulgenormalization though we find that including moderate normalization evolution marginally improves fits to the GWB data. We conclude that theM_BH–M_bulgerelation which best describes the available GWB and electromagnetic data sets has intrinsic scatter that evolves as $\epsilon \left(z\right)={\epsilon }_0+(0.56\pm 0.4){\log }_{10}(1+z)$ and normalization that evolves asα(z) = α₀(1 + z)^0.84±0.35. The results of this work imply that theM_BH–M_bulgerelation we see today is not universal throughout cosmic time and that a diversity of seeding models and growth mechanisms may be at play in the early stages of SMBH–galaxy evolution. 

Abstract We test the impact of an evolving supermassive black hole mass scaling relation (M_BH–M_bulge) on the predictions for the gravitational-wave background (GWB). The observed GWB amplitude is 2–3 times higher than predicted by astrophysically informed models, which suggests the need to revise the assumptions in those models. We compare a semi-analytic model’s ability to reproduce the observed GWB spectrum with a static versus evolving-amplitudeM_BH–M_bulgerelation. We additionally consider the influence of the choice of galaxy stellar mass function (GSMF) on the modeled GWB spectra. Our models are able to reproduce the GWB amplitude with either a large number density of massive galaxies or a positively evolvingM_BH–M_bulgeamplitude (i.e., theM_BH/M_bulgeratio was higher in the past). If we assume that theM_BH–M_bulgeamplitude does not evolve, our models require a GSMF that implies an undetected population of massive galaxies (M_⋆≥ 10¹¹M_⊙atz> 1). When theM_BH–M_bulgeamplitude is allowed to evolve, we can model the GWB spectrum with all fiducial values and anM_BH–M_bulgeamplitude that evolves asα(z) =α₀(1 +z)^1.04±0.5.