⚛️ quantum physics

An HHL-Based Quantum-Classical Solver for the Incompressible Navier-Stokes Equations with Approximate QST

यह शोध पत्र एक हाइब्रिड क्वांटम-क्लासिकल सॉल्वर प्रस्तुत करता है जो इनकंप्रेसिबल नेवियर-स्टोक्स समीकरणों को कुशलतापूर्वक हल करने के लिए हैरो-हॅसिमिड-लॉयड (HHL) एल्गोरिदम को चेबिशेव-आधारित अनुमानित क्वांटम स्टेट टोमोग्राफी के साथ एकीकृत करता है, और IBM के Qiskit फ्रेमवर्क का उपयोग करके लिड-ड्रिवन कैविटी और टेलर-ग्रीन वोर्टेक्स प्रवाह के सटीक सिमुलेशन के माध्यम से इस दृष्टिकोण को सफलतापूर्वक मान्य करता है।

मूल लेखक: Moshe Inger, Steven Frankel

प्रकाशित 2026-03-20

📖 6 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Moshe Inger, Steven Frankel

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

मुख्य चित्र: एक क्वांटम-क्लासिकल टीम की जुगलबंदी

कल्पना कीजिए कि आप यह अनुमान लगाने की कोशिश कर रहे हैं कि एक जटिल तरल (जैसे पंख के ऊपर से गुजरने वाली हवा या पाइप में बहता पानी) कैसे गति करेगा। इसे कंप्यूटेशनल फ्लूइड डायनेमिक्स (CFD) कहा जाता है। यह हजारों नर्तकों के एक भीड़भाड़ वाले कमरे में नाचने के रास्ते का अनुमान लगाने जैसा है।

इस भविष्यवाणी का सबसे कठिन हिस्सा यह पता लगाना नहीं है कि नर्तक कैसे हिलते हैं; बल्कि यह समझना है कि वे एक विशिष्ट "घनत्व" (असंपीड्यता/incompressibility) बनाए रखने के लिए एक-दूसरे पर कैसे दबाव डालते हैं। गणितीय शब्दों में, इसके लिए हर एक स्टेप पर पॉइसन इक्वेशन (Poisson Equation) नामक एक विशाल, जटिल पहेली को हल करना आवश्यक होता है।

समस्या: एक साधारण कंप्यूटर पर, इस पहेली को हल करना समुद्र तट पर रेत के हर एक कण को एक-एक करके गिनने जैसा है। यह कंप्यूटर के 90% समय को खा जाता है और सब कुछ धीमा कर देता है।

समाधान: इस पेपर के लेखकों ने एक हाइब्रिड टीम बनाई है। उन्होंने "नृत्य" (तरल की गति) को एक साधारण कंप्यूटर पर ही रखा, लेकिन "रेत गिनने" वाली पहेली (दबाव/pressure) को बहुत तेज़ी से हल करने के लिए एक क्वांटम कंप्यूटर को काम पर लगाया।

पात्रों का परिचय

क्लासिकल कंप्यूटर (मैनेजर):
यह वह मानक कंप्यूटर है जिसे हम रोज़ाना उपयोग करते हैं। यह तरल की गति (वेग/velocity) के भारी कार्यों को संभालता है और समग्र सिमुलेशन का प्रबंधन करता है। यह भरोसेमंद है लेकिन विशिष्ट "दबाव" वाले गणित में धीमा है।
क्वांटम कंप्यूटर (तेज़ विशेषज्ञ):
यह नया खिलाड़ी है। यह दबाव की पहेली को हल करने के लिए क्वांटम मैकेनिक्स के नियमों का उपयोग करता है। संख्याओं को एक-एक करके जाँचने के बजाय, यह एक साथ कई संभावनाओं की जाँच करता है। यह पेपर एक प्रसिद्ध एल्गोरिदम का उपयोग करता है जिसे HHL (इसके आविष्कारकों के नाम पर) कहा जाता है, जो इन विशिष्ट प्रकार की पहेलियों के लिए एक सुपर-फास्ट कैलकुलेटर की तरह काम करता है।
अनुवादक (चेबिशेव पॉलिनोमिअल्स - Chebyshev Polynomials):
यही पेचीदा हिस्सा है। क्वांटम कंप्यूटर पहेली को हल करता है और उत्तर एक "क्वांटम भाषा" (संभावनाओं के बादल) में देता है। यदि आप पूरे उत्तर को पढ़ने की कोशिश करेंगे, तो इसमें बहुत समय लगेगा (जो क्वांटम कंप्यूटर का उपयोग करने के उद्देश्य को ही विफल कर देगा)।
इसे ठीक करने के लिए, लेखकों ने चेबिशेव पॉलिनोमिअल्स पर आधारित एक चतुर अनुवादक का उपयोग किया। इसे एक गाने को सुनने और हर एक वाद्य यंत्र की आवाज़ को ट्रांसक्राइब करने के बजाय केवल मुख्य धुन के नोट्स लिखने जैसा समझें। यह उन्हें पूरे क्वांटम स्टेट को पढ़े बिना दबाव के क्षेत्र का "सार" जल्दी से प्राप्त करने की अनुमति देता है।

