Green's Functions from Sample-based Krylov Quantum Diagonalization: An Impurity Solver for Dynamical Mean-Field Theory

المادة التي تتصرف بطرق غريبة وعنيدة، مثل الموصلات الفائقة التي تنقل الكهرباء دون مقاومة أو المواد المغناطيسية التي تشغل محركات أقراصنا الصلبة، غالبًا ما تتحدى القواعد البسيطة التي نستخدمها لتفسير المواد الصلبة العادية. في هذه المواد، لا تتحرك الإلكترونات بشكل مستقل؛ بل تدفع وتجذب بعضها البعض بقوة شديدة تجعل سلوك الإلكترون الواحد مرتبطًا ارتباطًا وثيقًا بسلوك جيرانه. لفهم كيفية عمل هذه المواد، يجب على العلماء حساب "دالة غرين" (Green's function)، وهي أداة رياضية تعمل كخريطة لكيفية تحرك الإلكترون عبر هذه البيئة المزدحمة والفوضوية. تكشف هذه الخريطة عن مستويات طاقة المادة وكيفية استجابتها للضوء أو الكهرباء، لكن حسابها للأنظمة الكبيرة يعد مهمة هائلة. إذ ينمو عدد الترتيبات الممكنة للإلكترونات بشكل انفجاري مع حجم النظام لدرجة أن أقوى الحواسيب الفائقة تنفد ذاكرتها بسرعة، مما يجبر الباحثين على دراسة نسخ مصغرة ومبسطة فقط من المواد التي تهمهم بالفعل.

لقد اقترح فريق من الباحثين الآن طريقة جديدة لمعالجة هذه المشكلة باستخدام الحواسيب الكمومية، ليس من خلال محاولة محاكاة النظام بأكره دفعة واحدة، بل عبر أخذ سلسلة من اللقطات الصغيرة والذكية. وفي دراسة ركزت على نظام نموذجي يحاكي سلوك الإلكترونات في معدن، قام الباحثون بتكييف تقنية تسمى "القطرية كيرلوف القائمة على العينات" (sample-based Krylov quantum diagonalization). فبدلاً من أن يطلبوا من الحاسوب الكمومي الاحتفاظ بالحالة الكاملة والضخمة للنظام في ذاكرته، طلبوا منه تطوير النظام لفترات زمنية قصيرة جدًا ثم قياس النتائج. ومن خلال تكرار هذه العملية مرات عديدة وجمع النتائج الأكثر شيوعًا، تمكنوا من إعادة بناء الخريطة الكاملة لحركة الإلكترونات دون الحاجة أبدًا إلى تخزين الكمية المستحيلة من البيانات التي تتطلبها الطرق التقليدية. وقد صُمم هذا النهج ليعمل على الحواسيب الكمومية الصاخبة وغير المثالية المتاحة اليوم، والتي لا تستطيع بعد تشغيل الدوائر الطويلة والمعقدة المطلوبة من الطرق القديمة.

اختبر الباحثون طريقتهم على نموذج محدد يُعرف باسم "نموذج أندرسون أحادي الشائبة" (single-impurity Anderson model)، والذي يصف ذرة واحدة تتفاعل مع بحر من الإلكترونات. حيث أعدوا محاكاة تكون فيها هذه الذرة المركزية متصلة بسلسلة من اثني عشر موقعًا، تمثل البيئة المحيطة. في هذا الإعداد، تحمل الذرة المركزية قوة تنافر قوية تمنع إلكترونين من احتلال نفس الموقع، وهي سمة رئيسية تقود التحول بين المعدن الذي يوصل الكهرباء والعازل الذي لا يوصلها. استخدم الفريق حاسوبًا كلاسيكيًا لمحاكاة الجهاز الكمومي، مما سمح لهم بمقارنة طريقتهم الجديدة مقابل الحل الدقيق والمثالي لهذا النظام. ووجدوا أنه من خلال أخذ عينات من التطور الكمومي، يمكنهم إعادة بناء خريطة طاقة الإلكترون بدقة عالية باستخدام جزء ضئيل فقط من إجمالي الحالات الممكنة. بالنسبة لنظام مكون من اثني عشر موقعًا، قللت الطريقة من المساحة الحسابية المطلوبة بمقدار ثلاثة إلى أربعة مراتب عشرية، مما يعني أنها احتاجت للتعامل مع بضعة آلاف من الحالات فقط بدلًا من مئات الملايين التي كانت ستتطلبها عملية حسابية كاملة.

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

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

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

يقدم هذا النهج مسارًا واعدًا لنظرية المجال المتوسط الديناميكي (dynamical mean-field theory)، وهو إطار قوي يُستخدم لدراسة المواد ذات الارتباط القوي. في هذا الإطار، يتم رسم شبكة معقدة من الذرات وتحويلها إلى مشكلة شائبة واحدة، والتي يجب حلها بشكل متكرر لإيجاء خصائص المادة بأكملها. لطالما كان "حل الشائبة" (impurity solver) هو عنق الزجاجة، وهي الخطوة التي تحسب سلوك الإلكترون لذلك الموقع الواحد. وإذا أمكن تشغيل هذا الحل على حاسوب كمومي مع عدد أكبر من مواقع الحمام (bath sites) مما هو ممكن حاليًا، فسيسمح ذلك للعلماء بنمذجة المواد بدقة أكبر بكثير. ويقترح الباحثون أن طريقتهم يمكن أن تعمل كـ "حلل للشائبة"، مما يتيح دراسة أنظمة أكبر وأكثر واقعية من أي وقت مضى. ورغم أن الدراسة كانت محدودة بنموذج معين ولم تتضمن بعد الحلقة الكاملة ذاتية الاتساق للنظرية، إلا أن النتائج توفر إثباتًا قويًا لمفهوم إمكانية إجراء مثل هذه الحسابات.

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

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