Quantum Simulation of Stokes Flow via Schrödingerisation and Artificial Compressibility

तकनीकी सारांश: श्रोडिंगरइजेशन (Schrödingerisation) और कृत्रिम संपीड्यता (Artificial Compressibility) के माध्यम से स्टोक्स प्रवाह का क्वांटम सिमुलेशन

समस्या विवरण
असंपीड्य स्टोक्स प्रवाह (incompressible Stokes flow) का अनुकरण करना माइक्रोफ्लुइडिक्स, निम्न-रेनॉल्ड्स-संख्या हाइड्रोडायनामिक्स और पोरस मीडिया ट्रांसपोर्ट जैसे अनुप्रयोगों के लिए महत्वपूर्ण है। हालांकि, शास्त्रीय संख्यात्मक विधियां एक मौलिक बाधा का सामना करती हैं: स्टोक्स प्रणाली की सैडल-पॉइंट संरचना (saddle-point structure) के लिए असतत स्तर पर इनसंपीड्यता प्रतिबंध ($\nabla \cdot u = 0$) को सटीक रूप से संतुष्ट करना आवश्यक होता है। यह युग्मन (coupling) बड़े, दुर्बल-कंडीशन वाले रैखिक प्रणालियों (ill-conditioned linear systems) को जन्म देता है जिन्हें हल करना कम्प्यूटेशनल रूप से महंगा होता है, और जैसे-जैसे समस्या की विमा (dimensionality) बढ़ती है, संसाधन आवश्यकताएं अत्यधिक बढ़ जाती हैं। जबकि कृत्रिम संपीड्यता जैसी शास्त्रीय नियमितीकरण तकनीकें मौजूद हैं, वे उच्च-विमीय सेटिंग्स में कम्प्यूटेशनल रूप से निषेधात्मक बनी रहती हैं क्योंकि आवश्यक संसाधनों का घातीय स्केलिंग (exponential scaling) होता है।

कार्यप्रणाली
लेखक एक एकीकृत क्वांटम एल्गोरिद्मिक ढांचे का प्रस्ताव करते हैं जो कृत्रिम संपीड्यता को श्रोडिंगरइजेशन तकनीक के साथ एकीकृत करके इन चुनौतियों का समाधान करता है।

1. कृत्रिम संपीड्यता नियमितीकरण (Artificial Compressibility Regularization): सख्त इनसंपीड्यता स्थिति को एक छद्म-दबाव विकास समीकरण (pseudo-pressure evolution equation) पेश करके शिथिल किया जाता है। विशेष रूप से, निरंतरता समीकरण $\nabla \cdot u = 0$ को $\nabla \cdot u = \varepsilon p$ द्वारा प्रतिस्थापित किया जाता है, जहाँ $\varepsilon$ एक छोटा कृत्रिम संपीड्यता पैरामीटर है। यह मूल सैडल-पॉइंट समस्या को मोमेंटम समीकरण और एक परवलयिक दबाव विकास समीकरण (parabolic pressure evolution equation) से जुड़ी एक प्रारंभिक-मान समस्या (initial-value problem) में बदल देता है। लेखक त्रुटि अनुमान प्रदान करते हैं जो वास्तविक असंपीड्य समाधान के प्रति रैखिक अभिसरण ($O(\varepsilon)$) प्रदर्शित करते हैं।

2. स्टैगर्ड-ग्रिड विविक्तकरण (Staggered-Grid Discretization): संख्यात्मक स्थिरता और उचित वेग-दबाव युग्मन सुनिश्चित करने के लिए, लेखक स्टैगर्ड-ग्रिड फइन वॉल्यूम विधि का उपयोग करते हैं। यह दृष्टिकोण वेग घटकों और दबाव को अलग-अलग नोडल बिंदुओं पर रखता है, जिससे कोलोकेटेड ग्रिड (collocated grids) में सामान्यतः दिखने वाले अवांछित दोलनों (spurious oscillations) से बचा जा सकता है।

3. श्रोडिंगरइजेशन मैपिंग (Schrödingerisation Mapping): मुख्य नवाचार गैर-हर्मिटियन (non-Hermitian), विसरणी (dissipative) आंशिक अंतर समीकरणों (PDEs) को एक यूनिटरी क्वांटम सिस्टम पर मैप करना है। एक "वार्प्ड फेज ट्रांसफॉर्मेशन" ($w = e^{-q}u$) और सहायक चर $q$ पर फूरियर ट्रांसफॉर्म का उपयोग करके, लेखक मूल गतिकी को एक उच्च-विमीय हिल्बर्ट स्पेस में समाहित करते हैं। यह एम्बेडिंग विसरणी प्रणाली को एक श्रोडिंगर-प्रकार के समीकरण ($\partial_t \hat{v} = i H_{Stokes} \hat{v}$) में परिवर्तित कर देता है जिसमें एक हर्मिटियन हैमिल्टनियन $H_{Stokes}$ होता है।