उन्होंने इसका परीक्षण कैसे किया: दो "जिम्नास्टिक" रूटीन

यह साबित करने के लिए कि उनकी टीम काम करती है, उन्होंने दो प्रसिद्ध फ्लूइड डायनेमिक्स चुनौतियों पर इसका परीक्षण किया:

1. लिड-ड्रिवन कैविटी (एक "जेली का डिब्बा")

एक वर्गाकार बॉक्स की कल्पना करें जो पानी से भरा है। ऊपर का ढक्कन दाईं ओर फिसलता है, जिससे पानी उसके साथ खिंचता है, जिससे एक विशाल घूमता हुआ भंवर (vortex) बनता है।

परीक्षण: उन्होंने इसे 16x16 ग्रिड पर सिम्युलेट किया।
परिणाम: हाइब्रिड टीम ने सफलतापूर्वक घूमते हुए भंवर को फिर से बनाया। क्वांटम कंप्यूटर का उत्तर पूर्ण क्लासिकल उत्तर की तुलना में लगभग 92% सटीक था। त्रुटियां ज्यादातर कोनों में थीं जहाँ पानी बहुत स्थिर (शून्य के करीब) हो जाता है, जिसे सटीक रूप से मापना कठिन होता है।

2. टेलर-ग्रीन वोर्टेक्स (एक "परफेक्टली स्विर्लिंग डोनट")

यह एक गणितीय परीक्षण है जहाँ तरल एक सटीक, अनुमानित पैटर्न में घूमता है जो समय के साथ कम होता जाता है। चूंकि हम गणितीय रूप से सटीक उत्तर जानते हैं, इसलिए यह एक आदर्श परीक्षण है।

परिणाम: हाइब्रिड टीम यहाँ अविश्वसनीय रूप से सटीक थी! उन्होंने तरल की गति के लिए 99% सटीकता और दबाव के लिए 96% सटीकता हासिल की। इसने साबित कर दिया कि क्वांटम विधि केवल "अनुमान" नहीं लगाती है; यह वास्तव में भौतिकी की सही गणना करती है।

बाधाएं जिनका उन्हें सामना करना पड़ा

भले ही यह सफल रहा, लेकिन टीम को कुछ मुश्किलों का सामना करना पड़ा:

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

यह क्यों महत्वपूर्ण है

यह पेपर एक प्रूफ ऑफ कॉन्सेप्ट है। यह राइट ब्रदर्स की पहली उड़ान की तरह है: इसने यात्रियों को महासागर के पार नहीं पहुँचाया, लेकिन इसने साबित कर दिया कि एक मशीन उड़ सकती है।

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

निष्कर्ष:
हम अभी मौसम के पूर्वानुमान के लिए अपने सुपर कंप्यूटरों को क्वांटम कंप्यूटरों से बदलने के लिए तैयार नहीं हैं। लेकिन यह पेपर साबित करता है कि निकट भविष्य में, हम फ्लूइड सिमुलेशन के सबसे कठिन हिस्सों को तेज करने के लिए क्वांटम कंप्यूटरों को एक "टर्बोचार्जर" के रूप में उपयोग कर सकते हैं, जिससे बेहतर हवाई जहाज, कार और जलवायु मॉडल डिजाइन करने में मदद मिल सकती है।

यहाँ मोशे इंगर और स्टीवन फ्रैंकल के शोध पत्र "An H-H-L-Based Quantum-Classical Solver for the Incompressible Navier–Stokes Equations with Approximate QST" का विस्तृत तकनीकी सारांश दिया गया है।

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

