Mathematical Modeling of the Ni(OH)₂/NiOOH Reaction on a Nickel Rotating Disk Electrode: A Foundation for Studying Organic Electrooxidation

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

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

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

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

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

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