Thermal evolution of dark matter and gravitational-wave production in the early universe from a symplectic glueball model

तकनीकी सारांश: एक सिम्प्लेक्टिक ग्लूबॉल मॉडल से प्रारंभिक ब्रह्मांड में डार्क मैटर का तापीय विकास और गुरुत्वाकर्षण तरंग उत्पादन

समस्या और प्रेरणा
डार्क मैटर की प्रकृति भौतिकी की एक मौलिक अनसुलझी समस्या बनी हुई है। जबकि प्रयोगात्मक साक्ष्य इसके अस्तित्व का पुरजोर समर्थन करते हैं, इसकी कण पहचान अज्ञात है। ऐसे सैद्धांतिक मॉडल जो डार्क मैटर को एक दृढ़ता से युग्मित (strongly coupled), गैर-अबेलियन गेज सिद्धांत (क्वांटम क्रोमोडायनामिक्स (QCD) के समान) के सम्मिश्र अवस्था (composite state) के रूप में प्रस्तावित करते हैं, वे आकर्षक हैं क्योंकि वे बिना फाइन-ट्यूनिंग के स्वाभाविक रूप से द्रव्यमान पैमाने (mass scales) उत्पन्न कर सकते हैं और स्थिर, तटस्थ, तथा स्व-परस्पर क्रिया करने वाले उम्मीदवार प्रदान कर सकते हैं। पिछला साहित्य $SU(N)$गेज समूहों पर आधारित ऐसे "डार्क ग्लूबॉल" मॉडलों का व्यापक अध्ययन कर चुका है। हालांकि, अन्य गेज समूहों, विशेष रूप से सिम्प्लेक्टिक समूहों$Sp(N)$पर आधारित सिद्धांतों का व्यवहार कम खोजा गया है। यह कार्य कॉम्पैक्ट सिम्प्लेक्टिक समूह$Sp(2)$(जो$Sp(4, \mathbb{C})$का वास्तविक रूप है और$U(4)$ के साथ प्रतिच्छेद करता है) पर आधारित एक डार्क मैटर मॉडल की व्यवहार्यता की जांच करता है। इसका प्राथमिक उद्देश्य यह निर्धारित करना है कि क्या यह सिद्धांत एक प्रथम-क्रम (first-order) कन्फाइनमेंट/डीकन्फाइनमेंट चरण संक्रमण (phase transition) से गुजरता है, जो एक पता लगाने योग्य स्टोकेस्टिक गुरुत्वाकर्षण-तरंग पृष्ठभूमि (stochastic gravitational-wave background) उत्पन्न कर सकता है, और इसके समीकरण (equation of state - EoS) को गैर-परतदार (non-perturbatively) रूप से अभिलक्षणित करना है।

कार्यप्रणाली
लेखक एक यूक्लिडियन, आइसोट्रोपिक हाइपरक्यूबिक लैटिस (hypercubic lattice) पर $Sp(2)$ यांग-मिल्स सिद्धांत का गैर-परतदार अध्ययन करते हैं।