असंपीड्य नेवियर-स्टोक्स समीकरणों (incompressible Navier–Stokes equations) का संख्यात्मक समाधान कम्प्यूटेशनल फ्लूइड डायनेमिक्स (CFD) का एक आधार स्तंभ है। हालाँकि, मानक प्रोजेक्शन विधियों को विसरण-मुक्त स्थिति ( $\nabla \cdot u = 0$ ) को लागू करने के लिए प्रत्येक टाइम स्टेप पर दबाव (pressure) के लिए पॉइसन समीकरण (Poisson equation) को हल करने की आवश्यकता होती है।

बाधा (The Bottleneck): पॉइसन समीकरण के विव כךन (discretization) से रैखिक बीजगणितीय समीकरणों का एक बड़ा, विरल (sparse) तंत्र ($Ax = b$) प्राप्त होता है। इस प्रणाली को हल करने में आमतौर पर कुल निष्पादन समय का 90% तक खर्च होता है।
क्वांटम अवसर (The Quantum Opportunity): हैरो-हसिमिड-लॉयड (HHL) एल्गोरिदम विरल रैखिक प्रणालियों को हल करने के लिए एक सैद्धांतिक घातीय गति (exponential speedup) प्रदान करता है। हालाँकि, दो प्रमुख बाधाएँ इसके व्यावहारिक अनुप्रयोग को रोकती हैं:
1. रीडआउट सीमा (Readout Limitation): HHL समाधान के समान एक क्वांटम अवस्था $|x\rangle$ आउटपुट करता है। पूर्ण वेक्टर (statevector) निकालने के लिए घातीय संख्या में मापन की आवश्यकता होती है, जो गति (speedup) को निष्प्रभावी कर देता है।
2. स्टेट प्रिपरेशन (State Preparation): शास्त्रीय डेटा को क्वांटम अवस्था (एम्प्लीट्यूड एनकोडिंग) में लोड करना एक खुली चुनौती है।
3. हाइब्रिड बाधाएं (Hybrid Constraints): जबकि HHL $2^n$ डेटा बिंदुओं को $n$ क्विबिट्स में एनकोड करता है, डेटा को क्लासिकल और क्वांटम प्रोसेसर के बीच वापस भेजने की आवश्यकता वर्तमान में लाभ को सीमित करती है, क्योंकि आउटपुट के लिए अभी भी $2^n$ क्लासिकल बिट्स की आवश्यकता होती है।

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

लेखक एक हाइब्रिड क्वांटम-क्लासिकल इटरेटिव सॉल्वर का प्रस्ताव करते हैं जो 2D असंपीड्य नेवियर-स्टोक्स समीकरणों के लिए एक क्लासिकल प्रोजेक्शन स्कीम में HHL एल्गोरिदम को एकीकृत करता है।

A. क्लासिकल फ्रेमवर्क (Classical Framework)

शासी समीकरण (Governing Equations): सॉल्वर एक स्प्लिट-स्टेप प्रोजेक्शन विधि का उपयोग करता है। एक मध्यवर्ती वेग $u^*$ स्पष्ट रूप से (explicitly) गणना किया जाता है, जिसके बाद असंपीड्यता को लागू करने के लिए प्रेशर करेक्शन स्टेप किया जाता है।
विवीकरण (Discretization): डोमेन को 2D कार्टेशियन ग्रिड पर सेकंड-ऑर्डर सेंट्रल डिफरेंस स्कीम का उपयोग करके विवीकृत किया गया है। 2D ग्रिड को 1D वेक्टर में मैप किया गया है ताकि पॉइसन समीकरण ( $\nabla^2 p = \rho/\Delta t \nabla \cdot u^*$ ) के लिए एक विरल, ब्लॉक-ट्राइडायगोनल मैट्रिक्स $A$ बनाया जा सके।
कर्विलिनियर मैपिंग (Curvilinear Mapping): बाउंडरी लेयर्स (विशेष रूप से लिड-ड्रिवन कैविटी के लिए) को संभालने के लिए, एक हाइपरबोलिक स्ट्रेचिंग फंक्शन का उपयोग किया गया है ताकि उच्च ग्रेडिएंट वाले क्षेत्रों के पास ग्रिड पॉइंट्स को केंद्रित किया जा सके, जिससे रिज़ॉल्यूशन में सुधार हो।

B. क्वांटम सबरूटीन (HHL)

