An official website of the United States government
Context.High-precision pulsar timing is highly dependent on the precise and accurate modelling of any effects that can potentially impact the data. In particular, effects that contain stochastic elements contribute to some level of corruption and complexity in the analysis of pulsar-timing data. It has been shown that commonly used solar wind models do not accurately account for variability in the amplitude of the solar wind on both short and long timescales. Aims.In this study, we test and validate a new, cutting-edge solar wind modelling method included in theenterprisesoftware suite (widely used for pulsar noise analysis) through extended simulations. We use it to investigate temporal variability in LOFAR data. Our model testing scheme in itself provides an invaluable asset for pulsar timing array (PTA) experiments. Since, improperly accounting for the solar wind signature in pulsar data can induce false-positive signals, it is of fundamental importance to include in any such investigations. Methods.We employed a Bayesian approach utilising a continuously varying Gaussian process to model the solar wind. It uses a spherical approximation that modulates the electron density. This method, which we refer to as a solar wind Gaussian process (SWGP), has been integrated into existing noise analysis software, specificallyenterprise. Our Validation of this model was performed through simulations. We then conduct noise analysis on eight pulsars from the LOFAR dataset, with most pulsars having a time span of ∼11 years encompassing one full solar activity cycle. Furthermore, we derived the electron densities from the dispersion measure values obtained by the SWGP model. Results.Our analysis reveals a strong correlation between the electron density at 1 AU and the ecliptic latitude (ELAT) of the pulsar. Pulsars with |ELAT|< 3° exhibit significantly higher average electron densities. Furthermore, we observed distinct temporal patterns in the electron densities in different pulsars. In particular, pulsars within |ELAT|< 3° exhibit similar temporal variations, while the electron densities of those outside this range correlate with the solar activity cycle. Notably, some pulsars exhibit sensitivity to the solar wind up to 45° away from the Sun in LOFAR data. Conclusions.The continuous variability in electron density offered in this model represents a substantial improvement over previous models, that assume a single value for piece-wise bins of time. This advancement holds promise for solar wind modelling in future International Pulsar Timing Array (IPTA) data combinations.
more »
« less
Pulsar observations at low frequencies: applications to pulsar timing and solar wind models
ABSTRACT Efforts are underway to use high-precision timing of pulsars in order to detect low-frequency gravitational waves. A limit to this technique is the timing noise generated by dispersion in the plasma along the line of sight to the pulsar, including the solar wind. The effects due to the solar wind vary with time, influenced by the change in solar activity on different time-scales, ranging up to ∼11 yr for a solar cycle. The solar wind contribution depends strongly on the angle between the pulsar line of sight and the solar disc, and is a dominant effect at small separations. Although solar wind models to mitigate these effects do exist, they do not account for all the effects of the solar wind and its temporal changes. Since low-frequency pulsar observations are most sensitive to these dispersive delays, they are most suited to test the efficacy of these models and identify alternative approaches. Here, we investigate the efficacy of some solar wind models commonly used in pulsar timing using long-term, high-cadence data on six pulsars taken with the Long Wavelength Array, and compare them with an operational solar wind model. Our results show that stationary models of the solar wind correction are insufficient to achieve the timing noise desired by pulsar timing experiments, and we need to use non-stationary models, which are informed by other solar wind observations, to obtain accurate timing residuals.
more »
« less
- PAR ID:
- 10363164
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 511
- Issue:
- 3
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 3937-3950
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Context. Radio pulsars can be used for many studies, including the investigation of the ionized interstellar medium and the solar wind via their dispersive effects. These phenomena affect the high-precision timing of pulsars and are among the main sources of noise in experiments searching for low-frequency gravitational waves in pulsar data. Aims. In this paper, we compare the functionality and reliability of three commonly used schemes to measure temporal variations in interstellar propagation effects in pulsar timing data. Methods. We carried out extensive simulations at low observing frequencies (100–200 MHz) by injecting long-term correlated noise processes with power-law spectra and white noise, to evaluate the robustness, accuracy, and precision of the following three mitigation methods: epoch-wise (EW) measurements of interstellar dispersion; the DMX method of simultaneous, piece-wise fits to interstellar dispersion; and DM GP, which models dispersion variations through Gaussian processes using a Bayesian analysis method. We then evaluated how reliably the input signals were reconstructed and how the various methods reacted to the presence of achromatic long-period noise. Results. All the methods perform well, provided the achromatic long-period noise is modeled for DMX and DM GP. The most precise method is DM GP, followed by DMX and EW, while the most accurate is EW, followed by DMX and DM GP. We also tested different scenarios including simulations ofL-band times of arrival and realistic DM injection, with no significant variation in the obtained results. Conclusions. Given the nature of our simulations and our scope, we deem that EW is the most reliable method to study the Galactic ionized media. Follow-up works should be conducted to confirm this result via more realistic simulations. We note that DM GP and DMX seem to be the best-performing techniques in removing long-term correlated noise, and hence for gravitational wave studies. However, full simulations of pulsar timing array experiments are needed to support this interpretation.more » « less
-
ABSTRACT Using the first station of the Long Wavelength Array (LWA1), we examine polarized pulsar emission between 25 and 88 MHz. Polarized light from pulsars undergoes Faraday rotation as it passes through the magnetized interstellar medium. Observations from low-frequency telescopes are ideal for obtaining precise rotation measures (RMs) because the effect of Faraday rotation is proportional to the square of the observing wavelength. With these RMs, we obtained polarized pulse profiles to see how polarization changes in the 25–88 MHz range. The RMs were also used to derive values for the electron-density-weighted average Galactic magnetic field along the line of sight. We present RMs and polarization profiles of 15 pulsars acquired using data from LWA1. These results provide new insight into low-frequency polarization characteristics and pulsar emission heights, and complement measurements at higher frequencies.more » « less
-
Abstract Pulsar timing is a powerful tool that, by accounting for every rotation of a pulsar, precisely measures the spin frequency, spin frequency derivatives, astrometric position, binary parameters when applicable, properties of the interstellar medium, and potentially general relativistic effects. Typically, this process demands fairly stringent scheduling requirements for monitoring observations as well as deep domain knowledge to “phase connect” the timing data. We present an algorithm that automates the pulsar-timing process for binary pulsars, whose timing solutions have an additional level of complexity, although the algorithm works for isolated pulsars as well. Using the statisticalF-test and the quadratic dependence of the reducedχ2near a minimum, the global rotation count of a pulsar can be determined efficiently and systematically. We have used our algorithm to establish timing solutions for two newly discovered binary pulsars, PSRs J1748−2446aq and J1748−2446at, in the globular cluster Terzan 5, using ∼70 Green Bank Telescope observations from the last 13 yr.more » « less
-
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.more » « less
