پژوهش همکاری بین‌المللی [[ابزار طیف‌سنجی انرژی تاریک]] تحت عنوان «DESI 2024 VI: قیدهای کیهان‌شناختی از اندازه‌گیری‌های نوسانات آکوستیک باریونی»، نتایج علمی حاصل از تحلیل داده‌های سال اول (DR1) این آزمایش نسل چهارم را گزارش می‌کند ((همکاری بین‌المللی DESI، مقاله DESI 2024 VI: Cosmological Constraints from the Measurements of Baryon Acoustic Oscillations، آرشیو arXiv:2404.03002، سال ۲۰۲۴ میلادی)). این تحلیل با تکیه بر ۶ میلیون کهکشان و اختروش، توزیع ماده در کیهان را در بازه سرخ‌گرایی ${0.1 < z < 4.2}$ نقشه‌برداری کرده و به عنوان یک خط‌کش استاندارد برای بازسازی تاریخچه انبساط عالم به کار رفته است.
<یافته بنیادین پژوهش: هرچند داده‌های فاصله خط‌کش استاندارد DESI در چارچوب مدل کیهان‌شناختی استاندارد تخت $\Lambda\text{CDM}$ سازگار باقی می‌مانند، اما ترکیب آن‌ها با داده‌های تابش زمینه کیهانی و رصدهای ابرنواخترهای نوع Ia، ترجیح آماری معناداری (بین ۲.۵ تا ۳.۹ سیگما) را برای یک انرژی تاریک دینامیکی و متغیر با زمان ($w_0 > -1$ و $w_a < 0$) به جای ثابت کیهان‌شناسی اینشتین نشان می‌دهد.>
# مبانی فیزیکی نوسانات آکوستیک باریونی (BAO) به مثابه خط‌کش استاندارد #
در جهان اولیه و پیش از دوران بازترکیب (حدود ۳۸۰ هزار سال پس از مهبانگ)، فوتون‌ها و باریون‌ها به واسطه پراکندگی کامپتون و برهم‌کنش‌های الکترومغناطیسی، به صورت یک شاره درهم‌تنیده رفتار می‌کردند. امواج صوتی ناشی از برهم‌کنش فشار تابش و گرانش ماده تاریک درون این پلاسما با سرعت ${c_s(z)}$ منتشر می‌شدند:
$$
c_s(z) = \frac{c}{\sqrt{3 \left(1 + \frac{3\rho_b(z)}{4\rho_\gamma(z)}\right)}}
$$
با خنک شدن کیهان در سرخ‌گرایی ${z_d \approx 1060}$، باریون‌ها از کشش فوتون‌ها رها شدند («دوران درگ باریونی») و حرکت امواج صوتی متوقف شد. بیشترین مسافتی که این امواج صوتی در طول این مدت پیمودند، **«افق صوتی در دوران درگ»** «${r_d}$» نامیده می‌شود:
$$
r_d = \int_{z_d}^\infty \frac{c_s(z)}{H(z)} \, dz \approx 147.05 \left(\frac{\Omega_m h^2}{0.1432}\right)^{-0.23} \left(\frac{\Omega_b h^2}{0.02236}\right)^{-0.13} \left(\frac{N_{\text{eff}}}{3.04}\right)^{-0.1} \, \text{Mpc}
$$
این فاصله مشخص (${r_d \sim 150\text{ Mpc}}$) در توزیع کهکشان‌ها به صورت یک قله موضعی در تابع هم‌بستگی دوتایی و نوسانات تناوبی در طیف توان ماده منجمد شده است.
با اندازه‌گیری زاویه ظاهری جدایی کهکشان‌ها در امتداد عمود بر خط دید (${\Delta\theta = r_d / D_M(z)}$) و بازه سرخ‌گرایی در راستای خط دید (${\Delta z = r_d H(z) / c}$)، می‌توان دو کمیت هندسی کلیدی را مستقل از جزئیات اخترفیزیکی مقیاس خرد اندازه گرفت:
1. **فاصله هم‌حرکت عرضی:** ${D_M(z) / r_d}$
2. **فاصله هابل در امتداد خط دید:** ${D_H(z) / r_d = c / (H(z) r_d)}$
در مواردی که نسبت سیگنال به نوفه اجازه تفکیک جهت‌ها را نمی‌دهد، از میانگین زاویه‌ای ${D_V(z) / r_d}$ استفاده می‌شود:
$$
D_V(z) \equiv \left( z D_M(z)^2 D_H(z) \right)^{1/3}
$$
# داده‌های رصدی DESI و نمونه‌های ردگیر #
طیف‌سنج DESI نصب‌شده بر روی تلسکوپ ۴ متری مایال در رصدخانه کیت‌پیک آریزونا، با بهره‌گیری از ۵۰۰۰ فیبر نوری رباتیک در کانون تلسکوپ، طیف نوری اهداف را در گستره وسیعی از کیهان ثبت می‌کند. نمونه‌های استفاده‌شده در تحلیل سال اول (DR1) عبارت‌اند از:
* **کهکشان‌های درخشان (BGS):** تعداد ۳۰۰,۰۱۷ کهکشان در بازه سرخ‌گرایی ${0.1 < z < 0.4}$ با سرخ‌گرایی مؤثر ${z_{\text{eff}} = 0.295}$.
* **کهکشان‌های سرخ درخشان (LRG):** در دو بازه مجزا شامل ۵۰۶,۹۰۵ کهکشان در ${0.4 < z < 0.6}$ و ۷۷۱,۸۷۵ کهکشان در ${0.6 < z < 0.8}$.
* **کهکشان‌های خط نشری (ELG) و نمونه ترکیبی LRG+ELG:** در بازه ${0.8 < z < 1.1}$ تعداد ۱,۸۷۶,۱۶۴ کهکشان ترکیب شده و در بازه ${1.1 < z < 1.6}$ تعداد ۱,۴۱۵,۶۸۷ کهکشان ELG تحلیل شده است.
* **اختروش‌ها به عنوان ردگیر مستقیم (QSO):** تعداد ۸۵۶,۶۵۲ اختروش در بازه ${0.8 < z < 2.1}$ با ${z_{\text{eff}} = 1.491}$.
* **جنگل لایمن-آلفا (Lyman-$\alpha$ Forest):** بررسی خطوط جذبی هیدروژن خنثی در طیف ۷۰۹,۵۶۵ اختروش دوردست در بازه سرخ‌گرایی ${1.77 < z < 4.16}$ با سرخ‌گرایی مؤثر ${z_{\text{eff}} = 2.330}$.
دقت کلی اندازه‌گیری مقیاس ایزوتروپیک BAO در این پیمایش به ۰.۴۹ درصد رسیده است.
# روش تحلیل و پروتکل کورسازی داده‌ها (Blinding) #
برای جلوگیری از پیش‌داوری شناختی و سوگیری تأیید در اندازه‌گیری پارامترهای کیهانی، داده‌های DESI تحت یک فرآیند کورسازی دقیق تحلیل شدند:
* در کاتالوگ کهکشان‌ها و اختروش‌ها، سرخ‌گرایی‌ها با یک تابع ریاضی ناشناخته جابه‌جا شدند تا موقعیت قله BAO دستکاری شود.
* در تحلیل جنگل لایمن-آلفا، یک مؤلفه اضافی ساختگی به توابع هم‌بستگی افزوده شد.
* خط لوله‌های تحلیل، آزمون‌های اعتبارسنجی روی داده‌های شبیه‌سازی‌شده و ارزیابی سیستماتیک‌ها همگی در حالت نابینا نهایی شدند و تنها پس از قفل شدن کامل کدهای محاسباتی، فرآیند رمزگشایی صورت گرفت.
# نتایج در مدل استاندارد کیهان‌شناسی تخت (Flat $\Lambda\text{CDM}$) #
داده‌های حاصل از خط‌کش استاندارد BAO در DESI به تنهایی وابستگی مستقیم به حاصل‌ضرب ثابت هابل و افق صوتی (${H_0 r_d}$) و کسر چگالی ماده کل (${\Omega_m}$) دارند. تحلیل داده‌های DESI به تنهایی در مدل استاندارد تخت نتیجه زیر را به دست می‌دهد:
$$
\Omega_m = 0.295 \pm 0.015, \quad r_d h = (101.8 \pm 1.3) \, \text{Mpc}
$$
برای شکستن تباهیدگی میان ${H_0}$ و ${r_d}$ و دستیابی به مقدار عددی ثابت هابل، به اطلاعات کالیبراسیون اولیه نیاز است:
1. **ترکیب با سنتز هسته‌ای مهبانگ (BBN):** با استفاده از داده‌های فراوانی دوتریوم و تعیین چگالی باریونی (${\Omega_b h^2 = 0.02218 \pm 0.00055}$)، مقدار ثابت هابل به صورت زیر تعیین می‌شود:
$$
H_0 = (68.53 \pm 0.80) \, \text{km} \, \text{s}^{-1} \, \text{Mpc}^{-1}
$$
2. **ترکیب با مقیاس زاویه‌ای آکوستیک پلانک (${\theta_*}$):** ترکیب داده‌های DESI با قید هندسی ${\theta_*}$ تابش زمینه کیهانی، دقت را افزایش می‌دهد:
$$
H_0 = (68.52 \pm 0.62) \, \text{km} \, \text{s}^{-1} \, \text{Mpc}^{-1}
$$
3. **ترکیب کامل با پلانک و همگرایی گرانشی ACT:**
$$
\Omega_m = 0.3069 \pm 0.0050, \quad H_0 = (67.97 \pm 0.38) \, \text{km} \, \text{s}^{-1} \, \text{Mpc}^{-1}
$$
این مقادیر ثابت هابل همچنان تنش آماری حدود ۳.۷ تا ۴ سیگما را با اندازه‌گیری‌های نردبان فاصله مستقیم تلسکوپ فضایی هابل بر پایه قیفاووسی‌ها (پروژه SH0ES با ${H_0 = 73.04 \pm 1.04}$) حفظ می‌کنند.
# آزمون انحنای فضایی فضا-زمان ($\Omega_K$) #
با رها کردن فرض تخت بودن هندسه عالم و آزاد گذاشتن چگالی انحنا (${\Omega_K}$)، داده‌های DESI به دلیل پوشش چندین سرخ‌گرایی مختلف می‌توانند تحول آهنگ انبساط را به طور مستقل ردیابی کنند. قید حاصل از ترکیب DESI و تابش زمینه کیهانی تخت بودن فضا را با دقت بسیار بالایی تأیید می‌کند:
$$
\Omega_K = 0.0024 \pm 0.0016
$$
این نتیجه، تمایل اندکی را که در داده‌های اولیه پلانک به سمت انحنای بسته (${\Omega_K < 0}$) وجود داشت، با داده‌های مستقل مقیاس ساختارها برطرف می‌سازد.
# فیزیک انرژی تاریک: شواهد دینامیکی بودن معادله حالت #
در فیزیک نظری، انرژی تاریک با معادله حالت ${w(a) = P_{\text{de}} / \rho_{\text{de}}}$ توصیف می‌شود که برای ثابت کیهان‌شناسی همواره ${w = -1}$ است.
## ۱. مدل با معادله حالت ثابت ($w\text{CDM}$) ##
تحت فرض ثابت بودن ${w}$ در طول زمان، داده‌های DESI به تنهایی مقدار زیر را نتیجه می‌دهند:
$$
w = -0.99_{-0.13}^{+0.15}
$$
ترکیب DESI با پلانک و داده‌های ابرنواختر پانتئون‌پلاس به مقدار دقیق ${w = -0.997 \pm 0.025}$ می‌رسد که کاملاً با ثابت کیهان‌شناسی سازگار است.
## ۲. مدل دینامیکی با معادله حالت متغیر ($w_0 w_a\text{CDM}$) ##
برای آزمودن فرضیه تکامل زمانی انرژی تاریک، از فرمول‌بندی دوپارامتری شوالیه-پولارسکی-لیندر (CPL) استفاده می‌شود:
$$
w(a) = w_0 + w_a (1 - a)
$$
در این فرمول‌بندی، ${w_0}$ مقدار معادله حالت در زمان حال (${a = 1}$) و ${w_a}$ معرف نرخ تغییرات آن نسبت به زمان است (ثابت کیهان‌شناسی متناظر با ${w_0 = -1}$ و ${w_a = 0}$ است).
نتایج حاصل از ترکیب داده‌ها در این فضا نشان‌دهنده ترجیح سیستماتیک برای ناحیه ${w_0 > -1}$ و ${w_a < 0}$ است:
* **ترکیب DESI و پلانک:** اختلاف با مدل استاندارد در سطح ۲.۶ سیگما:
$$
w_0 = -0.45_{-0.21}^{+0.34}, \quad w_a = -1.79_{-1.00}^{+0.48}
$$
* **افزودن داده‌های ابرنواختری (SN Ia):**
* ترکیب با پانتئون‌پلاس (Pantheon+): ترجیح در سطح ۲.۵ سیگما با ${w_0 = -0.827 \pm 0.063}$ و ${w_a = -0.75_{-0.25}^{+0.29}}$
* ترکیب با یونیون۳ (Union3): ترجیح در سطح ۳.۵ سیگما با ${w_0 = -0.64 \pm 0.11}$ و ${w_a = -1.27_{-0.34}^{+0.40}}$
* ترکیب با داده‌های سال پنجم پیمایش انرژی تاریک (DESY5): ترجیح در سطح ۳.۹ سیگما با ${w_0 = -0.727 \pm 0.067}$ و ${w_a = -1.05_{-0.27}^{+0.31}}$
<تفسیر فیزیکی: رفتار ترجیح‌داده‌شده نشان می‌دهد که انرژی تاریک در سرخ‌گرایی‌های نزدیک (${z < 0.5}$) مانند یک میدان کوانتومی کوینتسنس با ${w > -1}$ عمل می‌کند و در سرخ‌گرایی‌های بالاتر به قلمرو انرژی فانتوم با ${w < -1}$ وارد می‌شود. هرچند آزمون‌های داده‌ورزی بیشتر برای رد خطاهای سیستماتیک احتمالی در داده‌های نورسنجی ابرنواخترها ضروری است، اما این نخستین نشانه رصدی قوی از نارسایی ثابت کیهان‌شناسی پایدار است.>
# قیدهای فیزیک ذرات: جرم نوترینوها و تعداد گونه‌های نسبیتی #
نوترینوهای پرجرم به دلیل پراکندگی سرعتی بالا، مانع از رمبش گرانشی ماده در مقیاس‌های کوچک‌تر از طول پویش آزاد خود می‌شوند و رشد ساختارهای کیهانی را سرکوب می‌کنند. علاوه بر این، در کیهان متأخر به ماده غیرنسبیتی تبدیل شده و بر آهنگ انبساط اثر می‌گذارند.
ترکیب داده‌های DESI و پلانک در مدل تخت $\Lambda\text{CDM}$ منجر به یکی از سخت‌گیرانه‌ترین کران‌های بالایی برای مجموع جرم نوترینوها شده است:
$$
\sum m_\nu < 0.072 \, \text{eV} \quad (95\% \text{ CL})
$$
این قید بالایی با در نظر گرفتن پیش‌شرط جرم‌های نامنفی حاصل شده است. در صورت اعمال قیدهای آزمایش‌های نوسان نوترینوهای زمینی:
* برای چینش جرم بهنجار (Normal Hierarchy با پیش‌شرط ${\sum m_\nu > 0.059\text{ eV}}$): قید بالا به ${\sum m_\nu < 0.113\text{ eV}}$ می‌رسد.
* برای چینش جرم معکوس (Inverted Hierarchy با پیش‌شرط ${\sum m_\nu > 0.10\text{ eV}}$): قید بالا ${\sum m_\nu < 0.145\text{ eV}}$ خواهد بود.
همچنین تعداد مؤثر گونه‌های تابشی نسبیتی اولیه برابر با ${N_{\text{eff}} = 3.10 \pm 0.17}$ به دست آمده است که سازگاری کاملی با پیش‌بینی مدل استاندارد فیزیک ذرات (${N_{\text{eff}} = 3.044}$) دارد.
# جایگاه علمی و چشم‌انداز آینده #
داده‌های سال اول DESI تنها حدود یک‌سوم از کل داده‌های نهایی این پروژه پنج‌ساله را تشکیل می‌دهند. تداوم این رصدها و تکمیل کاتالوگ با بیش از ۴۰ میلیون طیف نوری، مشخص خواهد کرد که آیا کشش مشاهده‌شده به سوی مدل‌های انرژی تاریک متغیر پایدار خواهد ماند و پارادایم کیهان‌شناسی استاندارد را دستخوش بازنگری بنیادی خواهد کرد، یا اینکه اثرات ناشناخته در روش‌های کالیبراسیون و داده‌برداری عامل این انحراف بوده‌اند.