4. क्वांटम सर्किट निर्माण: लेखक समय विकास ऑपरेटर $U_{Stokes}(\tau) = \exp(i H_{Stokes} \tau)$ के सिमुलेशन के लिए एक स्पष्ट क्वांटम सर्किट डिजाइन करते हैं।

   • हैमिल्टनियन को प्राथमिक घटकों में प्रथम-क्रम ली-ट्रोटर-सुजुकी (Lie-Trotter-Suzuki) फॉर्मूला का उपयोग करके विघटित किया जाता है।
   • विभेदक ऑपरेटरों (शिफ्ट और लैपलेसियन) को बाइनरी एनकोडिंग और शिफ्ट ऑपरेटरों ($S_+, S_-$) का उपयोग करके दर्शाया जाता है।
   • सर्किट सहायक $q$-स्पेस और मोमेंटम $\eta$-स्पेस के बीच स्विच करने के लिए क्वांटम फूरियर ट्रांसफॉर्म (QFT) का उपयोग करता है, जिससे गैर-यूनिटरी मूल गतिकी को यूनिटरी क्वांटम विकास के माध्यम से सिम्युलेट किया जा सके।
   • समाधान को इनवर्स QFT और सहायक चर सब्सपेस जहाँ $q > 0$ है, पर प्रोजेक्टिव मेजरमेंट के माध्यम से प्राप्त किया जाता है।

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

• एकीकृत ढांचा: यह शोध पत्र एक सुसंगत रूपरेखा प्रदान करने वाला पहला स्पष्ट क्वांटम एल्गोरिदम प्रस्तुत करता है जो कृत्रिम संपीड्यता के साथ श्रोडिंगरइजेशन को जोड़ता है, जो युग्मित वेग-दबाव प्रणाली के लिए एक सुव्यवस्थित निरूपण प्रदान करता है।
• स्पष्ट सर्किट डिजाइन: लेखक एक विस्तृत क्वांटम सर्किट कार्यान्वयन का वर्णन करते हैं, जिसमें हैमिल्टनियन का मौलिक गेट्स (CNOT, Hadamard, Phase, और multi-controlled rotation gates) में अपघटन शामिल है, और एक कठोर गेट काउंट विश्लेषण प्रदान करते हैं।
• जटिलता विश्लेषण: एक कठोर जटिलता विश्लेषण एक सैद्धांतिक क्वांटम लाभ स्थापित करता है। लेखक प्रदर्शित करते हैं कि उच्च-विमीय समस्याओं के लिए, प्रस्तावित विधि शास्त्रीय विधियों की तुलना में समस्या की विमा में घातीय गति (exponential speedup) प्राप्त करती है। विशेष रूप से, जबकि शास्त्रीय जटिलता कुछ मान्यताओं के तहत $O(\delta^{-(d+4)})$ के रूप में स्केल करती है, क्वांटम दृष्टिकोण $d$ पर इस घातीय निर्भरता से बचता है।
• संख्यात्मक सत्यापन: Qiskit प्लेटफॉर्म पर किए गए संख्यात्मक सिमुलेशन के माध्यम से एल्गोरिदम की वैधता और स्केलेबिलिटी की पुष्टि की गई है, जिसमें बेंचमार्क समस्याओं पर सटीकता परीक्षण और पैरामीटर अभिसरण विश्लेषण शामिल हैं।

परिणाम
Qiskit पर संख्यात्मक सिमुलेशन एल्गोरिदम की सटीकता की पुष्टि करते हैं। लेखकों ने एक विश्लेषणात्मक समाधान वाले युग्मित स्टोक्स सिस्टम पर इस पद्धति का परीक्षण किया, यह प्रदर्शित करते हुए कि जैसे-जैसे विविक्तकरण पैरामीटर ($n_x, n_q$) बढ़ते हैं और समय अंतराल ($dt$) तथा कृत्रिम संपीड्यता पैरामीटर ($\varepsilon$) घटते हैं, पॉइंटवाइज एरर कम होता जाता है। परिणाम विभेदक ऑपरेटर विविक्तकरण की शुद्धता और समग्र क्वांटम सर्किट डिजाइन की स्थिरता की पुष्टि करते हैं।

महत्व और दावे
यह शोध पत्र दावा करता है कि यह शास्त्रीय स्टोक्स फ्लो सिमुलेशन में निहित विमा संबंधी बाधा (dimensionality bottleneck) को दूर करके कम्प्यूटेशनल फ्लूइड डायनेमिक्स में एक महत्वपूर्ण सफलता प्रदान करता है। एक गैर-यूनिटरी, सैडल-पॉइंट समस्या को यूनिटरी क्वांटम विकास में बदलकर, यह विधि उच्च-विमीय तरल गतिकी की समस्याओं को हल करने के लिए क्वांटम कंप्यूटिंग की क्षमता का लाभ उठाती है जो वर्तमान में शास्त्रीय कंप्यूटरों के लिए अव्यवहार्य हैं। लेखक इस बात पर जोर देते हैं कि उनका दृष्टिकोण समस्या की विमा में "प्रमाणित घातीय गति" (provable exponential speedup) प्रदान करता है, जो इसे जटिल ज्यामिति और उच्च-विमीय हाइड्रोडायनामिक्स के भविष्य के अनुप्रयोगों के लिए एक व्यवहार्य उम्मीदवार बनाता है। यह कार्य ऐसे एल्गोरिदम को निकट-अवधि (near-term) क्वांटम हार्डवेयर पर पोर्ट करने और अंततः गैर-रेखीय नेवियर-स्टोक्स समीकरणों तक इस ढांचे का विस्तार करने की दिशा में एक आधारभूत कदम के रूप में स्थित है।