SwiftC : طيف الطاقة الزاوي القابل للاشتقاق بسرعة خارج حدود ليمبر
SwiftC : fast differentiable angular power spectra beyond Limber

المجلة: Journal of Cosmology and Astroparticle Physics، المجلد: 2026، العدد: 1
DOI: https://doi.org/10.1088/1475-7516/2026/01/043
تاريخ النشر: 2026-01-01
المؤلف: Zhenyun Du وآخرون
الموضوع الرئيسي: المجرات: التكوين، التطور، الظواهر

نظرة عامة

يقدم هذا القسم SwiftC ℓ، وهو خط أنابيب حسابي يعتمد على JAX مصمم لحساب طيف القوة الزاوي بسرعة ودقة تتجاوز تقريب ليمبر، وهو أمر حاسم لتفسير البيانات من المسوحات الواسعة من المرحلة الرابعة القادمة. باستخدام طريقة مبتكرة تعتمد على FFTLog، يحقق SwiftC ℓ دقة عالية من خلال عوامل نمو وانحياز تعتمد على k، مما يتناسب مع مجموعة متنوعة من الاستقصاءات الكونية مثل تجمع المجرات، والانحراف الضعيف، وتأثيرات الخلفية الكونية الميكروية (CMB). يظهر خط الأنابيب كفاءة ملحوظة، حيث يحسب 120 طيف قوة زاوي عبر 1,030 مضاعف في 5 مللي ثانية فقط على نواة GPU واحدة، وهو أسرع بحوالي 40 مرة من المدخل الرائد في تحدي N5K مع الحفاظ على دقة مماثلة.

في الختام، يتم وضع SwiftC ℓ كأداة حيوية لعلم الكونيات عالي الدقة، خاصة للتحليلات التي تتضمن عدم التوزيع غير الغاوسي الأولي وتأثير ساكس-ولف المتكامل، حيث يكون الكشف عن الإشارة حرجًا عند المضاعفات المنخفضة. تسهل قابلية تمييز خط الأنابيب أخذ العينات المعتمدة على التدرجات وتوقعات فيشر القوية، مما يجعله مناسبًا للتحليلات المستقبلية للهياكل واسعة النطاق. يُنصح المستخدمون بتحسين المعلمات، خاصة عدد النقاط في تحليل FFTLog، لضمان الدقة، خاصة في الصناديق الضيقة. بشكل عام، من المتوقع أن يعزز SwiftC ℓ بشكل كبير قدرات التحليل للمسوحات الواسعة القادمة، مما يوفر تقديرات سريعة وموثوقة للمعلمات الكونية.

مقدمة

ت outlines المقدمة ضرورة وجود خطوط أنابيب تحليل متقدمة في ضوء المسوحات القادمة للهيكل واسع النطاق من المرحلة الرابعة، بما في ذلك مسح فيرا سي. روبين لوراثة الفضاء والزمن (LSST)، ويوديك، وأداة الطيف الطاقي المظلم (DESI). تعد هذه المسوحات بمجموعات بيانات شاسعة تتطلب كل من الدقة العالية والكفاءة الحسابية في تحليلها. غالبًا ما تكون الانزياحات الحمراء للأجسام المرصودة غير دقيقة، مما يعقد حساب طيف القوة الزاوي المتوقع، وهو أمر أساسي لفهم الهياكل الكونية. تبسط الطرق التقليدية، مثل تقريب ليمبر، الحسابات لكنها تفشل عند المقاييس الزاوية الأكبر، حيث تحدث ظواهر فلكية كبيرة.

لمعالجة هذه التحديات، يقدم المؤلفون SwiftC ℓ، وهو خط أنابيب جديد وسريع وقابل للاشتقاق لحساب طيف القوة الزاوي الذي يتجاوز قيود تقريب ليمبر. يستخدم SwiftC ℓ طريقة تعتمد على FFTLog تلتقط بدقة اعتماد k واعتماد الانزياح الأحمر لطيف القوة. يدمج مجموعة متنوعة من الاستقصاءات الفلكية، بما في ذلك تجمع المجرات والانحراف الضعيف، ويستفيد من مكتبة Python-JAX للاشتقاق التلقائي، مما يجعله مناسبًا لأخذ العينات المعتمدة على التدرجات وتوقعات فيشر. يتابع البحث لتفصيل الإطار النظري، وطريقة FFTLog، ومقارنات الأداء مع الأكواد الموجودة، وينتهي بتطبيقات على بيانات مشابهة لـ LSST من خلال تحليلات إحصائية متنوعة.

النتائج

