Two-dimensional quantum lattice gas algorithm for anisotropic Burger-like equations

यहाँ निकोलो फोनियो, पियरे सागौट और ज्यूसेपी डी मोलफेटा के शोध पत्र "Two-dimensional quantum lattice gas algorithm for anisotropic Burger-like equations" का विस्तृत तकनीकी सारांश दिया गया है।

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

कंप्यूटेशनल फ्लूइड डायनेमिक्स (CFD) में क्वांटम कंप्यूटिंग (QC) का अनुप्रयोग, भौतिक सटीकता और संख्यात्मक स्थिरता (numerical stability) के साथ क्वांटम लाभ को संतुलित करने की चुनौतियों का सामना करता है। मौजूदा दृष्टिकोण अक्सर शास्त्रीय लैट्टिस गैस मॉडल (जैसे लैट्टिस बोल्ट्ज़मैन मेथड - LBM) को क्वांटम अवस्थाओं में एनकोड करने पर निर्भर करते हैं, जो गैर-यूनिटरी (non-unitary) कोलिजन ऑपरेटर्स से ग्रस्त हो सकते हैं, जिसके लिए जटिल मल्टी-टाइम-स्टेप कार्यान्वयन या लीनियराइजेशन तकनीकों की आवश्यकता होती है जो टर्बुलेंस (turbulence) के साथ संघर्ष करती हैं।

विशेष रूप से, लेखक येपेज़ क्वांटम लैट्टिस गैस (QLG) मॉडल (Phys. Rev. A, 2006) का पुनर्मूल्यांकन करते हैं। जबकि यह मॉडल यूनिटरी है और वैचारिक रूप से पूर्ण है, पिछले स्वरूपों में मैक्रोस्कोपिक श्यानता (viscosity) के व्युत्पन्न में एक सन्निकटन (approximation) शामिल था। इसके अलावा, इस मॉडल को दो आयामों (2D) में विस्तारित करना, जबकि विश्लेषणात्मक सुगमता बनाए रखी जाए और संसाधनों (qubits) के उपयोग को न्यूनतम किया जाए, एक खुला विषय बना हुआ है। लक्ष्य एक सुधारे हुए श्यानता (viscosity) को व्युत्पन्न करना, न्यूनतम संसाधनों (प्रति साइट 2 क्वबिट्स) के साथ मॉडल को 2D में सामान्य बनाना, और परिणामी मैक्रोस्कोपिक आंशिक अंतर समीकरणों (PDEs) को निर्धारित करना है।

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

A. सैद्धांतिक ढांचा: Q-D1Q2 मॉडल

लेखक 2 क्वबिट प्रति साइट वाले 1D क्वांटम लैट्टिस गैस (Q-D1Q2) का विश्लेषण करते हैं।

• अवस्था प्रतिनिधित्व (State Representation): प्रत्येक लैट्टिस साइट $x$ एक अवस्था $|\psi(x,t)\rangle = |q_0\rangle \otimes |q_1\rangle$ रखती है, जहाँ क्वबिट्स विविक्त वेगों (discrete velocities) के साथ चलने वाली कण आबादी का प्रतिनिधित्व करते हैं।
• एल्गोरिदम: विकास तीन चरणों से मिलकर बना है:
  1. आरंभीकरण (Initialization): शास्त्रीय वितरण फलनों $f_i$ के आधार पर क्वबिट्स को आरंभीकृत किया जाता है।
  2. कोलिजन (Collision): द्रव्यमान संरक्षण के लिए एक यूनिटरी ऑपरेटर $\hat{C}$ (यूलर कोण $\theta, \zeta, \xi$ द्वारा पैरामीटराइज्ड) लागू किया जाता है।
  3. स्ट्रीमिंग और मापन (Streaming & Measurement): वितरणों को अपडेट करने के लिए आबादी का मापन (अवस्था का पतन) किया जाता है, जिसके बाद शास्त्रीय स्ट्रीमिंग होती है।
• व्युत्पत्ति (Derivation): लेखक क्वांटम विकास से प्राप्त लैट्टिस बोल्ट्ज़मैन समीकरण पर एक चैपन-एनस्कोग विस्तार (Chapman-Enskog expansion) करते हैं। वे डिफ्यूसिव स्केलिंग रिजीम ($\delta x \sim \epsilon, \delta t \sim \epsilon^2$) को मानते हुए, स्पेस ($\delta x$) में दूसरे क्रम और समय ($\delta t$) में प्रथम क्रम तक परिमित अंतर समीकरण का विस्तार करते हैं।

B. श्यानता सुधार (Viscosity Correction)

