DOI: https://doi.org/10.1038/s41467-026-69231-4
PMID: https://pubmed.ncbi.nlm.nih.gov/41651854
تاريخ النشر: 2026-02-06
المؤلف: Jiaxing Peng وآخرون
الموضوع الرئيسي: مواد وتصنيع المكثفات الفائقة
نظرة عامة
في هذا البحث، تم إنشاء إطار نمذجة شامل للتحقيق في ديناميات تخزين الطاقة والشحن لأقطاب الكربون النانوية المسامية غير المتبلورة. يشير نموذج FEC المقترح إلى أن المسام الفائقة الصغر، بينما تظهر سعة أعلى، لديها موصلية أيونية أقل مقارنة بالمسام الصغيرة الأكبر. يُعزى هذا التباين في ديناميات الشحن إلى زيادة تكرار الاصطدامات بين الأيونات وجدران المسام داخل المسام الفائقة الصغر.
علاوة على ذلك، تم تطوير نموذج متعدد المقاييس للتنبؤ بمقاومة خلية سوبركاباسيتور الكربون المسامي، باستخدام محاكاة الديناميكا الجزيئية (MD). حقق هذا النموذج توافقًا كميًا مع القياسات الكهروكيميائية التجريبية، مما يساهم بفعالية في سد الفجوة بين المجهر الحسابي والمجهر التجريبي في سياق ديناميات الشحن. قد تكون المنهجيات المقدمة لنمذجة وتحليل الكربونات المسامية غير المتبلورة قابلة للتطبيق أيضًا في مجالات أخرى تتعلق بالوسائط المسامية، بما في ذلك التحفيز، وإزالة الملح الكهروكيميائية، وتخزين الهيدروجين، والتقاط الكربون، وفصل الغاز، واستعادة النفط.
مقدمة
في هذا القسم، يناقش المؤلفون أهمية مواد الكربون المسامي في تطبيقات تخزين الطاقة، لا سيما كأقطاب للسوبركاباسيتورات والبطاريات. تعتبر الخصائص الهيكلية للكربونات النانوية المسامية، مثل حجم المسام، والشكل، والترابط، والمساحة السطحية المحددة (SSA)، حاسمة لتعزيز قدرات تخزين الطاقة. تشير النتائج التجريبية إلى أن المسام ذات الحجم القريب من القطر الأيوني يمكن أن تحسن السعة، بينما تساهم المجالات الصغيرة الشبيهة بالجرافين في جدران المسام أيضًا بشكل إيجابي. كانت المحاكاة الجزيئية، لا سيما باستخدام نماذج المسام المثالية والديناميكا الجزيئية المروّدة (QMD)، أساسية في فهم السلوك السعوي على المستوى الذري. ومع ذلك، توجد قيود في النماذج الحالية، لا سيما فيما يتعلق بتجاهل تفاعلات فان der Waals والاعتماد على معلمات المحاكاة دون دمج البيانات التجريبية.
لمعالجة هذه القيود، طور المؤلفون نموذجًا واقعيًا للأقطاب غير المتبلورة مستندًا إلى تشتت الأشعة السينية بزاوية صغيرة (SAXS)، وامتصاص الغاز، وبيانات الكثافة للكربونات النانوية المسامية. تم استخدام هذا النموذج في محاكاة الديناميكا الجزيئية (MD) ذات الجهد الثابت مع إلكتروليت سائل أيوني، مما أسفر عن نتائج سعة تفاضلية تتماشى جيدًا مع القياسات التجريبية. تستخدم الدراسة أيضًا طريقة كرة فورونوي لتحليل نقل الأيونات داخل المسام ذات الأحجام المختلفة، كاشفة أن المسام الفائقة الصغر (<0.7 نانومتر) مسؤولة بشكل أساسي عن تخزين الشحن من خلال تبادل الأيونات، بينما تساهم المسام الصغيرة الأكبر (>0.7 نانومتر) بشكل ضئيل. كما يبرز المؤلفون أهمية ترابط المسام في تحديد السعة ويقترحون نماذج متعددة المقاييس لربط المحاكاة الجزيئية بالقياسات الماكروسكوبية، مما يعزز الفهم لديناميات الشحن في أقطاب الكربون المسامية.
طرق
في هذا القسم، يوضح المؤلفون الطرق المستخدمة في إعداد المواد وإجراء القياسات التجريبية ذات الصلة بدراستهم لأقطاب الكربون المسامية. تم الحصول على الكربون المسامي (YP-50F) من شركة Kuraray Chemical، بينما تم الحصول على مادة الربط من البولي تترافلوروإيثيلين (PTFE) والإلكتروليت من Aladdin و Lanzhou Yulu Fine Chemical Co. على التوالي. تم تحديد توزيع حجم المسام (PSD) والمساحة السطحية المحددة (SSA) للكربون المسامي باستخدام نظام Quantachrome Autosorb IQ، مع قياسات إيزوثرم امتصاص النيتروجين عند 77 كلفن بعد التنشيط عند 200 درجة مئوية. تم إجراء تجارب تشتت الأشعة السينية بزاوية صغيرة (SAXS) باستخدام مصدر أشعة سينية أحادية اللون لتحليل الخصائص الهيكلية لعينات الكربون.
شمل إعداد القطب الكربوني المسامي خلط مسحوق الكربون مع مادة الربط PTFE في الإيثانول، مما أسفر عن فيلم بسمك حوالي 150 ميكرومتر. تم إجراء محاكاة الديناميكا الجزيئية (MD) مع قطبين كربونيين مساميين مغمورين في إلكتروليت سائل أيوني، تحديدًا [EMIM][BF4]. استخدمت المحاكاة طريقة الجهد الثابت للحفاظ على الجهد بين الأقطاب، مع استخدام تنفيذ معزز بواسطة وحدة معالجة الرسوميات (GPU) من أجل الكفاءة. تم تحقيق توازن النظام عند 298 كلفن، وتم إجراء التحليلات على الحالة التوازنية بعد 150 نانوثانية من الاستقطاب. تم اشتقاق السعة التفاضلية من الشحن التوازني على ذرات القطب بالنسبة لجهد القطب، مما يوفر رؤى حول ديناميات الشحن وتوزيع الأيونات داخل الهيكل المسامي.
مناقشة
يقدم قسم المناقشة في ورقة البحث تحليلًا مفصلًا للنموذج الذري للكربونات المسامية غير المتبلورة، مع التركيز بشكل خاص على الكربون النانوي YP-50F. تم تحديد الهيكل ثلاثي الأبعاد باستخدام قياسات تشتت الأشعة السينية بزاوية صغيرة (SAXS)، كاشفًا عن قمتين كبيرتين للتشتت عند حوالي 3.4 و4.7 نانومتر$^{-1}$. كانت الطريقة الكلاسيكية للحقل العشوائي الغاوسي (cGRF) المستخدمة في البداية غير كافية، مما أدى إلى تطوير طريقة القمم المتعددة للتشتت GRF (mGRF)، التي نجحت في التقاط القمم المرصودة وإنتاج نموذج مصفوفة قطب ثلاثي الأبعاد مع مساحة سطح محددة (SSA) تبلغ 1573 م$^2$غ$^{-1}$، تتماشى بشكل وثيق مع القيم التجريبية. كما قامت الطريقة الهجينة للديناميكا الجزيئية العكسية ثلاثية الأبعاد (H3D-RMD) بمزيد من تحسين الهيكل الذري، مما أسفر عن SSA تبلغ 1694 م$^2$غ$^{-1}$، مما يحقق صحة النموذج مقابل قياسات امتصاص الغاز.
تستكشف الدراسة أيضًا قدرات تخزين الشحن لأقطاب الكربون النانوية المسامية من خلال محاكاة الديناميكا الجزيئية ذات الجهد الثابت، مما يظهر أن المسام الفائقة الصغر تظهر سعة تفاضلية جاذبية أعلى مقارنة بالمسام الصغيرة الكبيرة. يُعزى ذلك إلى التراكم التفضيلي للأنيونات في المسام الأصغر، مما يؤدي إلى قيمة إيجابية للجهد الصفري للشحن (PZC). تشير تحليلات الشحن المحلي للقطب (LEC) إلى أن المسام الفائقة الصغر تسهل تخزين شحن أكبر بسبب خصائصها الهيكلية، على الرغم من الموصلية الأيونية الأبطأ. تم نمذجة ديناميات الشحن باستخدام نموذج مقاومة متعدد المقاييس (MIM)، الذي يتنبأ بفعالية بنتائج طيف impedance الكهروكيميائي (EIS)، مما يربط المحاكاة الجزيئية بالبيانات التجريبية الماكروسكوبية. لا يعزز هذا الإطار الشامل للنمذجة فقط فهم تخزين الشحن في أقطاب الكربون النانوية المسامية، بل يقدم أيضًا منهجيات قابلة للتطبيق في مجالات متنوعة تتعلق بالوسائط المسامية.
DOI: https://doi.org/10.1038/s41467-026-69231-4
PMID: https://pubmed.ncbi.nlm.nih.gov/41651854
Publication Date: 2026-02-06
Author(s): Jiaxing Peng et al.
Primary Topic: Supercapacitor Materials and Fabrication
Overview
In this research, a comprehensive modeling framework was established to investigate the energy storage and charging dynamics of amorphous nanoporous carbon electrodes. The proposed FEC model indicates that ultramicropores, while exhibiting higher capacitance, have lower ionic conductivity compared to larger micropores. This discrepancy in charging dynamics is attributed to the increased frequency of collisions between ions and the pore walls within ultramicropores.
Furthermore, a multi-scale model was developed to predict the impedance of a porous carbon supercapacitor cell, utilizing molecular dynamics (MD) simulations. This model achieved quantitative agreement with experimental electrochemical measurements, effectively bridging the gap between computational microscopy and experimental macroscopy in the context of charging dynamics. The methodologies introduced for modeling and analyzing amorphous porous carbons may also be applicable to other domains involving porous media, including catalysis, electrochemical desalination, hydrogen storage, carbon capture, gas separation, and oil recovery.
Introduction
In this section, the authors discuss the significance of porous carbon materials in energy storage applications, particularly as electrodes for supercapacitors and batteries. The structural characteristics of nanoporous carbons, such as pore size, shape, connectivity, and specific surface area (SSA), are critical for enhancing energy storage capabilities. Experimental findings indicate that pores sized near the ionic diameter can improve capacitance, while smaller graphene-like domains in pore walls also contribute positively. Molecular simulations, particularly using ideal pore models and quenched molecular dynamics (QMD), have been instrumental in understanding the capacitive behavior at the atomic level. However, limitations exist in current models, particularly regarding the neglect of van der Waals interactions and the reliance on simulation parameters without incorporating experimental data.
To address these limitations, the authors developed a realistic amorphous electrode model informed by small angle x-ray scattering (SAXS), gas adsorption, and density data of nanoporous carbons. This model was utilized in constant-potential molecular dynamics (MD) simulations with an ionic liquid electrolyte, yielding differential capacitance results that align well with experimental measurements. The study further employs the Voronoi sphere method to analyze ion transport within different-sized pores, revealing that ultramicropores (<0.7 nm) are primarily responsible for charge storage through ion exchange, while larger micropores (>0.7 nm) contribute minimally. The authors also highlight the importance of pore connectivity in determining capacitance and propose multiscale models to link molecular simulations with macroscopic measurements, thereby enhancing the understanding of charging dynamics in porous carbon electrodes.
Methods
In this section, the authors detail the methods used for preparing materials and conducting experimental measurements relevant to their study of porous carbon electrodes. The porous carbon (YP-50F) was sourced from Kuraray Chemical, while the polytetrafluoroethylene (PTFE) binder and electrolyte were obtained from Aladdin and Lanzhou Yulu Fine Chemical Co., respectively. The porous carbon’s pore size distribution (PSD) and specific surface area (SSA) were characterized using a Quantachrome Autosorb IQ system, with nitrogen adsorption isotherms measured at 77 K after activation at 200 °C. Small-angle X-ray scattering (SAXS) experiments were performed using a monochromatic X-ray source to analyze the structural properties of the carbon samples.
The preparation of the porous carbon electrode involved mixing carbon powder with PTFE binder in ethanol, resulting in a film approximately 150 µm thick. Molecular dynamics (MD) simulations were conducted with two porous carbon electrodes immersed in an ionic liquid electrolyte, specifically [EMIM][BF4]. The simulations employed a constant-potential method to maintain voltage between the electrodes, utilizing a GPU-accelerated implementation for efficiency. The system was equilibrated at 298 K, and analyses were performed on the equilibrium state after 150 ns of polarization. The differential capacitance was derived from the equilibrium charge on the electrode atoms relative to the electrode potential, providing insights into the charging dynamics and ion distribution within the porous structure.
Discussion
The discussion section of the research paper presents a detailed analysis of the atomic model of amorphous porous carbons, specifically focusing on the nanoporous carbon YP-50F. The 3D structure was characterized using small-angle X-ray scattering (SAXS) measurements, revealing two significant scattering peaks at approximately 3.4 and 4.7 nm$^{-1}$. The classical Gaussian random field (cGRF) method initially used was insufficient, leading to the development of the multiple scattering peaks GRF (mGRF) method, which successfully captured the observed peaks and produced a 3D electrode matrix model with a specific surface area (SSA) of 1573 m$^2$g$^{-1}$, closely aligning with experimental values. The hybrid 3D reverse molecular dynamics (H3D-RMD) method further refined the atomic structure, yielding an SSA of 1694 m$^2$g$^{-1}$, validating the model against gas adsorption measurements.
The study also explores the charge storage capabilities of nanoporous carbon electrodes through constant-potential molecular dynamics simulations, demonstrating that ultramicropores exhibit higher gravimetric differential capacitance compared to large micropores. This is attributed to the preferential accumulation of anions in smaller pores, leading to a positive potential of zero charge (PZC) value. The local electrode charge (LEC) analysis indicates that ultramicropores facilitate a more significant charge storage due to their structural characteristics, despite slower ionic conductivity. The charging dynamics were modeled using a multi-scale impedance model (MIM), which effectively predicts the electrochemical impedance spectroscopy (EIS) results, bridging molecular simulations with macroscopic experimental data. This comprehensive modeling framework not only enhances the understanding of charge storage in nanoporous carbon electrodes but also offers methodologies applicable to various fields involving porous media.