Viscous Gubser flow with conserved charges to benchmark fluid simulations

1. समस्या विवरण (Problem Statement)

सापेक्षतावादी द्रवगतिकी (Relativistic hydrodynamics) भारी आयन टकरावों में उत्पन्न होने वाले क्वार्क-ग्लूऑन प्लाज्मा (QGP) के मॉडलिंग के लिए एक मानक उपकरण है। जबकि आदर्श और श्यान (viscous) द्रवगतिकी (जैसे कि ब्योर्केन और गुब्सर प्रवाह) के लिए बेंचमार्क परीक्षण मौजूद हैं, अधिकांश मौजूदा बेंचमार्क शून्य रासायनिक क्षमता (vanishing chemical potentials) (शून्य शुद्ध बैरियन संख्या, स्ट्रेंजनेस और विद्युत आवेश) मानकर चलते हैं।

हालाँकि, QCD चरण आरेख (phase diagram) को समझने, क्रिटिकल पॉइंट का पता लगाने और न्यूट्रॉन स्टार विलय (neutron star mergers) के मॉडलिंग के लिए ऐसे हाइड्रोडायनामिक कोड की आवश्यकता है जो परिमित रासायनिक क्षमता (finite chemical potentials) और संरक्षित आवेशों (conserved charges) (बैरियन संख्या $B$, स्ट्रेंजनेस $S$, विद्युत आवेश $Q$) को सटीक रूप से संभाल सकें। वर्तमान संख्यात्मक कोड (जैसे कि स्मूथड पार्टिकल हाइड्रोडायनामिक्स या ग्रिड-आधारित विधियों का उपयोग करने वाले) ऐसे कठोर यूनिट टेस्ट की कमी रखते हैं जो निम्नलिखित को संयोजित करते हों:

1. श्यान प्रभाव (shear viscosity)।
2. परिमित रासायनिक क्षमता।
3. संरक्षित आवेश धाराएं (conserved charge currents)।

ऐसे बेंचमार्क के बिना, उन कोडों की संख्यात्मक सटीकता को सत्यापित करना कठिन है जो बड़े रासायनिक क्षमता वाले सिस्टम का अनुकरण करने का प्रयास कर रहे हैं, विशेष रूप से समीकरण अवस्था (Equation of State - EoS), तापमान विकास और आवेश घनत्व के बीच युग्मन (coupling) के संबंध में।

2. कार्यप्रणाली (Methodology)

लेखकों ने संरक्षित आवेशों के साथ विस्कस गुब्सर फ्लो (Viscous Gubser Flow with Conserved Charges - VGCC) के लिए एक अर्ध-विश्लेषणात्मक समाधान (semi-analytical solution) विकसित किया है। इस कार्यप्रणाली में निम्नलिखित चरण शामिल हैं:

• समरूपता और निर्देशांक (Symmetry and Coordinates): वे गुब्सर प्रवाह की $SO(3)_q \otimes SU(1,1) \otimes Z_2$ समरूपता का उपयोग करते हैं। सिस्टम को मिल्न निर्देशांकों $(\tau, r, \phi, \eta)$ में वर्णित किया गया है और वेइल रीस्केलिंग (Weyl rescaling) और कोऑर्डिनेट मैपिंग के माध्यम से इसे डी सिटर स्पेस (de Sitter space) $(\rho, \theta, \phi, \eta)$ में परिवर्तित किया गया है। यह रूपांतरण द्रव वेग को एक स्थिर प्रोफाइल ($\hat{u}^\mu = (-1, 0, 0, 0)$) में बदलकर आंशिक अंतर समीकरणों (PDEs) को साधारण अंतर समीकरणों (ODEs) में बदल देता है।
• संरक्षण नियम (Conservation Laws):
  • ऊर्जा-संवेग (Energy-Momentum): इसे विस्कस हाइड्रोडायनामिक्स के इज़राइल-स्टियर्ट फॉर्मलिज्म (Israel-Stewart formalism) का उपयोग करके हल किया गया है, जिसमें शियर स्ट्रेस टेंसर $\pi^{\mu\nu}$ के लिए एक रिलैक्सेशन-प्रकार का समीकरण शामिल है।
  • संरक्षित धाराएं (Conserved Currents): लेखक सिद्ध करते हैं कि सटीक गुब्सर समरूपता में, डिफ्यूजन धाराएं (diffusion currents) स्वतः ही शून्य हो जाती हैं (परिशिष्ट A)। फलस्वरूप, आवेश घनत्वों ($n_B, n_S, n_Q$) का विकास ऊर्जा-संवेग टेंसर से अलग (decouple) हो जाता है और यह $D n_Y + n_Y \theta = 0$ के माध्यम से स्वतंत्र रूप से विकसित होता है।
