Cylindrical cosmological simulations with StePS

तकनीकी सारांश: StePS के साथ बेलनाकार ब्रह्मांडीय सिमुलेशन (Cylindrical Cosmological Simulations)

समस्या विवरण
मानक ब्रह्मांडीय $N$-body सिमुलेशन आमतौर पर एक आवधिक घन बॉक्स (periodic cubic box) के भीतर पदार्थ का विकास करते हैं, जो एक तीन-टोरस ($T^3 = S^1 \times S^1 \times S^1$) टोपोलॉजी लागू करता है। हालांकि यह संख्यात्मक रूप से सुविधाजनक है और सांख्यिकीय समरूपता (statistical homogeneity) को बनाए रखता है, यह टोपोलॉजी सांख्यिकीय समदैशिकता (statistical isotropy) को तोड़ देती है क्योंकि घन ज्यामिति केवल अष्टफलकीय समूह ($O_h$) के तहत अपरिवर्तित रहती है, न कि निरंतर घूर्णन समूह $SO(3)$ के तहत। फलस्वरूप, आवधिक बक्सों में कोणीय संवेग (angular momentum) सामान्यतः संरक्षित नहीं रहता है, और घन समरूपता गुरुत्वाकर्षण क्षेत्र में सूक्ष्म अनिसोट्रॉपी (anisotropies) उत्पन्न कर सकती है, जो बड़े पैमाने पर क्लस्टरिंग सांख्यिकी को पूर्वाग्रहित कर सकती है, विशेष रूप से तब जब सिमुलेशन वॉल्यूम रुचि के पैमानों की तुलना में काफी बड़ा न हो। इसके अतिरिक्त, मानक आवधिक सीमाएं उन प्रणालियों के सिमुलेशन के लिए उपयुक्त नहीं हैं जिनमें अंतर्निहित बेलनाकार समरूपताएं होती हैं, जैसे कि कॉस्मिक फिलामेंट्स या विशिष्ट अनिसोट्रोपिक ब्रह्मांडीय मॉडल।

कार्यप्रणाली (Methodology)
लेखक एक संकुचित सिमुलेशन ढांचे को पेश करते हैं जो $S^1 \times \mathbb{R}^2$ (स्लैब) टोपोलॉजिकल मैनिफोल्ड पर आधारित है। इस ज्यामिति में, सिमुलेशन डोमेन एक एकल अक्ष (z-अक्ष, जिसका आवर्त $L_z$ है) के अनुदिश आवधिक है, लेकिन रेडियल दिशा ($\mathbb{R}^2$) में खुला है।

• ज्यामितीय संकुचन (Geometric Compactification): लेखक गैर-अनंत $\mathbb{R}^2$ प्लेन को संकुचित करने के लिए रेडियल दिशा में एक 'इनवर्स स्टीरियोग्राफिक प्रोजेक्शन' का उपयोग करते हैं। यह अनंत रेडियल निर्देशांक को एक परिमित डोमेन में मैप करता है, जिससे एक "ज़ूम-इन" रणनीति संभव होती है जहाँ केंद्रीय अक्ष के पास स्थानिक और द्रव्यमान रिज़ॉल्यूशन उच्च होता है और बाहरी सीमा की ओर बढ़ते हुए धीरे-धीरे कम होता जाता है। यह बाहरी मुक्त सीमा पर किनारों के आर्टिफैक्ट्स (edge artifacts) को दबा देता है।
• गति के समीकरण (Equations of Motion): न्यूटनियन गति के समीकरणों को कोमूविंग कोऑर्डिनेट्स (comoving coordinates) में हल किया जाता है। गुरुत्वाकर्षण बल की गणना $S^1 \times \mathbb{R}^2$ टोपोलॉजी के अनुसार z-अक्ष के अनुदिश आवधिक छवियों (periodic images) को जोड़कर की जाती है। प्रत्यक्ष योग (direct summation) के धीमे अभिसरण को संभालने के लिए, लेखक बेलनाकार ज्यामिति के अनुकूल एक 'एवल्ड समेशन' (Ewald summation) विधि लागू करते हैं, जो बल को रियल-स्पेस और रेसिप्रोकल-स्पेस घटकों में विभाजित करती है।
• बल गणना एल्गोरिदम (Force Calculation Algorithms): StePS (STEreographically Projected cosmological Simulations) ढांचे के भीतर कार्यान्वयन दो बल गणना विधियों का समर्थन करता है:
  1. प्रत्यक्ष $O(N^2)$: यह GPU पर कुशलतापूर्वक मैप होता है, जो कम कण संख्या या छोटे बेलनाकार वॉल्यूम के लिए उपयुक्त है।
  2. ऑक्ट्री $O(N \log N)$: CPU क्लस्टर्स के लिए कार्यान्वित। ऑक्ट्री ग्रिड की क्यूबिक प्रकृति से उत्पन्न होने वाली व्यवस्थित त्रुटियों को कम करने के लिए, कोड प्रत्येक टाइमस्टेप पर डोमेन में रैंडम शिफ्ट और रोटेशन लागू करता है। एक महत्वपूर्ण अंशांकन चरण रेडियल बल त्रुटियों को ठीक करता है ताकि स्पूरियस रेडियल कोलैप्स या विस्तार को रोका जा सके, जिससे बैकग्राउंड डेंसिटी की स्थिरता सुनिश्चित होती है।