في هذا القسم، يقدم المؤلفون نتائج تحليلهم 3 × 2pt باستخدام خط أنابيب SwiftC ℓ، مقارنين أدائه مع المشاركين السابقين في التحدي على كل من بنى GPU وCPU. من الجدير بالذكر، أنه على عكس المشاركات الأخرى في التحدي، يستخدم SwiftC ℓ طريقة FFTLog عبر نطاق الزخم الزاوي بالكامل دون التحول إلى تقريب ليمبر. يسمح هذا الاختيار لـ SwiftC ℓ بالحفاظ على متطلبات دقة $\chi^2 < 1$، بينما تستهدف خطوط الأنابيب الأخرى متطلبات أكثر صرامة $\chi^2 < 0.2$ لـ $\ell < 200$. تكشف التحليلات أن SwiftC ℓ ينفذ التحليل الكامل في 0.006 ثوانٍ على GPU، وهو أسرع بحوالي 40 مرة من المدخل الرائد في التحدي، FKEM، و0.5 ثانية على CPU، وهو ما يعادل ضعف الوقت لـ FKEM. تظهر دقة أطياف القوة الزاوي عبر صناديق الانزياح الأحمر المختلفة أن SwiftC ℓ لا يفي فقط بمتطلبات الدقة، بل يتجاوزها، متفوقًا على تقريب ليمبر حتى عند قيم $\ell$ العالية. بالإضافة إلى ذلك، يوضح المؤلفون الإعداد الحسابي لعرض صناديق مختلفة، مشيرين إلى أن تحليل FFTLog يتطلب المزيد من النقاط لإعادة بناء دقيقة لوظائف النافذة الأرق. أوقات التشغيل لصناديق نصف العرض وربع العرض هي حوالي 40 مللي ثانية و10 مللي ثانية، على التوالي. تشمل النتائج أيضًا مقارنة استعادة المعلمات باستخدام أخذ عينات جماعية مع كل من بيانات FKEM وبيانات تقريب ليمبر. تشير النتائج إلى أن متجه بيانات FKEM يستعيد المعلمات بفعالية، بينما يقدم تقريب ليمبر بعض التحولات. استخدم أخذ العينات الجماعية على وحدة معالجة الرسوميات Nvidia Tesla A100 12 متجولًا مع 2000 نقطة لكل منها، مكتملًا في 162 ثانية، بينما استغرق أسلوب هاملتونيان مونت كارلو (HMC) 23 دقيقة و34 ثانية لسلسلة واحدة من 2000 نقطة.

المناقشة

في هذا القسم، يناقش المؤلفون الإطار النظري والأساليب الحسابية لحساب طيف القوة الزاوي باستخدام مجموعة متنوعة من الاستقصاءات الكونية ضمن خط أنابيب SwiftC ℓ. يتم تعريف طيف القوة الزاوي \( C_{ij}^{\ell} \) من حيث نوى خط الرؤية \( \Delta_{i,j}^{\ell}(k) \) وطيف القوة المادية \( P(k, \chi, \chi’) \)، مع دمج المساهمات من تجمع المجرات، والانحراف الضعيف، وCMB. يوضح المؤلفون الصيغ الرياضية لمساهمات مختلفة، بما في ذلك تجمع المجرات، وانحياز التكبير، وتشوهات الفضاء الانزياحي (RSD)، وعدم التوزيع غير الغاوسي المحلي الأولي، وتأثيرات CMB، مع التأكيد على مرونة نهجهم لاستيعاب كل من النماذج الخطية وغير الخطية.

تستخدم الاستراتيجية الحسابية طريقة FFTLog للتكامل العددي الفعال، مما يسمح بمعالجة عوامل النمو المعتمدة على k دون الحاجة إلى تقسيم طيف القوة. يتحقق المؤلفون من تنفيذهم مقابل المعايير من تحدي N5K، موضحين أن SwiftC ℓ يحقق توافقًا أقل من 1% مع الأكواد الموجودة بينما يكون أسرع بحوالي 40 مرة من أفضل مدخل. يستكشفون أيضًا أداء خط الأنابيب من خلال تحليلات MCMC وهمية وتوقعات فيشر، مما يبرز فائدته للدراسات المستقبلية للهياكل واسعة النطاق وقدرته على التعامل مع مساحات المعلمات المعقدة بكفاءة. يخلص المؤلفون إلى أن SwiftC ℓ هو أداة قيمة لتحليلات علم الكونيات عالية الدقة، خاصة في سياق المسوحات الواسعة القادمة.

Journal: Journal of Cosmology and Astroparticle Physics, Volume: 2026, Issue: 1
DOI: https://doi.org/10.1088/1475-7516/2026/01/043
Publication Date: 2026-01-01
Author(s): Zhenyun Du et al.
Primary Topic: Galaxies: Formation, Evolution, Phenomena

Overview

The section presents SwiftC ℓ, a JAX-based computational pipeline designed for the fast and accurate calculation of the angular power spectrum beyond the Limber approximation, crucial for interpreting data from upcoming stage IV wide-field surveys. Utilizing an innovative FFTLog-based method, SwiftC ℓ achieves high precision through k-dependent growth factors and biases, accommodating various cosmological probes such as galaxy clustering, weak lensing, and cosmic microwave background (CMB) effects. The pipeline demonstrates remarkable efficiency, computing 120 angular power spectra over 1,030 multipoles in just 5 milliseconds on a single GPU core, and is approximately 40 times faster than the leading entry in the N5K challenge while maintaining comparable accuracy.

