Search for: All records

Abstract We present timing solutions spanning nearly two decades for five redback (RB) systems found in globular clusters (GCs), created using a novel technique that effectively “isolates” the pulsar. By accurately measuring the time of passage through periastron (T₀) at points over the timing baseline, we use a piecewise-continuous binary model to get local solutions of the orbital variations that we pair with long-term orbital information to remove the orbital timing delays. The isolated pulse times of arrival can then be fit to describe the spin behavior of the millisecond pulsar (MSP). The results of our timing analyses via this method are consistent with those of conventional timing methods for binaries in GCs as demonstrated by analyses of NGC 6440D. We also investigate the observed orbital phase variations for these systems. Quasiperiodic oscillations in Terzan 5P’s orbit may be the result of changes to the gravitational quadruple moment of the companion as prescribed by the Applegate model. We find a striking correlation between the standard deviation of the phase variations as a fraction of a system’s orbit ( ${\sigma }_{\mathrm{\Delta }{T}_0}$ ) and the MSP’s spin frequency, as well as a potential correlation between ${\sigma }_{\mathrm{\Delta }{T}_0}$ and the binary’s projected semimajor axis. While long-term RB timing is fraught with large systematics, our work provides a needed alternative for studying systems with significant orbital variations, especially when high-cadence monitoring observations are unavailable. 

We present a 34 yr timing solution of the redback pulsar system Terzan 5A (Ter5A). Ter5A, also known as B1744−24A or J1748−2446A, has a 11.56 ms pulse period, a ~0.1 Msun dwarf companion star, and an orbital period of 1.82 hr. Ter5A displays highly variable eclipses and orbital perturbations. Using new timing techniques, we have determined a phase-connected timing solution for this system over 34 yr. This is the longest ever published for a redback pulsar. We find that the pulsar’s spin variability is much larger than most globular cluster pulsars. In fact, of the nine redback pulsars with published or in-preparation long-term timing solutions, Ter5A is by far the noisiest. We see no evidence of strong correlations between orbital and spin variability of the pulsar. We also find that long-term astrometric timing measurements are likely too contaminated by this variability to be usable, and therefore they require careful short-term timing to determine reasonable positions. Finally, we measure an orbital period contraction of  −2.5(3) x 10^-13, which is likely dominated by the general relativistic orbital decay of the system. The effects of the orbital variability due to the redback nature of the pulsar are not needed to explain the observed orbital period derivative, but they are constrained to less than ~30% of the observed value. 

We present the discovery of PSR J1947–1120, a new huntsman millisecond pulsar with a red giant companion star in a 10.3 day orbit. This pulsar was found via optical, X-ray, and radio follow-up of the previously unassociated gamma-ray source 4FGL J1947.6–1121. PSR J1947–1120 is the second confirmed pulsar in the huntsman class and establishes this as a bona fide subclass of millisecond pulsars. We use MESA models to show that huntsman pulsars can be naturally explained as neutron star binaries whose secondaries are currently in the “red bump” region of the red giant branch, temporarily underfilling their Roche lobes and hence halting mass transfer. Huntsman pulsars offer a new view of the formation of typical millisecond pulsars, allowing novel constraints on the efficiency of mass transfer and recycling at an intermediate stage in the process. 

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 Terzan 5 is a rich globular cluster within the galactic bulge containing 39 known millisecond pulsars, the largest known population of any globular cluster. These faint pulsars do not have sufficient signal-to-noise ratio (S/N) to measure reliable flux density or polarization information from individual observations in general. We combined over 5.2 days of archival data, at 1500 and 2000 MHz, taken with the Green Bank Telescope over the past 12 years. We created high-S/N profiles for 32 of the pulsars and determined precise rotation measures (RMs) for 28. We used the RMs, pulsar positions, and dispersion measures to map the projected parallel component of the Galactic magnetic field toward the cluster. The 〈 B ∣∣ 〉 shows a rough gradient of ∼6 nG arcsec −1 (∼160 nG pc −1 ) or, fractionally, a change of ∼20% in the R.A. direction across the cluster, implying Galactic magnetic field variability at sub-parsec scales. We also measured average flux densities S ν for the pulsars, ranging from ∼10 μ Jy to ∼2 mJy, and an average spectral index α = −1.35, where S ν ∝ ν α . This spectral index is flatter than most known pulsars, likely a selection effect due to the high frequencies used in pulsar searches to mitigate dispersion and scattering. We used flux densities from each observation to constrain the scintillation properties toward the cluster, finding strong refractive modulation on timescales of months. The inferred pulsar luminosity function is roughly power law, with slope ( d log N ) / ( d log L ) = − 1 at the high-luminosity end. At the low-luminosity end, there are incompleteness effects, implying that Terzan 5 contains many more pulsars. 

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