Let the Qudit Do the Jacobi: A Structured Quantum Algorithm for Spectral Decomposition

तकनीकी सारांश: स्पेक्ट्रल अपघटन के लिए एक संरचित क्वांटम एल्गोरिदम

समस्या विवरण
स्पेक्ट्रल अपघटन (Spectral decomposition) भौतिकी और इंजीनियरिंग में एक मौलिक कार्य है, जिसके लिए हैमिल्टनियन या यूनिटरी ऑपरेटरों के आइजनवैल्यू (eigenvalues) और आइजनवेक्टर (eigenvectors) का निर्धारण करना आवश्यक होता है। जबकि शास्त्रीय संख्यात्मक रैखिक बीजगणित (classical numerical linear algebra) सुदृढ़ विधियाँ प्रदान करता है, क्वांटम दृष्टिकोणों के सामने एक द्वंद्व है:

1. टोमोग्राफिक दृष्टिकोण: क्वांटम ऑपरेटर को शास्त्रीय सूचना में बदलने के लिए प्रक्रिया टोमोग्राफी (process tomography) के माध्यम से शास्त्रीय एल्गोरिदम लागू करना। यह संसाधन-गहन है और ऑपरेटर की क्वांटम प्रकृति को खो देता है।
2. सुसंगत (Coherent) दृष्टिकोण: ऑपरेटर को क्वांटम डोमेन में ही बनाए रखना। मौजूदा विधियाँ जैसे क्वांटम फेज एस्टिमेशन (QPE), नियंत्रित-यूनिटरी संचालन और गहरे सर्किट की आवश्यकता रखती हैं, जो प्रयोगात्मक रूप से अत्यधिक मांग वाले होते हैं। वेरिएशनल दृष्टिकोण (जैसे VQE) अक्सर विकर्णन (diagonalization) को एक वैश्विक अनुकूलन समस्या के रूप में देखते हैं, जिसमें शास्त्रीय पुनरावृत्त विधियों की तरह संरचित अभिसरण गारंटी (convergence guarantees) का अभाव होता है।

एक ऐसे संरचित, पुनरावृत्त क्वांटम एल्गोरिदम के लिए एक अंतराल है जो मैट्रिक्स पुनर्निर्माण या नियंत्रित-यूनिटरी संचालन के बिना, सीधे क्वांटम हार्डवेयर पर अज्ञात यूनिटरी ऑपरेटरों को विकेंद्रीकृत कर सके, विशेष रूप से एक ऐसा जो उच्च-आयामी प्रणालियों (qudits) की मूल क्षमताओं का लाभ उठा सके।

कार्यप्रणाली: जैकोबी क्वडिट एल्गोरिदम (JQA)
लेखक जैकोबी क्वडिट एल्गोरिदम (JQA) का प्रस्ताव करते हैं, जो शास्त्रीय जैकोबी विकेंद्रीकरण विधि का एक क्वडिट-नेटिव (qudit-native) कार्यान्वयन है। इसका मूल दर्शन एक अज्ञात यूनिटरी ऑपरेटर $M$ को, जो एक एकल $d$-आयामी क्वडिट पर कार्य करता है, स्थानीय समानता रूपांतरणों (local similarity transformations) के एक अनुक्रम के माध्यम से उसके ऑफ-डायगोनल तत्वों को क्रमिक रूप से समाप्त करना है, न कि मैट्रिक्स का पुनर्निर्माण करना।

• रोटेशन का अपघटन: सामान्य मैट्रिसेस के लिए शास्त्रीय जैकोबी एल्गोरिदम में, एक दो-पैरामीटर गिवेंस रोटेशन $U^G_{(p,q)}(\phi, \alpha)$ को $(p,q)$ सब-ब्लॉक पर लागू किया जाता है ताकि ऑफ-डायगोनल तत्वों को न्यूनतम किया जा सके। JQA इस एकल दो-पैरामीटर अनुकूलन को दो क्रमिक, एकल-पैरामीटर क्वांटम वेरिएशनल अनुकूलनों में बदल देता है:

  1. एक "वास्तविक" यूनिटरी रोटेशन $U^R_{(p,q)}(\phi_1) = \exp(i\phi_1 \hat{g}^R_{(p,q)})$।
  2. एक "काल्पनिक" यूनिटरी रोटेशन $U^I_{(p,q)}(\phi_2) = \exp(i\phi_2 \hat{g}^I_{(p,q)})$।
     यहाँ, $\hat{g}^R$ और $\hat{g}^I$, $su(d)$ ली बीजगणित (Lie algebra) के ऑफ-डायगोनल जनरेटर हैं।

