🔬 materials science

DFT calculations of magnetocrystalline anisotropy energy with fixed spin moment

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

المؤلفون الأصليون: Justyn Snarski-Adamski (Institute of Molecular Physics, Polish Academy of Sciences, Poznan, Poland), Joanna Marciniak (Institute of Molecular Physics, Polish Academy of Sciences, Poznan, Poland, Uppsa

نُشر 2026-03-11

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: Justyn Snarski-Adamski (Institute of Molecular Physics, Polish Academy of Sciences, Poznan, Poland), Joanna Marciniak (Institute of Molecular Physics, Polish Academy of Sciences, Poznan, Poland, Uppsala University, Uppsala, Sweden), Wojciech Marciniak (Institute of Molecular Physics, Polish Academy of Sciences, Poznan, Poland, Poznan University of Technology, Poznan, Poland), Justyna Rychły-Gruszecka (Institute of Molecular Physics, Polish Academy of Sciences, Poznan, Poland), Mirosław Werwiński (Institute of Molecular Physics, Polish Academy of Sciences, Poznan, Poland)

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

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

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

يقدم هذا البحث طريقة جديدة ذكية لحل هذا الارتباك، باستخدام طريقة تسمى عزم المغزل الثابت (FSM).

التشبيه: "قرص" الراديو

تخيل أن القوة المغناطيسية للمادة (عزم المغزل) تشبه مقبض التحكم في مستوى الصوت في راديو.

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

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

الاكتشاف الكبير: "الخريطة العالمية"

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

عندما ترسم القبضة (MAE) مقابل مستوى الصوت (عزم المغزل)، ترسم جميع النماذج المختلفة نفس المنحنى. هم فقط يهبطون في نقاط مختلفة على طول هذا المنحنى.

النموذج (أ) قد يقول إن النقطة الطبيعية هي عند مستوى صوت 30% بقبضة قدرها 5.
النموذج (ب) قد يقول إن النقطة الطبيعية هي عند مستوى صوت 35% بقبضة قدرها 4.

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

لماذا يهم هذا في الحياة الواقعية؟

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

العقبة

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

باختصار

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

ملخص تقني لورقة بحثية بعنوان: "حسابات نظرية الكثافة الوظيفية لطاقة تباين الخواص المغناطيسية البلورية باستخدام عزم المغزل الثابت"

1. بيان المشكلة

يعتمد تطوير الجيل الجديد من المغناطيسات الدائمة، وخاصة تلك ذات المحتوى المحدود من العناصر الأرضية النادرة، بشكل كبير على نظرية الكثافة الوظيفية (DFT) للتنبؤ بالخصائص الجوهرية مثل طاقة تباين الخواص المغناطيسية البلورية (MAE). ومع ذلك، يوجد تحدٍ كبير: فإن قيم الـ MAE التي يتم حسابها باستخدام جهود تبادل-ارتباط (XC) مختلفة (مثل LDA مقابل GGA) غالبًا ما تعطي نتائج متباينة ومتناقضة، وأحيانًا مختلفة نوعيًا.

على سبيل المثال، قد يتنبأ أحد الدوال بوجود تباين خواص أحادي المحور قوي (MAE موجب) مع دالة معينة، وبتباين خواص في المستوي (MAE سالب) مع دالة أخرى. هذا التضارب يخلق غموضًا لدى التجريبيين ويعيق التصميم الموثوق للمواد المغناطيسية الجديدة. علاوة على ذلك، تُجرى حسابات DFT القياسية عادةً عند حالة الاستقرار (0 كلفن)، بينما تتطلب التطبيقات العملية فهم الـ MAE عند درجة حرارة الغرفة وما فوقها، حيث يمكن أن تكون التبعيات معقدة (مثل السلوك الشبيه بالموجة الجيبية).

2. المنهجية

يقترح المؤلفون ويستخدمون طريقة عزم المغزل الثابت النسبية الكاملة (FR-FSM) لحل هذه التناقضات.

المفهوم الجوهري: بدلاً من حساب الـ MAE فقط عند عزم المغزل المغناطيسي المتزن ذاتيًا، تقوم طريقة FR-FSM بتقييد إجمالي عزم المغزل المغناطيسي ( $m_s$ ) بقيم ثابتة عبر سلسلة من القيم أثناء الحسابات المتسقة ذاتيًا.
التنفيذ:
- كود FPLO: الأداة الرئيسية المستخدمة هي كود FPLO (الإصدارات 14 و18)، والتي تطبق FR-FSM عن طريق إضافة مجال مغناطيسي مساعد ( $B_{FSM}$ ) إلى معادلة ديراك النسبية الكاملة. يعمل هذا المجال كمعامل لاجرانج لتقارب عزم المغزل نحو القيمة المستهدفة بشكل تكراري.
- المقارنة مع VASP: يشير البحث إلى أنه بينما يستخدم VASP نهج "عزوم المغزل المحلية المقيدة" (دالة الجزاء)، فإنه لا يثبت إجمالي عزم المغزل في حسابات SOC بنفس الفعالية التي يحققها FR-FSM في FPLO.
بروتوكول الحساب:
- يتم حساب الـ MAE كفرق طاقة بين محوري المغناطيسية السهل والصعب.
- أُجريت الحسابات باستخدام دوال XC متنوعة: LDA (von Barth-Hedin, Perdew-Zunger, Perdew-Wang)، وGGA (PBE)، ودوال التبادل فقط (exchange-only).
- استُخدمت طريقة التتراهيدرون (tetrahedron method) لتكامل منطقة بريلوين لضمان دقة عددية عالية.
- بالنسبة للسبائك، تم استخدام تقريب البلورة الافتراضية (VCA) لنمذجة تغيرات التركيب.

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

التوفيق بين جهود الـ XC: يوضح البحث أنه بينما تختلف قيم الـ MAE عند التوازن عبر الدوال المختلفة، فإن منحنيات MAE( $m_s$ ) التي تولدها FR-FSM تتداخل بشكل كبير بغض النظر عن جهد الـ XC المستخدم. النقاط المتزنة التي يتم الحصول عليها من الدوال المختلفة تقع ببساطة في مواضع مختلفة على طول هذا المنحنى العالمي.
تحديد القيم القصوى الافتراضية: تسمح طريقة FR-FSM للباحثين برسم خريطة كاملة لاعتماد MAE( $m_s$ )، مما يتيح تقدير الحد الأقصى الافتراضي للـ MAE الذي يمكن أن تحققه المادة إذا تم ضبط عزمها المغناطيسي للوصول إلى القيمة المثلى.
إطار تصميم السبائك: توفر الطريقة أداة لتحديد الإضافات المثالية للسبائك. من خلال تحليل خرائط MAE كدالة لكل من التركيب وعزم المغزل، يمكن للباحثين تحديد كيفية نقل النظام نحو العزم الأمثل ( $m_s$ ) لتحقيق أقصى قدر من التباين.
الارتباط الحراري: يقترح المؤلفون أنه يمكن ربط علاقة MAE( $m_s$ ) بـ MAE( $T$ ) (الاعتماد على درجة الحرارة)، مما يوفر وسيلة للتنبؤ بالسلوكيات الحرارية المعقدة (مثل تلك الملاحظة في Fe $_2$ B و Co $_2$ B) التي تغفل عنها حسابات الحالة الأرضية القياسية.

4. النتائج الرئيسية

التداخل العالمي (الشكل 1): بالنسبة لمواد مثل FePt، CeFe $_{12}$ $_{12}$ ، FeNi، FeB، Fe $_2$ $_{2}$ P، وMnBi، فإن قيم الـ MAE المحسوبة باستخدام LDA وGGA ودوال التبادل فقط تقع جميعها على منحنى MAE( $m_s$ $m_{s}$ ) واحد مستمر.
- ملاحظة: تتنبأ الدوال المختلفة بعزوم مغزل متزنة مختلفة (مرتبة كالتالي: BH < PW92 < PBE < Exchange-only)، مما يؤدي إلى قيم MAE متزنة مختلفة. ومع ذلك، فإنها جميعًا تأخذ عينة من نفس العلاقة الفيزيائية الأساسية.
- مثال: في CeFe $_{12}$ ، ينتقل الـ MAE المتزن من +2 MJ/m $^3$ (أحادي المحور) إلى -1 MJ/m $^3$ (في المستوي) اعتمادًا على الدالة المستخدمة. وتكشف FR-FSM أن هذا مجرد انتقال على طول المنحنى، وليس تناقضًا جوهريًا.
تحسين السبائك (الشكل 2 و3):
- في (Fe $_{1-x}$ Co $_x$ ) $_5$ SiB $_2$ و (Fe $_{1-x}$ Co $_x$ ) $_5$ PB $_2$ ، يؤدي استبدال Si بـ P إلى إزاحة خريطة MAE، مما يسمح باستراتيجيات استهداف للسبائك.
- تظهر الدراسة أن التغيرات في MAE الناتجة عن تمدد الحجم أو سبائك الكوبالت تتبع نفس اتجاه منحنى FSM، مما يؤكد أن عزم المغزل هو المحرك الأساسي لتغيرات الـ MAE.
ما وراء LDA/GGA (MnBi): بالنسبة للأنظمة ذات الارتباط القوي مثل MnBi، كانت منحنيات LDA/GGA القياسية غير كافية. أدت إضافة حد Hubbard $U$ (أي PBE+U) إلى تغيير منحنى MAE( $m_s$ ) بشكل كبير، مما أثبت أن FR-FSM هي أيضًا أداة ضرورية لاختبار تصحيحات ما وراء LDA.

5. الأهمية والقيود

الأهمية: يحول إطار FR-FSM حساب الـ MAE من تنبؤ بنقطة واحدة (غالبًا ما يكون غير موثوق بسبب اختيار الدالة) إلى علاقة مستقلة عن الدالة. فهو يوفر أداة قوية لـ:
1. التحقق من البيانات التجريبية مقابل الحدود النظرية.
2. توجيه التركيب الكيميائي للمغناطيسات الدائمة الجديدة للوصول إلى "قمة" منحنى الـ MAE.
3. فهم تبعات درجة الحرارة المعقدة.
القيود:
- التكلفة الحسابية: تتطلب FR-FSM حسابات ذات جهد كامل (full-potential) مع شبكات k-point كثيفة، مما يجعلها مكلفة حاسوبيًا.
- توافر البرمجيات: حاليًا، يدعم كودان فقط (FPLO و RSPt) طريقة عزم المغزل الثابت النسبية الكاملة. إن عدم تنفيذها في الأكواد واسعة الانتشار مثل VASP (الذي يستخدم طريقة تقييد مختلفة) يحد من اعتمادها على نطاق أوسع.
- التقريب: تفترض الطة أن إجمالي عزم المغزل متطابق لكل من المحور السهل والمحور الصعب، وهو تقريب صالح للأنظمة المدروسة ولكنه قد يؤدي إلى نتائج غير فيزيائية في الأنظمة ذات الاختلافات الكبيرة في العزم الناتجة عن التباين.

الخلاصة: تثبت هذه الورقة أن FR-FSM هي منهجية حاسمة للتصميم العقلاني للمغناطيسات الدائمة، حيث توفر وسيلة للتوفيق بين النتائج المتضاربة لـ DFT وتحديد الحدود العليا النظرية للأداء المغناطيسي.