एक प्रमुख कार्यप्रणाली चरण $J_1 - J_0$ (साम्यावस्था पर आबादी के सापेक्ष कोलिजन टर्म का व्युत्पन्न) की कठोर गणना है।

• पूर्ववर्ती कार्य (येपेज़) में इस पद का सन्निकटन किया गया था, जिससे श्यानता अभिव्यक्ति $\nu \propto \cot^2(\theta)$ प्राप्त हुई थी।
• लेखक साम्यावस्था आबादी और कोलिजन टर्म के लिए एक सटीक क्लोज्ड-फॉर्म समाधान व्युत्पन्न करते हैं, जिसके परिणामस्वरूप एक सुधारा हुआ श्यानता सूत्र प्राप्त होता है जिसमें लो-मैक (low-Mach) सीमा में उच्च-क्रम के पद शामिल होते हैं।

C. 2D सामान्यीकरण (Q-D2Q2)

मॉडल को एक वर्गाकार लैट्टिस का उपयोग करके 2D में विस्तारित किया गया है जिसमें प्रति साइट केवल 2 क्वबिट्स हैं (न्यूनतम संसाधन उपयोग)।

• स्ट्रीमिंग: आबादी $x$ और $y$ दोनों दिशाओं में स्ट्रीम होती है।
• कोलिजन: समान 1D कोलिजन ऑपरेटर स्थानीय रूप से लागू किया जाता है।
• विश्लेषण: लेखक यह देखने के लिए विभिन्न विविक्त वेग सेटों (जैसे वर्गाकार लैट्टिस, त्रिकोणीय लैट्टिस) के लिए मैक्रोस्कोपिक PDE व्युत्पन्न करते हैं कि विसंगति (anisotropy) और एडवेक्शन पद कैसे उभरते हैं।

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

1. श्यानता का विश्लेषणात्मक सुधार:
   लेखक 1D मॉडल में किनेमैटिक श्यानता $\nu$ के लिए एक सुधारा हुआ अभिव्यक्ति व्युत्पन्न करते हैं। नया सूत्र है:
   $$\nu = -\frac{\delta x^2}{\delta t} \frac{1}{2} \left( 1 - \frac{\sqrt{\alpha^2 + 1}}{\sin^2(\theta)} \right)$$
   जहाँ $\alpha = \cot(\theta)\cos(\zeta - \xi)$ है। यह पिछले $\cot^2(\theta)$ सन्निकटन को ठीक करता है, और $O(\epsilon^4)$ तक बर्गर्स-जैसे समीकरण की वैधता को बढ़ाता है।

2. न्यूनतम 2D सामान्यीकरण:
   वे एक Q-D2Q2 मॉडल प्रस्तावित करते हैं (2 आयाम, प्रति साइट 2 क्वबिट्स)। पारंपरिक 2D लैट्टिस गैसों (जैसे FHP) के विपरीत, जिन्हें 6 या 7 वेगों की आवश्यकता होती है, यह मॉडल केवल 2 का उपयोग करता है।

• परिणाम: मॉडल सफलतापूर्वक विसंगत बर्गर्स-जैसे डायनेमिक्स (anisotropic Burgers-like dynamics) को पुनरुत्पादित करता है।
• समीकरण: व्युत्पन्न मैक्रोस्कोपिक समीकरण है:
  $$\partial_t \rho + \mathbf{d} \cdot \nabla \rho - \nabla \cdot (\mathbf{D} \nabla \rho) = 0$$
  जहाँ $\mathbf{d}$ एक एडवेक्टिव वेक्टर है और $\mathbf{D}$ एक विसंगत डिफ्यूजन टेंसर (anisotropic diffusion tensor) है जो चुने गए वेग सेट पर निर्भर करता है।

3. स्थिरता के साथ अनिश्चित रूप से छोटी श्यानता:
   यह शोध पत्र प्रदर्शित करता है कि QLG मॉडल एकल-रिलैक्सेशन-टाइम LBM मॉडलों में पाई जाने वाली संख्यात्मक अस्थिरता (जो $\tau \to 0.5$ होने पर अस्थिर हो जाते हैं) के बिना अनिश्चित रूप से छोटी श्यानता प्राप्त कर सकता है। यह इसकी अंतर्निहित यूनिटरी संरचना और H-थ्योरम की संतुष्टि के कारण है।

4. परिणाम (Results)

