Nonlinear phonon dispersion in disordered solids and non-Debye vibrational spectra

ملخص تقني: تشتت الفونونات غير الخطي في المواد الصلبة المضطربة والأطياف الاهتزازية غير الديبية

بيان المشكلة
بينما يمكن وصف المواد الصلبة البلورية بشكل جيد باستخدام نموذج ديباي الكلاسيكي — حيث تتبع الفونونات الصوتية منخفضة التردد علاقة تشتت خطية $\omega(k) \propto k$ وتتبع كثافة الحالات الاهتزازية (VDoS) قانون $\omega^2 \propto D(\omega)$ — تُظهر المواد الصلبة المضطربة (الزجاج) انحرافات كبيرة تُعرف بالأنومالية غير الديبية (non-Debye anomalies). وأبرز هذه الأنوماليات هي "ذروة بوسون" (Boson Peak - BP)، وهي عبارة عن فائض في الأنماط الاهتزازية فوق توقعات ديباي. لا يزال أصل هذه الأنوماليات موضوعاً للنقاش؛ حيث يُقترح آليتان رئيستان: (1) الليونة غير الخطية في علاقة تشتت الفونونات $\omega(k)$ الناتجة عن الاضطراب الهيكلي، و(2) وجود أنماط اهتزازية شبه موضعية (quasi-localized) غير فونونية متأصلة في البنى المضطربة. وتتمثل الفجوة الحرجة في الفهم الحالي في عدم القدرة على الفصل الكمي وتحديد المساهمات النسبية لهاتين الآليتين في ظهور السلوك غير الديبي وذروة بوسون. علاوة على ذلك، فإن طبيعة طول المقياس الناتج عن الاضطراب والذي يحكم تشتت الفونونات وتوهين الموجات في المواد غير البلورية ليست مفهومة تماماً.

المنهجية
لمعالجة هذه التحديات، استخدم المؤلفون محاكاة حاسوبية واسعة النطاق لثلاث فئات متميزة من المواد الصلبة المضطربة ثلاثية الأبعاد:

1. الشبكات المرنة المضطربة: تتكون من نوابض هوكية مسترخية مع تباين في متوسط أعداد التنسيق ($z$)، تتراوح من عتبة ماكسويل الحرجة ($z_c=6$) إلى ما فوقها بكثير.
2. زجاجات القوة العكسية متعددة التشتت (PIPL): أنظمة تنافرية تخضع لتاريخ حراري متغير، يتم قياسه من خلال درجة الحرارة الأم $T_p$ التي يخرج عندها النظام من حالة التوازن.
3. الزجاجات الثنائية: تشمل مخاليط ثنائية تنافرية وزجاجات "الكرات اللزجة" (sticky-sphere) ذات قوى جاذبة قابلة للضبط عبر معامل $r_c$.

استخدمت الدراسة تقنية حاسوبية مبتكرة تسمى "طريقة الموجة المفروضة" (imposed-wave method) لحساب علاقة تشتت الفونونات $\omega(k)$ في الأنظمة المضطربة. وخلافاً للبلورات، حيث تكون الفونونات دوال ذاتية دقيقة لمصفوفة هسيان (Hessian)، فإن الاضطراب يمنع تحديد تردد $\omega$ فريد لكل متجه موجة $k$. تقوم طريقة الموجة المفروضة بتعيين تردد فريد لكل متجه موجة مقبول من خلال فرض مجال قوة شبيه بالموجة على الجسيمات وحساب إزاحة الاستجابة الخطية. وهذا يسمح باستخراج $\omega(k)$ والسرعة الجماعية المرتبطة بها.

كما استخدم المؤلفون طريقة كثير الحدود النواة (Kernel Polynomial Method - KPM) لحساب كثافة الحالات الاهتزازية $D(\omega)$ لأحجام أنظمة كبيرة جداً ($N \sim 10^6$)، مما مكن من تحليل الذيول منخفضة التردد وفصل المساهمات الفونونية وغير الفونونية. كما تم حساب مقياس اضطراب لا بُعدي، $\chi$، وهو نسبة تقلبات معامل القص إلى متوسطه عبر المقاييس الميزوسكوبية، لتوحيد وصف الاضطراب عبر النماذج المختلفة.

المساهمات والنتائج الرئيسية

