GEM: An implementation of the ghost-Gutzwiller approximation for simulating interacting quantum systems

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

وقد قدمت طريقة أخف وأسرع تُعرف باسم "تقريب غوتزويلر" (Gutzwiller approximation) مسارًا مختصرًا منذ فترة طويلة من خلال التركيز فقط على السلوك المتوسط للإلكترونات، وتجاهل التقلبات عالية الطاقة المعقدة. ومع ذلك، فإن هذا التبسيط يأتي على حساب شيء ما: فهو يفتقد التفاصيل الحاسمة حول طيف الطاقة الذي يحدد كيفية سلوك المواد فعليًا. ومن أجل سد هذه الفجوة، تم تطوير إطار نظري جديد يسمى "تقريب غوتزويلر الشبح" (ghost-Gutzwiller approximation). فهو يقدم "المدارات الشبحية" (ghost orbitals)—وهي حالات مساعدة وهمية تعمل كحاجز لاستعادة تفاصيل الطاقة العالية المفقودة دون التكلفة الكاملة للطريقة الثقيلة. وحتى الآن، كان هذا الوسط الواعد موجودًا فقط في النظرية وفي قصاصات برمجية متفرقة، مفتقرًا إلى حزمة برمجية موحدة ومفتوحة المصدر يمكن للباحثين استخدامها بسهما لاختبارها على نماذج واقعية.

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

تمت كتابة البرنامج بلغة بايثون (Python)، وهو مصمم للعمل بسلاسة مع الأدوات العلمية الحالية، وخاصة منظومة TRIQS، وهي إطار عمل قياسي للمحاكاة الكمومية. يتولى GEM المهام الشاقة المتمثلة في إعداد المعادلات وإدارة حلقة "الاستمرارية الذاتية" (self-consistency loop)، حيث تقوم المحاكاة بتحديث افتراضاتها بشكل متكرر حتى تستقر النتائج. يتضمن البرنامج "حلالًا" (solver) مدمجًا يستخدم "التشخيص الدقيق" (exact diagonalization)، وهي تقنية تجد حالات الطاقة الدقيقة لنظام كمومي صغير، ولكنه مبني أيضًا لاستقبال حلالات أخرى في المستقبل. تعني هذه النمطية أن الباحثين يمكنهم دمج طرق عددية مختلفة اعتمادًا على المشكلة المحددة التي يدرسونها، سواء كانت تتعلق بالحالات الأرضية عند درجة حرارة الصفر أو الديناميكا الحرارية عند درجات حرارة محددة.

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

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

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

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