DOI: https://doi.org/10.1103/rrz1-mdkg
تاريخ النشر: 2026-02-13
المؤلف: Serena Valtolina وآخرون
الموضوع الرئيسي: أبحاث النباضات والموجات الجاذبية
نظرة عامة
تبحث الدراسة في إمكانية استخدام ملاحظات تلسكوب فيرمي الكبير للمناطق من النباضات ذات الأشعة السينية لتأسيس مجموعة توقيت النباضات (PTA) لاكتشاف إشارات موجات الجاذبية (GW) عند ترددات النانو هيرتز. التركيز الأساسي هو على الخلفية العشوائية لموجات الجاذبية (GWB) المتوقعة من تراكب الانبعاثات من ثنائيات الثقوب السوداء الضخمة. أظهرت النتائج السابقة من مجموعات توقيت النباضات الراديوية إشارة GWB مع سعة إجهاد تبلغ حوالي $2-3 \times 10^{-15}$ عند تردد $1 \text{ yr}^{-1}$. بالمقابل، أسفر تحليل 12.5 سنة من بيانات فيرمي من 35 نباضًا عن حد أعلى يبلغ $1 \times 10^{-14}$ لسعة GWB.
تستخدم الدراسة طريقة احتمالية منتظمة لإعادة تحليل نفس مجموعة البيانات، مما يعزز نمذجة الارتباطات بين النباضات ويقلل من عدم اليقين في أشكال ملفات النبض. على الرغم من أن النهج القائم على الفوتون تلو الفوتون الذي تم اختباره سابقًا لم يأخذ في الاعتبار هذه الارتباطات، تشير المحاكاة إلى أنه أكثر قوة إحصائيًا لاستعادة GWB. تؤدي التحليلات المحسنة إلى حد أعلى لسعة إجهاد GWB يبلغ $1.2 \times 10^{-14}$، مما يتماشى مع النتائج السابقة، وبالتالي يعزز صحة الطريقة الجديدة مع توفير حد أعلى متسق لاكتشاف موجات الجاذبية.
مقدمة
في هذه المقدمة، يناقش المؤلفون المنهجية وأهمية توقيت النباضات، خاصة في سياق تجارب مجموعة توقيت النباضات (PTA) التي تهدف إلى اكتشاف موجات الجاذبية (GWs). يتنبأ نموذج التوقيت بأوقات الوصول (TOAs) لانبعاثات النباضات بناءً على معلمات تتعلق بدوران النباض وحركته بالنسبة للأرض. يمكن أن تنشأ بقايا التوقيت، وهي الفجوات بين TOAs المرصودة وتلك المتوقعة من النموذج، من عوامل مختلفة، بما في ذلك عدم استقرار النباضات الداخلية والتفاعلات مع الوسط بين النجمي. من الجدير بالذكر أن هذه البقايا قد تشفر أيضًا معلومات حول خلفية موجات الجاذبية العشوائية (GWB) الناتجة عن ثنائيات الثقوب السوداء الضخمة، مع توقعات للارتباطات من دالة الارتباط هيلينغز ودوانز (HD).
يبرز المؤلفون التحديات في توقيت النباضات بسبب الضوضاء من مصادر مختلفة، خاصة في الملاحظات الراديوية، حيث يمكن أن تؤدي التشتت والتبعثر إلى إخفاء إشارة GWB. كما يذكرون إمكانية الملاحظات ذات الأشعة السينية، التي، على الرغم من تدفق الفوتونات الأقل، تقدم نموذج ضوضاء أبسط بسبب غياب التفاعلات مع الوسط بين النجمي. تختتم المقدمة بالإشارة إلى النتائج الأخيرة من تعاونات PTA التي تشير إلى وجود دليل على عملية حمراء مشتركة تتماشى مع الارتباط HD، وتحدد هدف الورقة: تقديم أول بحث عن إشارة GW مرتبطة بـ HD باستخدام بيانات الفوتون من النباضات ذات الأشعة السينية، باستخدام طريقة جديدة تعتمد على صياغة منتظمة لاحتمالية PTA في فضاء فورييه. كما يتم الإشارة بإيجاز إلى هيكل الورقة، مما يدل على نهج شامل لاستنتاج GWB من ملاحظات الفوتون.
طرق
في هذا القسم، يقدم المؤلفون نظرة عامة على الإطار البايزي المستخدم لتحليل بقايا التوقيت في تجارب مجموعة توقيت النباضات الراديوية (PTA)، إلى جانب ملخص للمنهجيات والنتائج من الإصدار الأول للبيانات من GPTA. هذا الإطار ضروري لفهم المعالجة الإحصائية لبقايا التوقيت، والتي تعتبر حاسمة لاكتشاف موجات الجاذبية (GWs).
علاوة على ذلك، يقدم المؤلفون صياغة تم تطويرها في المرجع [25] التي تسهل استنتاج إشارة مرتبطة بموجات الجاذبية بناءً على ملاحظات فوتونات الأشعة السينية من النباضات. يتم التحقق من صحة هذه الطريقة مقابل بيانات PTA الراديوية، مما يضمن قوتها وقابليتها للتطبيق في سياق اكتشاف موجات الجاذبية. للحصول على اشتقاق شامل والتحقق من صحة هذا النهج، يتم توجيه القراء إلى المرجع [25].
نتائج
في هذا القسم، يقدم المؤلفون النتائج من الإصدار الأول للبيانات من مجموعة توقيت النباضات لموجات الجاذبية (GPTA)، مع التركيز على تحليل ملاحظات الفوتون من 35 نباضًا. باستخدام التحليل ذي الخطوتين الموضح في القسم II B، استخدم المؤلفون أداة SHOOGLE لاشتقاق توزيع الخلفية لثوابت فورييه لكل نباض، تحت معلمات هايبر ثابتة $\rho_0 = [\log_{10} A_0, \gamma_0, \log_{10} \kappa_0] = [-10.2, 5, -9.2]$. تم نمذجة ضوضاء توقيت النباض (TN) باستخدام كثافة طيفية للقوة على شكل قانون قوة ذي ذيل مسطح عبر أول عشرة صناديق تردد، بينما تم نمذجة الخلفية لموجات الجاذبية (GWB) عبر أدنى خمسة صناديق تردد.
تضمن التحليل حدين مختلفين لوزن الفوتون: (i) حد ثابت يبلغ $w > 0.05$، متسق مع الدراسات السابقة، و(ii) حد مثالي لكل نباض، تم تحديده لضمان أن مجموع الأوزان المربعة للفوتونات المدرجة يمثل 99% من الإجمالي. النتائج، الموضحة في رسم بياني للاحتمالية-الاحتمالية (الشكل 2)، تقارن استعادة سعة GWB من 200 مجموعة بيانات محاكاة. تمثل الخط الأحمر الصلب الاستعادات من بيانات الفوتون باستخدام طريقة احتمالية منتظمة، بينما يعكس الخط الأحمر المتقطع الاستعادات من أوقات الوصول (TOAs) المحسوبة بعد طي بيانات الفوتون. تشير المناطق الرمادية إلى فترات الثقة 1σ و3σ لنموذج استعادة غير متحيز، مما يوفر رؤية شاملة لنتائج التحليل.
مناقشة
في هذا القسم، يناقش المؤلفون المنهجية لتحليل بقايا التوقيت من ملاحظات النباضات ذات الأشعة السينية للبحث عن الخلفيات لموجات الجاذبية (GWB). يتم اشتقاق بقايا التوقيت، التي يشار إليها بـ $\delta t$، من أوقات الوصول المرصودة (TOAs) ويتم التعبير عنها كدالة لنموذج توقيت يتم تقييمه عند معلمات مرجعية. يستخدم المؤلفون تقريبًا منخفض الرتبة لنمذجة التأخيرات الإضافية الناتجة عن عمليات عشوائية مختلفة، بما في ذلك ضوضاء النباض الداخلية وتغيرات قياس التشتت، باستخدام أساس فورييه. يتم صياغة دالة الاحتمالية لبقايا التوقيت كإحدى التوزيعات الغاوسية، ويتم اشتقاق التوزيع الخلفي لمعلمات الضوضاء والمعلمات الفائقة من خلال الاستدلال البايزي.
تقدم الورقة نهج احتمالية فورييه منتظم يفصل التحليل إلى خطوتين: أولاً، تحليل النباضات الفردية للحصول على توزيعات خلفية لثوابت فورييه، وثانيًا، دمج هذه النتائج للبحث عن GWB عبر مجموعة النباضات. هذه الطريقة مفيدة لأنها يمكن أن تطبق على كل من بيانات الأشعة السينية والراديو، مما يسمح بتحليل أكثر قوة يأخذ في الاعتبار عدم اليقين في أشكال ملفات النبض ومعلمات الضوضاء. يبرز المؤلفون التطبيق الناجح لهذه الطريقة على مجموعات البيانات المحاكاة، مما يوضح فعاليتها في استعادة إشارات GWB وتوفير حدود عليا على سعات GWB، والتي تتماشى مع النتائج السابقة من تحليلات توقيت النباضات الراديوية.
DOI: https://doi.org/10.1103/rrz1-mdkg
Publication Date: 2026-02-13
Author(s): Serena Valtolina et al.
Primary Topic: Pulsars and Gravitational Waves Research
Overview
The research investigates the potential of using Fermi Large Area Telescope observations of gamma-ray pulsars to establish a pulsar timing array (PTA) for detecting gravitational wave (GW) signals at nanohertz frequencies. The primary focus is on the stochastic gravitational wave background (GWB) anticipated from the superposition of emissions from supermassive black hole binaries. Previous findings from radio PTAs suggested a GWB signal with a strain amplitude of approximately $2-3 \times 10^{-15}$ at a frequency of $1 \text{ yr}^{-1}$. In contrast, an analysis of 12.5 years of Fermi data from 35 pulsars yielded an upper limit of $1 \times 10^{-14}$ for the GWB amplitude.
The study employs a regularized likelihood method to reanalyze the same dataset, effectively modeling cross-pulsar correlations and marginalizing over uncertainties in pulse profile shapes. Although the photon-by-photon approach previously tested did not account for these correlations, simulations indicate that it is statistically more robust for GWB recovery. The improved analysis results in an upper limit for the GWB strain amplitude of $1.2 \times 10^{-14}$, which aligns with earlier findings, thereby reinforcing the validity of the new method while providing a consistent upper limit for gravitational wave detection.
Introduction
In this introduction, the authors discuss the methodology and significance of pulsar timing, particularly in the context of pulsar timing array (PTA) experiments aimed at detecting gravitational waves (GWs). The timing model predicts the times of arrival (TOAs) of pulsar emissions based on parameters related to the pulsar’s rotation and its motion relative to Earth. Timing residuals, which are the discrepancies between observed TOAs and those predicted by the model, can arise from various factors, including intrinsic pulsar instabilities and interactions with the interstellar medium. Notably, these residuals may also encode information about a stochastic gravitational wave background (GWB) generated by supermassive black hole binaries, with correlations predicted by the Hellings and Downs (HD) correlation function.
The authors highlight the challenges in pulsar timing due to noise from various sources, particularly in radio observations, where dispersion and scattering can obscure the GWB signal. They also mention the potential of gamma-ray observations, which, despite lower photon flux, offer a simpler noise model due to the absence of interstellar medium interactions. The introduction concludes by noting recent findings from PTA collaborations that suggest evidence for a common red process consistent with the HD correlation, and outlines the paper’s objective: to present the first search for an HD-correlated GW signal using photon data from gamma-ray pulsars, employing a novel method based on a regularized formulation of the PTA likelihood in Fourier space. The structure of the paper is also briefly outlined, indicating a comprehensive approach to GWB inference from photon observations.
Methods
In this section, the authors provide an overview of the Bayesian framework utilized for analyzing timing residuals in radio Pulsar Timing Array (PTA) experiments, alongside a summary of the methodologies and results from the initial data release of the GPTA. This framework is essential for understanding the statistical treatment of timing residuals, which are critical for detecting gravitational waves (GWs).
Furthermore, the authors introduce a formalism developed in Reference [25] that facilitates the inference of a GW-correlated signal based on observations of gamma-ray photons from pulsars. This method is validated against radio PTA data, ensuring its robustness and applicability in the context of gravitational wave detection. For a comprehensive derivation and validation of this approach, readers are directed to Reference [25].
Results
In this section, the authors present the results from the first data release of the Gravitational Wave Pulsar Timing Array (GPTA), focusing on the analysis of photon observations from 35 pulsars. Utilizing the two-step analysis outlined in Section II B, the authors employed the SHOOGLE tool to derive the posterior distribution of Fourier coefficients for each pulsar, under fixed hyperparameters $\rho_0 = [\log_{10} A_0, \gamma_0, \log_{10} \kappa_0] = [-10.2, 5, -9.2]$. The pulsar timing noise (TN) was modeled using a flat-tail power law-shaped power spectral density (PSD) across the first ten frequency bins, while the gravitational wave background (GWB) was modeled over the lowest five frequency bins.
The analysis incorporated two different photon weight thresholds: (i) a fixed threshold of $w > 0.05$, consistent with previous studies, and (ii) an optimal threshold for each pulsar, determined to ensure that the sum of squared weights for the included photons accounts for 99% of the total. The results, illustrated in a probability-probability plot (Figure 2), compare the recovery of GWB amplitude from 200 simulated datasets. The solid red line represents recoveries from the photon data using a regularized likelihood method, while the dashed red line reflects recoveries from time of arrivals (TOAs) computed after folding the photon data. The gray areas indicate the 1σ and 3σ confidence intervals for an unbiased recovery model, providing a comprehensive view of the analysis outcomes.
Discussion
In this section, the authors discuss the methodology for analyzing timing residuals from gamma-ray pulsar observations to search for gravitational wave backgrounds (GWB). The timing residuals, denoted as $\delta t$, are derived from observed times of arrival (TOAs) and are expressed as a function of a timing model evaluated at reference parameters. The authors employ a low-rank approximation to model additional delays caused by various stochastic processes, including pulsar intrinsic noise and dispersion measure variations, using a Fourier basis. The likelihood function for the timing residuals is formulated as a Gaussian distribution, and the posterior distribution for noise parameters and hyperparameters is derived through Bayesian inference.
The paper introduces a regularized Fourier likelihood approach that separates the analysis into two steps: first, analyzing individual pulsars to obtain posterior distributions for Fourier coefficients, and second, combining these results for a GWB search across the pulsar array. This method is advantageous as it can be applied to both gamma-ray and radio data, allowing for a more robust analysis that accounts for uncertainties in pulse profile shapes and noise parameters. The authors highlight the successful application of this method to simulated datasets, demonstrating its effectiveness in recovering GWB signals and providing upper limits on GWB amplitudes, which are consistent with previous findings from radio pulsar timing analyses.