Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRobust Estimation of Scattering in Pulsar Timing Analysis
92
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a robust approach to incorporating models for the time-variable broadening of the pulse profile due to scattering in the ionized interstellar medium into profile-domain pulsar timing analysis. We use this approach to simultaneously estimate temporal variations in both the dispersion measure (DM) and scattering, together with a model for the pulse profile that includes smooth evolution as a function of frequency, and the pulsars timing model. We show that fixing the scattering timescales when forming time-of-arrival estimates, as has been suggested in the context of traditional pulsar timing analysis, can significantly underestimate the uncertainties in both DM, and the arrival time of the pulse, leading to bias in the timing parameters. We apply our method using a new, publicly available, GPU accelerated code, both to simulations, and observations of the millisecond pulsar PSR J1643$-$1224. This pulsar is known to exhibit significant scattering variability compared to typical millisecond pulsars, and we find including low-frequency ($< 1$ GHz) data without a model for these scattering variations leads to significant periodic structure in the DM, and also biases the astrometric parameters at the $4sigma$ level, for example, changing proper motion in right ascension by $0.50 pm 0.12$. If low frequency observations are to be included when significant scattering variations are present, we conclude it is necessary to not just model those variations, but also to sample the parameters that describe the variations simultaneously with all other parameters in the model, a task for which profile domain pulsar timing is ideally suited.
rate research
Read More
We introduce a method for performing a robust Bayesian analysis of non-Gaussianity present in pulsar timing data, simultaneously with the pulsar timing model, and additional stochastic parameters such as those describing red spin noise and dispersion measure variations. The parameters used to define the presence of non-Gaussianity are zero for Gaussian processes, giving a simple method of defining the strength of non-Gaussian behaviour. We use simulations to show that assuming Gaussian statistics when the noise in the data is drawn from a non-Gaussian distribution can significantly increase the uncertainties associated with the pulsar timing model parameters. We then apply the method to the publicly available 15 year Parkes Pulsar Timing Array data release 1 dataset for the binary pulsar J0437$-$4715. In this analysis we present a significant detection of non-Gaussianity in the uncorrelated non-thermal noise, but we find that it does not yet impact the timing model or stochastic parameter estimates significantly compared to analysis performed assuming Gaussian statistics. The methods presented are, however, shown to be of immediate practical use for current European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA) datasets.
We extend profile domain pulsar timing to incorporate wide-band effects such as frequency-dependent profile evolution and broadband shape variation in the pulse profile. We also incorporate models for temporal variations in both pulse width and in the separation in phase of the main pulse and interpulse. We perform the analysis with both nested sampling and Hamiltonian Monte Carlo methods. In the latter case we introduce a new parameterisation of the posterior that is extremely efficient in the low signal-to-noise regime and can be readily applied to a wide range of scientific problems. We apply this methodology to a series of simulations, and to between seven and nine yr of observations for PSRs J1713$+$0747, J1744$-$1134, and J1909$-$3744 with frequency coverage that spans 700-3600MHz. We use a smooth model for profile evolution across the full frequency range, and compare smooth and piecewise models for the temporal variations in DM. We find the profile domain framework consistently results in improved timing precision compared to the standard analysis paradigm by as much as 40% for timing parameters. Incorporating smoothness in the DM variations into the model further improves timing precision by as much as 30%. For PSR J1713+0747 we also detect pulse shape variation uncorrelated between epochs, which we attribute to variation intrinsic to the pulsar at a level consistent with previously published analyses. Not accounting for this shape variation biases the measured arrival times at the level of $sim$30ns, the same order of magnitude as the expected shift due to gravitational-waves in the pulsar timing band.
A new Bayesian method for the analysis of folded pulsar timing data is presented that allows for the simultaneous evaluation of evolution in the pulse profile in either frequency or time, along with the timing model and additional stochastic processes such as red spin noise, or dispersion measure variations. We model the pulse profiles using `shapelets - a complete ortho-normal set of basis functions that allow us to recreate any physical profile shape. Any evolution in the profiles can then be described as either an arbitrary number of independent profiles, or using some functional form. We perform simulations to compare this approach with established methods for pulsar timing analysis, and to demonstrate model selection between different evolutionary scenarios using the Bayesian evidence. %s The simplicity of our method allows for many possible extensions, such as including models for correlated noise in the pulse profile, or broadening of the pulse profiles due to scattering. As such, while it is a marked departure from standard pulsar timing analysis methods, it has clear applications for both new and current datasets, such as those from the European Pulsar Timing Array (EPTA) and International Pulsar Timing Array (IPTA).
Pulsar timing is a technique that uses the highly stable spin periods of neutron stars to investigate a wide range of topics in physics and astrophysics. Pulsar timing arrays (PTAs) use sets of extremely well-timed pulsars as a Galaxy-scale detector with arms extending between Earth and each pulsar in the array. These challenging experiments look for correlated deviations in the pulsars timing that are caused by low-frequency gravitational waves (GWs) traversing our Galaxy. PTAs are particularly sensitive to GWs at nanohertz frequencies, which makes them complementary to other space- and ground-based detectors. In this chapter, we will describe the methodology behind pulsar timing; provide an overview of the potential uses of PTAs; and summarise where current PTA-based detection efforts stand. Most predictions expect PTAs to successfully detect a cosmological background of GWs emitted by supermassive black-hole binaries and also potentially detect continuous-wave emission from binary supermassive black holes, within the next several years.
The Vela pulsar is among a number of pulsars which show detectable optical pulsations. We performed optical observations of this pulsar in January and December 2009 with the Iqueye instrument mounted at the ESO 3.5 m New Technology Telescope. Our aim was to perform phase fitting of the Iqueye data, and to measure the optical pulse profile of the Vela pulsar at high time resolution, its absolute phase and rotational period. We calculated for the first time an independent optical timing solution and obtained the most detailed optical pulse profile available to date. Iqueye detected a distinct narrow component on the top of one of the two main optical peaks, which was not resolved in previous observations, and a third statistically significant optical peak not aligned with the radio one. The quality of the Iqueye data allowed us to determine the relative time of arrival of the radio-optical-gamma-ray peaks with an accuracy of a fraction of a millisecond. We compare the shape of the Iqueye pulse profile with that observed in other energy bands and discuss its complex multi-wavelength structure.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
