Solid-angle based nearest-neighbor algorithm adapted for systems with low coordination number

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

لقد طور فريق من الباحثين من جامعة أوتريخت وجامعة باريس ساكلي طريقة محسنة لحل هذه المشكلة المحددة. فقد ركزوا على تقنية موجودة تسمى خوارزمية أقرب جار القائمة على الزاوية الصلبة (solid-angle-based nearest-neighbor algorithm)، والتي تحدد الجيران من خلال التحقق من مقدار الرؤية المحيطة بالجسيم التي يحجبها رفاقه. وبينما تعمل هذه الطريقة بشكل جيد في الأنظمة الكثيفة والفوضوية، إلا أنها تميل إلى ارتكاب خطأ منهجي في الشبكات المفتعة والمنظمة حيث تكون الجسيمات قليلة ومتباعدة. وفي هذه البيئات منخفضة الكثافة، غالبًا ما تصل الخوارزمية الأصلية إلى مسافات بعيدة جدًا، فتحصي جسيمات تنتمي إلى طبقة الجيران التالية وكأنها جزء من الدائرة الأولى. وقد قدم الباحثون تصحيحًا هندسيًا لإصلاح هذا التكرار في العد دون إضافة أي إعدادات قابلة للضبط جديدة إلى الحساب. ونجحت طريقتهم المعدلة، التي يسمونها mSANN، في تحديد العدد الصحيح للجيران في الهياكل المعقدة التي تتراوح من أنماط خلية النحل إلى بلورات الماس، مما يوفر خريطة أكثر دقة للعالم المجهري.

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

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

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

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

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

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

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