• वेरिएशनल खोज: $\phi_1$ और $\phi_2$ पैरामीटर डेरिवेटिव-मुक्त एक-आयामी खोजों के माध्यम से प्रयोगात्मक रूप से निर्धारित किए जाते हैं। एल्गोरिदम एक स्थानीय लागत फलन (cost function) $\tau_{(p,q)}$ को कम करता है, जो $(p,q)$ सब-ब्लॉक में ऑफ-डायगोनल तत्वों के वर्ग परिमाण का योग है। इसे निम्न प्रकार से प्राप्त किया जाता है:

  • क्वडिट को कम्प्यूटेशनल बेसिस अवस्थाओं $|p\rangle$ और $|q\rangle$ में प्रारंभ करना।
  • रोटेशन लागू करना और मैट्रिक्स तत्व के परिमाणों के अनुमान के लिए एक एकल डायगोनल ऑब्जर्वेबल $B$ (गैर-विकेंद्रित स्पेक्ट्रम के साथ) को मापना।
  • लागत फलन के लिए आवश्यक मैट्रिक्स तत्व परिमाणों के पुनर्निर्माण के लिए पूर्ण स्टेट टोमोग्राफी के बिना इसे दोहराना।

• पुनरावृत्त योजना (Iterative Scheme): एल्गोरिदम एक चक्रीय पिवोटिंग रणनीति (cyclic pivoting strategy) का उपयोग करता है। प्रत्येक "स्वीप" (sweep) में, प्रत्येक इंडेक्स जोड़ी $(p,q)$ को क्रमिक रूप से संबोधित किया जाता है। यह प्रक्रिया तब तक दोहराई जाती है जब तक कि वैश्विक ऑफ-डायगोनल लागत फलन $\tau(M)$ एक थ्रेशोल्ड से नीचे नहीं आ जाता।

• आइजनवैल्यू निष्कर्षण: अभिसरण होने पर, ऑपरेटर लगभग विकेंद्रित (diagonal) हो जाता है। आइजनवैल्यू (eigenphases) को एक इंटरफेरोमेट्रिक प्रोटोकॉल का उपयोग करके निकाला जाता है। बेसिस अवस्थाओं के सुपरपोजिशन तैयार करके, विकेंद्रित ऑपरेटर को लागू करने और उसके बाद हैडामार्ड जैसी इंटरफेरेंस प्रक्रियाओं का उपयोग करके, आइजनवैल्यू के बीच के फेज अंतर का अनुमान लगाया जाता है। पूर्ण स्पेक्ट्रम को एक ग्लोबल फेज तक पुनर्निर्मित किया जाता है।

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

1. संरचित क्वांटम एल्गोरिदम: अज्ञात यूनिटरी ऑपरेटरों को एक एकल क्वडिट पर विकेंद्रीकृत करने के लिए एक संरचित, पुनरावृत्त एल्गोरिदम का प्रस्ताव, जो स्पष्ट मैट्रिक्स पुनर्निर्माण या नियंत्रित-यूनिटरी गेट्स से बचता है।
2. क्रमिक वेरिएशनल अनुकूलन: यह प्रदर्शन कि दो-पैरामीटर जैकोबी रोटेशन को दो क्रमिक एकल-पैरामीटर वेरिएशनल खोजों में विभाजित किया जा सकता है। यह मैपिंग स्वाभाविक रूप से प्रयोगात्मक रूप से सुलभ क्वडिट संचालन के अनुरूप है और शास्त्रीय विधि के अभिसरण व्यवहार को संरक्षित करती है।
3. इंटरफेरोमेट्रिक आइजनवैल्यू निष्कर्षण: पूर्ण ऑपरेटर के शास्त्रीय पोस्ट-प्रोसेसिंग के बिना, हस्तक्षेप (interference) का उपयोग करके पूरी तरह से क्वांटम डोमेन के भीतर आइजनवैल्यू निकालने के लिए एक प्रोटोकॉल का परिचय।
4. अनुभवजन्य स्केलिंग: संख्यात्मक सत्यापन कि एल्गोरिदम $O(d^2)$ के कुल पिवट ऑपरेशंस के साथ एक निर्धारित सटीकता के लिए $O(d^2)$ के रूप में स्केल करता है।

परिणाम
Haar-रैंडम यूनिटरी मैट्रिसेस के एनसेम्बल पर $d=30$ तक के आयामों के लिए संख्यात्मक सिमुलेशन किए गए।