In conclusion, SwiftC ℓ is positioned as a vital tool for high-precision cosmology, particularly for analyses involving primordial non-Gaussianity and the integrated Sachs-Wolfe effect, where signal detection is critical at low multipoles. The pipeline’s differentiability facilitates gradient-based sampling and robust Fisher forecasts, making it suitable for future large-scale structure analyses. Users are advised to optimize parameters, particularly the number of points in the FFTLog decomposition, to ensure accuracy, especially in narrow bins. Overall, SwiftC ℓ is expected to significantly enhance the analysis capabilities of upcoming wide-field surveys, providing rapid and reliable cosmological parameter estimations.

Introduction

The introduction outlines the necessity for advanced analysis pipelines in light of upcoming stage IV large-scale structure surveys, including the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST), Euclid, and the Dark Energy Spectroscopic Instrument (DESI). These surveys promise vast datasets that demand both high precision and computational efficiency in their analysis. The redshifts of observed objects are often imprecise, complicating the computation of the projected angular power spectrum, which is essential for understanding cosmic structures. Traditional methods, such as the Limber approximation, simplify calculations but fail at larger angular scales, where significant astrophysical phenomena occur.

To address these challenges, the authors introduce SwiftC ℓ, a novel, fast, and differentiable pipeline for computing the angular power spectrum that transcends the limitations of the Limber approximation. SwiftC ℓ employs an FFTLog-based method that accurately captures the k-dependence and redshift dependence of the power spectrum. It integrates a variety of astrophysical probes, including galaxy clustering and weak lensing, and leverages the Python-JAX library for automatic differentiation, making it suitable for gradient-based sampling and Fisher forecasts. The paper proceeds to detail the theoretical framework, the FFTLog method, and performance comparisons with existing codes, culminating in applications to LSST-like data through various statistical analyses.

Results

In this section, the authors present the results of their 3 × 2pt analysis using the SwiftC ℓ pipeline, comparing its performance against previous challenge entrants on both GPU and CPU architectures. Notably, unlike the other challenge entries, SwiftC ℓ employs the FFTLog method across the entire angular momentum range without switching to the Limber approximation. This choice allows SwiftC ℓ to maintain a $\chi^2 < 1$ accuracy requirement, while other pipelines target a stricter $\chi^2 < 0.2$ for $\ell < 200$. The analysis reveals that SwiftC ℓ executes the full analysis in 0.006 seconds on a GPU, approximately 40 times faster than the leading challenge entry, FKEM, and 0.5 seconds on a CPU, which is about twice as long as FKEM. The accuracy of the angular power spectra across various redshift bins demonstrates that SwiftC ℓ not only meets but exceeds the accuracy requirements, outperforming the Limber approximation even at high $\ell$ values. Additionally, the authors detail the computational setup for different bin widths, indicating that the FFTLog decomposition necessitates more points for accurate reconstruction of thinner window functions. The running times for half-width and quarter-width bins are approximately 40 ms and 10 ms, respectively. The results also include a comparison of parameter recovery using ensemble sampling with both FKEM and Limber approximation data vectors. The findings indicate that the FKEM data vector effectively recovers parameters, while the Limber approximation introduces some shifts. The ensemble sampling on the Nvidia Tesla A100 GPU node utilized 12 walkers with 2000 points each, completing in 162 seconds, while the Hamiltonian Monte Carlo (HMC) method took 23 minutes and 34 seconds for a single chain of 2000 points.

Discussion

In this section, the authors discuss the theoretical framework and computational methods for calculating the angular power spectrum using various cosmological probes within the SwiftC ℓ pipeline. The angular power spectrum \( C_{ij}^{\ell} \) is defined in terms of line-of-sight kernels \( \Delta_{i,j}^{\ell}(k) \) and the matter power spectrum \( P(k, \chi, \chi’) \), incorporating contributions from galaxy clustering, weak lensing, and the CMB. The authors detail the mathematical formulations for different contributions, including galaxy clustering, magnification bias, redshift space distortions (RSD), primordial local non-Gaussianity, and CMB effects, emphasizing the flexibility of their approach to accommodate both linear and non-linear models.

The computational strategy employs the FFTLog method for efficient numerical integration, allowing for the treatment of k-dependent growth factors without necessitating a split in the power spectrum. The authors validate their implementation against benchmarks from the N5K challenge, demonstrating that SwiftC ℓ achieves subpercent agreement with existing codes while being approximately 40 times faster than the best-performing entry. They also explore the performance of the pipeline through mock MCMC and Fisher forecasting analyses, highlighting its utility for future large-scale structure studies and its capability to handle complex parameter spaces efficiently. The authors conclude that SwiftC ℓ is a valuable tool for high-precision cosmological analyses, particularly in the context of upcoming wide-field surveys.