Self-Folding and Self-Scrolling Mechanisms of Edge-Deformed Graphene Sheets: A Molecular Dynamics Study
تكشف دراسة الديناميكا الجزيئية هذه أن سلوكيات الطي الذاتي واللف الذاتي لصفائح الغرافين المشوهة عند الحواف إلى طيات نانوية أو لفافات نانوية قائمة على الغرافين تتحدد بشكل حاسم من خلال زوايا التواء الحواف الأولية ودرجة الحرارة، مع وجود تكوينات محددة تؤدي إلى مورفولوجيات أحادية الطي، أو متعددة الطيات، أو ثنائية الحالة.
المؤلفون الأصليون:Marcelo Lopes Pereira Junior, Luiz Antonio Ribeiro Junior
تخيل ورقة من مادة رقيقة للغاية لدرجة أنها ليست سوى طبقة واحدة من الذرات، ومع ذلك فهي قوية بما يكفي لتحمل دراجة هوائية ومرنة بما يكفي لتُنسج في شكل نسيج. هذه هي الجرافين، وهي شكل من أشكال الكربون الذي استحوذ على خيال العلماء لقدرتها الفريدة على الانحناء والالتواء والتدحرج دون انكسار. وبينما نفكر غالبًا في مثل هذه المواد على أنها مسطحة وثابتة، فإن للطبيعة طريقة لطيها. ففي ظل الظروف المناسبة، يمكن لهذه الصفائح أن تلتف تلقائيًا في لولبات ضيقة أو تنطوي فوق بعضها البعض مثل قطعة من الورق، مما يخلق أشكالًا جديدة ذات خصائص كهربائية وميكانيكية مختلفة. إن فهم كيفية ولماذا تحدث هذه التحولات بدقة أمر بالغ الأهمية للمهندسين الذين يأملون في بناء آلات متناهية الصغر أو بطاريات متطورة من الكربون، لكن العملية تحدث بسرعة كبيرة وعلى مقياس صغير جدًا بحيث لا تستطيع حتى أقوى المجاهر مراقبتها في الوقت الفعلي.
ولحل هذا اللغز، لجأ الباحثان مارسيلو لوبس بيريرا جونيور ولويز أنطونيو ريبيرو جونيور إلى أداة قوية تسمى محاكاة الديناميكا الجزيئية. فبدلاً من بناء نماذج فيزيائية، استخدما حاسوبًا فائق القدرة لإنشاء نسخة افتراضية من صفيحة الجرافين، مما سمح لهما بمراقبة حركة كل ذرة كما لو كانا يشاهدان فيلمًا عالي السرعة. بدآ بصفائح تم لف حوافها مسبقًا، لمحاكاة إعداد تجريبي محدد حيث يتم غزل طبقات الكربون لتشكيل هياكل معقدة. ومن خلال ضبط زاوية هذه الالتواءات وتشغيل المحاكاة عند درجات حرارة مختلفة، استطاعا رؤية كيفية تفاعل الصفيحة بدقة عبر الزمن، وتتبع التغيرات في الطاقة التي تدفع المادة لتغيير شكلها.
اكتشف الباحثان أن مقدار الالتواء عند الحافة يعمل كمفتاح يحدد الشكل النهائي. فعندما تم التواء الحواف بمقدار متوسط، لم تنطوِ الصفيحة على نفسها؛ بل انطوت فوق بعضها البعض، مكونةً هيكلًا مرتبًا ومتراكمًا. حدث هذا الطي بسرعة هائلة، غالبًا في غض_ن أجزاء من تريليونات من الثانية. ومع ذلك، عندما زادت زاوية الالتواء عن حد معين، تغير السلوك تمامًا؛ حيث بدأت الصفيحة تلتف في شكل لولبي، يشبه المخطوطة أو جريدة ملفوفة. كان هذا التحول إلى شكل لولبي سريعًا أيضًا، ولكنه تطلب التواءً أوليًا أقوى للتغلب على الصلابة الطبيعية لذرات الكربون. ووجد الفريق أنه إذا كان الالتواء طفيفًا جدًا، فإن الصفيحة ببساطة تعود إلى حالتها المسطحة، غير قادرة على التغلب على صلابتها لتشكيل شكل جديد.
كان أحد أكثر النتائج إثارة للدهشة هو كيفية تأثير درجة الحرارة على هذه التغييرات. ففي درجات الحرارة المنخفضة، استطاعت الصفيحة تشكيل هياكل معقدة، مثل شكل مطوي مزدوج حيث تلتقي الحواف وتترابط معًا. ولكن مع تسخين المحاكاة، أدى الاهتزاز الشديد للذرات إلى تعطيل هذه الروابط الدقيقة. وفي بعض الحالات، بدلاً من أن تنطوي الصفيحة لتصبح طبقتين في درجة حرارة الغرفة، انهارت لتصبح طبقة واحدة عند التسخين، لأن الطاقة الحرارية كانت قوية جدًا بحيث لا تسمح للطي المزدوج بالتماسك. ومع ذلك، أثبتت الأشكال اللولبية أنها أكثر مرونة؛ فبمجرد أن تشكلت الجرافين في شكل لولبي، كانت تميل للبقاء على حالها رغم الحرارة.
حدثت النتيجة الأكثر غرابة ربما عندما قام الباحثون بلف حافتي الصفيحة بمقدار مختلف. في هذا السيناريو، لم تختر الجرافين مسارًا واحدًا فقط، بل خلقت هيكلًا هجينًا حيث انطى جانب بينما التف الجانب الآخر في شكل لولبي. هذا الوضع المزدوج، الذي يجمع بين قسم مطوي وقسم لولبي في نفس قطعة المادة، لم يسبق ملاحظته من قبل. أظهرت المحاكاة أن هذا الشكل المختلط كان مستقرًا ويمكن أن يتشكل بغض النظر عن درجة الحرارة، مما يشير إلى أن المادة لديها ميل طبيعي لإيجاد توازن بين هذين الشكلين المختلفين عندما تُدفع في اتجاهات متضاربة.
كشفت الدراسة أيضًا أن الترتيب المحدد للذرات على طول حافة الصفيحة لم يغير النتيجة. فسواء كان ترتيب الحافة في خط مستقيم أو في نمط متعرج، تصرفت الصفيحة بنفس الطريقة، أي بالانطواء أو الالتفاف بناءً على زاوية الالتواء ودرجة الحرارة فقط. وهذا يشير إلى أن الآلية التي تقود هذه التغييرات عالمية لهذا النوع من مادة الكربون، وتعتمد على الشكل العام والطاقة أكثر من اعتمادها على التفاصيل المجهرية للحافة. ومن خلال رسم خرائط لهذه المسارات، قدم الباحثون صورة واضحة لكيفية سلوك الجرافين المتأثر بتشوه الحواف، مما يوفر دليلاً للتجارب المستقبلية. لقد أظهروا أنه من خلال التحكم الدقيق في الالتواء عند الحواف، من الممكن توجيه المادة لتصبح نوعًا محددًا من الهياكل المطوية أو اللولبية، مما يفتح الباب لتصميم مواد قائمة على الكربون بأشكال مخصصة ودقيقة.
بيان المشكلة تمثل الطيات النانوية (GNFs) واللفات النانوية (GNSs) القائمة على الجرافين تعديلات هيكلية هامة لطبقة الجرافين الأحادية (SLG)، مما يؤدي إلى تغييرات في الخصائص الميكانيكية والإلكترونية البصرية. وبينما رصدت الدراسات التجريبية ظواهر الطي واللف الذاتي الناتجة عن تشوهات الحواف، واستخدمت الأعمال النظرية ميكانيكا الاستمرارية أو الديناميكا الجزيئية (MD) لدراسة هذه العمليات، إلا أن الأطر الزمنية المحددة، والحالات شبه المستقرة الوسيطة، والمورفولوجيات النهائية الناتجة عن الجرافين أحادي الطبقة المشوه عند الحواف لا تزال غير مفهومة تماماً. غالباً ما اعتمدت التحقيقات النظرية السابقة على عوامل خارجية (مثل أنابيب الكربون النانوية أو الأشرطة النانوية) لتحفيز التشوه، بدلاً من التركيز على آليات التشوه الذاتي المدفوعة حصرياً بهندسة الحواف.
المنهجية استخدم المؤلفون عمليات محاكاة الديناميكا الجزيئية التفاعلية الذرية الكاملة (MD) باستخدام حقل القوة ReaxFF، المنفذ ضمن حزمة برمجيات LAMMPS. استخدمت الدراسة مجموعة معاملات C/H/O التي طورها Chenoweth وآخرون.
النماذج: كانت التكوينات الأولية عبارة عن صفائح جرافين ذات حواف "أرمرتشير" (armchair) ملفوفة (AESG) بأبعاد 100 Å × 200 Å (7,776 ذرة). تضمنت هذه النماذج زوايا التواء محددة (ϕ و θ) طُبقت على الحواف لمحاكاة التجارب المختبرية لتشكيل الصفائح مزدوجة اللف (biscrolled). وشملت زوايا الالتواء التي تم استقصاؤها القيم: −3π,−2π,0,2π, و 3π.
ظروف المحاكاة: أُجريت عمليات المحاكاة في مجموعة NVT باستخدام منظم حرارة Nosé-Hoover. تم تغيير درجة الحرارة من 100 K إلى 1000 K بزيادات قدرها 100 K. تم دمج معادلات الحركة باستخدام خوارزمية velocity-Verlet بخطوة زمنية قدرها 0.25 fs لإجمالي مدة بلغت 250 ps.
التحليل المقارن: تضمنت الدراسة أيضاً محاكاة لجرافين ذي حواف "زيجزاج" (ZESG) لتحديد ما إذا كانت طبوغرافيا الحافة (armchair مقابل zigzag) تؤثر على آليات التشوه.
النتائج الرئيسية كشفت عمليات المحاكاة عن نتائج مورفولوجية متميزة بناءً على زوايا الالتواء الأولية ودرجة الحرارة:
زوايا الالتواء ≤π: كانت زوايا الالتواء الأقل من π غير كافية لتحفيز التشوه الذاتي. في هذه الحالات، لم تستطع تفاعلات فان دير فالس بين الطبقات التغلب على صلابة الانحناء لـ SLG، وعادت الهياكل لتسترخي إلى شكل الجرافين أحادي الطبقة التقليدي.
زوايا الالتواء >2π (الطي):
AESG(0, 2π): تطورت إلى مورفولوجيا مطوية لمرة واحدة ومكدسة بالكامل (GNF).
AESG(2π, 2π): تطورت إلى مورفولوجيا مطوية مرتين ومكدسة بالكامل. ومع ذلك، كان هذا الهيكل المحدد المطوي مرتين مستقراً فقط تحت درجة حرارة حرجة تبلغ حوالي 900 K. عند 1000 K، تغلبت التقلبات الحرارية على بئر طاقة الوضع المطلوب للحالة المطوية مرتين، مما أدى إلى ظهور مورفولوجيا مطوية لمرة واحدة بدلاً من ذلك.
الحركية (Kinetics): حدث الانتقال من الحواف الملتوية إلى الهياكل المطوية خلال فترة فائقة السرعة (حوالي 2 ps عند 300 K، وأقل من 1 ps عند 1000 K).
زوايا الالتواء >2π (اللف):
AESG(0, 3π) و AESG(3π, 3π): انتقلت هذه التكوينات بسرعة إلى لفات نانوية (GNSs) مستقرة بعد فترة استرخاء أولية.
ديناميكيات AESG(3π, 3π): تشكلت حالة "بسكرويل" (biscroll) شبه مستقرة في البداية نتيجة للتفاعل بين اللف الجزئي من الجانب الأيسر والأيمن. انهارت هذه الحالة في النهاية لتتحول إلى GNS واحدة. حدث تشكل الـ GNS النهائي خلال 76 ps.
الاستقرار: على عكس حالات الطي، حدث تشكل GNSs المستقرة من زوايا الالتواء العالية هذه حتى في درجات الحرارة العالية (حتى 1000 K).
تكون الحالة المزدوجة:
AESG(2π, 3π): تطور هذا النظام، الذي يمتلك زوايا التواء مختلفة على الحواف المتقابلة، إلى حالة مزدوجة فريدة تجمع بين الهياكل المطوية والملفوفة. ووجد أن هذه العملية مستقلة عن درجة الحرارة.
زوايا الالتواء السالبة:
AESG(-2π, 2π): تطور إلى هيكل منهار بحواف مطوية مرتين، استقر بفعل تكون روابط كربون-كربون تساهمية نتيجة التقارب الشديد بين الحواف.
AESG(-3π, 3π): ينتج حالة مزدوجة تجمع بين الهياكل المطوية والملفوفة، بشكل مشابه للنتائج المعروضة لحالة زاوية الالتواء الموجبة.
استقلالية طبوغرافيا الحافة: أظهرت المقارنات بين شبكات "أرمرتشير" (AESG) و"زيجزاج" (ZESG) أن طبوغرافيا الحافة لا تؤثر على آليات الطي أو اللف الأساسية؛ حيث كانت التطورات الزمنية لطاقة الوضع والمورفولوجيات النهائية متشابهة لكلا نوعي الشبكة.
الطاقة (Energetics): أكد تحليل طاقة الوضع أن التكوينات الملفوفة والمطوية مفضلة طاقياً. أظهرت منحنيات طاقة الوضع اتجاهاً تنازلياً أثناء تشكل GNSs، واستقرت بعد 40–100 ps اعتماداً على النظام.
الأهمية والادعاءات يزعم المؤلفون أن هذه الدراسة توفر رؤى حاسمة حول التشوه الذاتي لصفائح الجرافين المنحنية، وتحديداً معالجة الأطر الزمنية والحالات الوسيطة التي يصعب رصدها تجريبياً. تسلط الورقة الضوء على عدة نتائج جديدة:
سجلت أول رصد لـ هيكل منهار بحواف مطوية مرتين مستقر بفعل تكون روابط تساهمية في حالة AESG(-2π, 2π).
حددت حالة مزدوجة حيث تتواجد الهياكل المطوية والملفوفة معاً في المورفولوجيا النهائية (لوحظ في AESG(2π, 3π) و AESG(-3π, 3$\pi) وغيرها)، وهو تكوين لم يسبق الإبلاغ عنه.
أثبتت أن عملية التشوه الذاتي مستقلة عن درجة الحرارة بالنسبة للحالة المزدوجة واللف ذي زاوية الالتواء العالية، بينما المورفولوجيا المحددة للطي ثنائي الطي حساسة لدرجة الحرارة (مستقرة فقط تحت ~900 K).
أوضحت الدراسة أن طبوغرافيا الحافة (armchair مقابل zigzag) ليست عاملاً محدداً لآليات الطي أو اللف الذاتي.
يعمل هذا العمل على سد الفجوة بين الملاحظات التجريبية للجرافين المشوه عند الحواف والفهم النظري، مما يقدم رؤية ذرية مفصلة لتشكل GNFs و GNSs دون الحاجة إلى عوامل تحفيز خارجية.