Simulation of PBH formation in a matter-dominated universe

तकनीकी सारांश: पदार्थ-प्रधान ब्रह्मांड में PBH निर्माण का सिमुलेशन

समस्या विवरण
प्राइमॉर्डियल ब्लैक होल (PBH) का निर्माण ब्रह्मांड के प्रमुख घटक के समीकरण अवस्था (equation of state) के प्रति अत्यधिक संवेदनशील है। विकिरण-प्रधान युग के दौरान PBH निर्माण अपेक्षाकृत अच्छी तरह से समझा गया है, लेकिन एक प्रारंभिक पदार्थ-प्रधान युग के दौरान की प्रक्रिया कम निश्चित है। एक पदार्थ-प्रधान युग में, दबाव प्रवणता (pressure gradients) नगण्य होती है, जिससे अति-सघन क्षेत्र अधिक आसानी से ढह सकते हैं। पूर्व विश्लेणात्मक अनुमान अक्सर विPerturbative सन्निकटन या ज़ेल्डोविच सन्निकटन (Zel'dovich approximation) पर निर्भर करते हैं, जो अत्यधिक गैर-रैखिक शासन (highly nonlinear regime) में लागू नहीं हो सकते हैं। इस शासन की जांच करने में एक प्रमुख चुनौती गैर-गोलाकार समरूपता को संख्यात्मक सिमुलेशन में लागू करना है, विशेष रूप से डस्ट फ्लूइड (dust fluids) के मामले में जहाँ शेल-क्रॉसिंग सिंगुलैरिटी (shell-crossing singularities) ब्लैक होल क्षितिज के पूरी तरह से हल होने से पहले ही संख्यात्मक क्रैश का कारण बनती हैं।

कार्यप्रणाली
लेखक एक विस्तारमान पदार्थ-प्रधान पृष्ठभूमि में PBH निर्माण का अनुकरण करने के लिए बाउमगार्ट-शापिरो-शिबाता-नाकामुरा (BSSN) फॉर्मलिज्म का उपयोग करके पूर्णतः गैर-रैखिक संख्यात्मक सापेक्षता (numerical relativity) का उपयोग करते हैं।

• प्रारंभिक स्थितियाँ: सिमुलेशन वक्रता विचलन $\zeta$ के एक कार्यात्मक रूप का उपयोग करता है जिसमें दीर्घवृत्ताकार (ellipticity) शामिल है, जो त्रिकैक्सियल (triaxial) विन्यास बनाता है। प्रारंभिक डेटा एक डस्ट फ्लूइड के लिए दीर्घ-तरंग प्रवणता विस्तार समाधानों का उपयोग करके निर्मित किया गया है, जो एक आयाम पैरामीटर $\mu$ और दीर्घवृत्ताकार पैरामीटर ($e, p$) द्वारा अभिलक्षित है।
• गेज और निर्देशांक: एक नवीन समय स्लाइसिंग स्थिति को अपनाया गया है, जो क्षितिज प्रवेश के करीब आने पर ई-फोल्डिंग संख्या-आधारित गेज से एक सिंगुलैरिटी-अवॉयडिंग "1+log" शैली की स्लाइसिंग में परिवर्तित हो जाती है। सिमुलेशन बॉक्स के केंद्रीय क्षेत्र को बड़ा करने के लिए एक गैर-कार्तिय निर्देशांक प्रणाली का उपयोग किया गया है, जिससे वहां रिज़ॉल्यूशन में सुधार होता है जहाँ पतन (collapse) होता है।
• पदार्थ विवरण: दो अलग-अलग पदार्थ मॉडलों का परीक्षण किया गया है:
  1. डस्ट फ्लूइड: एक दबावहीन तरल विवरण। लेखक नोट करते हैं कि यह विवरण अनिवार्य रूप से शेल-क्रॉसिंग सिंगुलैरिटी की ओर ले जाता है जहाँ तरल घनत्व अनंत हो जाता है, जिससे संख्यात्मक कोड क्रैश हो जाता है।
  2. कोलिजनलेस पार्टिकल्स (Collisionless Particles): फ्लूइड सन्निकटन के टूटने से बचने के लिए, लेखक पूर्ण संख्यात्मक सापेक्षता के भीतर एक पार्टिकल-इन-सेल (PIC) दृष्टिकोण लागू करते हैं। वांछित घनत्व विषमता को पुनरुत्पादित करने के लिए कणों को द्रव्यमान दिया जाता है, और उनकी गति जियोडेसिक समीकरणों द्वारा नियंत्रित होती है। यह सिमुलेशन को उस बिंदु के आगे जारी रखने की अनुमति देता है जहाँ तरल तत्व क्रॉस करते हैं।

मुख्य परिणाम

• गोलाकार समरूपता सत्यापन: कोड को गोलाकार सममान डस्ट पतन के लिए लेमाइट्र-टोलमैन-बॉन्डी (LTB) समाधान के विरुद्ध सत्यापित किया गया था। संख्यात्मक परिणामों ने क्षितिज निर्माण तक के विश्लेणात्मक समाधान को सफलतापूर्वक पुनरुत्पादित किया और स्थानीय नग्न सिंगुलैरिटी (PBH निर्माण) और वैश्विक नग्न सिंगुलैरिटी के बीच के दहलीज (threshold) की सही पहचान की।
• डस्ट फ्लूइड सिमुलेशन (गैर-गोलाकार):
  • त्रिकैक्सियल पतन के लिए, डस्ट फ्लूइड सिमुलेशन शेल-क्रॉसिंग सिंगुलैरिटी के कारण क्रैश हो जाते हैं जब तक कि आयाम $\mu$ काफी बड़ा न हो।
  • क्षितिज निर्माण को देखने के लिए आवश्यक दहलीज आयाम, ज़ेल्डोविच सन्निकटन और हूप कंजेक्चर (hoop conjecture) पर आधारित पिछले विश्लेणात्मक अनुमानों द्वारा अनुमानित मान से लगभग दोगुना पाया गया है।
  • यह सुझाव देता है कि फ्लूइड सन्निकटन के वैध दायरे के भीतर, PBH निर्माण पहले की तुलना में कठिन है।
• कोलिजनलेस पार्टिकल सिमुलेशन:
  • कोलिजनलेस कणों के साथ समान प्रणाली को मॉडल करते समय, सिमुलेशन शेल क्रॉसिंग पर क्रैश नहीं होते हैं।
  • ब्लैक होल निर्माण (apparent horizon गठन द्वारा इंगित) $\mu \approx 0.045$ (दीर्घवृत्ताकार $e=0.2$ के लिए) जैसे निम्न प्रारंभिक आयामों के लिए भी देखा गया है।
  • यह दहलीज पिछले विश्लेणात्मक अनुमानों की तुलना में एक क्रम (order of magnitude) छोटा है और डस्ट फ्लूइड सिमुलेशन से प्राप्त दहलीज से काफी कम है।
  • हालांकि, पार्टिकल सिमुलेशन कणों के इंटरसेक्शन के बाद हैमिल्टोनियन बाधाओं (Hamiltonian constraint) का महत्वपूर्ण उल्लंघन प्रदर्शित करते हैं, जो यह दर्शाता है कि जबकि पतन आगे बढ़ता है, उच्च-घनत्व वाले विकास की मात्रात्मक सटीकता पार्टिकल विवरण और ज्यामिति के बीच बेमेल होने के कारण सीमित है।

महत्ता और दावे
यह शोध पत्र दावा करता है कि पदार्थ क्षेत्र के भौतिक विवरण पर निर्भर करते हुए, पदार्थ-प्रधान युग के दौरान PBH निर्माण मानदंड में एक बड़ा सैद्धांतिक अनिश्चितता मौजूद है:

1. फ्लूड सन्निकटन सीमा: यदि डस्ट फ्लूइड विवरण में शेल क्रॉसिंग से जुड़ी सिंगुलैरिटी के निर्माण को एक भौतिक दमन तंत्र (physical suppression mechanism) के रूप में व्याख्या किया जाता है (या यदि फ्लूइड सन्निकटन पतन से पहले टूट जाता है), तो PBH निर्माण के लिए दहलीज पिछले विश्लेणात्मक अनुमानों की तुलना में काफी अधिक है, जिससे PBH निर्माण अधिक कठिन हो जाता है।
2. पार्टिकल डायनेमिक्स: यदि कोलिजनलेस पार्टिकल विवरण को सिंगल-स्ट्रीम फ्लूइड सन्निकटन के वैध विस्तार के रूप में स्वीकार किया जाता है, तो PBH निर्माण उम्मीद से बहुत आसान (कम दहलीज) हो जाता है।
3. सैद्धांतिक अनिश्चितता: फ्लूइड और पार्टिकल परिणामों के बीच का अंतर यह उजागर करता है कि मानक विश्लेणात्मक मानदंड (ज़ेल्डोविच सन्निकटन) त्रिकैक्सियल पतन की गैर-रैखिक गतिशीलता को सटीक रूप से कैप्चर नहीं कर सकते हैं। लेखक निष्कर्ष निकालते हैं कि दहलीज मान उच्च घनत्व पर पदार्थ क्षेत्र के उपचार के प्रति अत्यधिक संवेदनशील है।

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