DOI: https://doi.org/10.3847/1538-4357/ae3a96
تاريخ النشر: 2026-02-17
المؤلف: Run-Chao Chen وآخرون
الموضوع الرئيسي: انفجارات أشعة غاما والسوبرنوفا
نظرة عامة
في الدراسات الأخيرة، تم ملاحظة تذبذبات متماسكة عند ترددات الكيلوهرتز في بعض انفجارات أشعة غاما (GRBs)، مما يشير إلى ديناميكيات شبه دورية محتملة داخل محركاتها المركزية. ومن الأمثلة البارزة GRB 230307A، الذي أظهر إشارة دورية متماسكة تعتمد على الطاقة تم تفسيرها على أنها توقيع دوران لمغناطيسار جديد الشكل ناتج عن اندماج ثنائي مضغوط. في ضوء هذه النتائج، قمنا بإجراء بحث شامل عن إشارات مشابهة عبر 532 انفجار أشعة غاما قصير باستخدام بيانات أحداث مؤرخة زمنياً من تجربة المصدر المتقطع والانفجاري (BATSE) على مرصد كومبتون لأشعة غاما.
كشفت تحليلاتنا عن حالة واحدة ذات دلالة إحصائية: GRB 960616 (إشارة BATSE 5502)، التي تتميز بحلقة انبعاث رئيسية تستمر حوالي 30 مللي ثانية، والتي أظهرت تعديلًا متماسكًا عند 1100 هرتز، وكان ذلك ملحوظًا بشكل خاص فوق 320 كيلو إلكترون فولت، مع سعة كسرية تبلغ حوالي 47%. من خلال محاكاة مونت كارلو المعتمدة على البيانات، وضعنا حدًا أعلى يبلغ حوالي 8% على نسبة الطاقة الكلية المنبعثة التي يمكن أن تُعزى إلى هذا التذبذب شبه الدوري (QPO). تعزز هذه النتيجة، جنبًا إلى جنب مع التكرارية الملحوظة التي تبلغ حوالي 0.91 مللي ثانية، الفرضية القائلة بأن التكرارات المللي ثانية قد تظهر خلال انفجارات أشعة غاما في سيناريوهات تتضمن الاندماجات.
مقدمة
تناقش مقدمة الورقة الأنظمة السلفية لانفجارات أشعة غاما (GRBs)، والتي تصنف إلى نوعين رئيسيين: النجوم الضخمة وأنظمة ثنائية مضغوطة. تخضع هذه الأنظمة لأحداث كارثية تؤدي إلى تشكيل محرك مركزي، عادةً ما يكون ثقبًا أسودًا عالي الاكتساب أو نجم نيوتروني يدور بسرعة (مغناطيسار مللي ثانية)، والذي يغذي نفاثة نسبية مسؤولة عن انبعاثات أشعة غاما. بينما تم إثبات العلاقة بين انفجارات أشعة غاما وأنظمتها السلفية من خلال ملاحظات متنوعة، لا تزال الطبيعة الدقيقة للمحركات المركزية غير مؤكدة.
قدمت الدراسات الأخيرة رؤى جديدة حول المحركات المركزية لانفجارات أشعة غاما من خلال الكشف عن تباين عالي التردد في الانبعاثات الفورية. ومن الجدير بالذكر، تم الإبلاغ عن تذبذب شبه دوري (QPO) عند حوالي 2.7 كيلو هرتز في انفجار أشعة غاما قصير، مما يشير إلى وجود نجم نيوتروني عالي الكتلة. علاوة على ذلك، تم تحديد تعديل دوري مستقر عند 909 هرتز في GRB 230307A، مما يربطه بمغناطيسار مللي ثانية ناشئ. مدفوعين بهذه النتائج، يهدف المؤلفون إلى تحليل بيانات أرشيفية من تجربة المصدر المتقطع والانفجاري (BATSE) للبحث عن إشارات دورية مشابهة في انفجارات أشعة غاما القصيرة، وخاصة تلك التي تقل مدتها عن ثانيتين. تستخدم الدراسة نهجًا منهجيًا للكشف عن إشارات دورية عالية التردد في بيانات الأحداث المؤرخة زمنياً لـ 532 انفجار أشعة غاما قصير، مع التركيز على الاعتماد الزمني والطاقة للتعديلات المحتملة. توضح الورقة هيكلها، موضحة الأدوات، وطرق تقليل البيانات، وتحليل التكرارية، مما يؤدي إلى مناقشة الآثار لفهم المحركات المركزية لانفجارات أشعة غاما القصيرة.
طرق
في هذا القسم، يوضح المؤلفون منهجيتهم للكشف عن التكرارية المتماسكة في بيانات الأحداث المجزأة باستخدام اختبار رايلي. يتم حساب قوة رايلي عند تردد تجريبي \( f \) كالتالي
\[
R(f) = \frac{1}{N} \left[ \left( \sum_{j=1}^{N} \cos(2\pi f t_j) \right)^2 + \left( \sum_{j=1}^{N} \sin(2\pi f t_j) \right)^2 \right],
\]
حيث \( N \) تشير إلى عدد الفوتونات التي تصل في الوقت \( t_j \). يمتد التحليل عبر ترددات تجريبية من 500 إلى 2500 هرتز، مع تسجيل الحد الأقصى لقوة رايلي لكل جزء، مما يسمح بتحديد الترددات المرشحة. يطبق المؤلفون نظرية القيم المتطرفة لتقييم دلالة القمم المكتشفة، نمذجة أقصى قوى رايلي باستخدام توزيع غامبل المميز بمعلمات \( \mu \) و \( \beta \).
للتحقق من الإشارات الدورية المحتملة، ينفذ المؤلفون نهجًا إحصائيًا من خطوتين. أولاً، يثبتون أن تعديلًا دوريًا حقيقيًا يجب أن يتجاوز توقع الضوضاء، الذي يتم تحديده من خلال تجارب خلفية مستقلة. ثانيًا، يقيمون الدلالة العالمية لأي قيم متطرفة من خلال مقارنتها بالتوقع الإحصائي عبر جميع الترددات التجريبية. يتم حساب احتمال الإنذار الكاذب (FAP) لتحديد احتمال ملاحظة تقلبات الضوضاء التي تساوي أو تتجاوز القوة الملحوظة، مع تطبيق حد علوي محافظ من بونفيروني إذا لم يتم ضمان الاستقلالية بين التجارب بالكامل. تمكن هذه المنهجية الشاملة من الكشف القوي عن الإشارات الدورية وسط ضوضاء عشوائية.
النتائج
في تحليل التكرارية ضمن مجموعة بيانات تتكون من 186,327 جزءًا من انفجارات أشعة غاما القصيرة (GRBs)، تم إجراء فحص أولي لاستبعاد القيم الشاذة ذات القوة العالية المتأثرة بتلوث الضوضاء الحمراء. ومن الجدير بالذكر، تم تحديد GRB 980310B (رقم الزناد 6634) كحالة شاذة بسبب ذروته القصيرة والعنيفة في منحنى الضوء، مما أنتج فائضًا كبيرًا من الضوضاء الحمراء. وبالتالي، تم استبعاد جميع أجزاء البيانات المرتبطة بهذا GRB من التحليل الإضافي، مما أدى إلى مجموعة بيانات مصقولة تتكون من 186,127 قوة مرشحة \( R_x \) وتردداتها المقابلة \( f_x \).
لتقييم الدلالة الإحصائية لهذه المرشحات، تم ملاءمة دالة كثافة الاحتمال غامبل مع المجموعة العالمية من \( R_x \)، مما أسس عتبة مرجعية \( R_{50} \) التي تفصل بين الأنظمة المهيمنة على الضوضاء وتلك المعززة بالإشارة. تم الاحتفاظ بالمرشحات التي لديها \( R_x > R_{50} \) لمزيد من التقييم. من بين المرشحات، تم تسليط الضوء على زنادين: الزناد 6412 (GRB 971004) أظهر ميزة ذات دلالة هامشية عند 625 هرتز، والتي اعتبرت غير مرجحة لتمثيل إشارة دورية حقيقية بسبب عزلتها في كاشف واحد وقناة طاقة واحدة. بالمقابل، عرض الزناد 5502 (GRB 960616) ميزة ذات دلالة عند 1083 هرتز، مع احتمال إنذار كاذب مصحح تجريبيًا يبلغ حوالي \( 5.54 \times 10^{-4} \). تم دعم هذه الميزة من خلال أجزاء متداخلة عبر عدة كواشف وقنوات طاقة، جميعها تتزامن مع حلقة الانبعاث الرئيسية، مما يشير إلى تعديل دوري حقيقي. ستركز القسم التالي على التحقق من هذه التكرارية في GRB 960616.
مناقشة
في هذا القسم، يناقش المؤلفون الأدوات وأنواع البيانات المستخدمة في تحليلهم للمتغيرات عالية الطاقة التي تم الكشف عنها بواسطة أداة BATSE على CGRO. يتكون نظام BATSE من ثمانية وحدات كاشف غير موجهة، كل منها تحتوي على كاشف منطقة كبيرة (LAD) وكاشف طيفي (SD)، مما يوفر تغطية شاملة للسماء للكشف عن انفجارات أشعة غاما (GRBs). تم تحسين LAD لاكتشاف الانفجارات، بينما يعزز SD الدقة الطيفية عبر نطاق طاقة واسع. يلتقط وضع الأحداث المؤرخة زمنياً (TTE) أحداث الفوتون الفردية بدقة زمنية عالية، مما يسمح بإجراء دراسات مفصلة لخصائص الانفجارات، بما في ذلك التكرارات المعتمدة على الطاقة.
استخرج المؤلفون بيانات TTE لـ 532 انفجار أشعة غاما قصير تم تسجيلها بين مايو 1991 ومايو 2000، مستخدمين خوارزمية كتل بايزيانية لتقسيم البيانات إلى فترات متجانسة إحصائيًا مناسبة لتحليل التوقيت عالي الدقة. أجروا بحثًا شاملًا عن الإشارات الدورية، مستخدمين تحويل الموجة الوزنية Z (WWZ) لإنشاء طيف زمني-ترددي، والذي كشف عن تذبذب كبير عند حوالي 1100 هرتز في بعض قنوات الطاقة. وُجد أن هذا التذبذب متماسك وقريب من الدوري، مع عامل جودة يشير إلى تعديل قوي للانبعاث الفوري. شمل التحليل أيضًا محاكاة مونت كارلو لتقييم دلالة التكرارية المكتشفة، مؤكدًا احتمال إنذار كاذب يبلغ حوالي \(2.69 \times 10^{-6}\)، مما يشير إلى كشف قوي لإشارة 1100 هرتز. بشكل عام، تشير النتائج إلى تعديل متماسك للغاية في انبعاث GRB، مع آثار لفهم العمليات الفيزيائية الأساسية التي تحرك هذه الأحداث عالية الطاقة.
DOI: https://doi.org/10.3847/1538-4357/ae3a96
Publication Date: 2026-02-17
Author(s): Run-Chao Chen et al.
Primary Topic: Gamma-ray bursts and supernovae
Overview
In recent studies, coherent oscillations at kilohertz frequencies have been observed in select gamma-ray bursts (GRBs), indicating potential quasi-periodic dynamics within their central engines. A notable instance is GRB 230307A, which displayed a coherent, energy-dependent periodic signal interpreted as a spin signature of a newly formed millisecond magnetar resulting from a compact binary merger. In light of these findings, we undertook an extensive search for similar signals across 532 short GRBs utilizing time-tagged event data from the Burst and Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory.
Our analysis revealed one statistically significant case: GRB 960616 (BATSE trigger 5502), characterized by a main emission episode lasting approximately 30 ms, which exhibited coherent modulation at 1100 Hz, particularly pronounced above 320 keV, with a fractional amplitude of about 47%. Through data-driven Monte Carlo simulations, we established an upper limit of approximately 8% on the fraction of total radiated energy that can be attributed to this quasi-periodic oscillation (QPO). This finding, alongside the observed periodicity of approximately 0.91 ms, bolsters the hypothesis that millisecond periodicities may emerge during GRBs in scenarios involving mergers.
Introduction
The introduction of the paper discusses the progenitor systems of gamma-ray bursts (GRBs), which are classified into two main types: massive star and compact binary systems. These systems undergo catastrophic events that lead to the formation of a central engine, typically a hyper-accreting black hole or a rapidly rotating neutron star (millisecond magnetar), which powers a relativistic jet responsible for gamma-ray emissions. While the connection between GRBs and their progenitor systems has been established through various observations, the exact nature of the central engines remains uncertain.
Recent studies have provided new insights into the central engines of GRBs through the detection of high-frequency variability in prompt emissions. Notably, a quasi-periodic oscillation (QPO) at approximately 2.7 kHz was reported in a short GRB, suggesting the presence of a hypermassive neutron star. Furthermore, a stable 909 Hz periodic modulation was identified in GRB 230307A, linking it to a nascent millisecond magnetar. Motivated by these findings, the authors aim to analyze archival data from the Burst and Transient Source Experiment (BATSE) to search for similar periodic signals in short GRBs, particularly those with durations under 2 seconds. The study employs a systematic approach to detect high-frequency periodic signals in the Time-Tagged Event data of 532 short GRBs, with a focus on the time and energy dependence of potential modulations. The paper outlines its structure, detailing the instrumentation, data reduction methods, and the analysis of periodicity, culminating in a discussion of the implications for understanding the central engines of short GRBs.
Methods
In this section, the authors outline their methodology for detecting coherent periodicity in segmented event data using the Rayleigh test. The Rayleigh power at a trial frequency \( f \) is calculated as
\[
R(f) = \frac{1}{N} \left[ \left( \sum_{j=1}^{N} \cos(2\pi f t_j) \right)^2 + \left( \sum_{j=1}^{N} \sin(2\pi f t_j) \right)^2 \right],
\]
where \( N \) denotes the number of photons arriving at time \( t_j \). The analysis spans trial frequencies from 500 to 2,500 Hz, with the maximum Rayleigh power recorded for each segment, allowing the identification of candidate frequencies. The authors apply extreme-value theory to assess the significance of detected peaks, modeling the maximal Rayleigh powers using a Gumbel distribution characterized by parameters \( \mu \) and \( \beta \).
To validate potential periodic signals, the authors implement a two-step statistical approach. First, they establish that a genuine periodic modulation must exceed the noise expectation, determined through independent background trials. Second, they evaluate the global significance of any extreme values by comparing them against the statistical expectation across all trial frequencies. The false-alarm probability (FAP) is calculated to quantify the likelihood of observing noise fluctuations that meet or exceed the observed power, with a conservative Bonferroni upper bound applied if independence among trials is not fully assured. This comprehensive methodology enables robust detection of periodic signals amidst stochastic noise.
Results
In the analysis of periodicity within a dataset of 186,327 segments of short gamma-ray bursts (GRBs), a preliminary screening was conducted to exclude high-power outliers influenced by red-noise contamination. Notably, GRB 980310B (trigger number 6634) was identified as an anomalous case due to its intense, short spike in the light curve, which produced significant red-noise excess. Consequently, all data segments associated with this GRB were excluded from further analysis, resulting in a refined dataset of 186,127 candidate powers \( R_x \) and their corresponding frequencies \( f_x \).
To assess the statistical significance of these candidates, a Gumbel probability density function was fitted to the global ensemble of \( R_x \), establishing a reference threshold \( R_{50} \) that delineates noise-dominated from signal-enhanced regimes. Candidates with \( R_x > R_{50} \) were retained for further evaluation. Among the candidates, two triggers were highlighted: trigger 6412 (GRB 971004) exhibited a marginally significant feature at 625 Hz, which was deemed unlikely to represent a genuine periodic signal due to its isolation in a single detector and energy channel. In contrast, trigger 5502 (GRB 960616) displayed a significant feature at 1083 Hz, with a trial-corrected false-alarm probability of approximately \( 5.54 \times 10^{-4} \). This feature was supported by overlapping segments across multiple detectors and energy channels, all coinciding with the main emission episode, suggesting a genuine periodic modulation. The subsequent section will focus on the validation and characterization of this periodicity in GRB 960616.
Discussion
In this section, the authors discuss the instrumentation and data types utilized in their analysis of high-energy transients detected by the BATSE instrument on the CGRO. The BATSE system comprises eight uncollimated detector modules, each featuring a Large Area Detector (LAD) and a Spectroscopy Detector (SD), which together provide comprehensive sky coverage for detecting gamma-ray bursts (GRBs). The LAD is optimized for burst detection, while the SD enhances spectral resolution across a broad energy range. The Time-Tagged Event (TTE) mode captures individual photon events with high temporal resolution, allowing for detailed studies of burst characteristics, including energy-dependent periodicities.
The authors extracted TTE data for 532 short GRBs recorded between May 1991 and May 2000, employing a Bayesian blocks algorithm to segment the data into statistically homogeneous intervals suitable for high-resolution timing analysis. They conducted a thorough search for periodic signals, utilizing the weighted wavelet Z-transform (WWZ) to generate time-frequency spectrograms, which revealed a significant oscillation at approximately 1100 Hz in certain energy channels. This oscillation was found to be coherent and nearly periodic, with a quality factor suggesting a strong modulation of the prompt emission. The analysis further included Monte Carlo simulations to assess the significance of the detected periodicity, confirming a false-alarm probability of approximately \(2.69 \times 10^{-6}\), indicating a robust detection of the 1100 Hz signal. Overall, the findings suggest a highly coherent modulation in the GRB emission, with implications for understanding the underlying physical processes driving these high-energy events.