Spectral functions on a quantum computer through system-environment interaction

بيان المشكلة
تُعد الدوال الطيفية، التي تُقاس عبر مطيافية الانبعاث الضوئي بزاوية محددة (ARPES)، بالغة الأهمية لتوضيح بنية النطاق للمواد. ومن الناحية النظرية، تقابل هذه الدوال في الحالة التوازنية الدوال الديناميكية أحادية النقطة. ورغم صعوبة حسابها كلاسيكيًا، إلا أن قياسها على الحواسيب الكمومية يواجه عقبات كبيرة. فالأساليب القياسية، التي تتضمن حساب الارتباطات الديناميكية $\langle c^\dagger_j(t)c_0(0)\rangle$ وإجراء تحويلات فورييه، تعاني من عبء كبير في أخذ العينات. وتحديدًا، نظرًا لأن المؤثرات $c_j + c^\dagger_j$ لمواقع مختلفة $j$ لا تتبادل، فلا يمكن قياس سوى ارتباط واحد لكل عملية تصوير (shot). ولإعادة بناء الدالة الطيفية لجميع قيم $N$ من الزخم، يتطلب ذلك عبئًا في أخذ العينات بمقدار $O(N)$. وهذا الأمر يمثل مشكلة خاصة في حواسيب المصيدة الأيونية (ion-trap)، التي تمتلك دقة بوابات واتصالًا عاليين ولكنها تعاني من بطء نسبي في أوقات التصوير، مما يجعل عبء $O(N)$ في عدد مرات التصوير عائقًا أمام التطبيقات ذات المنفعة الواسعة.

المنهجية
يقترح المؤلفون طريقة لقياس الدوال الطيفية $A(k, \omega)$ بشكل مباشر عن طريق نمذجة التفاعل الفيزيائي بين النظام والبيئة كما يحدث في تجارب ARPES. ويتضمن هذا النهج ما يلي:

1. الاقتران بين النظام والبيئة: يتم اقتران نظام الاهتمام (الهاملتوني $H_{sys}$) ببيئة مكونة من $N$ من المواقع (الهاملتوني $H_{env}$) عبر حد تفاعل $H_{int}$. الهاملتوني الكلي هو $H = H_{sys} + \epsilon H_{int} + \omega H_{env}$.
2. تحضير الحالة والتطور: يتم تهيئة النظام في حالة ذاتية $|E\rangle$ للهاملتوني $H_{sys}$، بينما تُهيأ البيئة في حالة الفراغ $|0\rangle$ (لـ $A_+$) أو حالة الامتلاء الكامل $|all\rangle$ (لـ $A_-$). يتطور النظام الكلي تحت تأثير $H$ لزمن ثابت $t$.
3. القياس: بدلاً من قياس الارتباطات الديناميكية، يقوم البروتوكول بقياس عدد إشغال الزخم $n(k)$ في البيئة. يتطلب هذا إجراء تحويل فورييه فيرميوني (FFT) على الكيوبتات البيئية في نهاية عملية التطور، متبوعًا بقياس مصفوفات Pauli $Z$ المحلية.
4. النتيجة النظرية: في حد الاقتران الضعيف ($\epsilon \to 0$) والزمن الطويل ($t \to \infty$)، يتناسب القيمة المتوقعة $\langle n(k) \rangle$ مع الدالة الطيفية $A(k, \omega)$ (المدمجة مع نواة تقترب من دالة دلتا مع $t \to \infty$).
5. تنفيذ FFT الفيرميوني: للتعامل مع أحجام الأنظمة التعسفية (تحديدًا $N=27$)، طور المؤلفون تفكيك بوابة فعال لـ FFT فيرميوني من رتبة-$n$ عامة. وقد استخدموا هيكلًا تكراريًا يتضمن عمليات "التداخل" (interleave) لجعل الأنماط الفيرميونية غير المتجاورة متجاورة من أجل التحويل، مع تحسينها باستخدام تقنيات تقليص الرسم البيفي (graph decimation) لجميع البنى المتصلة ببعضها (all-to-all).

المساهمات الرئيسية

• كفاءة أخذ العينات: المساهمة الأساسية هي تقليل عبء أخذ العينات. من خلال قياس عدد إشغال الزخم للبيئة، يحصل هذا الأسلوب على البيانات الطيفية لجميع قيم $N$ من الزخم $k$ في عملية تصوير واحدة، مما يقلل عدد مرات التصوير المطلوبة بمقدار $O(N)$ مقارنة بأساليب الارتباط الديناميكي القياسية.
• تحسين وقت التشغيل: يترجم هذا الاختزال في عدد مرات التصوير إلى تحسين بمقدار $O(N)$ في وقت التشغيل الإجمالي، وهو أمر بالغ الأهمية لمنصات مثل حواسيب المصيدة الأيونية حيث يكون وقت التصوير هو العامل المحدد.
• المتانة تجاه خطأ Trotter: يجادل المؤلفون بأن نهجهم أكثر متانة تجاه أخطاء التطور الزمني (Trotterization). فبينما قد تعطي الطرق القياسية قيمًا سالبة غير فيزيائية للدوال الطيفية بسبب الأخطاء، فإن الطريقة المقترحة تقيس عدد الإشغال، وهو كمية تظل موجبة حتى مع التطور غير المثالي.
• اختبار الأداء على الأجهزة: قام المؤلفون بتنفيذ FFT فيرميوني من الرتبة-3 على 27 كيوبت واختبروا استراتيجيات تجميع مختلفة (CZ محلية، سلالم CX، وتقليص الرسم البيفي) على أنظمة Quantinuum H2-1 و H2-2 الأيونية. وحددوا أن التجميع ذو المقياس $\log_2(N)$ يوفر أدنى طاقة وأقل معدل عدم دقة (infidelity).
• عرض الخوارزمية الكاملة: تم عرض الخوارزمية الكاملة على نظام مكون من 54 كيوبت (27 موقعًا للنظام + 27 موقعًا للبيئة) على جهاز Quantinuum H2-2، حيث نجحت في حساب الدالة الطيفية لهاملتوني فيرميوني حر أحادي البعد.

النتائج

• أداء الأجهزة: نجحت الطريقة على جهاز H2-2 في إعادة إنتاج شكل نطاق الطيف لنظام مكون من 27 موقعًا. وبينما تسبب ضجيج الأجهزة في تخميد عام للإشارة (أرضية ضجيج بالقرب من 0.5)، ظلت السمات الفيزيائية، مثل مواضع النطاقات وعروضها، غير متحيزة.
• تحليل الخطأ: بالمقارنة مع المحاكاة الدقيقة، أظهرت الطريقة المقترحة أنها تحقق نفس الدقة بعدد أقل بكيد من خطوات Trotter مقارنة بأسلوب الارتباط الديناميكي القياسي، لا سيما عند صغر قيمة الاقتران $\epsilon$.
• تحليل التكلفة: بالنسبة لتجربة الـ 27 موقعًا المحددة، تحملت الطريقة عبئًا في بوابات الكيوبت الثنائية (بسبب البيئة وFFT) ولكنها وفرت عاملًا قدره $\sim 27$ في عدد مرات التصوير. ويقدر المؤلفون أنه لولا هذه الطريقة، لكان إجمالي تكلفة وقت التشغيل أعلى بنحو 11 مرة لنفس الدقة.
• تأثيرات الاقتران المحدود: يشير التحليل النظري ومحاكاة الفيرميونات الحرة إلى أن الاقتران المحدود $\epsilon$ يؤدي إلى تعريض النطاق وظهور "نطاقات شبحية". ومع ذلك، بالنسبة لـ $\epsilon t \lesssim \pi$، تكون هذه التأثيرات مهملة، وتُقرب الطريقة الدالة الطيفية بأمانة.

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