• EoS के माध्यम से युग्मन (Coupling via Equation of State): हालांकि आवेश का विकास संवेग से अलग है, लेकिन तापमान ($T$) और रासायनिक क्षमता ($\mu_Y$) EoS के माध्यम से जुड़े हुए हैं। लेखक दो अलग-अलग कॉन्फॉर्मल EoS मॉडल का उपयोग करते हैं:
  1. EoS1: स्पष्ट $\mu$ निर्भरता के साथ एक मासलेस क्वार्क-ग्लूऑन प्लाज्मा मॉडल।
  2. EoS2: विशेष रूप से कोड बेंचमार्किंग के लिए उपयोग किया जाने वाला एक पैरामीटराइज्ड कॉन्फॉर्मल मॉडल।
     ये EoS मॉडल हाइड्रोडायनामिक चरों ($E, n_Y$) और ऊष्मप्रवैगिकी चरों ($T, \mu_Y$) के बीच मैपिंग की अनुमति देते हैं।
• फ्रीज-आउट हाइपरसरफेस (Freeze-out Hypersurface): लेखक एक स्थिर ऊर्जा घनत्व ($E = E_{FO}$) द्वारा परिभाषित फ्रीज-आउट हाइपरसरफेस को व्युत्पन्न करते हैं। वे इस हाइपरसरफेस के सामान्य वेक्टर्स (normal vectors) को विश्लेषणात्मक रूप से भी ज्ञात करते हैं, जो संख्यात्मक कोडों के परीक्षण के लिए एक नया मीट्रिक प्रदान करता है।

3. मुख्य योगदान (Key Contributions)

• प्रथम अर्ध-विश्लेषणात्मक VGCC समाधान: यह शोध पत्र संरक्षित आवेशों वाले सिस्टम के लिए तापमान, रासायनिक क्षमता और शियर स्ट्रेस को एक साथ विकसित करने वाले विस्कस गुब्सर फ्लो का पहला अर्ध-विश्लेषणात्मक समाधान प्रस्तुत करता है।
• शून्य डिफ्यूजन का प्रमाण: यह सिद्ध करना कि गुब्सर समरूपता डिफ्यूजन धाराओं को प्रतिबंधित करती है, जो आवेश घनत्वों के लिए सरल निरंतरता समीकरणों (continuity equations) को सरल बनाती है।
• संख्यात्मक कोड के लिए नया बेंचमार्क: लेखक फ्रीज-आउट हाइपरसरफेस प्रोफाइल और सामान्य वेक्टर्स प्रदान करते हैं, जो हाइड्रोडायनामिक कोड के लिए एक नया "यूनिट टेस्ट" है। यह न केवल बल्क विकास बल्कि जटिल हाइपरसरफेस को पुनर्गठित करने की कोड की क्षमता का भी परीक्षण करता है।
• ccake का बेंचमार्किंग: इस समाधान का उपयोग ccake को बेंचमार्क करने के लिए किया गया है, जो कि एक नया विकसित (2+1)D सापेक्षतावादी विस्कस हाइड्रोडायनामिक कोड है, जो BSQ संरक्षित आवेशों पर आधारित है।

4. परिणाम (Results)

