Efficient Operator Selection and Warm-Start Strategy for Excitations in Variational Quantum Eigensolvers

समस्या विवरण

वेरिएशनल क्वांटम आइजनसॉल्वर (VQE) क्वांटम अनेक-निकाय प्रणालियों (many-body systems) में ग्राउंड-स्टेट समस्याओं को हल करने के लिए एक प्रमुख हाइब्रिड एल्गोरिदम है, फिर भी इसे क्वांट적인 लाभ (quantum advantage) प्राप्त करने में महत्वपूर्ण बाधाओं का सामना करना पड़ता है। प्राथमिक चुनौतियाँ हैं:

• बैरन प्लेटो (Barren Plateaus): लॉस फंक्शन में घातांकीय रूप से लुप्त होते ग्रेडिएंट्स (exponentially vanishing gradients) अभिसरण (convergence) में बाधा डालते हैं।
• सर्किट डेप्थ (Circuit Depth): पारंपरिक एन्सैट्स (ansätze) जैसे कि यूनिटरी कपल्ड क्लस्टर सिंगल्स डबल्स (UCCSD), इलेक्ट्रॉनों की संख्या के साथ बहुपद (polynomially) रूप से स्केल करते हैं, जिसके परिणामस्वरूप सर्किट की गहराई वर्तमान शोर वाले मध्यवर्ती-पैमाने के क्वांटम (NISQ) उपकरणों के लिए बहुत अधिक हो जाती है।
• ऑपरेटर चयन ओवरहेड (Operator Selection Overhead): ADAPT-VQE जैसे अनुकूली (adaptive) रणनीतियाँ गहराई की समस्याओं को कम करने के लिए ऑपरेटरों को पुनरावृत्ति से जोड़ती हैं, लेकिन चयन प्रक्रिया के दौरान बड़ी संख्या में क्वांटम सर्किटों के क्रमिक मूल्यांकन के कारण उच्च कम्प्यूटेशनल लागत का सामना करती हैं।
• आरंभीकरण संवेदनशीलता (Initialization Sensitivity): मानक आरंभीकरण (जैसे, शून्य पैरामीटर) या मोलर-प्लेसेट सेकंड-ऑर्डर (MP2) थ्योरी जैसे विलेनात्मक (perturbative) तरीके अक्सर उप-इष्टतम (suboptimal) शुरुआती बिंदु प्रदान करते हैं, विशेष रूप से दृढ़ सहसंबंध (strongly correlated) वाले क्षेत्रों में।

कार्यप्रणाली

लेखक एक नवीन प्रोटोकॉल प्रस्तावित करते हैं जो कुशल, वॉर्म-स्टार्टेड (warm-started) VQE एन्सैट्स बनाने के लिए ExcitationSolve ऑप्टिमाइज़र और एनर्जी सॉर्टिंग (ES) को जोड़ता है।

1. ExcitationSolve और एनर्जी सॉर्टिंग (ES):

• ExcitationSolve: एक ऑप्टिमाइज़र जो प्रत्येक पैरामीटर के लिए कॉस्ट फंक्शन को पुनर्गठित करता है, जिससे ग्रेडिएंट डिसेंट के बिना सीधा अनुकूलन संभव होता है। यह ऊर्जा परिदृश्य (energy landscape) का कुछ बिंदुओं पर मूल्यांकन करके एक ऑपरेटर के लिए इष्टतम पैरामीटर मान की पहचान करता है।
• एनर्जी सॉर्टिंग रणनीति: ऑपरेटरों को क्रमिक रूप से चुनने के बजाय, ES एक संदर्भ अवस्था (reference state) के विरुद्ध सभी ऑपरेटरों के ऊर्जा प्रभाव ($\Delta E_j = E_{ref} - E_j$) का एक साथ मूल्यांकन करता है। वे ऑपरेटर जो एक सीमा $\epsilon_A$ से अधिक हैं, उन्हें एक ही बार में चुना जाता है।
• वॉर्म-स्टार्ट (Warm-Start): ExcitationSolve प्रत्येक चयनित ऑपरेटर के लिए इष्टतम पैरामीटर $\theta_j$ तुरंत प्रदान करता है, जिससे एन्सैट को शून्य के बजाय परिमित, गैर-शून्य मानों के साथ आरंभीकृत किया जा सकता है।

2. क्लासिकल प्री-प्रोसेसिंग और विश्लेषणात्मक व्युत्पत्ति (Analytical Derivation):

