Vacancy defects in square-triangle tilings and their implications for quasicrystals formed by square-shoulder particles

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

بيان المشكلة
بينما تُلاحظ أشباه البلورات (Quasicrystals) ذات بلاطات المربع والمثلث على نطاق واسع في أنظمة المادة اللينة، إلا أنها تحتوي بشكل مميز على كثافة عالية من العيوب النقطية، والتي تظهر عادةً كبلاطات سداسية إضافية. لا تزال الدور الديناميكي الحراري لهذه العيوب غير مفهوم بشكل جيد؛ حيث لا يزال من غير الواضح ما إذا كان تركيزها العالي ناتجاً عن احتجاز حركي (Kinetic trapping) أثناء التجميع الذاتي، أم أنه سمة متأصلة في الطور المتوازن الذي يساهم في استقرار شبه البلورة. يبحث هذا العمل في المساهمة الديناميكية الحرارية لهذه العيوب في استقرار شبه بلورة المربع والمثلث ذات التماثل الاثني عشر.

المنهجية
تستخدم الدراسة نهجاً يتكون من جزأين يجمع بين نموذج بلاط قائم على الشبكة ونموذج جسيمات مجهري:

1. محاكاة مونت كارلو القائمة على البلاطات:

   • النموذج: يستخدم المؤلفون نموذجاً مثالياً قائماً على البلاطات يتكون من مربعات، ومثلثات، ونوعين من بلاطات العيوب السداسية غير المنتظمة: "الدروع" (Shields) و"البيض" (Eggs) (والأخير كيرالي، حيث يوجد كـ $Egg_L$ و $Egg_R$). تنشأ هذه العيوب طبيعياً عند إزالة رأس (Vertex) من بلاط مربع-مثلث مثالي.
   • الخوارزمية: تستخدم خوارزمية مونت كارلو (MC) عينات من الفضاء التشكلي عبر تحديد "الحواف المتحركة" (Edges) (الحواف الموجودة على الحدود أو جزء من بلاطة عيب) وتدويرها بمقدار $\pi/6$. يسمح هذا بإعادة الترتيب المحلي وتغيير أنواع البلاطات دون إنشاء فجوات أو تداخلات.
   • المجموعة الإحصائية: تُجرى المحاكاة في مجموعة $N\gamma T$ (عدد ثابت من الرؤوس $N$، وتوتر خطي $\gamma$، ودرجة حرارة $T$). يتم تعيين عقوبة طاقة حرة $\epsilon_D$ لكل بلاطة عيب للتحكم في تركيز العيوب.
   • حساب الإنتروبيا: يتم حساب فرق الإنتروبيا التشكيلية ($\Delta S$) بين البلاطات المعيبة والخالية من العيوب باستخدام التكامل الديناميكي الحراري. تحلل الدراسة مقياس الإنتروبيا مع حجم النظام وتركيز العيوب، مع التمييز بين مساهمات العيوب الفردية وإنتروبيا الخلط التوليفي.

2. نموذج الجسيمات المجهري:

   • النظام: جسيمات "اللب-التاج" (Core-corona) تتفاعل عبر جهد ذو كتف مربع (نواة صلبة $\sigma$، وكتف تنافري ناعم بارتفاع $\epsilon$ وعرض $\delta=1.4\sigma$).
   • حسابات الطاقة الحرة: يحسب المؤلفون الطاقة الحرة الاهتزازية لشبه البلورة والأطوار الدورية المنافسة (البلورات المربعة والسداسية) باستخدام تكامل أينشتاين والتكامل الديناميكي الحراري.
   • تكلفة العيب: لتحديد تكلفة الطاقة الحرة الاهتزازية لإنشاء فراغ، قام المؤلفون بتكييف طرق مستخدمة للبلورات الدورية. لقد أخذوا في الاعتبار حقيقة أن فراغاً واحداً في شبه البلورة يتحول إلى زوج من العيوب المستقلة (دروع أو بيض). ولمنع انتشار العيوب أثناء حسابات الطاقة الحرة، تم تقييد حركة الجسيمات داخل خلايا "افتراضية" تشبه خلايا فورونوي (Voronoi) محددة بواسطة هندسة البلاط المحلية.
   • تركيز التوازن: يتم بناء طاقة غيبس الحرة الإجمالية من خلال الجمع بين الطاقة الحرة الاهتزازية (من نموذج الجسيمات) والإنتروبيا التشكيلية (من نموذج البلاط). تُشتق تركيزات العيوب عند التوازن عن طريق تقليل هذه الطاقة الحرة.

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