एल्गोरिदम: HHL एल्गोरिदम का उपयोग $Ap = b $को हल करने के लिए किया जाता है, जहाँ$ b $मध्यवर्ती वेग का डाइवर्जेंस है और$ p$ दबाव है।
कार्यान्वयन (Implementation):
- एनकोडिंग (Encoding): मैट्रिक्स $A$ को पाउली स्ट्रिंग्स ( $A = \sum c_k P_k$ ) में विघटित किया जाता है।
- हैमिल्टोनियन सिमुलेशन (Hamiltonian Simulation): टाइम इवोल्यूशन $e^{-iAt}$ को $R=150$ स्टेप्स के साथ फर्स्ट-ऑर्डर ली-ट्रोटर प्रोडक्ट फॉर्मूला का उपयोग करके अनुमानित किया जाता है।
- आइगेनवैल्यू इनवर्जन (Eigenvalue Inversion): इसे Qiskit के ExactReciprocal फंक्शन का उपयोग करके एक अनसेलिया क्विबिट पर कंट्रोल्ड रोटेशन के माध्यम से किया जाता है।
- स्टेट प्रिपरेशन (State Preparation): लेखक इनपुट वेक्टर $|b\rangle$ को लोड करने के लिए एक आदर्श ओरेकल (ideal oracle) मानकर चलते हैं, जिसे वे वर्तमान सीमा के रूप में स्वीकार करते हैं।

C. अनुमानित क्वांटम स्टेट टोमोग्राफी (Approximate QST)

पूर्ण स्टेट रिकंस्ट्रक्शन की आवश्यकता के बिना रीडआउट बाधा को दूर करने के लिए, लेखक एक चेबिशेव बहुपद-आधारित (Chebyshev polynomial-based) QST दृष्टिकोण का उपयोग करते हैं:

स्पेक्ट्रल प्रोजेक्शन (Spectral Projection): प्रत्येक एम्प्लीट्यूड को मापने के बजाय, क्वांटम अवस्था $|p\rangle$ को $m$ चेबिशेव बहुपदों के एक ट्रंकेटेड बेसिस पर प्रोजेक्ट किया जाता है।
हैडामार्ड टेस्ट (Hadamard Test): बहुपदों के गुणांकों (coefficients) का अनुमान लगाने के लिए हैडामार्ड टेस्ट का उपयोग किया जाता है।
दक्षता (Efficiency): यह रीडआउट जटिलता को घातीय से घटाकर $O(md) $कर देता है, जहाँ$ m $बहुपदों की संख्या है और$ d$ आयाम है।
कैलिब्रेशन (Calibration): चूंकि HHL एक नॉर्मलाइज्ड अवस्था आउटपुट करता है, इसलिए लेखक भौतिक परिमाण (physical magnitude) को पुनः प्राप्त करने के लिए क्वांटम प्रेशर के $L_2$ नॉर्म को क्लासिकल रेफरेंस के साथ मैच करके प्रत्येक टाइम स्टेप पर एक एनर्जी कैलिब्रेशन करते हैं।

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

एंड-टू-एंड हाइब्रिड सॉल्वर: नेवियर-स्टोक्स समीकरणों को हल करने वाले एक पूर्णतः एकीकृत क्वांटम-क्लासिकल लूप का पहला प्रदर्शन, जहाँ HHL सबरूटीन क्लासिकल पॉइसन सॉल्वर को प्रतिस्थापित करता है।
अनुमानित QST एकीकरण: क्वांटम अवस्था से फ्लूइड डायनेजिक्स डेटा निकालने के लिए चेबिशेव-आधारित QST को सफलतापूर्वक लागू करना, जो घातीय मापन की आवश्यकता को दरकिनार करता है।
बेंचमार्क वैलिडेशन: IBM के Qiskit फ्रेमवर्क का उपयोग करके मानक CFD बेंचमार्क (Lid-Driven Cavity और Taylor-Green Vortex) के विरुद्ध सॉल्वर का सत्यापन।
ग्रिड मैपिंग विश्लेषण: कर्विलिनियर ग्रिड स्ट्रेचिंग (जो QST सटीकता में मदद करता है) और मैट्रिक्स कंडीशनिंग (जो HLR सटीकता को प्रभावित करता है) के बीच ट्रेड-ऑफ की जांच।

4. परिणाम (Results)

सॉल्वर का परीक्षण $16 \times 16$ ग्रिड और 8 क्लॉक क्विबिट्स ( $n_c=8$ ) के साथ दो बेंचमार्क समस्याओं पर किया गया।

