Lattice dynamics of the charge density wave compounds TaTe$_4$ and NbTe$_4$ and their evolution across solid solutions

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

تُظهر مركبات التتراكالكوچينيد المعدنية الانتقالية شبه أحادية البعد (TMTs)، وتحديداً TaTe$_4$ و NbTe$_4$، أطواراً إلكترونية معقدة، بما في ذلك موجات كثافة الشحنة (CDWs)، والحالات الطوبولوجية، والموصلية الفائقة تحت الضغط. وعلى الرغم من تشابههما البنيوي، إلا أنهما يُظهران اختلافات كبيرة في خصائص النقل وخصائص موجة كثافة الشحنة (CDW):

• TaTe$_4$ يُظهر موجة كثافة شحنة متوافقة (commensurate) عند درجة حرارة الغرفة ومقاومة مغناطيسية كبيرة.
• NbTe$_4$ يُظهر طبيعة موجة كثافة شحنة مثيرة للجدل (غير متوافقة مقابل غير متوافقة التجزئة/discommensurate) ومقاومة مغناطيسية أضعف.

بينما تشير دراسات البنية الإلكترونية إلى وجود عدم استقرار مدفوع بالفونونات (أنماط فونون رخوة) بدلاً من تعشيش سطح فيرمي (Fermi surface nesting) كمحرك أساسي لموجة كثافة الشحنة، إلا أن التوصيف التجريبي المباشر لديناميكا الشبكة عبر هذه المواد ومحاليلها الصلبة لا يزال محدوداً. إن فهم الأصل المجهري لموجة كثافة الشحنة والاختلافات بين TaTe$_4$ و NbTe$_4$ يتطلب تحقيقاً منهجياً لكيفية تطور الأنماط الاهتزازية عند ضبط الشبكة الفرعية للمعدن كيميائياً.

2. المنهجية

تستخدم الدراسة نهجاً نظرياً وتجريبياً مدمجاً:

• تحضير العينات: تم نمو بلورات أحادية من المحلول الصلب Ta$_{1-x}$Nb$_x$Te$_4$ ($x = 0.0, 0.2, 0.4, 0.6, 0.8, 1.0$) باستخدام طريقة التدفق الذاتي (self-flux method).
• التوصيف التجريبي:
  • حيود الأشعة السينية (XRD): استُخدم لتأكيد نقاء الطور وتتبع تطور معاملات الشبكة عبر سلسلة التطعيم.
  • مطيافية رامان (Raman Spectroscopy): أُجريت عند درجة حرارة الغرفة باستخدام إثارات ليزر بطول موجي 532 نانومتر، 638 نانومتر، و738 نانومتر لقياس الأنماط الاهتزازية.
• الحسابات النظرية:
  • نظرية الدالة الكثافية (DFT): أُجريت باستخدام برنامج VASP مع دالات GGA-PBEsol وإمكانات Projector Augmented-Wave (PAW).
  • حسابات الفونون: أُجريت باستخدام طريقة الفروق المحدودة (Phonopy) في خلية فائقة من نوع $(2\times2\times2)$ لرسم ديناميكا حدود المنطقة.
  • اللاخطية (Anharmonicity): أُجريت محاكاة الديناميكا الجزيئية (MD) عند 300 كلفن لأخذ التفاعلات اللاخطية في الاعتبار وإعادة تطبيع طيف الفونون، ومقارنة النتائج عند 0 كلفن و300 كلفن.

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

• تعيين الأنماط المنهجي: قدمت تعييناً متسقاً للأنماط النشطة في رامان لمركبي TaTe$_4$ و NbTe$_4$ من خلال ربط أطياف رامان بتوقعات نظرية DFT، مما حل التباسات تسميات التماثل ($A_{1g}$, $B_{1g}$, $B_{2g}$, $E_g$).
• تحديد عدم استقرار الفونون: أكدت الحسابات النظرية وجود عدم استقرار فونوني عند متجه موجة كثافة الشحنة (بالقرب من نقطة A) في كلا المركبين الأصليين، مما يدعم آلية مدفوعة بالشبكة لتكوين موجة كثافة الشحنة.
• اكتشاف سلوك التطور المزدوج: كشفت عن تباين متميز في كيفية تطور الأنماط الاهتزازية عبر المحلول الصلب:
  • الأنماط التي يهيمن عليها التيلوريوم (Te): تتطور بسلاسة في التردد.
  • الأنماط التي يهيمن عليها المعدن: تُظهر ترددات "مثبتة" (pinned) مع إعادة توزيع الكثافة.

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

أ. عدم استقرار الفونون وآلية موجة كثافة الشحنة (CDW)