• طبيعة العيوب: إن إزالة جسيم واحد (رأس) من بلاط مربع-مثلث لا يخلق مجرد فراغ بسيط، بل يؤدي إلى إعادة ترتيب محلي ينتج عنه زوج من البلاطات السداسية: درع واحد وبيضة واحدة (سواء كانت يسارية أو يمينية).
• كسب الإنتروبيا التشكيلية:
  • يؤدي إدخال العيوب إلى كسب كبير في الإنتروبيا التشكيلية. في الحد الديناميكي الحراري، وُجد أن الإنتروبيا التشكيلية لكل رأس لنظام بتركيز عيوب صفر ($\beta\epsilon_D = 0$) هي $0.554(1) k_B$، وهي أعلى بنحو خمس مرات من بلاط المربع والمثلث الخالي من العيوب ($0.120 k_B$).
  • والأهم من ذلك، أن كسب الإنتروبيا ليس مجرد مجموع إنتروبيات العيوب الفردية. فبينما تمتلك عيوب الدرع والبيض الفردية ثوابت إنتروبيا سالبة مقارنة بعيب نقطي بسيط في بلورة دورية (مما يشير إلى أنها أقل تفضيلاً بمفردها)، فإن الخلط التوليفي لأنواع العيوب الثلاثة (Shield, $Egg_L, Egg_R$) يولد حداً إضافياً كبيراً من الإنتروبيا. هذا الخلط الإنتروبي يرفع إجمالي كسب الإنتروبيا إلى قيمة موجبة صافية، مما يعمل على تثبيت وجود العيوب.
• تركيزات العيوب عند التوازن:
  • بتطبيق هذه النتائج على نموذج جسيمات الكتف المربع، يتوقع المؤلفون تركيزات عيوب عند التوازن مرتفعة بشكل غير طبيعي مقارنة بالبلورات الدورية.
  • عند درجات الحرارة المنخفضة ($k_B T/\epsilon = 0.1$)، تكون تركيزات العيوب ضئيلة ($\sim 10^{-8}$). ومع ذلك، مع ارتفاع درجة الحرارة، تنخفض تكلفة الطاقة الحرة الاهتزازية، وتهيمن الإنتروبيا على العملية.
  • عند $k_B T/\epsilon = 0.3$، يصل إجمالي تركيز العيوب إلى حوالي $1\%$ ($\sim 0.01$ لكل موقع شبكي)، وهو أعلى بكثير من تركيز الفراغات النموذجي ($\sim 10^{-4}$) في بلورات الكرات الصلبة.
• استقرار الطور:
  • بدون تضمين الإنتروبيا التشكيلية، لا يكون شبه بلورة المربع والمثلث مستقراً ديناميكياً حرارياً مقارنة بوجود تعايش بين البلورات المربعة والسداسية.
  • عند تضمين الإنتروبيا التشكيلية للعيوب، تظهر منطقة مستقرة لشبه البلورة ذات الـ 12 طية. وبالتالي، فإن التركيز العالي للعيوب هو سمة ديناميكية حرارية متأصلة مطلوبة لاستقرار هذا الطور. كما أن وجود العيوب يزيح قليلاً الكيمياء الكلية للتعايش والضغوط، ولكنه يوسع نطاق الكثافة المتوازن لشبه البلورة بنحو $7\%$.

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