• सिमुलेशन सेटअप: सिद्धांत को $Sp(2)$के डिफाइनिंग रिप्रेजेंटेशन में लिंक वेरिएबल्स$U_\mu(x)$ के साथ विल्सन एक्शन (Wilson action) द्वारा परिभाषित किया गया है। सिमुलेशन परिमित तापमान $T = 1/(a_{lat} N_\tau)$ पर किए जाते हैं, जहाँ $a_{lat}$ लैटिस स्पेसिंग है और $N_\tau$ यूक्लिडियन समय दिशा में विस्तार है।
• एल्गोरिदम: लेखक एक मार्कोव-चेन मोंटे कार्लो (MCMC) एल्गोरिदम का उपयोग करते हैं जो $Sp(2)$मैट्रिसेस के भीतर एम्बेडेड$SU(2)$ सबग्रुप्स में स्थानीय हीट-बाथ (heat-bath) और ओवर-रिलैक्सेशन अपडेट को जोड़ता है। एर्गोडिसिटी सुनिश्चित करने और ऑटोकोरिलेशन को कम करने के लिए, अपडेट्स को $Sp(2)$मैट्रिसेस के भीतर एम्बेडेड$SU(2)$ सबग्रुप्स में लागू किया जाता है।
• स्केल सेटिंग: लैटिस स्केल को विभिन्न $N_\tau$ मानों के लिए डीकन्फाइनमेंट ट्रांजिशन के लिए क्रिटिकल कपलिंग $\beta_c$ निर्धारित करके गैर-परतदार रूप से सेट किया जाता है। यह पॉलीकोव लूप (Polyakov loop) के चौथे-क्रम के बाइंडर कमुलेंट (Binder cumulant) के विश्लेषण और फाइनाइट-साइज़ स्केलिंग के माध्यम से प्राप्त किया जाता है। परिणामी $\beta$ बनाम $a_{lat}$ संबंध को व्यापक कपलिंग रेंज को कवर करने के लिए एक सुधारे गए लैटिस स्कीम में टू-लूप पर्टर्बेटिव अभिव्यक्ति के साथ इंटरपोलेट और मैच किया जाता है।
• समीकरण (EoS): थर्मोडायनामिक मात्राओं (दबाव $p$, ऊर्जा घनत्व $\epsilon$, एंट्रॉपी घनत्व $s$, और ऊर्जा-संवेग टेंसर का ट्रेस $\Delta$) की गणना "इंटीग्रल मेथड" का उपयोग करके की जाती है। इसमें परिमित तापमान पर औसत प्लाकेट ट्रेस और शून्य तापमान पर प्लाकेट ट्रेस के बीच के अंतर को कपलिंग $\beta$ के सापेक्ष एकीकृत किया जाता है।
• कंटिन्यूम एक्सट्रपलेशन: परिणाम कई लैटिस स्पेसिंग ($N_\tau = 5, 6, 7$) पर प्राप्त किए जाते हैं और $1/N_\tau^2$ के फलन के रूप में डेटा को फिट करके कंटिन्यूम लिमिट ($a_{lat} \to 0$) तक एक्सट्रपलेट किया जाता है। फाइनाइट-वॉल्यूम प्रभावों को बड़े स्थानिक आयतन ($N_s/N_\tau \gtrsim 4$) का उपयोग करके दबा दिया जाता है।

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

1. चरण संक्रमण की प्रकृति: यह अध्ययन पुष्टि करता है कि $Sp(2)$यांग-मिल्स सिद्धांत एक परिमित क्रिटिकल तापमान$T_c$ पर एक प्रथम-क्रम डीकन्फाइनमेंट चरण संक्रमण से गुजरता है। यह बाइंडर कमुलेंट कर्व्स के विभिन्न वॉल्यूम के लिए क्रॉसिंग और एक गैर-शून्य लेटेंट हीट (latent heat) की उपस्थिति से सिद्ध होता है।
2. थर्मोडायनामिक मात्राएँ: लेखक $0.7 \le T/T_c \le 2.75$ की सीमा में $Sp(2)$ सिद्धांत के EoS के पहले कंटिन्यूम-एक्सट्रपलेटेड परिणाम प्रदान करते हैं।
   • ऊर्जा-संवेग टेंसर का ट्रेस $\Delta/T^4$ $T_c$ के पास एक पीक प्रदर्शित करता है, जो प्रथम-क्रम संक्रमण की विशेषता है।
   • दबाव, ऊर्जा घनत्व और एंट्रॉपी घनत्व की गणना की गई है और वे स्टेफ़न-बोल्ट्ज़मैन सीमा की ओर धीरे-धीरे बढ़ते हुए पाए गए हैं, जो डीकन्फाइंड चरण में भी महत्वपूर्ण अंतःक्रियाओं को इंगित करता है।
   • ध्वनि की गति $c_s^2$ व्युत्पन्न की गई है, जो $SU(3)$ यांग-मिल्स सिद्धांत के गुणात्मक और मात्रात्मक रूप से समान व्यवहार को दर्शाती है।
