DOI: https://doi.org/10.1038/s41524-023-01184-4
تاريخ النشر: 2024-01-05
المؤلف: Marko Melander وآخرون
الموضوع الرئيسي: المحفزات الكهربائية لتحويل الطاقة
نظرة عامة
تناقش هذه القسم أهمية الأنظمة الكهروكيميائية في تحويل الطاقة النظيفة والتحديات المرتبطة بفهم كيميائها المعقدة. يبرز دور الجهد الثابت ومحاكاة المجموعة الكبرى (GCE) في استكشاف العمليات الكهروكيميائية كدالة لجهد القطب. ينتقد المؤلفون الحسابات الحالية لـ GCE على مستوى نظرية الكثافة الوظيفية (DFT)، التي تثبت مستوى فيرمي داخل خلية المحاكاة، مشيرين إلى أن هذه الطريقة غير كافية لنمذجة التفاعلات في المجال الخارجي والخلية ثنائية القطب المتحيزة.
للتغلب على هذه القيود، يقترح المؤلفون طريقة الجهد الداخلي الثابت (CIP)، التي تستخدم الجهد الداخلي المحلي للقطب كمعامل حراري لجهد القطب، بدلاً من مستوى فيرمي العالمي. يوضحون من خلال الدراسات العددية والتحليلية أنه بينما تعطي طرق CIP و GCE-DFT لمستوى فيرمي نتائج مكافئة للأقطاب المعدنية والتفاعلات في المجال الداخلي، فإن طريقة CIP أكثر تنوعًا وقابلية للتطبيق على الأنظمة في المجال الخارجي والخلية ثنائية القطب. يضع هذا التقدم طريقة CIP كطريقة قوية وفعالة لـ GCE-DFT لمحاكاة الواجهات الكهروكيميائية من المبادئ الأولى.
مقدمة
تناقش مقدمة ورقة البحث أهمية الأنظمة الكهروكيميائية في تطبيقات متنوعة، بما في ذلك تحويل الطاقة، والتخليق الكهربائي، والاستشعار. إحدى الميزات الرئيسية لهذه الأنظمة هي القدرة على التلاعب بالحرارة الديناميكية والتفاعلات عبر جهد القطب. يتم التحكم في جهد القطب العامل من خلال اتصالات خارجية بمصدر جهد، مما يسمح بالتعديل المباشر للجهد الكهروكيميائي للإلكترونات داخل مادة القطب. يمكن تشبيه هذا الإعداد بإطار ميكانيكا إحصائية، حيث يمثل القطب العامل واتصالاته مراحل متميزة، مع التعبير عن التوازن الكهروكيميائي كـ $\mu_{\text{bath}}^e = \mu_M^e$ عند التوازن الديناميكي الحراري.
يتم وصف الجهد الكهروكيميائي للإلكترونات في مرحلة ما بالمعادلة $\mu_i^e = \mu_{0,i}^e + zF\phi_i$، حيث $\mu_{0,i}^e$ هو الجهد الكيميائي في الحالة القياسية، و $a_i^e$ هو النشاط، و $\phi_i$ هو الجهد الداخلي، و $z$ هو شحنة الإلكترون، و $F$ هو ثابت فاراداي. يؤثر الجهد الداخلي على الجهد المطبق على القطب، والجهد الكهروكيميائي في القطب العامل مرتبط مباشرة بمستوى فيرمي الخاص به. بينما تعتبر تجارب الجهد الثابت شائعة، فإن تحقيق جهد ثابت في المحاكاة الذرية أثبت أنه تحدٍ، مما أدى إلى تطوير طرق محاكاة متنوعة مصنفة إلى أساليب هيكلية إلكترونية (GCE-DFT)، وطرق مجال القوة الكلاسيكية، وطرق الحقل المحدود القابلة للتطبيق على كل من DFT والجهود الكلاسيكية.
طرق
في هذه الدراسة، تم إجراء المحاكاة باستخدام كود GPAW، مع تمثيل دوال الموجة على شبكة فضاء حقيقي تتميز بتباعد قدره 0.18 Å. تم استخدام الوظيفة Perdew-Burke-Ernzerhof (PBE) بشكل متسق عبر جميع الحسابات، التي كانت مستقطبة. لنمذجة المذيب المستمر، تم استخدام نموذج المذيب SCMVD، مع تضمين نصف قطر فان دير فالز قدره 2.0 Å للذهب لتعكس بدقة السعة التفاضلية.
شملت المحاكاة حسابات جهد ثابت تم تنفيذها في كل من مستوى فيرمي القياسي ووضع الجهد الداخلي الثابت المعدل، باستخدام تنفيذ SJM المعدل. كانت النماذج الهيكلية دورية في اتجاهات xy بينما ظلت غير دورية في اتجاه z. تم تطبيق شروط الحدود ديريشليت لضمان أن الجهد الكهروستاتيكي $\phi = 0$ عند حد z للحالات الشحن المحددة، مما يتيح إنشاء قطب مرجعي داخل خلية المحاكاة، كما تم توضيحه في الطرق التكميلية.
نتائج
في هذا القسم، يناقش المؤلفون تنفيذ طريقة الجهد الداخلي الثابت (CIP) ضمن محاكاة نظرية الكثافة الوظيفية المقيدة العامة (GCE-DFT). يتم استخدام الجهد الداخلي كمتغير أساسي للتحكم في الجهد المطبق، مما يتناقض مع الطرق السابقة التي ركزت على مستوى فيرمي. يقدم المؤلفون تعريفًا رسميًا للجهد المطلق، $U_M^{abs} = \Delta M_S \phi – \mu_{M_e} + K$، حيث تمثل $\Delta M_S \phi$ فرق الجهد الكهروستاتيكي بين القطب والمحلول، و $\mu_{M_e}$ هو خاصية مادية جوهرية. يسمح اختيار مقياس المرجع، وخاصة الجهد الداخلي، بالتحكم المباشر في جهد القطب دون الإشارة صراحةً إلى مستوى فيرمي.
تعدل طريقة CIP نموذج الجل المذيب للحفاظ على جهد داخلي ثابت، مما يسهل المحاكاة الدقيقة للأنظمة الكهروكيميائية. من خلال تثبيت الجهد الداخلي $\phi_M$ بالنسبة للجهد الداخلي للمحلول $\phi_S$، الذي يتم تعيينه على الصفر على مسافة من القطب، توفر الطريقة تعريفًا واضحًا لجهد القطب. تضمن هذه الطريقة أن يتم تمثيل الجهد المطبق تجريبيًا بدقة في المحاكاة، مما يسمح باستكشاف عمليات كهروكيميائية متنوعة، بما في ذلك التفاعلات في المجال الداخلي والخارجي، بالإضافة إلى إعدادات الخلايا ثنائية القطب المتحيزة. تؤكد النتائج على أهمية الجهود الداخلية في تحديد سلوك القطب في السياقات الكهروكيميائية.
مناقشة
تتناول قسم المناقشة في ورقة البحث التقدم في محاكاة نظرية الكثافة الوظيفية المقيدة العامة (GCE-DFT)، مع التركيز بشكل خاص على إدخال طريقة DFT للجهد الداخلي الثابت (CIP). اعتمدت الطرق التقليدية لـ GCE-DFT على تثبيت الجهد الكهروكيميائي ($\mu_M^e$)، الذي يتناسب مع مستوى فيرمي، لنمذجة الأنظمة الكهروكيميائية. ومع ذلك، تشير النتائج الأخيرة إلى أن الجهد الكيميائي للحمام ($\mu_{\text{bath}}^e$) هو المعامل الأكثر صلة لتعريف جهد القطب، حيث يتم التحكم فيه تجريبيًا. تعالج طريقة CIP-DFT هذا من خلال السماح بالتحكم المباشر في $\mu_{\text{bath}}^e$، مما يعزز دقة المحاكاة، خاصة في الحالات التي تشمل التفاعلات في المجال الخارجي وإعدادات الخلايا ثنائية القطب حيث تفشل الطرق التقليدية.
تقدم الورقة تحليلًا شاملاً يوضح أن طريقة CIP تعطي نتائج قابلة للمقارنة مع تلك الخاصة بطرق مستوى فيرمي الثابت للتفاعلات في المجال الداخلي، بينما تحسن بشكل كبير من معالجة التفاعلات في المجال الخارجي. تسهل طريقة CIP-DFT أيضًا نمذجة الأنظمة ثنائية القطب المتحيزة دون الحاجة إلى مدارات موضعية أو طرق دالة غرين، والتي تكون عادةً صعبة في محاكاة DFT. يفتح هذا التنوع آفاقًا جديدة لمحاكاة العمليات الكهروكيميائية وديناميات نقل الشحن غير المتوازنة، مما يجعل CIP-DFT أداة واعدة للدراسات الحسابية المستقبلية في الكيمياء الكهربية. بشكل عام، تؤكد النتائج على أهمية تعريف جهود الأقطاب بدقة في النماذج الحسابية لتتوافق مع الظروف التجريبية، مما يعزز موثوقية التنبؤات النظرية في الأنظمة الكهروكيميائية.
DOI: https://doi.org/10.1038/s41524-023-01184-4
Publication Date: 2024-01-05
Author(s): Marko Melander et al.
Primary Topic: Electrocatalysts for Energy Conversion
Overview
The section discusses the significance of electrochemical systems in clean energy conversion and the challenges associated with understanding their complex chemistry. It highlights the role of constant potential and grand canonical ensemble (GCE) simulations in exploring electrochemical processes as a function of electrode potential. The authors critique the current GCE calculations at the density functional theory (DFT) level, which fix the Fermi level within the simulation cell, noting that this approach is inadequate for modeling outer sphere reactions and biased two-electrode cells.
To overcome this limitation, the authors propose a constant inner potential (CIP) method, which utilizes the local electrode inner potential as the thermodynamic parameter for electrode potential, rather than the global Fermi level. They demonstrate through numerical and analytical studies that while the CIP and Fermi level GCE-DFT methods yield equivalent results for metallic electrodes and inner-sphere reactions, the CIP method is more versatile and applicable to outer-sphere and two-electrode systems. This advancement positions the CIP approach as a robust and efficient GCE-DFT method for simulating electrochemical interfaces from first principles.
Introduction
The introduction of the research paper discusses the significance of electrochemical systems in various applications, including energy conversion, electrosynthesis, and sensing. A key feature of these systems is the ability to manipulate reaction thermodynamics and kinetics via electrode potential. The working electrode’s potential is controlled through external connections to a voltage source, allowing for direct modulation of the electrochemical potential of electrons within the electrode material. This setup can be likened to a statistical mechanics framework, where the working electrode and its connections represent distinct phases, with the electrochemical balance expressed as $\mu_{\text{bath}}^e = \mu_M^e$ at thermodynamic equilibrium.
The electrochemical potential of electrons in a phase is described by the equation $\mu_i^e = \mu_{0,i}^e + zF\phi_i$, where $\mu_{0,i}^e$ is the standard state chemical potential, $a_i^e$ is the activity, $\phi_i$ is the inner potential, $z$ is the electron charge, and $F$ is the Faraday constant. The inner potential influences the applied electrode potential, and the electrochemical potential in the working electrode is directly related to its Fermi level. While constant potential experiments are common, achieving constant potential in atomistic simulations has proven challenging, leading to the development of various simulation methods categorized into electronic structure approaches (GCE-DFT), classical force field methods, and finite-field methods applicable to both DFT and classical potentials.
Methods
In this study, simulations were conducted using the GPAW code, with wave functions represented on a real-space grid featuring a spacing of 0.18 Å. The Perdew-Burke-Ernzerhof (PBE) functional was consistently employed across all calculations, which were spin-polarized. To model the continuum solvent, the SCMVD dielectric continuum model was utilized, incorporating a van der Waals radius of 2.0 Å for gold to accurately reflect the differential capacitance.
The simulations included constant potential calculations executed in both standard Fermi level and modified constant inner potential modes, utilizing a revised SJM implementation. The structural models were periodic in the xy-directions while remaining non-periodic in the z-direction. Dirichlet boundary conditions were applied to ensure that the electrostatic potential $\phi = 0$ at the z-boundary for the specified charge states, thereby establishing a reference electrode within the simulation cell, as elaborated in the Supplementary Methods.
Results
In this section, the authors discuss the implementation of the Constant Inner Potential (CIP) approach within Generalized Constrained Electrochemical Density Functional Theory (GCE-DFT) simulations. The inner potential is utilized as a fundamental variable to control the applied potential, contrasting with previous methods that focused on the Fermi level. The authors present a formal definition of the absolute potential, $U_M^{abs} = \Delta M_S \phi – \mu_{M_e} + K$, where $\Delta M_S \phi$ represents the electrostatic potential difference between the electrode and the solution, and $\mu_{M_e}$ is an intrinsic material property. The choice of reference scale, particularly the inner potential, allows for direct control of the electrode potential without explicitly referencing the Fermi level.
The CIP method modifies the solvent jellium model to maintain a constant inner potential, facilitating accurate simulations of electrochemical systems. By fixing the inner potential $\phi_M$ relative to the solution’s inner potential $\phi_S$, which is set to zero at a distance from the electrode, the method provides a clear definition of the electrode potential. This approach ensures that the experimentally applied potential is represented accurately in simulations, allowing for the exploration of various electrochemical processes, including inner-sphere and outer-sphere reactions, as well as biased two-electrode setups. The findings underscore the significance of inner potentials in determining electrode behavior in electrochemical contexts.
Discussion
The discussion section of the research paper elaborates on the advancements in Grand Canonical Ensemble Density Functional Theory (GCE-DFT) simulations, particularly focusing on the introduction of the Constant Inner Potential (CIP) DFT method. Traditional GCE-DFT approaches have relied on fixing the electrochemical potential ($\mu_M^e$), which is proportional to the Fermi level, to model electrochemical systems. However, recent findings suggest that the bath chemical potential ($\mu_{\text{bath}}^e$) is the more relevant parameter for defining electrode potential, as it is experimentally controlled. The CIP-DFT method addresses this by allowing direct control over $\mu_{\text{bath}}^e$, thereby enhancing the accuracy of simulations, especially in cases involving outer-sphere reactions and two-electrode setups where traditional methods fall short.
The paper presents a comprehensive analysis demonstrating that the CIP approach yields results comparable to those of fixed Fermi level methods for inner-sphere reactions while significantly improving the treatment of outer-sphere reactions. The CIP-DFT method also facilitates the modeling of biased two-electrode systems without the need for localized orbitals or Green’s function methods, which are typically challenging in DFT simulations. This versatility opens new avenues for simulating electrochemical processes and non-equilibrium charge transfer dynamics, making CIP-DFT a promising tool for future computational studies in electrochemistry. Overall, the findings emphasize the importance of accurately defining electrode potentials in computational models to align with experimental conditions, thereby enhancing the reliability of theoretical predictions in electrochemical systems.