1. اكتشاف طول مقياس ميزوسكوبي ناتج عن الاضطراب ($\xi$):
   أثبت المؤلفون أن الليونة غير الخطية في تشتت الفونونات في المواد الصلبة المضطربة لا تحكمها المسافة البينية الذرية الدقيقة $a_0$ (كما في البلورات)، بل يحكمها طول مقياس ميزوسكوبي ناشئ أكبر بكثير $\xi \gg a_0$. وتتميز علاقة التشتت بالشكل $\omega(k) = f(k\xi)ck$، حيث تتبع أول رتبة من اللاخطية العلاقة $\omega(k) \approx ck - c\Upsilon \xi^2 k^3$ (مع $\Upsilon = 1/96$).

• علاقات القياس: وُجد أن $\xi$ يتدرج مع درجة الاضطراب: $\xi \sim 1/\sqrt{z-z_c}$ للشبكات المرنة بالقرب من انتقال عدم الانحباس (unjamming transition)، وتزديراد $\xi$ مع زيادة $T_p$ و $r_c$ في الزجاجات.
• مقياس اضطراب موحد: من الأهمية بمكان أن المؤلفين أظهروا أن النسبة $\xi/a_0$ تتناسب خطياً مع مقياس الاضطراب اللا بُعدي $\chi$ ($\xi/a_0 \sim \chi$)، مما يوحد تأثيرات بارامترات التحكم في الاضطراب المتنوعة.

2. الارتباط بتوهين الموجات:
   أظهر أن نفس الطول الميزوسكوبي $\xi$ يتحكم في معدلات توهين تشتت رايلي للموجات $\Gamma(k)$. وقد اشتق المؤلفون علاقة $\Gamma(k) \sim \xi^2 k^4$ لقيم $k$ الصغيرة، وهي متسقة مع نظرية المرونة غير المتجانسة (HET) وعلاقة رايلي-كليمينز. وهذا يثبت وجود رابط مباشر بين التشتت غير الخطي وتوهين الموجات، وكلاهما محكوم بنفس الطول الميزوسكوبي الناتج عن الاضطراب.

3. التفكيك الكمي للأطياف غير الديبية:
   قام البحث بفصل مساهمات تشتت الفونونات غير الخطي (الفونونية) والاهتزازات غير الفونونية في كثافة الحالات الاهتزازية $D(\omega)$ كمياً.

• المساهمة الفونونية: باستخدام $\omega(k)$ المقاسة، حسب المؤلفون كثافة الحالات الفونونية $D_{ph}(\omega)$، والتي تنحرف عن قانون ديباي كالتالي: $D_{ph}(\omega) \approx A_D \omega^2 + A_w \omega^4 + \dots$. ويعتمد المعامل $A_w$ على $\xi$.
• المساهمة غير الفونونية: يوجد توزيع بلا فجوة للأنماط غير الفونونية $D_G(\omega) \approx A_g \omega^4$.
• الأهمية النسبية: تحدد النسبة $A_g/A_w$ أي آلية تهيمن على بدء السلوك غير الديبي. وتظهر النتائج أن هذه النسبة حساسة للغاية لدرجة الاضطراب. في الزجاجات شديدة الاضطراب (غير المستقرة)، تهيمن الأنماط غير الفونونية ($A_g/A_w > 1$). ومع ذلك، في الزجاجات المختبرية المستقرة والواقعية (منخفضة الاضطراب)، غالباً ما تهيمن المساهمة الفونونية الناتجة عن التشتت غير الخطي ($A_g/A_w < 1$). وهذا يتحدى التفسيرات السابقة التي تعزو ذيل $\omega^4$ حصرياً للأنماط غير الفونونية.

4. تكوين ذروة بوسون:
   من خلال تحليل الكثافة التراكمية للحالات عند تردد ذروة بوسون $\omega_{BP}$، عرّف المؤلفون نسبة $R(\omega_{BP})$ التي تمثل كسر الاهتزازات غير الفونونية من إجمالي الأنماط الزائدة.

• النتيجة: وُجد أن $R(\omega_{BP})$ يمثل كسرًا معتبرًا من الواحد الصحيح (يتراوح بين $\sim 0.2$ إلى $\sim 0.5$) عبر نطاق واسع من الاضطراب الزجاجي.
• الاستنتاج: ذروة بوسون ليست مشكلة حصرياً من الفونونات الملطّنة (كما تقترح بعض النظريات الحديثة) ولا من الأنماط غير الفونونية وحدها. بدلاً من ذلك، تنشأ من تهجين كبير بين تلطن الفونونات غير الخطي والتهجين مع الإثارات شبه الموضعية غير الفونونية.

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