• अभिसरण (Convergence): एल्गोरिदम ऑफ-डायगोनल लागत फलन में निरंतर कमी प्रदर्शित करता है। यह शुरुआती चरणों में रैखिक कमी दिखाता है और डायगोनल फिक्स्ड पॉइंट के पास द्विघाती (quadratic) अभिसरण शासन की ओर त्वरित होता है, जो शास्त्रीय जैकोबी विधि की नकल करता है।
• स्वीप गणना: $d=10$ से $30$तक आयाम बढ़ने पर भी एक लक्षित सटीकता (जैसे,$\epsilon = 10^{-2}$या$10^{-4}$) तक पहुँचने के लिए आवश्यक जैकोबी स्वीप्स की संख्या लगभग स्थिर (6–7 स्वीप्स) रहती है।
• जटिलता: चूंकि प्रत्येक स्वीप में $O(d^2)$ पिवोट होते हैं और स्वीप की संख्या स्थिर है, इसलिए कुल ऑपरेशंस $O(d^2)$ के रूप में स्केल करते हैं।
• तुलना: क्रमिक एकल-पैरामीटर कार्यान्वयन (JQA) को मानक दो-पैरामीटर यूनिटरी जैकोबी एल्गोरिदम की तुलना में तुलनीय और परीक्षण किए गए मामलों में थोड़े कम स्वीप्स की आवश्यकता होती है।

महत्व और दावे
यह शोध पत्र JQA को शास्त्रीय विकेंद्रीकरण समस्याओं को तेजी से हल करने के तरीके के रूप में नहीं, बल्कि अज्ञात क्वांटम ऑपरेटरों या हैमिल्टनियन के प्रत्यक्ष विकेंद्रीकरण के उपकरण के रूप में प्रस्तुत करता है, जो प्रयोगात्मक रूप से सुलभ हैं लेकिन शास्त्रीय रूप से ज्ञात नहीं हैं।

लेखक इस कार्य के महत्व के रूप में निम्नलिखित दावे करते हैं:

• शास्त्रीय और क्वांटम के बीच सेतु: वैश्विक वेरिएशनल अनुकूलन या गहरे सुसंगत सर्किटों से आगे बढ़ते हुए, शास्त्रीय पुनरावृत्त मैट्रिक्स एल्गोरिदम और उनके क्वांटम वास्तविकताओं के बीच एक संरचित सेतु प्रदान करना।
• क्विड-नेटिव उपयुक्तता: यह उजागर करना कि एल्गोरिदम स्वाभाविक रूप से क्वडिट आर्किटेक्चर के लिए उपयुक्त है, जहाँ $SU(d)$के$SU(2)$ उपसमूहों को सीधे लागू किया जा सकता है, जिससे क्यूबिट-आधारित कार्यान्वयन के लिए आवश्यक कंपाइलेशन ओवरहेड से बचा जा सके।
• प्रयोगात्मक व्यवहार्यता: यह विधि केवल सिंगल-क्विड रोटेशन और एक एकल ऑब्जर्वेबल के मापन पर निर्भर करती है, जो एंसिला (ancilla) क्यूबिट्स या नियंत्रित-यूनिटरी ऑपरेशंस की आवश्यकता को समाप्त करती है, जिससे यह निकट-अवधि (near-term) और प्रारंभिक फॉल्ट-टोलरेंट उपकरणों के लिए उपयुक्त हो जाती है।
• समानांतरता (Parallelizability): एल्गोरिदम की संरचित प्रकृति विजातीय (disjoint) सबस्पेस पर कम्यूटिंग रोटेशन के संभावित समानांतर निष्पादन की अनुमति देती है, जो वितरित क्वांटम कंप्यूटिंग के लिए प्रासंगिक है।

लेखक संशोधित क्रमिक प्रक्रिया के लिए अभिसरण के सैद्धांतिक प्रमाण के संबंध में विनम्र रहते हैं, यह नोट करते हुए कि हालांकि एक औपचारिक प्रमाण उपलब्ध नहीं है, लेकिन शास्त्रीय एल्गोरिदम के साथ संरचनात्मक समानता और संख्यात्मक साक्ष्य दृढ़ता से अभिसरण गुणों के संरक्षण का सुझाव देते हैं। वे यह भी स्वीकार करते हैं कि प्रयोगात्मक प्लेटफार्मों में प्रतिबंधित कनेक्टिविटी पूर्ण विकेंद्रीकरण को रोकने के लिए नियंत्रण-सिद्धांत संश्लेषण (control-theoretic synthesis) की आवश्यकता हो सकती है।