Negative Thermal Expansion in Cubic Ice: A Collective Quantum Effect of the hydrogen-bond network

ملخص تقني: التمدد الحراري السالب في الجليد المكعب: تأثير كمي جماعي لشبكة الروابط الهيدروجينية

بيان المشكلة
يُلاحظ التمدد الحراري السالب (NTE)، وهو ظاهرة انكماش المواد الصلبة عند تسخينها، في العديد من المواد ذات الهياكل المفتوحة. في الجليد السداسي ($I_h$)، يحدث التمدد الحراري السالب تحت درجة حرارية تبلغ حوالي 70 كلفن. وبينما تفشل نماذج الديناميكا الشبكية الكلاسيكية ومحاكاة الديناميكا الجزيئية (MD) في إعادة إنتاج هذه الظاهرة الشاذة، تنجح محاكاة الديناميكا الجزيئية للمسار المتكامل (PIMD) التي تتضمن التأثيرات الكمية النووية (NQE). يشير هذا إلى أن التمدد الحراري السالب في الجليد هو ظاهرة كمية نووية. ومع ذلك، لا يزال هناك غموض حرج: هل ينشأ هذا التأثير من التقلبات الكمية المحلية لروابط (O–H) الفردية، أم أنه خاصية جماعية لشبكة الروابط الهيدروجينية؟ يتطلب التمييز بين هذين الآليتين نظاماً يحافظ على التنسيق رباعي الأوجه المحلي مع تغيير طوبولوجيا التراكم بعيدة المدى. وقبل هذا العمل، كان من الصعب إجراء مثل هذه المقارنة بسبب انتشار اضطراب التراكم في الجليد المكعب ($I_c$)، مما يؤدي إلى تشويه الديناميكا الشبكية.

المنهجية
عالج المؤلفون هذه المسألة من خلال الجمع بين حيود النيوترونات بالمسحوق والمحاكاة المتقدمة على الجليد المكعب ($I_c$) الخالي من اضطراب التراكم.

• تحضير العينة: تم تحضير $I_c$ نقي عبر نزع الغاز التوتوبوتاتيكي (topotactic degassing) من هيدرات الهيدرات $C_2$. تنتج هذه الطريقة طوراً مكعباً مستقراً (metastable) بوجود شبكة أكسجين ذات بنية ماسية ($Fd\bar{3}m$) واضطراب تراكم ضئيل، وهي مستقرة حتى حوالي 200 كلفن.
• القياسات التجريبية: جُمعت أنماط حيود النيوترونات بين 50 كلفن و200 كلفن باستخدام مقياس الحيود D20. استُخرجت معاملات الشبكة عبر تنقية "لو بيل" (Le Bail refinement)، كما حُددت معاملات الإزاحة الحرارية متباينة الخواص عبر تنقية "ريتفيلد" (Rietveld refinement).
• المحاكاة: أُجريت محاكاة الديناميكا الجزيئية الكلاسيكية (MD) والديناميكا الجزيئية للمسار المتكامل (PIMD) باستخدام جهد MB-pol متعدد الأجسام. عالجت محاكاة PIMD التأثيرات الكمية النووية (NQE) صراحةً من خلال تمثيلات البوليمر الحلقي.
• التحليل النظري:
  • التوزيع البروتوني الكمي: حلل الباحثون نصف قطر الدوران ($R_g$) للبوليمرات البروتونية الحلقية لقياس مدى تباين توزيع البروتون الكمي.
  • الجهود الفعالة: أنشأ المؤلفون نموذج تضمين (embedding model) لعزل مساهمة ديناميكيات البروتون المنفرد في جهد الأكسجين-أكسجين البيني الفعال ($V_{OO}$)، ومقارنة البيئات الكلاسيكية والكمية.
  • تحليل الفونونات: حُسبت معاملات غرونايسن (Grüneisen parameters) المحددة للأنماط عبر تشتت الفونون التوافقي والفونونات الكمية غير التوافقية (باستخدام دالات الارتباط في PIMD) لتحديد الأنماط الاهتزازية المحددة التي تقود تغير الحجم.

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

