The Density Profile of Dynamical Halos

तकनीकी सारांश: गतिशील प्रभामंडलों का घनत्व प्रोफाइल (The Density Profile of Dynamical Halos)

समस्या विवरण
डार्क मैटर प्रभामंडलों (halos) के मूलभूत गुणों, विशेष रूप से उनके घनत्व प्रोफाइल को पारंपरिक रूप से एक संदर्भ त्रिज्या जैसे $R_{200m}$ के भीतर कुल द्रव्यमान वितरण पर विश्लेषणात्मक फिट (जैसे, NFW, Einasto) द्वारा अभिलक्षित किया जाता है। हालाँकि, यह दृष्टिकोण दो अलग-अलग गतिशील आबादी को मिला देता है: वे कण जो गुरुत्वाकर्षण से बंधे (bound) हैं और प्रभामंडल के भीतर कक्षाओं में घूम रहे हैं, और वे कण जो अंदर आ रहे हैं (infalling) लेकिन अभी तक बंधे नहीं हैं। हाल के प्रस्तावों में एक "गतिशील प्रभामंडल" (dynamical halo) को सख्ती से घूमने वाले कणों के संग्रह के रूप में परिभाषित करने का सुझाव दिया गया है। जबकि यह परिभाषा आंतरिक संरचना को निरंतर अभिवृद्धि (accretion) के साथ मिश्रित होने से बचाती है, घूमने वाले प्रोफाइलों के संरचनात्मक गुणों (औसत और प्रकीर्णन/scatter) का पूर्णतः लक्षण वर्णन नहीं किया गया था। विशेष रूप से, यह स्पष्ट नहीं था कि निश्चित घूमने वाले द्रव्यमान पर कौन से संरचनात्मक स्वतंत्रता स्तर (degrees of freedom) शेष रहते हैं, अंतर्निहित प्रकीर्णन का परिमाण क्या है, और क्या निर्माण इतिहास के साथ सहसंबंध वास्तविक रूप से घूमने वाले घटक के गुण हैं या मिश्रित आबादी के कारण उत्पन्न कृत्रिम प्रभाव (artifacts) हैं।

कार्यप्रणाली (Methodology)
लेखकों ने एक कोल्ड डार्क मैटर (CDM) सिमुलेशन ($1024^3$ कणों के साथ $1 \, h^{-1}$ Gpc बॉक्स) का उपयोग किया, जिसे GADGET-2 द्वारा चलाया गया था। उन्होंने Rockstar हेलो फाइंडर का उपयोग करके एक डायनामिकल हेलो कैटलॉग तैयार किया, फिर कणों को उनके अभिवृद्धि समय ($a_{acc}$) और समय-औसत रेडियल वेग के आधार पर "घूमने वाले" (orbiting) या "अंदर आने वाले" (infalling) के रूप में टैग करने के लिए García et al. के एल्गोरिदम को लागू किया।

• घूमने वाला (Orbiting) परिभाषा: वे कण जिनका अभिवृद्धि समय प्रारंभिक है और जिनका समय-औसत रेडियल वेग शून्य के करीब है, उन्हें 'घूमने वाला' टैग किया गया है; वे कण जिनका अभिवृद्धि समय विलंबित है और जिनका रेडियल वेग ऋणात्मक है, उन्हें 'अंदर आने वाला' (infalling) टैग किया गया है।
• द्रव्यमान बिनिंग (Mass Binning): प्रभामंडलों को नौ संकीर्ण द्रव्यमान बिनों ($\log M_{orb} \in [13.40, 15.00]$) में विभाजित किया गया।
• पैरामीट्रिक मॉडल: घूमने वाले घनत्व प्रोफाइलों को Salazar et al. के मॉडल का उपयोग करके फिट किया गया, जिसे इस प्रकार परिभाषित किया गया है:
  $$\rho_{orb}(x|M_{orb}) = A \left(\frac{x}{\epsilon}\right)^{-\alpha(x|M_{orb})} \exp\left(-\frac{x^2}{2}\right)$$
  जहाँ $x \equiv r/r_h$, $r_h$ प्रभामंडल त्रिज्या (एक्सपोनेंशियल ट्रंकेशन स्केल) है, और ढलान $\alpha(x)$ शून्य से एक स्पर्शोन्मुख मान $\alpha_\infty$ तक परिवर्तित होती है। पैरामीटर $\epsilon$ को 0.037 पर स्थिर रखा गया था।
• फिटिंग रणनीति: व्यक्तिगत प्रभामंडलों को सर्वोत्तम-फिट $r_h$ और $\alpha_\infty$ निर्धारित करने के लिए फिट किया गया। इसके बाद लेखकों ने $\alpha_\infty$ और $a_{RF}$ (सापेक्ष निर्माण समय चर) के साथ सहसंबंध की जांच की।

मुख्य योगदान और परिणाम

1. एकल प्रभावी स्वतंत्रता स्तर (Single Effective Degree of Freedom):
   यह अध्ययन प्रदर्शित करता है कि एक निश्चित घूमने वाले द्रव्यमान पर, प्रोफाइल के दो मुक्त पैरामीटर ($r_h$ और $\alpha_\infty$) मजबूती से सह-संबंधित हैं। अधिक विस्तृत प्रभामंडल (बड़ा $r_h$) अधिक तीव्र आंतरिक ढलान (बड़ा $\alpha_\infty$) रखते हैं। यह सहसंबंध एक व्यक्तिगत गतिशील प्रभामंडल के प्रोफाइल को एक एकल मुक्त पैरामीटर, प्रभामंडल त्रिज्या $r_h$, द्वारा अभिलक्षित करने की अनुमति देता है। संबंध इस प्रकार है:
   $$\Delta\alpha_\infty = 0.05 + \ln R$$
   जहाँ $R = r_h / r_{h,st}(M_{orb})$ है और $\Delta\alpha_\infty$ स्टैक्ड ढलान से विचलन है।