• كشفت حسابات DFT عند 0 كلفن عن أنماط فونون تخيلية (غير مستقرة) بالقرب من نقطة A لكل من TaTe$_4$ ($\omega \approx -61.2$ سم$^{-1}$) و NbTe$_4$ ($\omega \approx -71.3$ سم$^{-1}$).
• تتوافق حالات عدم الاستقرار هذه مع التريميرية (trimerization) لذرات المعادن الانتقالية على طول السلاسل، وهو ما يتوافق مع البنية الفائقة لموجة كثافة الشحنة المرصودة ($P4/ncc$).
• أظهرت محاكاة الديناميكا الجزيئية (MD) عند 300 كلفن أن التفاعلات اللاخطية تعمل على استقرار الشبكة، مما يزيل الترددات التخيلية، وهذا يتماشى مع الملاحظة التجريبية لأنماط رامان المستقرة عند درجة حرارة الغرفة.

ب. تطور أنماط رامان في المحاليل الصلبة (Ta$_{1-x}$Nb$_x$Te$_4$)

حددت الدراسة سلوكين متميزين لأنماط التماثل $E_g$ مع زيادة تركيز النيوبيوم ($x$):

1. أنماط $E_g$ منخفضة التردد (يهيمن عليها التيلوريوم - Te):

   • السلوك: تتحول هذه الأنماط بسلاسة ورتابة (monotonically) في التردد من حد TaTe$_4$ إلى حد NbTe$_4$.
   • الكثافة: تظل ثابتة نسبياً عبر السلسلة.
   • الأصل الفيزيائي: هذه الاهتزازات مدفوعة بشكل أساسي بذرات التيلوريوم (Te) وتتضمن حركات بين السلاسل أو حركات السلسلة الصلبة. وهي تستكشف المتوسط الماكروسكوبي لمحتوى المعدن.

2. أنماط $E_g$ عالية التردد (يهيمن عليها المعدن):

   • السلوك: تُظهر هذه الأنماط سلوكاً شاذّاً. فبدلاً من التحول بسلاسة، تظهر قمتان متميزتان عند التركيزات المتوسطة:
     • إحدى القمم تظل "مثبتة" بالقرب من تردد TaTe$_4$ النقي ($\sim 170$ سم$^{-1}$).
     • والأخرى تظل "مثبتة" بالقرب من تردد NbTe$_4$ النقي ($\sim 213$ سم$^{-1}$).
   • الكثافة: تتناسب كثافة كل قمة مثبتة مباشرة مع تركيز المعدن المقابل (Ta أو Nb). مع زيادة $x$، تتلاشى قمة الـ Ta المثبتة وتختفي، بينما تنمو قمة الـ Nb المثبتة.
   • الأصل الفيزيائي: تهيمن هذه الأنماط على حركة ذرات المعادن الانتقالية (Ta/Nb) داخل السلسلة. يشير طبيعتها "المثبتة" إلى طابع قصير المدى، مما يعني أن تردد الاهتزاز يتحدد من خلال البيئة الكيميائية المحلية المباشرة (المعدن المحدد) بدلاً من التركيب المتوسط العالمي.

ج. الاتساق البنيوي

• أكدت بيانات XRD حدوث انكماش مستمر في معاملات وحدة الخلية ($a$ و $c$) مع استبدال التا (Ta) بالنيوبيوم (Nb)، مما ينفي حدوث انفصال طوري ماكروسكوبي. وبالتالي، فإن سلوك رامان الشاذ هو أمر جوهري للديناميكا الشبكية المحلية، وليس نتيجة لانفصال الأطوار.

5. الأهمية

• آلية تكوين موجة كثافة الشحونة (CDW): تدعم النتائج بقوة الآلية المدفوعة بالشبكة (الفونون) لانتقال موجة كثافة الشحونة في هذه المواد، حيث ترتبط أنماط الفونون غير المستقرة مباشرة بالتريميرية الهيكلية.
• الديناميكا المحلية مقابل العالمية: تسلط الدراسة الضوء على أنه بينما تتطور البنية البلورية عالمياً (تغير مستمر في معاملات الشبكة)، فإن أنماط معينة عالية التردد تحتفظ بـ "ذاكرة" للهوية الذرية المحلية. هذا يشير إلى أن تشوه الشبكة المرتبط بموجة كثافة الشحونة حساس للغاية للبيئة المحلية للمعدن الانتقالي.
• مسبار لتباينات موجة كثافة الشحنة: يشير السلوك المتميز لنمط $E_g$ عالي التردد الذي يهيمن عليه المعدن إلى أنه يعمل كمسبار حساس لخصائص موجة كثافة الشحونة المختلفة (المتوافقة مقابل غير المتوافقة التجزئة) الملحوظة في TaTe$_4$ و NbTe$_4$. إن التعديل المباشر لهذا النمط للشبكة المعدنية وتكوين التريمير يوحي بوجود صلة محتملة بين هذه الاهتزازات المحددة واستقرار أو بنية نطاقات موجة كثافة الشحنة.
• الأثر المنهجي: توضح أن التوصيف المنهجي للفونون عبر المحاليل الصلبة يعد أداة قوية لفك الاشتباك بين الآليات الإلكترونية والشبكية المتنافسة في المواد الكمومية منخفضة الأبعاد.