Matrix Product Operator Encodings of the Magnus Expansion and Dyson Series

तकनीकी सारांश: मैग्नस एक्सपेंशन और डायसन सीरीज़ के मैट्रिक्स प्रोडक्ट ऑपरेटर एनकोडिंग्स

समस्या विवरण
समय-निर्भर हैमिल्टोनियन, $H(t)$, द्वारा नियंत्रित एक-आयामी क्वांटम लैटिस मॉडल के वास्तविक समय के विकास (real-time evolution) का अनुकरण करना, कम्प्यूटेशनल मेनी-बॉडी फिजिक्स में एक महत्वपूर्ण चुनौती बनी हुई है। हालांकि समय-विकास ऑपरेटर $U(t, t_0)$ को औपचारिक रूप से एक टाइम-ऑर्डर्ड एक्सपोनेंशियल के माध्यम से परिभाषित किया जाता है, लेकिन इंटरैक्टिंग सिस्टम्स के लिए इसका विश्लेषणात्मक मूल्यांकन सामान्यतः असंभव है। मानक संख्यात्मक दृष्टिकोण अक्सर समय-विकास को छोटे चरणों में विभाजित करने पर निर्भर करते हैं जहाँ हैमिल्टोनियन को स्थिर माना जाता है (ट्रोटर-सुजुकी अपघटन) या टेन्सर नेटवर्क विधियों जैसे कि टाइम-डिपेंडेंट वेरिएशनल प्रिंसिपल (TDVP) या कायलर-आधारित (Krylov-based) दृष्टिकोणों के माध्यम से विकास ऑपरेटर का अनुमान लगाते हैं। हालाँकि, तेजी से उतार-चढ़ाव वाले हैमिल्टोनियन के लिए, इन विधियों को अत्यंत छोटे समय चरणों की आवश्यकता होती है, जिससे उच्च कम्प्यूटेशनल लागत आती है। इसके अतिरिक्त, समय-स्वतंत्र हैमिल्टोनियन के लिए टेलर सीरीज़ पर आधारित मौजूदा उच्च-क्रम टेन्सर नेटवर्क विधियाँ, बिना अत्यधिक बॉन्ड डायमेंशन (bond dimensions) के बोझ या थर्मोडायनामिक लिमिट में साइज-एक्सटेंसिविटी (size-extensivity) खोए, सीधे तौर पर समय-निर्भर मामलों में सामान्यीकृत नहीं होती हैं।

कार्यप्रणाली (Methodology)
लेखक एक मैट्रिक्स प्रोडक्ट ऑपरेटर (MPO) एनकोडिंग प्रस्तुत करते हैं जो विशेष रूप से समय-निर्भर हैमिल्टोनियन के लिए मैग्नस एक्सपेंशन (Magnus expansion) और डायसन सीरीज़ (Dyson series) दोनों के लिए तैयार की गई है। यह कार्य पिछले शोध (संदर्भ 22) पर आधारित है जिसने समय-स्वतंत्र हैमिल्टोनियन के टेलर सीरीज़ को साइज-एक्सटेंसिव MPOs के रूप में सफलतापूर्वक एनकोड किया था।

कार्यप्रणाली दो मुख्य मार्गों में आगे बढ़ती है:

1. MPO के रूप में मैग्नस एक्सपेंशन: समय-विकास ऑपरेटर को $U(t, t_0) = e^{\Omega(t, t_0)}$ के रूप में व्यक्त किया जाता है, जहाँ $\Omega$ विभिन्न समयों पर हैमिल्टोनियन के कम्यूटेटर्स (commutators) का विस्तार है। लेखक प्रदर्शित करते हैं कि प्रथम-डिग्री MPOs (जो स्थानीय हैमिल्टोनियन का प्रतिनिधित्व करते हैं) के कम्यूटेटर्स को स्वयं प्रथम-डिग्री MPOs के रूप में निर्मित किया जा सकता है। यह मैग्नस ऑपरेटर $\Omega$ को प्रथम-डिग्री MPOs के रैखिक संयोजन के रूप में बनाने की अनुमति देता है। इसके बाद, $e^\Omega$ की गणना स्थापित टेलर सीरीज़ MPO निर्माण का उपयोग करके की जाती है।
2. MPO के रूप में डायसन सीरीज़: लेखक डायसन सीरीज़, $U(t, t_0) = \mathcal{T} \exp(-i \int H(t') dt')$ को सीधे एक MPO के रूप में एनकोड करने का एक अधिक प्रत्यक्ष प्रस्ताव देते हैं। वे समय-निर्भर हैमकिलोनियन को ड्राइविंग फंक्शन्स $f_a(t)$ द्वारा भारित (weighted) समय-स्वतंत्र स्थानीय ऑपरेटरों $H^{(a)}$ के योग में विघटित करते हैं।
   • रीवायरिंग (Rewiring): $H(t)$ के MPO प्रतिनिधित्व की "रीवायरिंग" एक प्रमुख नवाचार है। एक एकल अपर-ट्राइएंगुलर संरचना के बजाय, लेखक प्रत्येक ड्राइविंग चैनल $a$ के लिए विशिष्ट वर्चुअल लेवल्स (जिन्हें $3a$ लेबल किया गया है) पेश करते हैं। यह MPO को ड्राइविंग फंक्शन्स के विभिन्न क्रमों (permutations) के बीच अंतर करने की अनुमति देता है।
   • एक्सटेंसिव MPO में रूपांतरण: टेलर सीरीज़ निर्माण के समान, लेखक प्रथम-डिग्री MPO (जो सिस्टम साइज के साथ रैखिक रूप से स्केल करता है) को एक साइज-एक्सटेंसिव MPO में परिवर्तित करते हैं। इसमें नए $3a$ लेवल्स से वापस बेस लेवल (1) पर तीरों (arrows) को रूट करना और उपयुक्त समय-क्रमित इंटीग्रल्स (time-ordered integrals), जिन्हें $[f_{a_1} \dots f_{a_k}]$ द्वारा दर्शाया गया है, से गुणा करना शामिल है।
   • उच्च-क्रम निर्माण: $N$-वें ऑर्डर डायसन MPOs का निर्माण करने के लिए एक एल्गोरिद्मिक प्रक्रिया (एल्गोरिदम 2) प्रदान की गई है, जो हैमिल्टोनियन पदों के गैर-विजातीय उत्पादों (non-disjoint products) की पहचान करती है और उन्हें सही समय-क्रमित इंटीग्रल्स द्वारा भारित करती है।

संपीड़न तकनीकें (Compression Techniques)
उच्च-क्रम एक्सपेंशन में निहित तेजी से बढ़ते बॉन्ड डायमेंशन को प्रबंधित करने के लिए, पेपर दो संपीड़न रणनीतियों का विवरण देता है:

• सटीक कॉलम संपीड़न (Exact Column Compression): वे लेवल्स जो ऑपरेटर अनुप्रयोग के समान "इतिहास" (अर्थात, '1' इंडेक्स को हटाने के बाद जिनके लेबल समान हैं) साझा करते हैं, उन्हें समतुल्य के रूप में पहचाना जाता है। उनके संगत पंक्तियों (rows) को जोड़ा जाता है, और अनावश्यक कॉलमों को हटा दिया जाता है।
• अनुमानित रो संपीड़न (Approximate Row Compression): यह तकनीक उन लेवल्स की पहचान करती है जिनका अनुसरण समान ऑपरेटर स्ट्रिंग्स द्वारा किया जाता है लेकिन अलग-अलग गुणांकों के साथ। एक राइट-हाफ ऑपरेटर स्पेस के लिए आधार (basis) का निर्माण करके और एक रैंक-रिवीलिंग QR डीकंपोज़िशन करके, उन लेवल्स को जो (क्रम $N$ तक) वास्तव में नए ऑपरेटर पेश नहीं करते हैं, रखे गए लेवल्स के रैखिक संयोजन के रूप में व्यक्त किया जाता है। यह चरण अनुमानित है लेकिन केवल $N$ से उच्च क्रमों पर त्रुटियां पेश करता है, जिससे $N$-वें क्रम के विस्तार की सटीकता सुरक्षित रहती है जबकि बॉन्ड डायमेंशन में काफी कमी आती है।

परिणाम और बेंचमार्क
लेखक समय-मॉड्यूलेटेड हैमिल्टोनियन के साथ स्पिन-1/2 चेन्स पर डायसन MPO निर्माण का बेंचमार्किंग करते हैं:

• परिमित सिस्टम (Finite Systems): एक 8-साइट सिस्टम पर क्वासी-एग्जैक्ट सिमुलेशन (उच्च-क्रम वर्नर इंटीग्रेशन का उपयोग करके) के विरुद्ध तुलना दिखाती है कि त्रुटि $dt^N$ के रूप में $O(dt^N)$ के रूप में स्केल करती है, जो $N$-वें ऑर्डर डायसन MPO के लिए सैद्धांतिक अभिसरण क्रम (convergence order) की पुष्टि करती है।
• अनंत सिस्टम (Infinite Systems): समय-मॉड्यूलेटेड अनिसोट्रॉपी वाले एक अनंत हीजेसबर्ग XXZ चेन पर मानक TDVP सिमुलेशन (अत्यंत छोटे समय चरण $dt=10^{-7}$ के साथ) के विरुद्ध तुलना की गई। डायसन MPO दृष्टिकोण ने काफी बड़े समय चरणों के साथ तुलनीय सटीकता प्राप्त की।
• रनटाइम दक्षता: एक रनटाइम-सटीकता ट्रेड-ऑफ विश्लेषण प्रदर्शित करता है कि प्रति टाइम स्टेप बढ़ी हुई सटीकता एक अधिक अनुकूल कम्प्यूटेशनल लागत की ओर ले जाती है। एक निश्चित लक्षित सटीकता के लिए, डायसन MPO विधि को मानक TDVP दृष्टिकोणों की तुलना में कम रनटाइम की आवश्यकता होती है।

महत्व और दावे
यह पेपर एक व्यवस्थित ढांचे के रूप में समय-निर्भर क्वांटम डायनेमिक्स के अनुकरण का दावा करता है जो है:

• मनमाना सटीक (Arbitrarily Accurate): समय चरण के साथ त्रुटि स्केलिंग को कम करने के लिए ट्रंकेशन के क्रम को बढ़ाया जा सकता है।
• साइज-एक्सटेंसिव (Size-Extensive): निर्माण सीधे थर्मोडायनामिक लिमिट में मान्य है, जो हैमिल्टोनियन की सहज जोड़ (naive addition) के विपरीत है।
• बहुमुखी (Versatile): यह लघु- और दीर्घ-रेंज इंटरैक्शन को संभालता है और परिमित एवं अनंत दोनों सिस्टम्स पर लागू होता है।
• क्वांटम सिमुलेशन के लिए लागू: लेखक सुझाव देते हैं कि यह विधि समय-निर्भर हैमिल्टोनियन के लिए कुशल सर्किट के बेंचमार्किंग और अनुकूलन (optimization) के लिए एक प्राकृतिक ढांचा प्रदान करती है, जो कुशल सर्किट के संकलन (compilation) का मार्गदर्शन कर सकती है।

यह कार्य डायसन सीरीज़ MPO को मैग्नस दृष्टिकोण (दो-चरणीय प्रक्रिया) की तुलना में समय-निर्भर सिस्टम के लिए एक अधिक स्वाभाविक और कुशल सूत्रीकरण के रूप में स्थापित करता है, हालांकि मैग्नस दृष्टिकोण फ्लुओकेट डायनेमिक्स (Floquet dynamics) के अध्ययन के लिए एक शक्तिशाली उपकरण बना हुआ है। लेखक नोट करते हैं कि भविष्य का कार्य उच्च आयामों के लिए प्रोजेक्टेड एंटैंगल्ड पेयर ऑपरेटर्स (PEPOs) के माध्यम से सामान्यीकरण और आवर्ती रूप से संचालित (periodically driven) सिस्टम्स पर अनुप्रयोगों की खोज कर सकता है।