تقترح هذه الورقة مخططاً للتوازي للخوارزميات الكمومية التبديدية يستفيد من الموضعية الهندسية لتنفيذ مؤثرات القفز في آن واحد، مما يؤدي إلى تقليل عمق الدارة بشكل أسي وتحسين عملية تطبيق هذه الأساليب بشكل كبير للحواسيب الكمومية القريبة من المدى الحالي.
في السعي لبناء حواسيب كمومية مفيدة، يبحث العلماء باستمرار عن طرق لجعل هذه الآلات الهشة تنجز المزيد من العمل بموارد أقل. أحد الأساليب الواعدة يستعير استراتيجية من الحوسبة الكلاسيكية تُعرف بطريقة "ماركوف تشين مونت كارلو" (Markov Chain Monte Carlo)، وهي تقنية تُستخدم لحل المشكلات المعقدة من خلال محاكاة مسارات عشوائية عبر مشهد شاسع من الاحتمالات. وفي العالم الكمومي، تطورت هذه الفكرة إلى فئة من الطرق تسمى الخوارزميات التبديدية (dissipative algorithms). فبدلاً من محاولة إجبار النظام الكمومي على البقاء معزولاً بشكل مثالي، تسمح هذه الطرق للنظام بالتفاعل عمدًا مع محيطه، تمامًا مثل فنجان قهوة ساخن يبرد ليتناسب مع درجة حرارة الغرفة. ومن خلال تصميم هذا التفاعل بعناية، يستقر النظام طبيعيًا في الحالة المنشودة، مثل حالة الطاقة الأدنى لتكوين مادة ما، وهو ما غالبًا ما يكون الهدف من المحاكاة الكمومية. ومع ذلك، فإن تشغيل هذه الخوارزميات على أجهزة حقيقية كان يمثل عقبة رئيسية؛ إذ تتطلب العملية محاكاة سلسلة من الانتقالات المحددة، وإنماؤها واحدًا تلو الآخر في تسلسل صارم يجعل الدوائر الحاسوبية عميقة وبطيئة للغاية، مما يتجاوز في كثير من الأحيان قدرات الأجهزة الحالية.
لقد وجد فريق من الباحثين في جامعة ييل، وجامعة تورنتو، ومختبر المحيط الهادئ الشمال الغربي الغربي طريقة لتسريع هذه العملية بشكل كبير من خلال تغيير كيفية تنفيذ هذه الانتقالات. في عملهم، عالجوا مشكلة "عمق الدائرة" (circuit depth)، وهو ما يعبر أساسًا عن عدد الخطوات التي يجب أن يتخذها الحاسوب الكمومي على التوالي لإكمال مهمة ما. كان النهج التقليدي لهذه الخوارزميات التبديدية يتضمن تطبيق كل انتقال بشكل متسلسل، والانتظار حتى ينتهي الواحد قبل بدء التالي، مما خلق سلسلة طويلة وضيقة من العمليات التي تستغرق وقتًا طويلاً جدًا للاكتمال. أدرك الباحثون أنه نظرًا لأن التفاعلات في العديد من الأنظمة الفيزيائية هي تفاعلات محلية — أي أن الجسيم يؤثر في الغالب على جيرانه المباشرين بدلاً من البعيدين — فإنه يمكنهم تجميع هذه الانتقالات معًا. ومن خلال إثبات أن هذه الانتقالات يمكن حصرها في مناطق صغيرة ومنفصلة من المعالج الكمومي، أظهروا أنه يمكن تنفيذ العديد منها في الوقت ذاته.
أثبت الفريق أنه من خلال تشغيل هذه العمليات بالتوازي، يمكنهم تقليل الوقت المطلوب لكل خطوة من الحسابات بشكل أسي. وقد اختبروا هذه الطريقة الجديدة على نظام محاكى مكون من مائة كيوبت مرتبة في خط أحادي الأبعاد، وهو إعداد شائع لدراسة المواد المغناطيسية. في هذا الاختبار المحدد، قلل نهجهم المتوازي من عمق الدائرة المطلوبة بمعامل قدره ثلاثة وخمسون مقارنة بالطريقة التسلسلية القياسية. وهذا يعد اكتشافًا مهمًا لأنه يشير إلى أن الخوارزميات التي كانت في السابق عميقة للغاية بحيث يصعب تشغيلها على الأجهزة الكمومية القريبة من الواقع، يمكن الآن تنفيذها بسهولة أكبر. لم يكتفِ الباحثون باقتراح هذه الفكرة نظريًا فحسب، بل قدموا برهانًا رياضيًا صارمًا يوضح أن تشغيل هذه الانتقالات المحلية بالتوازي لا يخل بدقة النتيجة النهائية؛ إذ لا يزال النظام يستقر في الحالة الصحيحة بنفس الموثوقية التي يتمتع بها الإصدار الأبطأ والمتسلسل، ولكنه يصل إليها بشكل أسرع بكثير.
يعالج هذا العمل مقايضة حرجة حدت من الاستخدام العملي للخوارمايات الكمومية التبديدية. ففي السابق، كان على العلماء الاختيار بين استخدام انتقال واحد، وهو سريع لكل خطوة ولكنه يستغرق وقتًا طويلاً جدًا للوصول إلى حل، أو استخدام انتقالات عديدة في وقت واحد، والتي تتقارب بسرعة ولكنها تتطلب دائرة عميقة جدًا لدرجة يستحيل بناؤها. وتكسر الطريقة الجديدة هذا المأزق؛ فمن خلال توطين التفاعلات وتشغيلها بالتوازي، حقق الباحثون أفضل ما في العالمين: زمن تقارب سريع مع عمق دائرة يمكن التحكم فيه. وقد أكدت عمليات المحاكاة التي أجروها أن طريقة الانتقال الواحد ستتطلب عمقًا كبيرًا لدرجة يستحيل تنفيذها فعليًا، بينما يهبط النهج المتوازي الجديد بالمتطلبات إلى مستوى قابل للتنفيذ بالنسبة للحواسيب الكمومية الأولية المقاومة للأخطاء.
إن تداعيات هذا الاكتشاف فورية في مجال المحاكاة الكمومية. فمن خلال جعل هذه الخوارزميات أكثر عملية، فتح الباحثون الباب لمحاكاة الظواهر الفيزيائية المعقدة، مثل كيفية سلوك المواد عند درجات حرجات مختلفة أو كيفية وصولها إلى حالات الاستقرار الأساسية، على الأجهزة التي يتم تطويرها حاليًا. تعتمد الدراسة على التجارب العددية والبراهين الرياضية بدلاً من اختبارات الأجهزة المادية، لكن النتائج واضحة وقوية ضمن نطاق نماذجهم. لا يدعي العمل أنه حل كل المشكلات في الحوسبة الكمومية، ولكنه يوفر مسارًا ملموسًا وقابلًا للتوسع لواحد من أكثر فئات الخوارزميات الكمومية واعدًا. إنه يوضح أنه من خلال فهم الطبيعة المحلية للتفاعلات الكمومية، يمكن للعلماء إعادة هيكلة حساباتهم لتناسب القيود الفيزيائية للآلات التي يبنونها، محولين بذلك الإمكانية النظرية إلى واقع عملي.
ملخص تقني: موازاة الخوارزميات الكمومية التبديدية
بيان المشكلة تقدم الخوارزميات الكمومية التبديدية إطار عمل لإعداد الحالات الكمومية (مثل حالات غيبس الحرارية أو حالات الحالة الأرضية) عبر محاكاة الأنظمة الكمومية المفتوحة التي تحكمها صيغة ليندبلاديان (Lindbladian). هذه الأساليب مستوحاة من تقنيات ماركوف تشين مونت كارلو (MCMC) الكلاسيكية، ولكنها تواجه عقبات تنفيذية كبيرة على الأجهزة الكمومية القريبة من المدى الحالي.
تتمثل العقبة الرئيسية المحددة في هذا العمل في عمق الدائرة (circuit depth) المطلوب لمحاكاة الليندبلاديان. يتم دفع ديناميكيات الليندبلاديان بواسطة مجموع من "مؤثرات القفز" (Ka). في المقال المنهجية الحالية:
التبديد ذو القفزة الواحدة [6]: يستخدم مؤثر قفز واحد لتقليل تكلفة الكيوبت، ولكنه يعاني من زمن خلط كبير، مما يؤدي إلى تدهور في مقياس العمق الإجمالي.
التبديد الكتلي المتسلسل [7]: يطبق مؤثرات قفز متعددة لتحقيق خلط سريع (لوغاريتمي في حجم النظام). ومع ذلك، نظرًا لتطبيق هذه المؤثرات بشكل متسلسل لمشاركة كيوبت مساعد واحد، فإن عمق الدائرة لكل خطوة زمنية يزداد خطيًا مع عدد مؤثرات القفز. وهذا يخلق مقايضة حيث يتم تعويض الانخفاض في زمن الخلط بزيادة العمق لكل خطوة، مما يؤدي إلى عمق دائرة إجمالي غير مثالي (تحديدًا O~(n5) للأنظمة أحادية البعد، مقارنة بـ O~(n9) لطرق القفزة الواحدة).
المنهجية يقترح المؤلفون مخطط التبديد الكتلي المتوازي الذي يستفيد من المحلية الهندسية للهاميلتونيّات الفيزيائية لتنفيذ مؤثرات القفز في وقت واحد. تعتمد المنهجية على ثلاث ركائز:
المحلية عبر حدود ليب-روبنسون (Lieb-Robinson Bounds): يثبت البحث أنه بالنسبة للهاميلتونيّات المحلية هندسيًا، فإن تطور هايزنبرغ لمؤثر محلي Aa(t) يكون محصورًا فعليًا داخل "مخروط ضوئي". ومن خلال استخدام حدود ليب-روبنسون، يثبت المؤلفون أن كل مؤثر قفز عالمي Ka يمكن تقريبه بواسطة مؤثر محليKaΩ يعمل فقط على منطقة Ω بنصف قطر r.
من الضروري ملاحظة أن نصف القطر r المطلوب لتحقيق دقة ϵ يتناسب فقط لوغاريتميًا مع حجم النظام n والزمن t (r=Ω(log(nt/ϵ))).
تضمن هذه المحلية تنفيذ مؤثرات القفز التي تعمل في مناطق منفصلة (يفصل بينها مسافة >2r) بشكل متوازٍ دون خطأ ملموس.
التنفيذ المتوازي: بدلاً من تطبيق O(n) من مؤثرات القفز بشكل متسلسل، يقوم الخوارزمية بتقسيمها إلى O(rD) من الطبقات (حيث D هو البعد المكاني). ضمن كل طبقة، يتم تطبيق مؤثرات القفز المنفصلة في وقت واحد باستخدام كيوبتات مساعدة متميزة.
يستخدم البحث صيغة تروتر (Trotter formula) من الدرجة الأولى مدمجة مع تمدد ستينسبرج (Stinesprang dilation) لمحاكاة الديناميكيات المحلية.
يُفترض أن دالة المرشح f(t) التي تحدد مؤثرات القفز تضمحل بسرعة (الافتراض 1) لضمان تقارب حدود التكامل المطلوبة للمحلية.
تحليل التعقيد: يميز المؤلفون بين عمق الاستعلام (طبقات استدعاء الأوراكل) وعمق الدائرة (طبقات البوابات الإجمالية). من خلال الموازاة، ينخفض عمق الدائرة لكل خطوة زمنية من الخطية في عدد مؤثرات القفز (المتسلسل) إلى اللوغاريتم المتعدد في n (المتوازي).
المساهمات الرئيسية
الإثبات النظري للمحلية: يقدم البحث إثباتًا صارمًا (النظرية 1 و2) بأن توطين مؤثرات القفز إلى نصف قطر لوغاريتمي يحافظ على زمن الخلط والحالة المستقرة لديناميكيات الليندبلاديان الأصلية ضمن خطأ يمكن التحكم فيه ϵ.
مخطط الموازاة: يقدم استراتيجية موازاة محددة تستغل الفصل المكاني لمؤثرات القفز المحلية، مما يقلل عمق الدائرة لكل خطوة زمنية بشكل كبير مقارنة بالتبديد الكتلي المتسلسل.
تقدير الموارد: على عكس الطرق العامة القائمة على الأوراكل، يقوم هذا العمل ببناء مؤثرات القفز صراحةً من هاميلتونيّ معطى، مما يوفر تقديرات لتعقيد البوابات من البداية إلى النهاية مصممة لتطبيقات محددة (مثل إعداد الحالة الأرضية).
التحسين التقاربي: يظهر التحليل أنه بينما يوفر التبديد الكتلي المتسلسل تحسينًا رباعي القوة (quartic improvement) على طرق القفزة الواحدة، فإن الطريقة المتوازية المقترحة تحقق تحسينًا رباعي القوة إضافيًا، مما يقلل إجمالي عمق الدائرة من O~(n5) إلى O~(n) للأنظمة أحادية البعد (الجدول 1).
النتائج
التحقق العددي: تحقق المؤلفون من النهج على نموذج إيزينج للمجال العرضي أحادي البعد (1D TFIM) مع 100 كيوبت.
تقليل العمق: أظهرت التجارب انخفاضًا بمقدار 53 ضعفًا في عمق الاستعلام مقارنة بالأساس المرجعي للتبديد الكتلي المتسلسل.
الجدوى: تؤكد الدراسة أن طريقة القفزة الواحدة غير مجدية للأنظمة الكبيرة بسبب أزمنة الخلط المفرطة، بينما يقلل النهج المتوازي بشكل كبير من متطلبات الموارد، مما يجعل الخوارزميات التبديدية أكثر قابلية للتطبيق على الأجهزة الكمومية القريبة من المدى والقدرة على تحمل الأخطاء في مراحلها الأولى.
الأهمية والادعاءات يدعي البحث أنه من خلال معالجة عقبة التنفيذ المتسلسل عبر التوطين الهندسي، يمكن للخوارزميات الكمومية التبديدية تحقيق تحسينات متعددة الحدود في إجمالي عمق الدائرة. وهذا يجعل إعداد الحالات الحرارية والأرضية أكثر عملية للأجهزة ذات الموارد المحدودة.
يؤكد المؤلفون أن عملهم يختلف عن دراسات التوطين السابقة (مثل [10]) من حيث:
التعميم لما وراء أخذ عينات غيبس ليشمل فئة أوسع من الخوارزميات التبديدية.
التركيز تحديدًا على قدرة الموازاة والمقياس للأنظمة الكبيرة، بدلاً من مجرد الصمود أمام الضوضاء في المدى القريب.
تقديم حدود عليا نظرية لعمق الدائرة عبر طرق التربيع وتفكيك الأوراكل، بدلاً من الاعتماد على التجميع المتغير (variational compilation) الذي يفتقر إلى الضمانات النظرية للأنظمة الكبيرة.
في نهاية المطاف، يجادل العمل بأن المقايضة باستخدام المزيد قليلاً من الكيوبتات المساعدة (لتمكين الموازاة) تفوق في أهميتها التخفيض الكبير في عمق الدائرة، وهو مقياس حاسم للأجهزة المبكرة ذات القدرة على تحمل الأخطاء حيث تتراكم أخطاء البوابات مع عمق الدائرة.