• हार्ट्री-फॉक (HF) संदर्भ अवस्था से शुरू होने वाले प्रारंभिक चयन चरण के लिए, लेखक प्रदर्शित करते हैं कि डबल एक्साइटेशन का चयन पूरी तरह से क्लासिकल रूप से किया जा सकता है।
• इलेक्ट्रॉनिक स्ट्रक्चर हैमिल्टनियन को टू-स्टेट सबस्पेस (HF स्टेट + एक डबल एक्साइटेशन) में कम करके, वे केवल एक-इलेक्ट्रॉन और दो-इलेक्ट्रॉन इंटीग्रल्स का उपयोग करके अधिकतम ऊर्जा कमी ($\Delta E_{min}$) और इष्टतम प्रारंभिक पैरामीटर ($\theta_{ES}$) की गणना करने के लिए विश्लेषणात्मक सूत्र (समीकरण 4 और 5) व्युत्पन्न करते हैं। यह ऑपरेटर चयन के दौरान क्वांटम संसाधनों की आवश्यकता को समाप्त करता है।

3. OVP-CEOs के साथ एकीकरण:

• इस पद्धति को वन-वेरिएशनल-पैरामीटर कपल एक्सचेंज ऑपरेटर्स (OVP-CEOs) के लिए अनुकूलित किया गया है, जो मानक क्वबिट एक्साइटेशन की तुलना में CNOT गेट काउंट (13 से 9) और सर्किट डेप्थ (11 से 7) को कम करते हैं।
• रेडंडेंसी (जहाँ OVP-CEOs के $+$ और $-$ दोनों वेरिएंट HF स्टेट पर समान रूप से कार्य करते हैं) और ऑर्बिटल मिक्सिंग के मुद्दों को संबोधित करने के लिए, एक अतिरिक्त चयन चरण लागू किया गया है। यह चरण विकसित होती अवस्था पर प्रत्येक ऑपरेटर जोड़ी के प्रभाव के आधार पर इष्टतम वेरिएंट ($+$ या $-$) को निर्धारित करता है, जिससे सर्किट डेप्थ को दोगुना किए बिना अभिसरण सुनिश्चित होता है।

4. पर्टर्बेशन थ्योरी के साथ सैद्धांतिक संबंध:

• लेखक अपने दृष्टिकोण को सेकंड-ऑर्डर एपस्टीन-नेस्बेट (EN2) पर्टर्बेशन थ्योरी से सैद्धांतिक रूप से जोड़ते हैं। वे दिखाते हैं कि EN2 एम्पलीट्यूड्स, एक सहसंबंध माप ($b/a$) के संदर्भ में ES एम्पलीट्यूड्स के प्रथम-क्रम टेलर विस्तार (first-order Taylor expansion) के अनुरूप हैं।
• यह स्थापित करता है कि ES और EN2 कमजोर सहसंबंध वाले क्षेत्रों में लगभग समान हैं, लेकिन ES दृढ़ सहसंबंध वाले क्षेत्रों में भी मजबूत रहता है जहाँ EN2 (और MP2) विफल हो सकते हैं।

मुख्य योगदान

• सिंगल-स्वीप ऑपरेटर चयन: एक प्रोटोकॉल जो ऑपरेटर पूल पर एक एकल स्वीप के माध्यम से एक अनुमानित ग्राउंड स्टेट एन्सैट बनाता है, जिससे चयन के लिए आवश्यक क्वांटम सर्किट मूल्यांकन की संख्या में भारी कमी आती है।
• क्लासिकल वॉर्म-स्टार्ट: विश्लेषणात्मक सूत्रों की व्युत्पत्ति जो प्रासंगिक डबल एक्साइटेशन के चयन और उनके पैरामीटर्स के आरंभीकरण को केवल क्लासिकल कंप्यूटेशन का उपयोग करके सक्षम बनाती है, जो क्वांटम ऑप्टिमाइज़र के लिए एक उच्च-गुणवत्ता वाला शुरुआती बिंदु प्रदान करती है।
• क्वाड्रेटिक स्पीडअप: अनुभवजन्य प्रदर्शन (empirical demonstration) कि प्रस्तावित विधि (ExcitationSolve + ES) पूल के आकार के संबंध में मानक ADAPT-VQE (ग्रेडिएंट डिसेंट + ग्रेडिएंट सिलेक्शन) और बिना सॉर्टिंग वाले एडेप्टिव ExcitationSolve की तुलना में क्वाड्रेटिक अभिसरण गति प्राप्त करती है।
• OVP-CEO एकीकरण: एक सफल अनुकूलन जो सर्किट डेप्थ को कम करने के साथ-साथ ऑपरेटर चयन के लिए प्रबंधनीय कम्प्यूटेशनल ओवरहेड को संतुलित करता है।