3. लेटेंट हीट: एक परिमित लेटेंट हीट $L_h$ निर्धारित की गई है। कंटिन्यूम-एक्सट्रपलेटेड मान $L_h/T_c^4 \simeq 1.69(12)$ पाया गया है। यह मान $SU(3)$($\simeq 1.4$) की तुलना में थोड़ा अधिक है, लेकिन गेज समूह के आयाम ($Sp(2)$के लिए 10 बनाम$SU(3)$ के लिए 8) के आधार पर स्केलिंग अपेक्षाओं के अनुरूप है।
4. कन्फाइंड फेज मॉडलिंग: निम्न-तापमान कन्फाइंड चरण में, EoS को गैर-पारस्परिक क्रिया करने वाले विशाल ग्लूबॉल्स (massive glueballs) के गैस द्वारा अच्छी तरह से वर्णित किया गया है। सबसे हल्का ग्लूबॉल एक स्केलर ($J^{PC}=0^{++}$) है जिसका द्रव्यमान $m \approx 3.58 \sqrt{\sigma}$ है। परिणाम $T_c$ के पास स्टेट्स के घनत्व (density of states) में घातांकीय वृद्धि दिखाते हैं, जो एक हेगेड्रोन-जैसे (Hagedorn-like) स्पेक्ट्रम के अनुरूप है।
5. गुरुत्वाकर्षण-तरंग उत्पादन: लेखक प्रथम-क्रम चरण संक्रमण द्वारा उत्पन्न गुरुत्वाकर्षण-तरंग (GW) स्पेक्ट्रम का अनुमान लगाते हैं।
   • गणना की गई लेटेंट हीट और एक अनुमानित इंटरफेस टेंशन का उपयोग करते हुए, वे ट्रांजिशन स्ट्रेंथ पैरामीटर $\alpha \simeq 0.62(4)$ निकालते हैं, जो संक्रमण को "मध्यवर्ती" (intermediate) से "मजबूत" (strong) के रूप में वर्गीकृत करता है।
   • GW स्पेक्ट्रम को बबल कोलिजन (bubble collisions), साउंड वेव्स और हाइड्रोडायनामिक टर्बुलेंस के योगदान को ध्यान में रखते हुए मॉडल किया गया है।
   • अनुमान है कि प्रमुख योगदान साउंड वेव्स से आएगा, जिसकी पीक फ्रीक्वेंसी मिली-हर्ट्ज़ (mHz) रेंज में है। यह संकेत भविष्य के अंतरिक्ष-आधारित डिटेक्टरों जैसे LISA, DECIGO, और BBO की संवेदनशीलता सीमा के भीतर होने की संभावना रखता है, यदि चरण संक्रमण सैकड़ों GeV के तापमान पर होता है।

डार्क मैटर व्याख्या और महत्व
पेपर सबसे हल्के स्केलर ग्लूबॉल ($0^{++}$) को एक डार्क मैटर उम्मीदवार के रूप में व्याख्यायित करता है।

