Matrix Product Operator Encodings of the Magnus Expansion and Dyson Series

ملخص تقني: ترميز مؤثرات ناتج ضرب المصفوفات لتوسيع ماغنوس وسلسلة دايسون

بيان المشكلة
لا تزال محاكاة التطور الزمني الفعلي للنماذج الشبكية الكمومية أحادية البعد التي تحكمها هاملتونيانات تعتمد على الزمن، $H(t)$، تشكل تحدياً كبيراً في الفيزياء الحاسوبية متعددة الأجسام. وبينما يتم تعريف مؤثر التطور الزمني $U(t, t_0)$ رسمياً عبر أس مرتب زمنياً، فإن التقييم التحليلي مستح مستحيل عموماً للأنظمة المتفاعلة. تعتمد النهج العددية القياسية غالباً على تقسيم التطور الزمني إلى خطوات صغيرة حيث يتم التعامل مع الهاملتونياني ككمية ثابتة (تفكيك تروتر-سوزوكي) أو تقريب مؤثر التطور عبر طرق الشبكة الموترية مثل مبدأ التباين الزمني المعتمد على الوقت (TDVP) أو النهج القائمة على كرولوف (Krylov). ومع ذلك، بالنسبة للهاملتونياناتات سريعة التقلب، تتطلب هذه الطرق خطوات زمنية صغيرة للغاية، مما يؤدي إلى تكاليف حاسوبية باهظة. علاوة على ذلك، فإن طرق الشبكة الموترية عالية الرتبة الموجودة، مثل تلك القائمة على متسلسلة تايلور للهاملتونيانات المستقلة عن الزمن، لا يمكن تعميمها بسهولة على الحالات المعتمدة على الزمن دون التسبب في أبعاد رابطة (bond dimensions) باهظة أو فقدان الامتداد الحجمي (size-extensivity) في الحد الديناميكي الحراري.

المنهجية
يقدم المؤلفون ترميزاً لمؤثر ناتج ضرب المصفوفات (MPO) لكل من توسيع ماغنوس وسلسلة دايسون، مصمماً خصيصاً للهاملتونيانات المعتمدة على الزمن. يبني هذا العمل على بحث سابق (المرجع 22) الذي نجح في ترميز متسلسلة تايلور للهاملتونيانات المستقلة عن الزمن كمؤثرات MPO ذات امتداد حجمي.

تتبع المنهجية مسارين رئيسيين:

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. يقومون بتفكيك الهاملتونياني المعتمد على الزمن إلى مجموع من المؤثرات المحلية المستقلة عن الزمن $H^{(a)}$ الموزونة بدوال دفع $f_a(t)$.
   • إعادة التوصيل (Rewiring): الابتكار الرئيسي هو "إعادة توصيل" تمثيل MPO للهاملتونياني $H(t)$. فبدلاً من هيكل واحد مثلثي علوي، يقدم المؤلفون مستويات افتراضية متميزة (يُرمز لها بـ $3a$) لكل قناة دفع $a$. وهذا يسمح لـ MPO بالتمييز بين التباديل المختلفة لدوال الدفع.
   • التحويل إلى MPO ذي امتداد حجمي: على غرار بناء متسلسلة تايلور، يقوم المؤلفون بتحويل MPO من الدرجة الأولى (الذي يتناسب خطياً مع حجم النظام) إلى MPO ذي امتداد حجمي. يتضمن ذلك إعادة توجيه الأسهم من مستويات $3a$ الجديدة للعودة إلى المستوى الأساسي (1) والضرب في التكاملات المرتبة زمنياً لدوال الدفع المناسبة، والمشار إليها بـ $[f_{a_1} \dots f_{a_k}]$.
   • البناء عالي الرتبة: يتم توفير إجراء خوارزمي (الخوارزمية 2) لبناء MPOs لـ دايسون من الرتبة $N$ عبر تحديد وجمع المساهمات من المنتجات غير المنفصلة لطلبات الهاملتونياني، الموزونة بالتكاملات المرتبة زمنياً الصحيحة.

تقنيات الضغط
لإدارة البعد الرابط المتزايد بسرعة والمتأصل في التوسيعات عالية الرتبة، يفصل البحث استراتيجيتين للضغط:

• الضغط العمودي الدقيق: يتم تحديد المستويات في MPO التي تشترك في نفس "التاريخ" لتطبيق المؤثر (أي أن تسمياتها متطابقة بعد إزالة مؤشرات '1' التي تمثل الحدود التي لم تبدأ بعد) كعناصر متكافئة. يتم جمع صفوفها المقابلة، وإزالة الأعمدة الزائدة.
• الضغط الصفي التقريبي: تحدد هذه التقنية المستويات التي تتبعها نفس سلاسل المؤثرات ولكن بمعاملات مختلفة. من خلال بناء أساس لفضاء المؤثر النصف الأيمن وإجراء تحليل تفكيك QR الكاشف للرتبة، يتم التعبير عن المستويات التي لا تقدم مؤثرات جديدة حقاً (حتى الرتبة $N$) كتركيبات خطية من المستويات المحتفظ بها. هذه الخطوة تقريبية ولكنها لا تُدخل أخطاء إلا في الرتب الأعلى من $N$، مما يحافظ على دقة توسيع الرتبة $N$ مع تقليل البعد الرابط بشكل كبير.

النتائج والاختبارات المعيارية
قام المؤلفون باختبار بناء MPO لـ دايسون على سلاسل سبين-1/2 مع هاملتونيانات معدلة زمنياً:

• الأنظمة المحدودة: أظهرت المقارنات مع عمليات المحاكاة شبه الدقيقة (باستخدام تكامل فيرنر عالي الرتبة) على نظام مكون من 8 مواقع أن الخطأ يتناسب مع $O(dt^N)$ لـ MPO دايسون من الرتبة $N$، مما يؤكد رتبة التقارب النظرية.
• الأنظمة اللانهائية: تمت مقارنة محاكاة MPO لـ دايسون مع محاكاة TDVP قياسية باستخدام خطوات زمنية صغيرة للغاية ($dt=10^{-7}$) على سلسلة XXZ لـ هيزنبرغ لانهائية مع تباين معدل زمنياً. حقق نهج MPO لـ دايسون دقة مماثلة مع خطوات زمنية أكبر بكثير.
• كفاءة وقت التشغيل: يظهر تحليل المقايضة بين وقت التشغيل والدقة أن الدقة المعززة لكل خطوة زمنية تترجم إلى تكلفة حاسوبية أكثر ملاءمة. بالنسبة لدقة مستهدفة ثابتة، يتطلب طريقة MPO لـ دايسون وقتاً أقل من تشغيل نهج TDVP القياسي.

الأهمية والادعاءات
يدعي البحث تقديم إطار عمل منهجي لمحاكاة الديناميكا الكمومية المعتمدة على الزمن يتميز بكونه:

• دقيقاً بشكل تعسفي: يمكن زيادة رتبة القطع لتقليل مقياس الخطأ مع الخطوة الزمنية.
• ذا امتداد حجمي: البناء صالح مباشرة في الحد الديناميكي الحراري، على عكس الإضافات الساذجة لقوى الهاملتونياني.
• متعدد الاستخدامات: يتعامل مع التفاعلات قصيرة وطويلة المدى ويطبق على كل من الأنظمة المحدودة واللانهائية.
• قابل للتطبيق في المحاكاة الكمومية: يشير المؤلفون إلى أن الطريقة توفر إطاراً طبيعياً لمعايرة وتحسين الدوائر الكمومية للهاملتونيانات المعتمدة على الزمن، مما قد يوجه تجميع الدوائر الفعالة.

يضع هذا العمل MPO لـ دايسون كصياغة أكثر طبيعية وكفاءة للأنظمة المعتمدة على الزمن مقارنة بنهج ماغنوس ذي الخطوتين، رغم أن الأخير يظل أداة قوية لدراسة ديناميكا فلوكيه (Floquet). ويشير المؤلفون إلى أن العمل المستقبلي يمكن أن يستكشف التعميمات إلى الأبعاد الأعلى عبر مؤثرات الأزواج المتشابكة المسقطة (PEPOs) والتطبيقات في الأنظمة المدفوعة دورياً.