Green's function theory of magnetism in Bi$_2$CuO$_4$: anisotropic Heisenberg XYZ model

رقصة المغناطيس: مقدمة تمهيدية

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

لفهم كيفية حركة هذه الدورات، يستخدم العلماء أداة رياضية تُسمى نظرية دالة غرين (Green's function theory). فكر في هذا كطريقة متطورة للغاية للتنبؤ بـ "الموسيقى" التي تعزفها المادة عندما تنقر عليها. تماماً كما يهتز وتر الجيتار بنغمات محددة، تهتز البلورة المغناطيسية بترددات محددة عند تعرضها للاضطراب. تُسمى هذه الاهتزازات موجات الدوران (spin waves) أو الماجنونات (magnons). وأحياناً، تحتوي المادة على "فجوة" في موسيقاها—أي صمت حيث لا يمكن عزف النغمات المنخفضة. هذه الفجوة حاسمة لأنها تخبرنا بمدى صلابة النظام المغناطيسي وكيف ستتفاعل المادة مع مجال مغناطيسي خارجي. يهتم العلماء بهذا الأمر لأن هذه المواد المغناطيسية هي اللبنات الأساسية للتقنيات المستقبلية، مثل ذاكرات الكمبيوتر فائقة السرعة التي يمكن التحكم فيها بالكهرباء.

قصة بلورة النحاس الملتوية

الآن، دعونا ننتقل إلى مختبر ر. أو. كوزيان وإي. إي. كراسوفسكي، اللذين قررا استقصاء بلورة غريبة تُسمى Bi₂CuO₄ (أكسيد البزموت والنحاس). هذه المادة تشبه صندوق الألغاز؛ فهي تشبه إلى حد ما الموصلات الفائقة ذات الحرارة العالية الشهيرة، لكن سلوكها مختلف تماماً. في الداخل، تشكل ذرات النحاس منصات مسطحة مربعة الشكل (plaquettes) تتراكم فوق بعضها مثل برج ملتوي. من المفترض أن تكون الدورات في ذرات النحاس هذه غير متوافقة في الاتجاه، ولكن لسنوات، لم يتفق العلماء على الاتجاه الذي تشير إليه بالضبط أو كيفية تفاعلها.

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

بنى الباحثون نموذجاً يُسمى نموذج XYZ، والذي يأخذ في الاعتبار هذه القوى الاتجاهية المختلفة. ثم طبقوا هذا النموذج على Bi₂CuO₄ لمعرفة ما إذا كان بإمكانه تفسير "الموسيقى" التي تصدرها البلورة. لقد بحثوا في تشتت استثارة الدوران (spin excitation dispersion)، وهو مجرد طريقة منمقة لوصف كيفية تغير طاقة موجات الدوران أثناء تحركها عبر البلورة.

إليكم ما وجدوه. أولاً، أكدوا أن الدورات في Bi₂CuO₄ مرتبة بالفعل في نمط ملتوي محدد. أظهر نموذجهم أن "موسيقى" البلورة تتكون من فرعين متميزين من الموجات. والأهم من ذلك، اكتشفوا أن انحيازاً مغناطيسياً ضئيلاً وضعيفاً في مستوى البلورة (التباين داخل المستوى - in-plane anisotropy) يخلق "فجوة" صغيرة في الفرع الأدنى من هذه الموجات. وهذا يعني أنه عند مستويات الطاقة المنخفضة جداً، ترفض البلورة الاهتزاز؛ فهي تظل صامتة حتى تدفعها بقوة كافية.

حسب المؤلفان أن هذه الفجوة تعتمد على قوة المجال المغناطيسي الخارجي. ومع زيادة المجال، صغرت الفجوة أكثر فأكثر حتى اختفت تماماً عند مجال حرج ($B_c$) يبلغ حوالي 0.4 تسلا. يمثل هذا الاختفاء حدثاً دراماتيكياً يُسمى انتقال قلب الدوران (spin-flop transition)، حيث يغير النظام المغناطيسي للبلورة بأكمله اتجاهه فجأة ليتماشى مع المجال الجديد.

من خلال مطابقة حساباتهم مع البيانات التجريبية من تشتت النيوترونات وقياسات الرنين، تمكن المؤلفون من تحديد القيم الدقيقة للتفاعلات المغناطيسية. وجدوا أن معامل التباين خارج المستوى يبلغ حوالي 0.013، والمعامل الضئيل داخل المستوى يبلغ 4.7 ⋅ 10⁻⁶. سمحت هذه الأرقام لهم بالتنبؤ بـ درجة حرارة نيل (Néel temperature - $T_N$)—وهي درجة الحرارة التي تفقد عندها البلورة نظامها المغناطيسي وتتحول إلى حساء فوضوي. أعطت حساباتهم قيمة قدرها 52 كلفن، وهي قريبة جداً من النطاق التجريبي الذي يتراوح بين 42-47 كلفن.

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

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