• स्थिरता: $N \ge 3$ वाले $SU(N)$सिद्धांतों के विपरीत,$Sp(2)$समूह नकारात्मक चार्ज-कंजुगेशन ($C=-1$) वाली अवस्थाओं का समर्थन नहीं करता है। फलस्वरूप, उन$SU(N)$मॉडलों में प्रस्तावित तंत्र जहाँ एक$C=-1$अवस्था को डिस्क्रीट सिमेट्री द्वारा स्थिर किया जाता है (जिससे$C=+1$अवस्थाएं क्षय होती हैं जबकि$C=-1$डार्क मैटर बनी रहती है), यहाँ **लागू नहीं होता** है।$Sp(2)$ मॉडल में स्थिरता क्षय के लिए उपलब्ध हल्के राज्यों की अनुपस्थिति और संभावित रूप से उच्च-आयामी ऑपरेटरों (जैसे ग्रेविटॉन) के माध्यम से दमित क्षय चैनलों पर निर्भर करती है।
• रेलिक अबंडेंस (Relic Abundance): लेखक डार्क सेक्टर के तापीय विकास पर चर्चा करते हैं, जिसमें "कैनिबलिज्म" (cannibalism) परिदृश्य (3-to-2 स्व-विलय/self-annihilation) शामिल है जो अवशेष घनत्व (relic density) को बदल सकता है। वे डार्क मैटर की प्रचुरता के लिए ऑर्डर-ऑफ-मैग्निट्यूड अनुमान प्रदान करते हैं, यह नोट करते हुए कि परिणाम रीहीटिंग तापमान और दृश्य (visible) एवं डार्क सेक्टर के तापमान अनुपात के प्रति संवेदनशील हैं।
• महत्व: पेपर दावा करता है कि $Sp(2)$ सिद्धांत एक एकल आयामी पैरामीटर पर निर्भर एक व्यवहार्य, "किफायती" (economic) डार्क मैटर मॉडल है। इसका प्राथमिक महत्व $Sp(2)$ गेज सिद्धांत के लिए EoS और लेटेंट हीट का पहला कठोर, गैर-परतदार लैटिस निर्धारण प्रदान करने में निहित है। यह स्थापित करता है कि ऐसे सिद्धांत मजबूत प्रथम-क्रम संक्रमणों का समर्थन कर सकते हैं जो अवलोकन योग्य गुरुत्वाकर्षण तरंगें उत्पन्न करने में सक्षम हैं। लेखक इस बात पर जोर देते हैं कि हालांकि उनके EoS के लैटिस परिणाम मजबूत हैं, GW स्पेक्ट्रा और डार्क मैटर अवशेष घनत्व के लिए भविष्यवाणियां गैर-इक्विलिब्रियम मापदंडों (जैसे बबल वॉल वेलोसिटी और न्यूक्लिएशन रेट) के अर्ध-मात्रात्मक अनुमानों पर निर्भर करती हैं, जो उनके वर्तमान इक्विलिब्रियम लैटिस सेटअप की सीमा से बाहर हैं।

सीमाएं और भविष्य का कार्य
लेखक स्पष्ट रूप से बताते हैं कि उनके GW और डार्क मैटर प्रचुरता के अनुमान बेंचमार्क उदाहरण हैं न कि प्रथम-सिद्धांत मात्रात्मक भविष्यवाणियां, क्योंकि वे उन गतिशील मात्राओं (बबल वॉल वेलोसिटी, इंटरफेस टेंशन, न्यूक्लिएशन रेट) के अनुमानों पर निर्भर करते हैं जिनके लिए रियल-टाइम सिमुलेशन या प्रभावी क्षेत्र सिद्धांत (effective field theory) दृष्टिकोण की आवश्यकता होती है जो उनके वर्तमान इक्विलिब्रियम अध्ययन के दायरे से बाहर हैं। भविष्य के कार्य में स्केल सेटिंग को परिष्कृत करना, तापमान सीमा का विस्तार करना और इंटरफेस टेंशन तथा बबल न्यूक्लिएशन रेट के लिए समर्पित गैर-परतदार गणना करना शामिल हो सकता है। इसके अतिरिक्त, $Sp(2)$में$C=-1$अवस्थाओं की अनुपस्थिति कुछ स्थिरता तंत्रों को रोकती है जो$SU(N)$ मॉडलों में उपलब्ध हैं, जो इस मॉडल की फेनोमेनोलॉजी को अलग करता है।