• प्रारंभिक स्थितियाँ (Initial Conditions): प्रारंभिक स्थितियाँ एक "ग्लास" (glass) कण वितरण का उपयोग करके उत्पन्न की जाती हैं, जिसे $S^1 \times \mathbb{R}^2$ टोपोलॉजी में प्रतिकारक गुरुत्वाकर्षण (repulsive gravity) के माध्यम से बनाया गया है। विक्षोभ (perturbations) को ज़ेल्डोविच सन्निकटन (Zel'dovich approximation) या द्वितीय-क्रम लैग्रेंजियन पर्टरबेशन थ्योरी (2LPT) का उपयोग करके अंकित किया जाता है, जो बेलनाकार ज्यामिति और ज़ूम-इन सेटअप के गैर-समान सैंपलिंग के लिए अनुकूलित है। गैर-तुच्छ विंडो फंक्शन (window function) को संभालने के लिए पावर स्पेक्ट्रम का अनुमान लगाने के लिए Feldman–Kaiser–Peacock (FKP) एस्टीमेटर का उपयोग किया जाता है।

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

1. नवीन टोपोलॉजी कार्यान्वयन: यह शोध पत्र $S^1 \times \mathbb{R}^2$ मैनिफोल्ड में ब्रह्मांडीय $N$-body सिमुलेशन के पहले कार्यान्वयन को प्रस्तुत करता है, जो StePS कोड को इसकी पिछली $T^3$ और गोलाकार $\mathbb{R}^3$ क्षमताओं से आगे बढ़ाता है।
2. संरक्षण नियम (Conservation Laws): लेखक प्रदर्शित करते हैं कि इस टोपोलॉजी में, लैग्रेंजियन सिलेंडर अक्ष के अनुदिश अनुवाद (translation) और ट्रांसवर्स प्लेन में घूर्णन ($SO(2)$) के तहत अपरिवर्तित रहता है। फलस्वरूप, z-अक्ष के अनुदिश रैखिक संवेग और उसी अक्ष के अनुदिश कोणीय संवेग संरक्षित रहते हैं, जो$T^3$ टोपोलॉजी की समरूपता-भंग (symmetry-breaking) समस्याओं को संबोधित करता है।
3. एल्गोरिद्मिक अनुकूलन: यह कार्य बल गणना (बेलनाकार एवल्ड समेशन), प्रारंभिक स्थिति निर्माण (बेलनाकार ग्लास), और हेलो फाइंडिंग (एक कस्टम गोलाकार-ओवरडेंसिटी फाइंडर जो $S^1$ आवधिकता और रेडियल रिज़ॉल्यूशन भिन्नता के लिए अनुकूलित है) के लिए आवश्यक संशोधनों का विवरण देता है।
4. सत्यापन (Validation): विधि को मान्य करने के लिए एक पूर्ण $\Lambda$CDM सिमुलेशन चलाया गया है, जो मानक ब्रह्मांडीय अपेक्षाओं के अनुरूप एक कॉस्मिक वेब के निर्माण को प्रदर्शित करता है।

परिणाम

• सिमुलेशन निष्पादन: एक सिमुलेशन $L_z = 1.0$ Gpc की ऊंचाई और $R_{sim} = 500$ Mpc की त्रिज्या वाले सिलेंडर में $2.4 \times 10^7$ कणों के साथ चलाया गया, जिसमें Planck 2018 $\Lambda$CDM पैरामीटर्स का उपयोग किया गया।
• संरचना निर्माण: सिमुलेशन सफलतापूर्वक गैर-रैखिक कॉस्मिक वेब (रिक्त स्थान, दीवारें, फिलामेंट्स और हेलो) को पुनरुत्पादित करता है। ज़ूम-इन रणनीति केंद्रीय अक्ष के पास कम द्रव्यमान वाले हेलो के उच्च-फिडेलिटी रिज़ॉल्यूशन की अनुमति देती है, जबकि बड़े वॉल्यूम अधिक रेडियस पर अधिक विशाल हेलो को हल करते हैं।
• पावर स्पेक्ट्रा: अनुमानित मैटर पावर स्पेक्ट्रा बड़े पैमाने पर रैखिक CAMB भविष्यवाणियों के साथ और छोटे पैमाने पर गैर-रैखिक Mead et al. (2021) "halofit" एमुलेटर के साथ अच्छी तरह से मेल खाते हैं, जो संरचना विकास को पकड़ने में विधि की सटीकता की पुष्टि करते है।
• हेलो गुण: लेखकों ने $z=0$ पर $\sim 1.1 \times 10^5$ डार्क मैटर हेलो की पहचान की और प्रत्येक के लिए द्रव्यमान, एकाग्रता (concentration) और कोणीय संवेग सहित 68 पैरामीटरों की गणना की।
• अनिसोट्रॉपी परीक्षण: गैलेक्सी स्पिन अनिसोट्रॉपी के हालिया दावों से प्रेरित होकर, लेखकों ने आवधिक $z$-अक्ष के सापेक्ष हेलो स्पिन अक्षों के वितरण का परीक्षण किया। कोलमोगोरोव-स्मिरनोव (Kolmogorov-Smirnov) परीक्षणों का उपयोग करते हुए, उन्होंने पाया कि $L_z = 1.0$ Gpc के लिए, स्पिन वितरण सांख्यिकीय समदैशिकता के अनुरूप है। वे निष्कर्ष निकालते हैं कि यह टोपोलॉजी, इस पैमाने पर, गैलेक्सी स्पिन में रिपोर्ट की गई विषमताओं की व्याख्या स्वाभाविक रूप से नहीं करती है।

महत्व और दावे
यह शोध दावा करता है कि $S^1 \times \mathbb{R}^2$ टोपोलॉजी मानक क्यूबिक सिमुलेशन के लिए एक विशिष्ट विकल्प प्रदान करती है, विशेष रूप से बेलनाकार समरूपता वाले सिस्टम (जैसे फिलामेंट्स) के अध्ययन के लिए या उन मॉडलों के लिए जहाँ एक विशिष्ट अक्ष के अनुदिश कोणीय संवेग का संरक्षण भौतिक रूप से प्रासंगिक है। लेखक इस बात पर जोर देते हैं कि हालांकि यह विधि विशिष्ट व्यवस्थित त्रुटियां (systematics) पेश करती है—जैसे कि रेडियली भिन्न रिज़ॉल्यूशन और चयन प्रभाव—यह एक स्व-संगत ढांचा प्रदान करती है जो $T^3$ बॉक्स की क्यूबिक अनिसोट्रॉपी से मुक्त है।

लेखक अपने दावों में विनम्र हैं, उनका कहना है कि वे सामान्य प्रयोजन के ब्रह्मांड विज्ञान के लिए स्थापित क्यूबिक $T^3$ सिमुलेशन को बदलने की उम्मीद नहीं करते हैं। इसके बजाय, यह विशिष्ट वैज्ञानिक प्रश्नों के लिए एक पूरक उपकरण के रूप में कार्य करता है जहाँ क्यूबिक समरूपता एक सीमित कारक है। वे इस बात पर भी प्रकाश डालते हैं कि टोपोलॉजी का एनिसोट्रोपिक बल क्षेत्र आर्टिफैक्ट्स उत्पन्न कर सकता है यदि आवधिक लंबाई बहुत कम ($L_z \lesssim 200$ Mpc) हो, जो कि भविष्य के विस्तृत लक्षण वर्णन का विषय है। इन ज्यामितिओं के और अधिक अन्वेषण को सुगम बनाने के लिए कोड (StePS v2) ओपन-सोर्स के रूप में जारी किया गया है।