परिणाम

इस पद्धति का परीक्षण नॉइलेस स्टेट-वेक्टर सिम्युलेटर्स का उपयोग करके 12 अणुओं (4 से 20 क्विबिट्स तक, जैसे H2, LiH, CH4, C2) पर बेंचमार्क किया गया।

• दक्षता (Efficiency): बड़े ऑपरेटर पूल वाले अणुओं (जैसे 804 ऑपरेटरों वाला CH4) के लिए, ExcitationSolve + ES विधि ने शुद्ध ExcitationSolve की तुलना में कम्प्यूटेशनल संसाधनों को $10^{-2}$ के कारक से कम कर दिया।
• अभिसरण (Convergence): यह दृष्टिकोण केवल प्रासंगिक ऑपरेटरों वाले कॉम्पैक्ट एन्सैट्स बनाता है। LiH के लिए, ES-निर्मित एन्सैट ने फिक्स्ड UCCSD एन्सैट के 96 ऑपरेटरों की तुलना में केवल 34 ऑपरेटरों का उपयोग किया, जिससे तेज़ अभिसरण हुआ और ऑप्टिमिज़ेशन प्लेटो (plateaus) से बचा जा सका।
• आरंभीकरण की गुणवत्ता: दृढ़ सहसंबंध वाले सिस्टम में (बॉन्ड लेंथ बढ़ाकर सिम्युलेट किया गया), ES ने MP2 और EN2 की तुलना में बेहतर शुरुआती अनुमान प्रदान किया। जहाँ MP2 लगातार सबसे खराब प्रदर्शन करता रहा, वहीं ES ने उस मजबूती को बनाए रखा जहाँ MP2 और EN2 सटीक समाधान से विचलित हो गए थे।
• उच्च-क्रम एक्साइटेशन (Higher-Order Excitations): पद्धति ने LiH के लिए ट्रिपल एक्साइटेशन को सफलतापूर्वक शामिल किया, जिससे केवल डबल एक्साइटेशन की तुलना में दो ऑर्डर बेहतर सटीकता प्राप्त हुई, और इसमें मूल्यांकन लागत में केवल मामूली वृद्धि हुई।
• OVP-CEO प्रदर्शन: प्रस्तावित चयन मानदंड के साथ OVP-CEOs का उपयोग करने से मानक एक्साइटेशन के समान ग्राउंड स्टेट एनर्जी तक अभिसरण हुआ, लेकिन कम सर्किट डेप्थ (कम CNOTs) के साथ, हालांकि पूल के बड़े आकार के कारण ऊर्जा मूल्यांकन की संख्या में थोड़ी वृद्धि हुई।

महत्व और दावे

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

• स्केलेबिलिटी (Scalability): ऑपरेटर चयन की कम्प्यूटेशनल जटिलता को कम करके और क्लासिकल प्री-प्रोसेसिंग के माध्यम से "वॉर्म-स्टार्ट" को सक्षम करके, यह पद्धति बड़े अणुओं (हाल के साहित्य की विशिष्ट 12-क्विबिट सीमा से परे) के सिमुलेशन को अधिक व्यवहार्य बनाती है।
• NISQ प्रासंगिकता: OVP-CEOs का एकीकरण NISQ हार्डवेयर के लिए एक अनुकूल ट्रेड-ऑफ प्रदान करता है, जहाँ सर्किट डेप्थ प्राथमिक सीमित कारक है।
• सैद्धांतिक अंतर्दृष्टि: यह पेपर हेयुरिस्टिक ऑपरेटर चयन और स्थापित पर्टर्बेशन थ्योरी के बीच एक मौलिक संबंध स्थापित करता है, यह दिखाते हुए कि ES केवल एक कुशल हेयुरिस्टिक नहीं है बल्कि EN2 थ्योरी का एक मजबूत सामान्यीकरण है।
• व्यापक प्रयोज्यता: हालाँकि यह VQE के लिए तैयार किया गया है, लेखक सुझाव देते हैं कि क्लासिकल प्री-प्रोसेसिंग और वॉर्म-स्टार्ट रणनीति अन्य कम्प्यूटेशनल रूप से महंगे तरीकों (जैसे क्वांटम मोंटे कार्लो, QSE) और हैमिल्टोनियन वेरिएशनल विधियों को भी लाभ पहुँचा सकती है जो उच्च-गुणवत्ता वाले प्रारंभिक अवस्थाओं पर निर्भर करते हैं।

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