• विकास गतिकी (Evolution Dynamics):
  • परिमित रासायनिक क्षमताओं का समावेश $\mu=0$ मामले की तुलना में सिस्टम को साम्यावस्था (equilibrium) से काफी दूर धकेल देता है।
  • तापमान प्रोफाइल (Temperature Profiles): $\mu/T$ के बड़े प्रारंभिक अनुपातों के लिए, तापमान प्रोफाइल कोर से दूर "शोल्डर्स" (shoulders) विकसित करते हैं। यह परिमित रासायनिक क्षमता पर बड़े विस्कस सुधारों से उत्पन्न एंट्रॉपी उत्पादन के कारण होता है, जिसे गिब्स-ड्यूहेम संबंध (Gibbs-Duhem relation) के माध्यम से लुप्त होते $n_Y \mu_Y$ पद की भरपाई करनी पड़ती है।
  • फ्रीज-आउट (Freeze-out): $\mu \neq 0$ पर, फ्रीज-आउट सतह न तो आइसोथर्मल (isothermal) है और न ही आइसेंट्रोपिक (isentropic)। कोर में एंट्रॉपी घनत्व पूंछ (tails) की तुलना में कम होता है, जो $\mu=0$ मामले के बिल्कुल विपरीत है।
• कोड बेंचमार्किंग (ccake बनाम अर्ध-विश्लेषणात्मक):
  • सहमति (Agreement): ऊर्जा घनत्व, संख्या घनत्व और द्रव वेग के लिए ccake के संख्यात्मक समाधान, अर्ध-विश्लेषणात्मक समाधान के साथ उत्कृष्ट सहमति दिखाते हैं।
  • विसंगतियां (Discrepancies): सबसे बड़ी विसंगतियां ऑफ-डायगोनल शियर स्ट्रेस घटकों ($\pi_{xy}$) में और देर के समय ($\tau > 1.6$ fm/c) रेडियल क्षेत्र $1-2$ fm में होती हैं। इन विसंगतियों का कारण सिस्टम का साम्यावस्था से बहुत दूर जाना है, जिससे कौज़ैलिटी (causality) का उल्लंघन होता है और SPH कणों में ऋणात्मक एंट्रॉपी उत्पन्न होती है, जिसके कारण सिमुलेशन क्रैश हो जाता है।
  • फ्रीज-आउट सतह: सिमुलेशन के अंत में क्रैश होने के बावजूद, ccake अर्ध-विश्लेषणात्मक रूप से व्युत्पन्न फ्रीज-आउट हाइपरसरफेस और सामान्य वेक्टर्स को सटीक रूप से पुनरुत्पादित करता है। यह पुष्टि करता है कि हाइड्रोडायनामिक शासन के भीतर कोड वैध है।

5. महत्व (Significance)

• भविष्य के सिमुलेशन का सत्यापन: यह कार्य RHIC के बीम एनर्जी स्कैन (BES) कार्यक्रम और भविष्य के प्रयोगों (जैसे FAIR, NICA) के लिए हाइड्रोडायनामिक कोडों को मान्य करने के लिए एक महत्वपूर्ण "गोल्ड स्टैंडर्ड" प्रदान करता है, जहाँ परिमित बैरियन घनत्व मुख्य केंद्र है।
• गैर-साम्यावस्था प्रभावों की समझ: परिणाम यह उजागर करते हैं कि परिमित रासायनिक क्षमताएं QGP के स्पेसटाइम विकास को कैसे बदलती हैं, विशेष रूप से गैर-आइसेंट्रोपिक फ्रीज-आउट को प्रेरित करके और तापमान प्रोफाइल को उन तरीकों से संशोधित करके जो शून्य-घनत्व सिमुलेशन में नहीं देखे जाते।
• क्रिटिकल पॉइंट की खोज के लिए उपकरण: संरक्षित आवेशों वाले सिस्टमों के सटीक अनुकरण को सक्षम करके, यह बेंचमार्क QCD क्रिटिकल पॉइंट की खोज का समर्थन करता है, जहाँ श्यानता (viscosity), रासायनिक क्षमता और EoS के बीच का परस्पर क्रिया अत्यंत महत्वपूर्ण है।
• ओपन सोर्स उपलब्धता: लेखक समाधान और बेंचमार्क डेटा का पायथन कार्यान्वयन प्रदान करते हैं, जिससे नए हाइड्रोडायनामिक मॉडलों के परीक्षण के लिए समुदाय द्वारा तत्काल अपनाना सुलभ हो सके।

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