Search for: All records

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 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 We present the detection of 661 known pulsars observed with the Australian SKA Pathfinder (ASKAP) telescope at 888 MHz as part of the Rapid ASKAP Continuum Survey (RACS). Detections were made through astrometric coincidence and we estimate the false alarm rate of our sample to be ∼0.5%. Using archival data at 400 and 1400 MHz, we estimate the power-law spectral indices for the pulsars in our sample and find that the mean spectral index is −1.78 ± 0.6. However, we also find that a single power law is inadequate for modeling all the observed spectra. With the addition of flux densities between 150 MHz and 3 GHz from various imaging surveys, we find that up to 40% of our sample show deviations from a simple power-law model. Using StokesVmeasurements from the RACS data, we measured the circular polarization fraction for 9% of our sample and find that the mean polarization fraction is ∼10% (consistent between detections and upper limits). Using the dispersion-measure-derived distance, we estimate the pseudo-luminosity of the pulsars and do not find any strong evidence for a correlation with the pulsars’ intrinsic properties. 

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. 

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 Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 ± 0.8 mm s^–1yr^–1, which deviates from the acceleration predicted by equilibrium Milky Way (MW) models at a 4σlevel. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the MW interacting with orbiting dwarf galaxies (∼1 mm s^–1yr^–1), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in Gaia DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 60 kyr, a companion mass of 0.3M_⊙(indicative of a white dwarf or low-mass star), and a semimajor axis of 1900 au. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios.