Characteristic Decomposition for Relativistic Numerical Simulations: I. Hydrodynamics

कल्पना कीजिए कि आप यह भविष्यवाणी करने की कोशिश कर रहे हैं कि एक तूफान किसी शहर के ऊपर से कैसे गुजरेगा। इसे सटीक रूप से करने के लिए, आप केवल एक स्थान पर हवा को नहीं देखते; आपको यह समझने की आवश्यकता है कि वायु दाब, तापमान और गति की विभिन्न "लहरें" कैसे परस्पर क्रिया करती हैं और माध्यम से यात्रा करती हैं। भौतिकी की दुनिया में, विशेष रूप से जब चीजें प्रकाश की गति के करीब चलती हैं या ब्लैक होल के तीव्र गुरुत्वाकर्षण के प्रभाव में होती हैं, तो वैज्ञानिक इन अराजक प्रवाहों को व्यक्तिगत, अनुमानित लहरों में तोड़ने के लिए जटिल गणित का उपयोग करते हैं। इस प्रक्रिया को "कैरक्टेरिस्टिक डीकंपोजिशन" (characteristic decomposition) कहा जाता है। इसे एक शेफ द्वारा एक जटिल स्टू (stew) को उसके व्यक्तिगत अवयवों—नमक, काली मिर्च, गाजर और बीफ—में अलग करने जैसा समझें—ताकि वह यह समझ सके कि प्रत्येक घटक वास्तव में कैसे व्यवहार करता है। इस अलगाव के बिना, ब्रह्मांडीय घटनाओं (जैसे न्यूट्रॉन सितारों का टकराव) के कंप्यूटर सिमुलेशन ऐसा है जैसे आप आंखों पर पट्टी बांधकर स्टू बनाने की कोशिश कर रहे हों; आपको एक परिणाम तो मिल सकता है, लेकिन वह सटीक नहीं होगा, और आप सबसे नाटकीय विस्फोटों को भी मिस कर सकते हैं।

दशकों से, वैज्ञानिक सामान्य तरल पदार्थों और यहाँ तक कि बिना चुंबकीय क्षेत्रों वाले, सापेक्षतावादी गति (प्रकाश की गति के करीब) पर चलने वाले तरल पदार्थों के लिए यह "अवयव पृथक्करण" करने में सक्षम रहे हैं। हालाँकि, जब आप इसमें शक्तिशाली चुंबकीय क्षेत्र जोड़ देते हैं—जिसे जनरल रिलेटिविस्टिक मैग्नेटोहाइड्रोडायनामिक्स (GRMHD) के रूप में जाना जाता है—तो यह रेसिपी एक उलझे हुए ढेर में बदल जाती है। गणित इतना जटिल हो जाता है कि सबसे सटीक कंप्यूटर उपकरण, जो इस बात पर निर्भर करते हैं कि ये लहरें बिल्कुल कैसे यात्रा करती हैं, उनका उपयोग नहीं किया जा सकता। सॉल ए. ट्यूकोल्स्की (Saul A. Teukolsky) का यह शोध पत्र एक चतुर नया गणितीय तरीका पेश करता है जिसे "क्वासी-इनवर्टिबल ट्रांसफॉर्मेशन" (quasi-invertible transformations) कहा जाता है। पूरे उलझे हुए गुत्थी को एक साथ सुलझाने के बजाय, लेखक यह दिखाते हैं कि कैसे एक सरल, ज्ञात समस्या (जहाँ तरल अपने सापेक्ष स्थिर है) से शुरू करके, इसे वास्तविक सिमुलेशन में उपयोग किए जाने वाले जटिल, चलते हुए चुंबकीय संस्करण में सुचारू रूप से बदला जा सकता है। यह शोध पत्र सापेक्षतावादी हाइड्रोडायनामिक्स (बिना चुंबकीय क्षेत्र वाले तरल पदार्थ) पर इस विधि का सफलतापूर्वक प्रदर्शन करता है, यह सिद्ध करते हुए कि यह तरीका काम करता है और ज्ञात परिणामों को बहुत सरल तरीके से पुनः प्राप्त करता है। हालाँकि, चुंबकीय क्षेत्रों वाला पूर्ण समाधान दूसरे शोध पत्र के लिए सुरक्षित रखा गया है, यह कार्य एक आवश्यक आधार तैयार करता है, जो यह सुझाव देता है कि हम अंततः अत्यधिक सटीक रूप से सबसे हिंसक ब्रह्मांडीय झटकों (shocks) को संभाल सकते हैं।

कॉस्मिक किचन: हमें लहरों को अलग करने की आवश्यकता क्यों है

यह समझने के लिए कि यह शोध पत्र क्यों महत्वपूर्ण है, आइए हम उस "किचन" की ओर वापस चलें जहाँ ये सिमुलेशन होते हैं। आधुनिक भौतिकी में, हम तारों के जन्म से लेकर ब्लैक होल के टकराव तक सब कुछ सिम्युलेट करते हैं। ये घटनाएँ उन समीकरणों द्वारा संचालित होती हैं जो बताते हैं कि पदार्थ और ऊर्जा कैसे प्रवाहित होते हैं। जब चीजें तेज चलती हैं या गुरुत्वाकर्षण मजबूत होता है, तो ये समीकरण "हाइपरबोलिक" (hyperbolic) हो जाते हैं। सरल शब्दों में, इसका अर्थ है कि जानकारी सिस्टम के माध्यम से विशिष्ट, स्वतंत्र लहरों में यात्रा करती है, ठीक वैसे ही जैसे तालाब में पत्थर फेंकने पर फैलती लहरों की तरह।

इनमें से प्रत्येक लहर की अपनी गति और दिशा होती है। किसी टक्कर या शॉक वेव को सटीक रूप से सिम्युलेट करने के लिए, कंप्यूटर को यह जानने की आवश्यकता होती है कि ये लहरें वास्तव में कैसे व्यवहार करती हैं। यदि कंप्यूटर यह अनुमान लगाने की कोशिश करता है कि वे कैसे परस्पर क्रिया करती हैं (बिना उनकी व्यक्तिगत गति को जाने), तो यह ऐसी त्रुटियाँ पैदा करता है जो सिमुलेशन को विफल कर सकती हैं या बेतुके परिणाम दे सकती हैं। "कैरक्टेरिस्टिक डीकंपोजिशन" इन व्यक्तिगत लहरों और उनकी गति को खोजने की गणितीय प्रक्रिया है। यह एक कार को यह अनुमान लगाकर चलाने की तुलना में, कि सड़क कहाँ जा रही है, बनाम एक सटीक जीपीएस मैप होने के बीच का अंतर है।

लंबे समय तक, वैज्ञानिकों के पास सरल तरल पदार्थों (हाइड्रोडायनामिक्स) और यहाँ तक कि बिना चुंबकीय क्षेत्रों वाले तेज गति से चलने वाले तरल पदार्थों के लिए एक आदर्श जीपीएस मैप था। लेकिन जैसे ही आप इसमें चुंबकीय क्षेत्र (मैग्नेटोहाइड्रडायनामिक्स या MHD) जोड़ते हैं और जनरल रिलेटिविटी के चरम गुरुत्वाकर्षण को शामिल करते हैं, वह मैप गायब हो जाता है। समीकरण इतने जटिल हो जाते हैं कि कोई भी आवश्यक "लहरों" (eigenvalues और eigenvectors) को नहीं खोज पाता, जो सबसे सटीक कंप्यूटर एल्गोरिदम चलाने के लिए आवश्यक हैं। इसने वैज्ञानिकों को "अनुमानित" तरीकों का उपयोग करने के लिए मजबूर किया—जैसे धुंधले मैप के साथ गाड़ी चलाना—जो कुछ चीजों के लिए ठीक काम करता है लेकिन तब विफल हो जाता है जब आपको यह भविष्यवाणी करने की आवश्यकता होती है कि एक तारा कब फटेगा या एक ब्लैक होल पदार्थ के एक टुकड़े को कब निगलेगा।

नया तरीका: गणित के लिए एक जादुई दर्पण

यह शोध पत्र, जो दो भागों की श्रृंखला का पहला भाग है, इस पहेली को हल करने का एक नया तरीका पेश करता है। लेखक, सॉल ए. ट्यूकोल्स्की, "क्वासी-इनवर्टिबल ट्रांसफॉर्मेशन" पर आधारित एक विधि प्रस्तावित करते हैं। यह सुनने में कठिन लग सकता है, लेकिन विचार आश्चर्यजनक रूप से सरल है।

कल्पित कीजिए कि आप एक भूलभुलैया (maze) को हल करने की कोशिश कर रहे हैं। भूलभुलैया कागज के एक टुकड़े पर बनी है जो मुड़ा हुआ, चिपचिपा और स्याही से ढका हुआ है (यह जटिल GRMHD समीकरण है)। इसे सीधे हल करने की कोशिश करना एक दुःस्वप्न है। लेकिन, आप जानते हैं कि यदि आप कागज को सीधा कर सकें और स्याही को धो सकें, तो भूलभुलैया बिल्कुल एक सरल, सीधे पथ की तरह दिखेगी जिसे आपने हजार बार हल किया है (यह "कोमूविंग फ्रेम" है, जहाँ तरल स्थिर है)।

समस्या यह है कि आप कागज को केवल सीधा नहीं कर सकते; गणित पूर्ण "अनडू" (undo) बटन की अनुमति नहीं देता क्योंकि चर (variables) सीमित होते हैं (जैसे कि एक तरल की गति और घनत्व आपस में जुड़े होते हैं)। आमतौर पर, यदि आप एक रूपांतरण को पूरी तरह से उल्टा नहीं कर सकते, तो आप फंस जाते हैं। लेकिन ट्यूकोल्स्की ने एक "क्वासी-इनवर्टिबल" ट्रिक की खोज की। यह एक जादुई दर्पण होने जैसा है जो छवि को पूरी तरह से प्रतिबिंबित नहीं करता है, लेकिन इसे इतना प्रतिबिंबित करता है कि आप भूलभुलैया के माध्यम से चल सकें।

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

हमने क्या सीखा और आगे क्या है

यह शोध पत्र स्पष्ट रूप से दिखाता है कि यह नई विधि काम करती है। यह सिद्ध करता है कि "कोमूविंग फ्रेम" (जहाँ तरल स्थिर है) से शुरू करके और इन क्वासी-इनवर्टिबल ट्रांसफॉर्मेशन का उपयोग करके, हम तरल पदार्थों के लिए कैरक्टेरिस्टिक डीकंपोजिशन को इस तरह से प्राप्त कर सकते हैं जो सटीक और गणनात्मक रूप से अनुकूल (computationally friendly) दोनों है। लेखक एक नया परिणाम भी पेश करते हैं: ऐसे तरल पदार्थों के लिए एक डीकंपोजिशन जहाँ रासायनिक संरचना बदलती है (जैसे परमाणु प्रतिक्रियाओं में), जो यह दर्शाता है कि यह विधि इन अतिरिक्त जटिलताओं को संभालने के लिए पर्याप्त लचीली है।

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

संक्षेप में, यह शोध पत्र यह दावा नहीं करता कि इसने पूरे ब्रह्मांडीय चुंबकीय तूफानों के रहस्य को सुलझा लिया है। इसके बजाय, इसने चाबी और ताला बनाया है। इसने सिद्ध कर दिया है कि "क्वासी-इनवर्टिबल" दरवाजा मौजूद है और यह सरल तरल पदार्थों के लिए आसानी से खुलता है। अब जब दरवाजा खुल गया है, तो दूसरा शोध पत्र ब्रह्मांड की सबसे हिंसक, चुंबकीय घटनाओं के लिए इन सिमुलेशन की पूरी शक्ति को अनलॉक करने के लिए इसके माध्यम से जाएगा। फिलहाल, हम जानते हैं कि अधिक सटीक, शॉक-कैप्चरिंग सिमुलेशन का मार्ग अब अनसुलझे गणित द्वारा अवरुद्ध नहीं है।