DOI: https://doi.org/10.3847/1538-4357/ae3d95
تاريخ النشر: 2026-02-26
المؤلف: H. Thankful Cromartie وآخرون
الموضوع الرئيسي: أبحاث النباضات والموجات الجاذبية
نظرة عامة
تقدم هذه القسم تحليلًا شاملاً لخمسة عشر عامًا من ملاحظات التوقيت لنجم النبضات الملي ثانية PSR J1231-1411، باستخدام بيانات من تلسكوبات الراديو نانسي وجرين بانك. يهدف البحث إلى تحديد كتلة النجم النابض من خلال الكشف عن تأخير شابيرو النسبي، باستخدام كل من الطرق التكرارية والبايزية. تسفر النتائج عن تقديرات لكتلة الرفيق والنجم النابض بقيم $m_c = 0.23^{+0.09}_{-0.06} M_\odot$ و $m_p = 1.87^{+1.11}_{-0.67} M_\odot$ (فترة ثقة 68.3%)، مع ميل مداري أكثر دقة بقيمة $i = 79.80^{+3.47}_{-4.70}^\circ$. من خلال تقييد الحد الأقصى لكتلة النجم النابض إلى $3 M_\odot$، يتم تحسين تقدير كتلة النجم النابض إلى $m_p = 1.71^{+0.70}_{-0.56} M_\odot$.
يتضمن التحليل ملاحظات جديدة من تلسكوب جرين بانك إلى جانب بيانات أرشيفية، باستخدام طرق سلسلة ماركوف مونت كارلو (MCMC) لاستنتاج قيم الكتلة والميل. النتائج متوافقة مع التوقعات النظرية وتبرز التحديات في قياس تأخير شابيرو في الأنظمة ذات الرفقاء منخفضي الكتلة والمدارات غير الحادة. يؤكد البحث على الحاجة إلى توقيت عالي الدقة لتحديد كتل النجوم النابضة بدقة، خاصة في الحالات التي تُعلم فيها هذه القياسات القيود على معادلة حالة النجم النيوتروني الداخلي. من المتوقع أن تعزز الحملات الرصدية المستقبلية مع المرافق من الجيل التالي دقة قيود الكتلة لمثل هذه النجوم النابضة.
مقدمة
تسلط المقدمة الضوء على أهمية الملاحظات متعددة الرسائل، وخاصة في نطاقات الموجات الكهرومغناطيسية وموجات الجاذبية، لتحديد معادلة الحالة (EoS) لمادة ذات كثافة فوق النووية. تشمل المنهجيات الرئيسية قياس تأخير شابيرو النسبي في نجوم النبضات الملي ثانية (MSPs)، مما أتاح تحديد كتل دقيقة لأثقل النجوم النيوترونية، وبالتالي تقييد معادلة الحالة الداخلية لها. بالإضافة إلى ذلك، يسمح نمذجة منحنى الضوء بالأشعة السينية باستخدام مستكشف تكوين النجم النيوتروني (NICER) التابع لناسا بقياسات متزامنة للكتلة ($M$) ونصف القطر ($R$)، مع تأثيرات النسبية العامة التي تؤثر على الانبعاثات المرصودة من الأشعة السينية بناءً على كثافة النجوم.
تدمج دراسة بارزة قام بها رايماكرز وآخرون (2021) بيانات رصدية متنوعة، بما في ذلك قياسات NICER وXMM-Newton، وملاحظات تأخير شابيرو الراديو، وأحداث موجات الجاذبية، لتقييد نسبة الكتلة إلى نصف القطر للنجوم النيوترونية. يذكرون أن نجمًا نيوترونيًا بكتلة $1.4 \, M_\odot$ له نصف قطر يقارب $12.33^{+0.76}_{-0.81}$ كم أو $12.18^{+0.56}_{-0.79}$ كم، اعتمادًا على معلمات EoS المستخدمة. يقدم المؤلفون تحليل توقيت النجوم النابضة لـ MSP J1231-1411، الذي، على الرغم من دقة محدودة بسبب هندسة مدار النظام، يوفر قياسات مهمة لتأخير شابيرو التي تعمل كأولويات محافظة لنمذجة منحنى الضوء لـ NICER. تساعد هذه الطريقة في حل التداخل في نسبة الكتلة إلى نصف القطر والتفريق بين النماذج المتنافسة لهندسة النقاط الساخنة. ستفصل الأقسام التالية من الورقة الطرق الرصدية، وتحليل البيانات، والنتائج المتعلقة بتقدير كتلة النجم النابض.
طرق
في هذا القسم، يوضح المؤلفون الطرق المستخدمة لقياس تأثير تأخير شابيرو في ثنائيات النجوم النابضة الملي ثانية (MSP)، وهو ظاهرة تنبأت بها النسبية العامة. يحدث تأخير شابيرو عندما يمر MSP خلف رفيقه النجمي، مما يتسبب في تأخير في وقت وصول نبضاته بسبب تشويه الزمكان للرفيق. يكون هذا التأخير أكثر وضوحًا عند مرحلة مدارية تبلغ 0.25 (التقاطع العلوي) ويمكن قياسه من خلال بقايا التوقيت، التي تظهر عادةً قمة تبلغ حوالي 10 ميكروثانية. من خلال تحليل مدى وشكل هذا التأخير، يستنتج المؤلفون المعلمات بعد كبلر، تحديدًا النطاق ($r$) والشكل ($s$)، والتي تعتبر حاسمة لتحديد كتل النجم النابض ($m_p$) ورفيقه ($m_c$) تحت فرضية النسبية العامة.
تشمل المنهجية استخدام برنامج توقيت النجوم النابضة PINT لتناسب أوقات الوصول الفرعية (TOAs) من الملاحظات في تلسكوبات نانسي وجرين بانك. يستخدم المؤلفون كل من تقنية التناسب التقليدية $\chi^2$ وتقنيات سلسلة ماركوف مونت كارلو (MCMC)، مع تقدير الضوضاء البايزية الذي يتم باستخدام برنامج ENTERPRISE. يستخدم التحليل تقويم DE440 الكوكبي من مختبر الدفع النفاث التابع لناسا ومقياس الزمن TDB لضمان الدقة في القياسات. أشارت النتائج الأولية من الملاحظات السابقة إلى نظام ذو ميل مرتفع مع MSP منخفض الكتلة نسبيًا، مما يبرز أهمية قياسات الكتلة الدقيقة لتقييد معادلة الحالة (EoS) للنجوم النيوترونية.
نتائج
في هذا القسم، يحلل المؤلفون تعقيدات مجموعة البيانات المتعلقة بـ PSR J1231-1411، مشيرين إلى التحديات مثل التلألؤ القوي، وتكوينات الخلفية غير المتجانسة، والضوضاء الكبيرة التي تؤثر على دقة التوقيت. يؤكدون على أهمية استخدام الاستدلال البايزي بدلاً من الانحدار الخطي البسيط بسبب وجود التغايرات بين المعلمات، خاصة في سياق قياس تأخير شابيرو، الذي يكون حساسًا لزاوية الميل المداري. يستخدم المؤلفون ثلاث طرق متميزة للتحليل: مخطط شبكة $\chi^2$ مع معلمات تأخير شابيرو ثابتة، وتناسب بايزي مع أولويات فلكية، وتناسب عام بالكامل يأخذ عينات من الضوضاء ونماذج التوقيت في وقت واحد.
تشير النتائج من تجارب MCMC إلى أن كتلة الرفيق ($m_c$) مقيدة بـ $0.23^{+0.09}_{-0.06} \, M_\odot$ مع ميل قدره $79.80^{+3.47}_{-4.70}^\circ$، وكتلة النجم النابض ($m_p$) تبلغ $1.87^{+1.11}_{-0.67} \, M_\odot$. تعطي التناسبات غير الخطية البايزية بالكامل نتائج مشابهة، مع $m_c = 0.21^{+0.06}_{-0.05} \, M_\odot$ وميل $81.06^{+3.36}_{-3.51}^\circ$، مما يشير إلى تحسين القيود على هذه المعلمات. تشير النتائج إلى أنه بينما تكون الطريقة الثابتة للضوضاء فعالة، فإن الطريقة غير الخطية تظهر وعدًا للتحليلات المستقبلية، على الرغم من متطلبات حسابية أعلى. بشكل عام، يبرز البحث أهمية تقنيات النمذجة المتقدمة في تحديد المعلمات الفلكية بدقة في أنظمة النجوم النابضة.
مناقشة
في هذا القسم، يناقش المؤلفون التحليل والنتائج المتعلقة بالنجم النابض PSR J1231-1411، الذي تم اكتشافه خلال بحث عن النجوم النابضة الراديوية التي تستهدف مصادر أشعة غاما الساطعة غير المرتبطة بفيرمي-لات. يتميز النجم النابض بأنه نجم نابض ملي ثانية (MSP) نموذجي بمدة 3.68 مللي ثانية في نظام ثنائي مع رفيق قزم أبيض. يوضح المؤلفون الحملات الرصدية الواسعة التي أجريت باستخدام تلسكوب جرين بانك (GBT) وتلسكوب نانسي الراديو (NRT)، مشيرين إلى التحديات التي تطرحها التلألؤ الانكساري وأهمية دقة التوقيت في قياس معلمات النجم النابض. يذكرون نتائج مهمة، بما في ذلك الكشف الأول عن تدهور فترة المدار ($\dot{P}_b$) وقياسات المنظور، والتي تعتبر حاسمة لتحديد كتل النجم النابض ورفيقه.
يؤكد المؤلفون على تعقيد تحليل توقيت النجوم النابضة، مشيرين إلى أن التحولات المختلفة للنموذج تؤثر بشكل كبير على الضوضاء البيضاء والحمراء الداخلية المقاسة، وكذلك على الكتل المستنتجة. يبرزون الطبيعة غير العادية لـ PSR J1231-1411 في سياق قياسات الكتلة الدقيقة، حيث أن مداره قريب من الدائرة وقربه يسمح بالكشف عن التأثيرات المدارية طويلة الأمد. تشير النتائج إلى قيود متواضعة على كتلة الرفيق ($m_c = 0.24^{+0.12}_{-0.07} M_\odot$)، وزاوية الميل ($i = 79.11^{+3.96}_{-5.69}^\circ$)، وكتلة النجم النابض ($m_p = 2.00^{+1.53}_{-0.76} M_\odot$). يناقش المؤلفون أيضًا آثار نتائجهم على الملاحظات المستقبلية وإمكانية تحسين دقة التوقيت مع تلسكوبات الراديو من الجيل التالي، والتي يمكن أن تعزز فهم النجوم النابضة وأنظمتها الثنائية.
DOI: https://doi.org/10.3847/1538-4357/ae3d95
Publication Date: 2026-02-26
Author(s): H. Thankful Cromartie et al.
Primary Topic: Pulsars and Gravitational Waves Research
Overview
This section presents a comprehensive analysis of 15 years of timing observations for the millisecond pulsar PSR J1231-1411, utilizing data from the Nançay and Green Bank radio telescopes. The study aims to constrain the pulsar’s mass through the detection of relativistic Shapiro delay, employing both frequentist and Bayesian methods. The findings yield companion and pulsar mass estimates of $m_c = 0.23^{+0.09}_{-0.06} M_\odot$ and $m_p = 1.87^{+1.11}_{-0.67} M_\odot$ (68.3% confidence interval), with a more precise orbital inclination of $i = 79.80^{+3.47}_{-4.70}^\circ$. By restricting the maximum pulsar mass to $3 M_\odot$, the pulsar mass estimate is refined to $m_p = 1.71^{+0.70}_{-0.56} M_\odot$.
The analysis incorporates new Green Bank Telescope observations alongside archival data, employing Markov Chain Monte Carlo (MCMC) methods to derive mass and inclination values. The results are consistent with theoretical predictions and highlight the challenges of measuring Shapiro delay in systems with low-mass companions and non-edge-on orbits. The study emphasizes the need for high-precision timing to accurately determine pulsar masses, particularly in cases where these measurements inform constraints on the neutron star interior equation of state. Future observational campaigns with next-generation facilities are anticipated to enhance the precision of mass constraints for such pulsars.
Introduction
The introduction highlights the significance of multi-messenger observations, particularly in electromagnetic and gravitational wave bands, for determining the equation of state (EoS) of supranuclear-density matter. Key methodologies include the measurement of relativistic Shapiro delay in radio millisecond pulsars (MSPs), which has enabled precise mass determinations of the heaviest neutron stars, thereby constraining their interior EoS. Additionally, X-ray light curve modeling using NASA’s Neutron Star Interior Composition Explorer (NICER) allows for simultaneous mass ($M$) and radius ($R$) measurements, with general relativistic effects influencing the observed X-ray emissions based on stellar compactness.
A notable study by Raaijmakers et al. (2021) integrates various observational data, including NICER and XMM-Newton measurements, radio Shapiro delay observations, and gravitational wave events, to constrain the mass-to-radius ratio of neutron stars. They report that a neutron star with a mass of $1.4 \, M_\odot$ has a radius of approximately $12.33^{+0.76}_{-0.81}$ km or $12.18^{+0.56}_{-0.79}$ km, depending on the EoS parameterization used. The authors present a pulsar timing analysis of MSP J1231-1411, which, despite limited precision due to the system’s orbital geometry, provides significant Shapiro delay measurements that serve as conservative priors for NICER’s light curve modeling. This approach aids in resolving the degeneracy in the mass-to-radius ratio and differentiating among competing models of hot spot geometry. Subsequent sections of the paper will detail the observational methods, data analysis, and findings related to the pulsar’s mass estimation.
Methods
In this section, the authors outline the methods used to measure the Shapiro delay effect in millisecond pulsar (MSP) binaries, a phenomenon predicted by general relativity. The Shapiro delay occurs when an MSP passes behind its stellar companion, causing a delay in the arrival time of its pulsations due to the companion’s spacetime distortion. This delay is most pronounced at an orbital phase of 0.25 (superior conjunction) and can be quantified through timing residuals, which typically show a spike of approximately 10 microseconds. By analyzing the extent and shape of this delay, the authors derive post-Keplerian parameters, specifically the range ($r$) and shape ($s$), which are crucial for determining the masses of the pulsar ($m_p$) and its companion ($m_c$) under the assumption of general relativity.
The methodology includes the use of the PINT pulsar timing software for fitting sub-banded times of arrival (TOAs) from observations at the Nançay and Green Bank telescopes. The authors employ both traditional $\chi^2$ fitting and Markov Chain Monte Carlo (MCMC) techniques, with Bayesian noise estimation conducted using the ENTERPRISE software. The analysis utilizes the DE440 planetary ephemeris from NASA’s Jet Propulsion Laboratory and the TDB timescale to ensure precision in the measurements. The initial findings from earlier observations indicated a highly inclined system with a relatively low-mass MSP, highlighting the importance of precise mass measurements for constraining the equation of state (EoS) of neutron stars.
Results
In this section, the authors analyze the complexities of the data set related to PSR J1231-1411, highlighting challenges such as strong scintillation, heterogeneous backend compositions, and significant noise affecting timing precision. They emphasize the importance of using Bayesian inference over simple linear regression due to the presence of covariances among parameters, particularly in the context of measuring the Shapiro delay, which is sensitive to the orbital inclination angle. The authors employ three distinct methods for analysis: a $\chi^2$ gridding scheme with fixed Shapiro delay parameters, a Bayesian fit with astrophysical priors, and a fully-generalized fit that simultaneously samples noise and timing models.
The results from the MCMC trials indicate that the companion mass ($m_c$) is constrained to $0.23^{+0.09}_{-0.06} \, M_\odot$ with an inclination of $79.80^{+3.47}_{-4.70}^\circ$, and a pulsar mass ($m_p$) of $1.87^{+1.11}_{-0.67} \, M_\odot$. The fully-Bayesian nonlinear MCMC fits yield similar results, with $m_c = 0.21^{+0.06}_{-0.05} \, M_\odot$ and inclination $81.06^{+3.36}_{-3.51}^\circ$, suggesting improved constraints on these parameters. The findings indicate that while the fixed-noise approach is effective, the nonlinear method shows promise for future analyses, albeit with higher computational demands. Overall, the study underscores the significance of advanced modeling techniques in accurately determining astrophysical parameters in pulsar systems.
Discussion
In this section, the authors discuss the analysis and findings related to the pulsar PSR J1231-1411, which was discovered during a radio pulsar search targeting bright unassociated Fermi-LAT γ-ray sources. The pulsar is characterized as a typical 3.68-ms millisecond pulsar (MSP) in a binary system with a white dwarf companion. The authors detail extensive observational campaigns conducted with the Green Bank Telescope (GBT) and the Nançay Radio Telescope (NRT), highlighting the challenges posed by diffractive scintillation and the importance of timing precision in measuring the pulsar’s parameters. They report significant findings, including the first detection of the orbital period decay ($\dot{P}_b$) and parallax measurements, which are crucial for determining the masses of the pulsar and its companion.
The authors emphasize the complexity of pulsar timing analysis, noting that various model permutations significantly affect the measured white and intrinsic red noise, as well as the derived masses. They highlight the unusual nature of PSR J1231-1411 in the context of precision mass measurements, as its orbit is nearly circular and its proximity allows for the detection of long-term orbital effects. The results indicate modest constraints on the companion mass ($m_c = 0.24^{+0.12}_{-0.07} M_\odot$), inclination angle ($i = 79.11^{+3.96}_{-5.69}^\circ$), and pulsar mass ($m_p = 2.00^{+1.53}_{-0.76} M_\odot$). The authors also discuss the implications of their findings for future observations and the potential for improved timing precision with next-generation radio telescopes, which could enhance the understanding of MSPs and their binary systems.