पॉइसन समीकरण सॉल्वर:
- क्लासिकल समाधानों की तुलना में ~2.3% (1D) और ~2.5% (2D) का औसत सापेक्ष त्रुटि (Average Relative Error - ARE) प्राप्त किया।
- सटीकता क्लॉक क्विबिट्स की संख्या के प्रति संवेदनशील थी, जिसमें $n_c=8$ पर महत्वपूर्ण सुधार देखा गया।
लिड-ड्रिवन कैविटी फ्लो ($Re=100$):
- वेलोसिटी फील्ड: हाइब्रिड सॉल्वर ने 8.17% के औसत ARE के साथ ग्लोबल वोर्टेक्स डायनेमिक्स को कैप्चर किया।
- प्रेशर ग्रेडिएंट: कैविटी के निचले कोनों में उच्च त्रुटियां (416% तक) देखी गईं जहाँ दबाव शून्य के करीब होता है। इसका कारण चेबिशेव एप्रोक्सिमेशन की वह कठिनाई है जहाँ एक सीमित बेसिस के साथ उच्च-परिमाण और शून्य के करीब मानों को एक साथ कैप्चर करना मुश्किल होता है।
- अभिसरण (Convergence): सेंटरलाइन्स के साथ वेलोसिटी प्रोफाइल घिया (Ghia) बेंचमार्क मानक की ओर अभिसरित हुए।
टेलर-ग्रीन वोर्टेक्स (TGV):
- सत्यापन के लिए इस परीक्षण में एक सटीक विश्लेषणात्मक समाधान का उपयोग किया गया।
- परिणाम: सॉल्वर ने वेलोसिटी के लिए ~1% और प्रेशर के लिए ~4% के औसत ARE के साथ उच्च सटीकता प्राप्त की।
- परिणाम पुष्टि करते हैं कि हाइब्रिड दृष्टिकोण विस्कस डिसिपेशन द्वारा संचालित वोर्टेक्स के घातीय क्षय (exponential decay) को सटीक रूप से मॉडल कर सकता है।

5. महत्व और भविष्य की दिशाएं (Significance and Future Directions)

प्रूफ ऑफ कॉन्सेप्ट: यह कार्य प्रदर्शित करता है कि क्वांटम सबरूटीन्स को बिना किसी विनाशकारी विचलन (catastrophic divergence) के, वर्तमान सिमुलेशन सीमाओं के बावजूद, इटरेटिव CFD वर्कफ़्लो में सफलतापूर्वक एम्बेड किया जा सकता है।
रीडआउट बाधा को संबोधित करना: चेबिशेव QST का उपयोग क्वांटम अवस्थाओं से उपयोगी भौतिक डेटा निकालने के लिए एक व्यवहार्य मार्ग प्रदान करता है, जो फ्लूइड डायनेमिक्स में व्यावहारिक क्वांटम लाभ के लिए एक महत्वपूर्ण कदम है।
सीमाएं और चुनौतियां:
- स्टेट प्रिपरेशन: इनपुट लोडिंग के लिए एक आदर्श ओरेकल पर निर्भरता एक सैद्धांतिक बाधा बनी हुई है।
- स्केलिंग: क्विबिट बाधाओं के कारण वर्तमान सिमुलेशन छोटे ग्रिड ( $16 \times 16$ ) तक सीमित हैं।
- बाउंड्री पर त्रुटि: जब एक ही डोमेन में उच्च-परिमाण मान मौजूद हों, तो QST पद्धति शून्य के करीब मानों वाले क्षेत्रों के साथ संघर्ष करती है।
भविष्य का कार्य:
- 3D नेवियर-स्टोक्स समीकरणों तक विस्तार करना।
- आदर्श स्टेट प्रिपरेशन को कुशल एल्गोरिदम (जैसे टेंसर नेटवर्क या संशोधित HHL) से बदलना।
- दोष-सहिष्णु (fault-tolerant) HHL से वेरिएशनल क्वांटम लीनियर सॉल्वर्स (VQLS) की ओर संक्रमण, जो NISQ हार्डवेयर के अनुकूल है।
- ऐसे पूर्ण क्वांटम स्कीमा विकसित करना जो क्लासिकल एडवेक्शन-डिफ्यूजन स्टेप्स को हैमिल्टोनियन इवोल्यूशन से बदल सकें।

निष्कर्षतः, यह शोध पत्र क्वांटम CFD के लिए एक आधारभूत ढांचा स्थापित करता है, यह सिद्ध करता है कि हाइब्रिड आर्किटेक्चर जटिल फ्लूइड डायनेमिक्स समस्याओं को हल कर सकते हैं और भविष्य के उच्च-रेनॉल्ड्स नंबर सिमुलेशन में क्वांटम कंप्यूटिंग को एकीकृत करने के लिए एक मजबूत मार्ग प्रदान करते हैं।