Point-group symmetry analysis of many-electron wavefunctions on a quantum computer

तकनीकी सारांश: क्वांटम कंप्यूटर पर अनेक-इलेक्ट्रॉन तरंग फलनों (wavefunctions) का बिंदु-समूह समरूपता विश्लेषण

समस्या विवरण
बिंदु-समूह (Point groups), जो आणविक प्रणालियों में स्थानिक समरूपता संचालनों (spatial symmetry operations) का प्रतिनिधित्व करते हैं, आणविक कक्षकों और स्पेक्ट्रोस्कोपी के विश्लेषण के लिए रसायन विज्ञान में मौलिक उपकरण हैं। जबकि क्वांटम एल्गोरिदम इन समरूपताओं का लाभ उठाने के लिए प्रस्तावित किए गए हैं—जैसे कि क्यूबिट संख्या को कम करने या सर्किट एंसेट्स (ansätze) को सरल बनाने के लिए—यथार्थवादी अनेक-इलेक्ट्रॉन तरंग फलनों के कार्यान्वयन और विस्तृत समरूपता विश्लेषण के व्यावहारिक तरीके अभी भी काफी हद तक अनछुए हैं। विशेष रूप से, निकट-अवधि (near-term) क्वांटम हार्डवेयर पर समरूपता संचालनों को मूर्त रूप से लागू करने और समरूपता गुणों (जैसे कि अपरिमेय निरूपण भार/irreducible representation weights) के विश्लेषण के लिए विधियों का अभाव है।

कार्यप्रणाली
लेखक अनेक-इलेक्ट्रॉन अवस्थाओं के बिंदु-समूह समरूपता का विश्लेषण करने के लिए एक अनैंसिल-मुक्त (ancilla-free), क्वांटम-शास्त्रीय हाइब्रिड विधि प्रस्तावित करते हैं, जो अबबेलियन (abelian) और गैर-अबेलियन (non-abelian) दोनों समूहों के लिए लागू है। विधि का मुख्य भाग दिए गए तरंग फलन $|\Psi\rangle$ के भीतर अपरिमेय निरूपणों ($\Gamma$) के "भार" ($w_\Gamma$) की गणना करना है। इस भार को $\Gamma$ के उप-स्थान (subspace) पर $|\Psi\rangle$ के प्रक्षेप (projection) के वर्ग गुणांक के रूप में परिभाषित किया गया है।

सैद्धांतिक ढांचा प्रोजेक्शन ऑपरेटर सूत्र पर आधारित है:
$$w_\Gamma = \frac{d_\Gamma}{|G|} \sum_C \chi^*_\Gamma(C) \sum_{g \in C} \langle\Psi|\hat{g}|\Psi\rangle$$
जहाँ $d_\Gamma$ निरूपण का आयाम है, $|G|$ समूह का क्रम है, $\chi_\Gamma(C)$ संयुग्मी वर्ग (conjugacy class) $C$ का अभिलक्षण (character) है, और $\hat{g}$ समरूपता संचालक है।

क्वांटम कंप्यूटर पर अनैंसिल क्यूबिट्स के बिना $\langle\Psi|\hat{g}|\Psi\rangle$ का मूल्यांकन करने के लिए, लेखक निरूपण मैट्रिसेस के विकर्णीकरण (diagonalization) का उपयोग करते हैं। समरूपता-अनुकूलित कक्षकों (जैसे कि हार्ट्री-फॉक) के लिए, समरूपता संचालक $U(g)$ को कक्षकीय रोटेशन (orbital rotations) में विभाजित किया जाता है। अपसरित कक्षकीय सेटों (degenerate orbital sets) के लिए निरूपण मैट्रिक्स $\bar{D}_\gamma(g)$ के विकर्णीकरण द्वारा, यूनिटरी $U(g)$ को एक आधार रोटेशन $\tilde{U}(g)$ और एक विकर्णीय चरण ऑपरेटर (diagonal phase operator) के गुणनफल के रूप में पुनर्गठित किया जाता है। अपेक्षित मान (expectation value) तब प्राप्त किया जाता है जब अवस्था $|\tilde{\Psi}(g)\rangle = \tilde{U}(g)|\Psi\rangle$ को तैयार किया जाता है, कंप्यूटेशनल आधार में मापा जाता है, और उपयुक्त चरण कारकों के साथ संयोजित किया जाता है। यह दृष्टिकोण हैडामार्ड टेस्ट (Hadamard tests) या लीनियर कॉम्बिनेशन ऑफ यूनिटरीज (LCU) तकनीकों की तुलना में गहरे सर्किट और अनैंसिल आवश्यकताओं से बचता है।

अध्ययन दो मुख्य सत्यापन रणनीतियों का उपयोग करता है:

1. शास्त्रीय संख्यात्मक सिमुलेशन: ffsim और pyscf पुस्तकालयों का उपयोग करते हुए, विधि को बेंजीन ($D_{6h}$) और स्टैगर्डड फेरोसीन ($D_{5d}$) के लिए लागू किया गया है। लेखक एकल स्लेटर निर्धारक (हार्ट्री-फॉक ग्राउंड और सिंगल-एक्साइटेड अवस्थाओं) और सहसंबद्ध तरंग फलनों (यूनिटरी क्लस्टर जैस्ट्रौ और लोकल यूनिटरी क्लस्टर जैस्ट्रौ फलन) का विश्लेषण करते हैं।
2. हार्डवेयर प्रदर्शन: विधि को IBM के ibm_kawasaki सुपरकंडक्टिंग क्वांटम डिवाइस (32 क्यूबिट तक) पर निष्पादित किया गया। लक्षित प्रणाली $D_{2h}$ समरूपता (जो $D_{6h}$ का एक उप-समूह है) वाली बेंजीन है। परीक्षण तरंग फलनों (ग्राउंड और फर्स्ट एक्साइटेड अवस्था) को डेंसिटी मैट्रिक्स रेनोर्मालाइजेशन ग्रुप (DMRG) गणनाओं का उपयोग करके तैयार किया गया, जिसे टेंसर-नेटवर्क तकनीकों के माध्यम से एक ब्रिक-वॉल एंसेट (brick-wall ansatz) में संकुचित किया गया, और फिर क्वांटम सर्किटों पर मैप किया गया। शोर (noise) को संबोधित करने के लिए, दो त्रुटि शमन (error mitigation) तकनीकों की तुलना की गई: गेट-फोल्डिंग के माध्यम से ज़ीरो नॉइज़ एक्सट्रपलेशन (ZNE) और एक क्वासी-प्रोबेबिलिस्टिक एरर मिटिगेशन दृष्टिकोण (QESEM)।

मुख्य परिणाम

• शास्त्रीय सिमुलेशन:
  • बेंजीन और फेरोसीन के हार्ट्री-फॉक ग्राउंड स्टेट्स के लिए, विधि पूर्णतः सममित अपरिमेय निरूपण ($A_{1g}$) के भार को लगभग 1.0 के रूप में सही ढंग से पहचानती है, जो क्लोज्ड-शेल विन्यास के अनुरूप है।
  • सिंगल-एक्साइटेड (SE) अवस्थाओं के लिए, विधि तरंग फलनों को कई अपरिमेय निरूपणों में सफलतापूर्वक विभाजित करती है। उदाहरण के लिए, बेंजीन में एक विशिष्ट SE अवस्था को $\approx 0.005 B_{1u} + 0.495 B_{2u} + 0.5 E_{1u}$ में विभाजित किया गया। सभी स्वतंत्र SE विन्यासों पर भार को जोड़ने पर, समूह सिद्धांत न्यूनीकरण सूत्रों के अनुरूप अपरिवर्तनीय परिणाम प्राप्त हुए।
  • सहसंबद्ध LUCJ तरंग फलनों के विश्लेषण में, विधि ने दिखाया कि ऊर्जा त्रुटि को न्यूनतम करने वाली अनुकूलन प्रक्रियाएं (बिना रेगुलेराइजेशन के) समरूपता भार को महत्वपूर्ण रूप से खराब कर सकती हैं (सही $A_{1g}$ भार से विचलन), जबकि रेगुलेराइजेशन ने समरूपता गुणों को बनाए रखने में मदद की, हालांकि इसमें ऊर्जा के अलग ट्रेड-ऑफ थे। यह सुझाव देता है कि यह मीट्रिक परीक्षण तरंग फलन की गुणवत्ता का आकलन करने के लिए उपयोगी है।
• हार्डवेयर प्रदर्शन:
  • ibm_kawasaki डिवाइस पर, कच्चे मापन (raw measurements) ने दिखाया कि क्यूबिट की संख्या बढ़ने के साथ समरूपता भार में महत्वपूर्ण गिरावट आती है (उदाहरण के लिए, बिना मिटिगेशन के 32 क्यूबिट के लिए ग्राउंड स्टेट के लिए $A_g$ भार गिरकर ~0.54 हो गया)।
  • त्रुटि शमन (error mitigation) लागू करना महत्वपूर्ण था। जबकि ZNE ने परिणामों में सुधार किया, क्वासी-प्रोबेबिलिस्टिक QESEM दृष्टिकोण ने उच्चतम निष्ठा (fidelity) प्रदान की, जिसने कुछ प्रतिशत त्रुटि के भीतर सही अपरिमेय निरूपण भार को पुनरुत्पादित किया (उदाहरण के लिए, 32 क्यूबिट पर ग्राउंड स्टेट के लिए $0.968 \pm 0.010$, जबकि सटीक मान $0.999$ है)।
  • प्रयोग ने पुष्टि की कि प्रस्तावित अनैंसिल-मुक्त विधि, टेंसर-नेटवर्क अवस्था तैयारी और उन्नत त्रुटि शमन के साथ मिलकर, वर्तमान हार्डवेयर पर समरूपता गुणों को सफलतापूर्वक पुनरुत्पादित कर सकती है।

महत्व और दावे
यह शोध पत्र सैद्धांतिक समरूपता प्रक्षेप (symmetry projection) और क्वांटम हार्डवेयर पर इसके व्यावहारिक कार्यान्वयन के बीच के अंतर को पाटने का दावा करता है। प्रस्तावित विधि निम्नलिखित के लिए एक व्यावहारिक उपकरण के रूप में कार्य करती है:

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

लेखक इस बात पर जोर देते हैं कि हालांकि हार्डवेयर प्रदर्शन में एक अबेलियन उप-समूह ($D_{2h}$) और पाउली-आधारित मूल्यांकन का उपयोग किया गया था, लेकिन अंतर्निहित कार्यप्रणाली सामान्य है और गैर-अबेलियन समूहों और मनमाने आधार फलनों (arbitrary basis functions) के लिए लागू है। यह कार्य सुझाव देता है कि जटिल सामग्रियों के भविष्य के क्वांटम सिमुलेशन के लिए समरूपता विश्लेषण एक व्यवहार्य और आवश्यक घटक है।