Search for: All records

Abstract Traditional pulsar timing techniques involve averaging large numbers of single pulses to obtain a high signal-to-noise ratio (S/N) profile, which is cross-correlated with a template to measure a time of arrival (TOA). However, the morphology of individual single pulses varies greatly due to pulse jitter, and pulses of different fluence contribute differently to the S/N of the pulse average. Our study proposes a method that accounts for these variations by identifying a range of “states” and timing each separately. We selected two 1 hr observations of PSR J2145−0750, each in a different frequency band, with the Green Bank Telescope. We normalized the pulse amplitudes using the dynamic spectrum to account for scintillation effects and probed different excision algorithms to reduce radio-frequency interference. We then measured four pulse parameters (amplitude, position, width, and energy) to classify the single pulses using automated clustering algorithms. For each cluster, we calculated an average pulse profile and template and used both to obtain a TOA and TOA error. Finally, we computed the weighted average TOA and TOA error, the latter as a metric of the total timing precision for the epoch. The TOA is shifted relative to the one obtained without clustering, and we can estimate the shift with this weighting using the same data. For the 820 and 1400 MHz bands, we obtained TOA uncertainties of 0.057 and 0.46μs, compared to 0.066 and 0.74μs when no clustering is applied. We conclude that tailoring this method could help improve the timing precision for certain bright pulsars in NANOGrav’s data set. 

Abstract Changes in a pulsar’s scintillation characteristics over time can be a useful way to probe the properties of the interstellar medium. Monitoring the scintillation bandwidth can be useful for correcting for scattering delay in pulse arrival time, which could be a significant source of error for precision pulsar timing projects. To better understand the magnitude and trend of scintillation bandwidth changes over long timescales (months–years), we performed long-term monitoring of PSR J0332+5434 with the Allen Telescope Array at radio frequencies ranging from 900 to 1956 MHz. In total, we collected more than 300 observations of PSR J0332+5434 over the course of 10 months, with observations taken daily to weekly. We find a mean scintillation bandwidth of 2.6 MHz (scattering delay of 0.06μs) with a standard deviation of 1.3 MHz (0.03μs) at 1140 MHz and a mean of 14.7 MHz (0.011μs) with a standard deviation of 9.1 MHz (0.006μs) at 1668 MHz. We see significant change over time on a timescale of ∼200 days, suggesting that, for pulsars that are timed with submicrosecond precision, correcting for similar changes over time will be vital. We find no evidence for periodicity in the scattering delays. We describe a novel method for measuring the scaling index and compare it to established methods. We find the mean scintillation bandwidth frequency scaling goes asν^{3.1 ± 1.0}. 

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

Abstract Low-frequency solar radio emission is sourced in the solar corona, with sub-100 MHz radio emission largely originating from the ∼10⁵K plasma around 2 optical radii. However, the region of emission has yet to be constrained at 35–45 MHz due to both instrumentation limitations and the rarity of astronomical events, such as total solar eclipses, which allow for direct observational approaches. In this work, we present the results from a student-led project to commission a low-frequency radio telescope array situated in the path of totality of the 2024 total solar eclipse in an effort to probe the middle corona. The Deployable Low-Band Ionosphere and Transient Experiment (DLITE) is a low-frequency radio array comprised of four dipole antennas, optimized to observe at 35–45 MHz, and capable of resolving the brightest radio sources in the sky. We constructed a DLITE station in Observatory Park, a dark-sky park in Montville, Ohio. Results of observations during the total solar eclipse demonstrate that DLITE stations can be quickly deployed for observations and provide constraints on the radius of solar emission at our center observing frequency of 42 MHz. In this work, we outline the construction of DLITE Ohio and the solar observation results from the total solar eclipse that transversed North America in 2024 April. 

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 Evidence for a low-frequency gravitational-wave background using pulsar timing arrays has generated recent interest in its underlying contributing sources. However, multiple investigations have seen that the significance of the evidence does not change with the choice of pulsar modeling techniques, but the resulting parameters from the gravitational wave searches do. PSR J1455−3330 is one of the longest-observed pulsars in the array monitored by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), but showed no evidence for long-timescale red noise, either intrinsic or the common signal found among many pulsars in the array. In this work, we argue that NANOGrav’s piecewise-constant function used to model variations in radio-frequency-dependent dispersive delay should not be used for this pulsar, and a much simpler physical model of a fixed solar wind density plus a linear or quadratic trend in dispersion measure is preferred. When the original model is replaced, (i) the pulsar’s timing parallax signal changes from an upper limit to a significant detection, (ii) red noise becomes significant, and (iii) the red noise is consistent with the common signal found for the other pulsars. Neither of these signals is radio-frequency dependent. While the same physical motivation will not apply to many of the pulsars currently used in pulsar timing arrays, we argue for careful, physically motivated timing and noise modeling of pulsars used in precision timing experiments. 

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 Accurate pulsar astrometric estimates are essential to almost all high-precision pulsar timing experiments. Traditional pulsar timing techniques refine these estimates by including them as free parameters when fitting a model to observed pulse time-of-arrival measurements. However, reliable submilliarcsecond astrometric estimations require years of observations. Even then, power from red noise can be inadvertently absorbed into astrometric parameter fits. This effect biases the resulting estimates and reduces the sensitivity to red noise processes, including gravitational waves (GWs). In this work, we seek to mitigate these shortcomings by using pulsar astrometric estimates derived from very long baseline interferometry (VLBI) as priors for the timing fit. First, we used VLBI and timing astrometric estimates of 18 millisecond pulsars to calibrate a rotation between the reference frames used in timing and VLBI, with a precision of ∼0.7 mas. Through this frame tie, we combined timing- and VLBI-based probabilities to obtain a maximum-posterior astrometric solution. We found offsets between our results and the timing-based astrometric solutions, which, if real, would lead to the absorption of spectral power at the frequencies of interest for single-source GW searches. However, we do not find significant power absorption due to astrometric fitting at the low-frequency domain of the GW background. 

