dorsal/arxiv
View SchemaTiming Gamma-ray Pulsars using Gibbs Sampling
| Authors | Colin J. Clark, Serena Valtolina, Lars Nieder, Rutger van Haasteren |
|---|---|
| Categories | |
| ArXiv ID | 2601.07592vv1 |
| URL | https://arxiv.org/abs/2601.07592 |
| License | http://creativecommons.org/licenses/by/4.0/ |
Abstract
Timing analyses of gamma-ray pulsars in the Fermi Large Area Telescope data set can provide sensitive probes of many astrophysical processes, including timing noise in young pulsars, orbital period variations in redback binaries, and the stochastic gravitational wave background (GWB). These goals can require careful accounting of stochastic noise processes, but existing methods developed to achieve this in radio pulsar timing analyses cannot be immediately applied to the discrete gamma-ray arrival time data. To address this, we have developed a new method for timing gamma-ray pulsars, in which the timing model fit is transformed into a weighted least squares problem by randomly assigning each photon to an individual Gaussian component of a template pulse profile. These random assignments are then numerically marginalised over through a Gibbs sampling scheme. This method allows for efficient estimation of timing and noise model parameters, while taking into account uncertainties in the pulse profile shape. We simulated Fermi-LAT data sets for gamma-ray pulsars with power-law timing noise processes, showing that this method provides robust estimates of timing noise parameters. We also describe a Gaussian-process model for orbital period variations in black-widow and redback binary systems that can be fit using this new timing method. We demonstrate this method on the black-widow binary millisecond pulsar B1957+20, where the orbital period varies significantly over the LAT data, but which provides one of the most stringent gamma-ray upper limits on the GWB.
{
"annotation_id": "7813d7b4-04e2-4f54-9892-501426738f20",
"date_created": "2026-02-17T05:53:12.608000Z",
"date_modified": "2026-02-17T05:53:12.608000Z",
"file_hash": "bad9a4e0e1af825db14214b30c8cb83c0ffb4d0abf24626b6b78441295d9fc57",
"private": false,
"record": {
"abstract": "Timing analyses of gamma-ray pulsars in the Fermi Large Area Telescope data set can provide sensitive probes of many astrophysical processes, including timing noise in young pulsars, orbital period variations in redback binaries, and the stochastic gravitational wave background (GWB). These goals can require careful accounting of stochastic noise processes, but existing methods developed to achieve this in radio pulsar timing analyses cannot be immediately applied to the discrete gamma-ray arrival time data. To address this, we have developed a new method for timing gamma-ray pulsars, in which the timing model fit is transformed into a weighted least squares problem by randomly assigning each photon to an individual Gaussian component of a template pulse profile. These random assignments are then numerically marginalised over through a Gibbs sampling scheme. This method allows for efficient estimation of timing and noise model parameters, while taking into account uncertainties in the pulse profile shape. We simulated Fermi-LAT data sets for gamma-ray pulsars with power-law timing noise processes, showing that this method provides robust estimates of timing noise parameters. We also describe a Gaussian-process model for orbital period variations in black-widow and redback binary systems that can be fit using this new timing method. We demonstrate this method on the black-widow binary millisecond pulsar B1957+20, where the orbital period varies significantly over the LAT data, but which provides one of the most stringent gamma-ray upper limits on the GWB.",
"arxiv_id": "2601.07592",
"authors": [
"Colin J. Clark",
"Serena Valtolina",
"Lars Nieder",
"Rutger van Haasteren"
],
"categories": [
"astro-ph.HE"
],
"license": "http://creativecommons.org/licenses/by/4.0/",
"title": "Timing Gamma-ray Pulsars using Gibbs Sampling",
"url": "https://arxiv.org/abs/2601.07592",
"version": "v1"
},
"schema_id": "dorsal/arxiv",
"source": {
"execution_id": "2a023388-abfd-4af6-9947-d6753ba6d16e",
"id": "arXiv Dataset",
"type": "Model",
"variant": "snapshot-2026-01-17",
"version": "0.1.0"
},
"user_id": 1000002
}