Landau-de Gennes numerical simulation of nematic liquid crystals utilizing radial basis functions

إليك ملخص تقني مفصل للورقة البحثية بعنوان "المحاكاة العددية لنموذج لاندا-دي جين للبلورات السائلة النيماتية باستخدام الدوال الشعاعية الأساسية" من إعداد جين-شنغ وو وإيفان إي. سماليوك.

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

تعتمد المحاكاة العددية التقليدية للبلورات السائلة النيماتية (LCs) باستخدام نظرية لاندا-دي جين (LdG) عادةً على طرق الفروق المحدودة القائمة على الشبكات المكعبة. وبين أنها فعالة في الأشكال الهندسية البسيطة، إلا أن هذه الطرق تواجه قيوداً كبيرة عند نمذجة الأنظمة المعقدة:

• الجمود الهندسي: تعاني الشبكات المكعبة في التوافق مع أشكال الحدود التعسفية أو المعقدة (مثل الجسيمات الغروية ذات الطوبولوجيا غير التقليدية مثل الأجسام ذات المقابض - handlebodies) دون إدخال تماثلات اصطناعية أو الحاجة إلى توليد شبكة معقدة.
• عدم الكفاءة الحسابية: يتطلب حل الميزات دقيقة المقياس، مثل العيوب المفردة (الخلوات - disclinations)، شبكات عالية الدقة موحدة عبر النطاق بأكمله، مما يؤدي إلى تكاليف حسابية واستهلاك مفرط للذاكرة.
• التطبيق المحدود: بينما تم تطبيق الدوال الشعاعية الأساسية (RBFs) بنجاح في مجالات مثل علم الزلازل وديناميكا السوائل، إلا أن تطبيقها على النمذجة الموترية لـ LdG للبلورات السائلة (تحديداً لأنظمة Q-tensor) لم يكن قد تم إثباته قبل هذا العمل.

2. المنهجية

طوّر المؤلفون إطار عمل عددي خالٍ من الشبكة يعتمد على طريقة الفروق المحدودة المولدة بالدوال الشعاعية الأساسية (RBFFD) لحل معادلات لاندا-دي جين.

• المعادلات الحاكمة: يعمل النظام على تقليل إجمالي دالة الطاقة الحرة ($F_{total}$)، والتي تشمل:

  • طاقة الحجم الحرة: توسع لاندا-دي جين (المصطلحات الحرارية $A, B, C$) والطاقة المرنة (تقريب الثابت الواحد $L$).
  • طاقة السطح: شروط تثبيت السطح ذات التدهور المخروطي عند الواجهات الغروية.
  • التقليل (Minimization): يتم إيجاد الهياكل المتوازنة عبر الانحدار المتدرج ($\frac{dQ}{dt} \propto -\frac{\delta F}{\delta Q}$) باستخدام مزيج من خوارزميات Adaptive Runge-Kutta (ARK23) وFast Inertia Relaxation Engine (FIRE).

• تنفيذ RBFFD:

  • الدوال الأساسية: يتم تقريب مكونات الـ Q-tensor كمزيج خطي من الدوال الشعاعية الأساسية (تحديداً الدوال الغاوسية $\phi(r) = e^{-\epsilon|r-r_i|^2}$).
  • تقريب المشتق: يتم حساب معاملات الفروق المحدودة ($d_j$) عن طريق حل نظام خطي حيث يتم تطبيق المؤثر على الدوال الأساسية. ولضمان الاستقرار والدقة العالية (من الدرجة الرابعة إلى الثامنة)، تدمج الطريقة حدوداً متعددة إضافية تصل إلى الدرجة الثانية.
  • الطبيعة الخالية من الشبكة: توضع العقد بحرية في الفضاء، غير مقيدة بهندسة الشبكة، مما يسمح بأشكال نطاق تعسفية.

• تحسين العقد التكيفي:

  • ابتكار محوري هو مخطط تكيفي حيث يتم ضبط كثافة العقد ديناميكياً بناءً على التشوه المرن المحلي.
  • تُعامل العقد كجسيمات ذات "شحنات" تتناسب مع المسافة المثالية للفصل ($r_0 \propto (\partial_i \partial_j Q_{ab})^{-1/2}$).
  • تعمل قوى التنافر الشبيهة بالكهرباء الساكنة على نقل العقد نحو المناطق ذات التدرج العالي (أنوية العيوب)، مما يحقق دقة تترا-وحتى 1.5-7 نانومتر بالقرب من العيوب مع الحفاظ على تباعد أكثر تباعداً في أماكن أخرى.