1. شمولية التمدد الحراري السالب: تكشف البيانات التجريبية أن $I_c$ الخالي من اضطراب التراكم يظهر حداً أقصى للكثافة بالقرب من 70 كلفن، وهو ما يطابق سلوك الجليد السداسي ($I_h$) بشكل وثيق. وعلى الرغم من اختلاف تسلسلات التراكم طويلة المدى، فإن كلا المتعدد الأشكال يظهران ملفات تعريف كثافة متطابقة تقريباً. يثبت هذا أن التمدد الحراري السالب هو خاصية جوهرية لشبكة الروابط الهيدروجينية رباعية الأوجه المفتوحة، وغير حساس بشكل كبير لطوبولوجيا التراكم.
2. الأصل الكمي: تتنبأ محاكاة MD الكلاسيكية بزيادة رتيبة في الكثافة عند التبريد لكلا المتعددين، وتفشل في التقاط الحد الأقصى للكثافة. فقط محاكاة PIMD، التي تتضمن التأثيرات الكمية النووية (NQE)، تعيد إنتاج حد الكثافة الأقصى تجريبياً بدقة نسبية خطأ تبلغ حوالي 0.2%.
3. تباين الخواص والاستجابة الشبكية: يظهر تحليل المسارات الكمية للبروتونات تبايناً ملحوظاً: البروتونات أكثر عدم تمركز بمقدار 1.5 مرة في الاتجاهات المستعرضة مقارنة بالاتجاهات الطولية لروابط (O–H) التساهمية. والأهم من ذلك، أن درجة الحرارة التي يصل فيها معدل نسبة التباين المستعرض إلى الطولي ($R_{g,z}/R_{g,x}$) إلى ذروته تتزامن مع درجة حرارة الحد الأقصى للكثافة. وهذا يربط انكماش الشبكة العياني مباشرة بالحالة الكمية متباينة الخواص للبروتونات.
4. التأثيرات الجماعية مقابل المحلية: تُظهر حسابات جهد (O–O) الفعال ($V_{OO}$) أن معاملة البروتون المضمن فقط كجسيم كمي (ضمن بيئة كلاسيكية) غير كافية لإعادة إنتاج الكثافة الشاذة. يتطلب الأمر تحقيقاً دقيقاً معاملة كل من البروتون وبيئته المحيطة كجسيمات كمية. وهذا يؤكد أن التمدد الحراري السالب هو تأثير كمي جماعي، وليس مجرد تصحيح موضعي لنقطة الصفر.
5. الآلية عبر الأنماط المستعرضة: يحدد تحليل الفونونات الأنماط الصوتية المستعرضة (TA) والأنماط البصرية المستمرة (TO) منخفضة التردد كدوافع رئيسية للتمدد الحراري السالب. تمتلك هذه الأنماط أكثر معاملات غرونايسن سالبية. ينشأ التمدد الحراري السالب من التسكين الحراري المتتابع لهذه الأنماط ذات معامل غرونايسن السالب وفقاً لإحصاءات بوز-أينشتاين. عند درجات الحرارة المنخفضة، تهيمن أنماط TA؛ ومع ارتفاع درجة الحرارة، تساهم أنماط TO منخفضة المستوى. فوق 70 كلفن تقريباً، يتجاوز الإثارة الحرارية للأنماط ذات التردد الأعلى ذات معاملات غرونايسن الموجبة المساهمات الانكماشية، مما يؤدي إلى تمدد حراري موجب قياسي.

الأهمية والادعاءات
يؤكد البحث أن التمدد الحراري السالب في جليد I هو ظاهرة كمية نووية جماعية محكومة بديناميكيات شبكة الروابط الهيدروجينية بدلاً من تقلبات الروابط المحلية. ويدعي المؤلفون أن تشابه سلوك التمدد الحراري السالب بين $I_c$ و $I_h$ يثبت أن الآلية مرتبطة بطوبولوجيا الشبكة رباعية الأوجه المشتركة. وتوضح الدراسة أن معاملات غرونايسن السالبة هي شرط ضروري ولكنه غير كافٍ للتمدد الحراري السالب؛ إذ يعد الوزن الإحصائي الكمي لفروع الفونونات (إحصاءات بوز-أينشتاين) أمراً ضرورياً، مما يفسر فشل النماذج الكلاسيكية. علاوة على ذلك، يوفر هذا العمل نظاماً مرجعياً بسيطاً وشفافاً كيميائياً (الجليد $I_c$ الخالي من اضطراب التراكم) لفهم الشذوذات الديناميكية الحرارية المدفوعة كمياً في الشبكات رباعية الأوجه المفتوحة الأخرى.