Abstract Pulsar timing array experiments have recently uncovered evidence for a nanohertz gravitational wave background by precisely timing an ensemble of millisecond pulsars. The next significant milestones for these experiments include characterizing the detected background with greater precision, identifying its source(s), and detecting continuous gravitational waves from individual supermassive black hole binaries. To achieve these objectives, generating accurate and precise times of arrival of pulses from pulsar observations is crucial. Incorrect polarization calibration of the observed pulsar profiles may introduce errors in the measured times of arrival. Further, previous studies have demonstrated that robust polarization calibration of pulsar profiles can reduce noise in the pulsar timing data and improve timing solutions. In this paper, we investigate and compare the impact of different polarization calibration methods on pulsar timing precision using three distinct calibration techniques: the Ideal Feed Assumption (IFA), Measurement Equation Modeling (MEM), and Measurement Equation Template Matching (METM). Three NANOGrav pulsars—PSRs J1643−1224, J1744−1134, and J1909−3744—observed with the 800 MHz and 1.5 GHz receivers at the Green Bank Telescope (GBT) are utilized for our analysis. Our findings reveal that all three calibration methods enhance timing precision compared to scenarios where no polarization calibration is performed. Additionally, among the three calibration methods, the IFA approach generally provides the best results for timing analysis of pulsars observed with the GBT receiver system. We attribute the comparatively poorer performance of the MEM and METM methods to potential instabilities in the reference noise diode coupled to the receiver and temporal variations in the profile of the reference pulsar, respectively. 

Abstract Free-floating objects (FFOs) in interstellar space—rogue planets, brown dwarfs, and large asteroids that are not gravitationally bound to any star—are expected to be ubiquitous throughout the Milky Way. Recent microlensing surveys have discovered several free-floating planets that are not bound to any known stellar systems. Additionally, three interstellar objects, namely 1I/’Oumuamua, 2I/Borisov, and 3I/ATLAS, have been detected passing through our solar system on hyperbolic trajectories. In this work, we search for FFOs on hyperbolic orbits that pass near millisecond pulsars (MSPs), where their gravitational influence can induce detectable perturbations in pulse arrival times. Using the NANOGrav 15 yr narrow band dataset, which contains high-precision timing data for 68 MSPs, we conduct a search for such hyperbolic scattering events between FFOs and pulsars. Although no statistically significant events were detected, this nondetection enables us to place upper limits (ULs) on the number density of FFOs as a function of their mass within our local region of the Galaxy. For example, the UL on the number density for Jupiter-mass FFOs (∼10^−2.5–10^−3.5M_⊙) obtained from different pulsars ranges from 5.25 × 10⁶pc⁻³to 5.37 × 10⁹pc⁻³, while the UL calculated by combining results from all the pulsars is 6.03 × 10⁵pc⁻³. These results represent the first constraints on FFO population derived from pulsar timing data.