• श्यानता का संख्यात्मक सत्यापन: 1D में सिमुलेशन ने सुधारा हुआ श्यानता सूत्र की पुष्टि की। प्रयोगात्मक श्यानता, पिछले सन्निकटन (Eq. 19) की तुलना में, सैद्धांतिक भविष्यवाणी (Eq. 22) से काफी बेहतर मेल खाती है, विशेष रूप से लो-विस्कोसिटी रिजीम में।
• शॉक फॉर्मेशन और स्थिरता:
  • मॉडल ने बर्गर्स समीकरण में शॉक फॉर्मेशन को सफलतापूर्वक सिम्युलेट किया।
  • ट्रेड-ऑफ (Trade-off): बहुत कम श्यानता (तीव्र शॉक) प्राप्त करने के लिए थर्मललाइजेशन (रिलैक्सेशन) के लिए लंबे सिमुलेशन समय और/या उच्च ग्रिड रिज़ॉल्यूशन की आवश्यकता होती है, लेकिन स्कीम अनकंडीशनली स्टेबल (unconditionally stable) रहती है।
• 2D विसंगत व्यवहार (Anisotropic Behavior):
  • विभिन्न वेग सेटों ने अलग-अलग मैक्रोस्कोपिक व्यवहार उत्पन्न किए।
  • सममित सेट (जैसे $c_0 = -c_1$) के परिणामस्वरूप 2D पर प्रोजेक्टेड 1D-जैसे डायनेमिक्स प्राप्त हुए।
  • असममित सेट (जैसे $c_0$ क्षैतिज, $c_1$ ऊर्ध्वाधर) ने क्रॉस-डेरिवेटिव पद ($\partial_{xy}\rho$) और अतिरिक्त एडवेक्टिव पद पेश किए, जो न्यूनतम 2D मॉडल की विसंगत प्रकृति की पुष्टि करते हैं।
• शास्त्रीय विधियों के साथ तुलना: व्युत्पन्न 2D PDE के लिए QLG सिमुलेशन और फाइनाइट डिफरेंस मेथड (FDM) सॉल्वर के बीच संख्यात्मक तुलना उत्कृष्ट गुणात्मक और मात्रात्मक सहमति ($L_2$ norm error analysis) दिखाती है, जो व्युत्पन्न PDE को मान्य करती है।

5. महत्व और भविष्य के परिप्रेक्ष्य (Significance and Future Perspectives)

• क्वांटम-नेटिव सिमुलेशन: यह कार्य "क्वांटम-नेटिव" फ्लूइड सिमुलेशन के लिए एक मार्ग स्थापित करता है जो केवल शास्त्रीय एल्गोरिदम को एनकोड नहीं करते हैं, बल्कि यूनिटरी क्वांटम डायनेमिक्स का उपयोग करते हैं ताकि विसरित घटनाओं (श्यानता) को स्वाभाविक रूप से रिवर्सिबल माइक्रोस्कोपिक नियमों के भीतर समाहित किया जा सके।
• संसाधन दक्षता: केवल प्रति साइट 2 क्वबिट्स के साथ 2D विसंगत डायनेमिक्स प्राप्त करके, यह मॉडल शास्त्रीय लैट्टिस गैस ऑटोमेटा (जैसे FHP) और मानक LBM के लिए एक संसाधन-कुशल विकल्प प्रदान करता है, जिन्हें समान भौतिकी को पकड़ने के लिए अधिक अवस्थाओं की आवश्यकता होती है।
• नेवियर-स्टोक्स से जुड़ाव: हालांकि वर्तमान 2D मॉडल बर्गर्स-जैसे डायनेमिक्स (एक स्केलर संरक्षण) को कैप्चर करता है, लेखक पहचानते हैं कि तीसरा वेग (3 क्वबिट्स) जोड़ने से मोमेंटम संरक्षण सक्षम हो सकता है, जो संभावित रूप से एक पूर्ण नेवियर-स्टोक्स सॉल्वर की ओर ले जा सकता है।
• हार्डवेयर प्रासंगिकता: न्यूनतम क्वबिट आवश्यकता और अनकंडीशनल स्टेबिलिटी इस दृष्टिकोण को विशेष रूप से निकट-अवधि क्वांटम हार्डवेयर (NISQ डिवाइसेस) पर कार्यान्वयन के लिए आशाजनक बनाती है।

निष्कर्षतः, यह शोध पत्र क्वांटम लैट्टिस गैस एल्गोरिदम के सैद्धांतिक आधार को परिष्कृत करता है, श्यानता व्युत्पन्न के एक लंबे समय से चले आ रहे सन्निकटन को ठीक करता है, और सफलतापूर्वक एक न्यूनतम, संसाधन-कुशल 2D विस्तार प्रदर्शित करता है जो जटिल गैर-रेखीय फ्लूइड डायनेमिक्स का अनुकरण करने में सक्षम है।