• التحقق التجريبي:

  • تم التحقق من صحة المحاكاة مقابل البيانات التجريبية باستخدام مجهر الفلورة المستحث بـ 3 فوتونات (3PEF-PM) ومجهر الضوء المستقطب التقليدي (POM).
  • شملت العينات بلورات سائلة محصورة في جسيمات غروية ذات أشكال "handlebody" (من الجنس 1 و5) وقطرات داخل مصفوفات بوليمرية معقدة.

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

• أول تنفيذ لـ LdG-RBF الموتري: يقدم هذا العمل أول تطبيق ناجح لـ RBFFD للنمذجة الموترية لـ LdG للبلورات السائلة النيماتية، متجاوزاً نماذج Frank-Oseen الاتجاهية السابقة.
• التحسين التكيفي الخالي من الشبكة: يسمح إدخال مخطط تحسين عقد تكيفي قائم على الكهرباء الساكنة بحل العيوب المفردة تلقائياً دون الحاجة لتوليد شبكة يدوية أو استخدام طرق الشبكات متعددة المستويات.
• المرونة الهندسية: تتعامل الطريقة بنجاح مع الطوبولوجيا المعقدة (مثل الأجسام ذات المقابض من الجنس 5) وشروط الحدود التعسفية التي يصعب نمذجتها بالشبكات المكعبة.
• الكفاءة الحسابية: تظهر الطريقة تحسينات كبيرة في السرعة وانخفاضاً في استهلاك الذاكن مقارنة بالطرق الشبكية التقليدية.

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

• دقة العيوب: التقطت المحاكاة بدقة القيود الطوبولوجية للهياكل المعقدة. على سبيل المثال:
  • الجنس-1 (Torus): تنبأت بشكل صحيح بوجود عيبين من النوع $+1$ وعيبين من النوع $-1$ (إجمالي عدد اللفات 0).
  • الجسم ذو المقابض (الجنس-5): تنبأت بإجمالي عدد لفات $-8$(3 عيوب موجبة و11 عيباً سالباً)، وهو ما يتوافق مع نظرية كرة الشعر المعممة ($\chi = 2-2g$).
  • نجح المخطط التكيفي في حل انقسام العيوب الصحيحة إلى أزواج من العيوب نصف الصحيحة بالقرب من السطح.
• المتعددات الأقطاب المرنة: أعادت الطريقة إنتاج أنماط التشويه متعدد الأقطاب المعروفة حول الجسيمات الغروية الكروية مع ظروف تثبيت مختلفة:
  • التثبيت العمودي ($\theta_e = 0^\circ$) والمسطح ($\theta_e = 90^\circ$): تشوهات رباعية الأقطاب ($l=2$).
  • التثبيت المخروطي ($\theta_e = 45^\circ$): تشوهات سداسية عشرية الأقطاب ($l=4$).
• المقارنة المرجعية: أظهرت المقارنات مع محاكاة الشبكة المنتظمة ما يلي:
  • الدقة: توافق ممتاز في تشوهات مجال الاتجاه، ومواقع العيوب، وعزوم متعددات الأقطب ($b_l$).
  • الكفاءة: تطلبت طريقة RBFFD حوالي 50% وقتاً أقل في الحساب لنفس عدد التكرارات عند استخدام مجموعة عقد ثابتة.
  • اقتصاد العقد: حققت RBFFD نفس التمثيل الهندسي بعدد أقل بكثير من العقد (على سبيل المثال، 800 ألف عقدة مقابل أكثر من 1.7 مليون عقدة لشبكة مكعبة $120^3$) من خلال إلغاء العقد غير النشطة داخل الجسيمات الصلبة.
• التوازي: الطريقة قابلة جداً للتوازي (GPU/CPU)، حيث أن حساب معاملات التفاضل لكل عقدة مستقل عن غيره.

5. الأهمية

تؤسس هذه الدراسة لـ RBFFD كأداة قوية ومتعددة الاستخدامات لمحاكاة أنظمة المادة اللينة ذات الهندسات المعقدة.

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

في الختام، يثبت وو وسماليوك أن طريقة RBFFD توفر توازناً متفوقاً بين المرونة الهندسية، والكفاءة الحسابية، والدقة العالية من الدرجة، لا سيما للأنظمة التي تتضمن عيوباً مفردة وحدوداً معقدة.