2. हेलो त्रिज्या में प्रकीर्णन (Scatter in Halo Radius):

   • एक निश्चित घूमने वाले द्रव्यमान पर, हेलो त्रिज्या का वितरण लगभग लॉग-नॉर्मल (log-normal) है। $\ln r_h$ में अंतर्निहित प्रकीर्णन $\sigma_{\ln r_h | M_{orb}} \approx 0.16$ (16%) है।
   • लेखकों ने अभिवृद्धि इतिहास पर इस प्रकीर्णन की निर्भरता की जांच की। उन्होंने कण अभिवृद्धि समय वितरण के 60वें पर्सेंटाइल के आधार पर एक सापेक्ष निर्माण समय $a_{RF}$ को परिभाषित किया।
   • $r_h$ और $a_{RF}$ के बीच एक सहसंबंध मौजूद है: देर से बनने वाले प्रभामंडल अधिक सघन (छोटा $r_h$) होते हैं।
   • इस सहसंबंध को ध्यान में रखने से निश्चित द्रव्यमान और निर्माण समय पर हेलो त्रिज्या में प्रकीर्णन घटकर $\sigma_{\ln r_h | M_{orb}, a_{RF}} \approx 0.11$ (11%) रह जाता है। लेखक नोट करते हैं कि यह कमी मामूली है, जो यह सुझाव देती है कि केवल सापेक्ष निर्माण समय ही प्रकीर्णन को पूरी तरह से समझाने के लिए पर्याप्त नहीं है।

3. अभिवृद्धि दर (Accretion Rate) के साथ तुलना:
   लेखकों ने हेलो त्रिज्या और सापेक्ष निर्माण समय के बीच सहसंबंध की तुलना हेलो त्रिज्या और प्रति इकाई द्रव्यमान अभिवृद्धि दर के बीच सहसंबंध के विरुद्ध की। उन्होंने पाया कि रुझान गुणात्मक रूप से समान हैं, जिसमें अभिवृद्धि दर सहसंबंध थोड़े कड़े प्रतिबंध दिखाता है लेकिन इसमें अधिक स्पष्ट गैर-गौसियन पूंछ (non-Gaussian tails) हैं। उन्होंने चेतावनी दी कि वे इस विशिष्ट विश्लेषण में गतिशील प्रभामंडलों के लिए अभिवृद्धि दरों का स्व-सुसंगत अनुमान नहीं लगा सके।

4. मॉडल प्रदर्शन:
   घूमने वाले प्रोफाइलों के लिए पैरामीट्रिक फिट, कुल द्रव्यमान प्रोफाइल पर लागू NFW फिट के तुल्य अवशेष (residuals) प्रदान करते हैं (औसत अवशेष $\approx -0.02$, प्रकीर्णन $\sim 0.15$)। यह Salazar et al. मॉडल की वैधता की पुष्टि करता है।

महत्व और दावे
यह शोध पत्र दावा करता है कि यह गतिशील रूप से परिभाषित प्रभामंडलों के लिए अभिवृद्धि-इतिहास निर्भरता के प्रथम लक्षण वर्णन को प्रदान करता है। घूमने वाले घटक को अलग करके, लेखक तर्क देते हैं कि मापे गए सहसंबंध (जैसे, संरचना और निर्माण समय के बीच सहसंबंध) बंधे हुए तंत्र के वास्तविक गुण हैं, न कि अंदर आने वाली सामग्री के मिश्रण के कारण उत्पन्न होने वाले परिभाषा संबंधी कृत्रिम प्रभाव।

इसका प्राथमिक महत्व एक गतिशील प्रभामंडल के द्रव्यमान और निर्माण समय को उसके संरचनात्मक गुणों ($r_h$ और $\alpha_\infty$) से "पढ़ने" की क्षमता में निहित है। मानक गोलाकार ओवरडेंसिटी (spherical overdensity) प्रभामंडलों के विपरीत, जहाँ त्रिज्या परिभाषात्मक रूप से संलग्न द्रव्यमान से जुड़ी होती है, गतिशील प्रभामंडल यह प्रकट करते हैं कि एक प्रभामंडल का विस्तार और संरचना कक्षीय गतिशीलता (orbital dynamics) से कैसे उत्पन्न होती है। यह अंशांकन (calibration) फॉरवर्ड मॉडलिंग के लिए एक आवश्यक आधार रेखा के रूप में कार्य करता है, जो अधिक सटीक गतिशील अपघटन (dynamical decompositions) को सक्षम बनाता है और लेंसिंग एवं क्लस्टरिंग डेटा से अनुमानित द्रव्यमान वितरण की व्याख्या में सुधार करता है। लेखक इस बात पर जोर देते हैं कि हालांकि उनका मॉडल कुल द्रव्यमान वितरण के लिए अनुभवजन्य फिट (empirical fits) का स्थान नहीं लेता है, लेकिन यह अंतर्नि้น गतिशील संरचना को समझने के लिए आवश्यक है।