DOI: https://doi.org/10.1038/s41467-025-58871-7
PMID: https://pubmed.ncbi.nlm.nih.gov/40234425
تاريخ النشر: 2025-04-15
المؤلف: Menglin Sun وآخرون
الموضوع الرئيسي: المحفزات الكهربائية لتحويل الطاقة
نظرة عامة
نقل الإلكترون المرتبط بالبروتون (PCET) هو عملية أساسية في التحفيز الكهربائي، وخاصة في تطبيقات تحويل الطاقة مثل تفاعلات تطور الهيدروجين (HER). تتناول هذه الدراسة التحديات المتعلقة بإدماج تأثيرات الكم النووي (NQEs) في النماذج النظرية لـ PCET تحت ظروف كانونية كبيرة. يقدم المؤلفون إطارًا حسابيًا موحدًا يستخدم أخذ عينات كانونية كبيرة واسعة النطاق وحقول قوة تعلم آلي مصممة خصيصًا للبيئات الكهروكيميائية. تكشف نتائجهم أن تأثيرات الكم النووي تؤثر بشكل كبير على محاكاة PCET، مما يظهر أن الطبيعة الموجية للبروتونات تسمح لها بالتنقل عبر الحواجز الطاقية الكلاسيكية، مما يؤدي إلى تقليل ملحوظ في طاقة التنشيط مقارنةً بالأساليب الكلاسيكية.
تمتد تداعيات هذا البحث إلى ما هو أبعد من HER، مما يشير إلى أن الرؤى حول NQEs يمكن أن تعزز فهمنا لمختلف تفاعلات PCET عبر عمليات تحويل الطاقة المختلفة. تعتبر القدرة على نمذجة ديناميات نقل البروتون بدقة عند واجهة القطب/المذيب أمرًا حاسمًا لتقدم التحفيز الكهربائي وتطوير تقنيات الطاقة الخضراء المبتكرة. يبرز هذا العمل أهمية المحاكاة على النطاق الذري في توضيح الحركيات والديناميكا الحرارية لـ PCET، وبالتالي يساهم في المجال الأوسع للكيمياء الكهربائية والفيزياء الكيميائية.
الطرق
تحدد هذه القسم منهجية خوارزمية مونت كارلو الهجينة الكانونية الكبرى (GC-HMC) المطورة لمحاكاة عمليات نقل الإلكترون المرتبط بالبروتون (PCET) في بيئات القطب/المذيب المعقدة. يتم تعريف دالة تقسيم مجموعة GC، مع إدماج متغيرات مثل العدد الإجمالي للإلكترونات ($N_e$)، والطاقة الكامنة ($E(R, N_e)$)، والجهد الكهروكيميائي ($\mu_e$). تستخدم الخوارزمية نهج مونت كارلو هاملتوني المقيد (HMC) لتعزيز كفاءة أخذ العينات، مما يسمح بإزاحات متعددة الجسيمات مع الحفاظ على حالة جهد ثابتة من خلال خوارزمية متروبوليس. تعتبر هذه الطريقة مفيدة بشكل خاص للأنظمة المفتوحة حيث يمكن أن يتقلب عدد الجسيمات، مما يجعلها أكثر ملاءمة من محاكاة الديناميكا الجزيئية التقليدية (MD).
يتم تفصيل طريقة HMC المقيدة بشكل أكبر، مع التأكيد على عمليتها المكونة من ثلاث خطوات: تطبيق تحويل إحداثي لمراقبة إحداثي التفاعل (RC)، ودمج RC ضمن كثافة الاحتمال، وأخذ عينات درجات الحرية باستخدام مخططات متروبوليس وHMC. بالإضافة إلى ذلك، يناقش القسم إدماج التأثيرات الكمومية من خلال طريقة مونت كارلو التكامل المساري (PIMC)، التي تعالج الطبيعة الكمومية للبروتونات في خطوات PCET. يتم تقديم دالة تقسيم مجموعة GC الكمومية، مع تسليط الضوء على تشابهها مع الصيغة الكلاسيكية، ويتم اشتقاق نسبة القبول للحركات التجريبية في الحالة الكمومية. بشكل عام، تهدف المنهجية إلى توفير إطار شامل لدراسة الخصائص الديناميكية الحرارية والحركية لعمليات PCET مع استيعاب كل من التأثيرات الكلاسيكية والكمومية.
النتائج
في هذا القسم، يقدم المؤلفون نتائج إحصائية من أخذ عيناتهم الكانونية الكبرى (GC)، مع التركيز على التغيرات في الجهد والشحنة خلال خطوات نقل الإلكترون المرتبط بالبروتون (PCET). يحددون علاقة بين جهد التخفيض \( U \) بالنسبة للقطب الهيدروجيني القياسي ودالة العمل \( W(R, N_e) \) للنظام، والتي تتأثر بمعامل الجهد الكهروكيميائي \( \mu_e \). تشير النتائج إلى أن المتوسط الجماعي لدالة العمل \( \langle W(R, N_e) \rangle_{GC} \) يتماشى مع \( -\mu_e \)، مما يثبت صحة نهج أخذ العينات GC الخاص بهم. تتماشى التقلبات في دالة العمل مع النماذج النظرية السابقة، ويظهر العدد الإجمالي للإلكترونات \( N_e \) توزيعًا طبيعيًا، مما يؤكد الالتزام بتوزيع مجموعة GC.
تكشف النتائج أيضًا أنه مع تقدم التفاعل من الحالة الأولية (IS) إلى الحالة النهائية (FS) في كل من آليتي فولمر وهيوروفيسكي، هناك زيادة تدريجية في العدد الإجمالي للإلكترونات، وهو أمر ضروري لتسهيل تفاعلات التخفيض. تدعم هذه الاتجاهات حسابات مفصلة تظهر علاقة مباشرة بين الجهود المطبقة الأكثر سلبية وزيادة أعداد الإلكترونات، مما يبرز استجابة النظام الإلكترونية للظروف المخفضة المتغيرة. يؤكد المؤلفون على أهمية محاكاة الأنظمة المفتوحة الكهروكيميائية بدقة، مشيرين إلى أن المخططات التقليدية للجهد الثابت أو الشحنة قد لا تلتقط الديناميات لعمليات PCET بشكل كافٍ.
المناقشة
في هذه الدراسة، نقدم نموذج DP-$N_e$، وهو طاقة تعلم آلي متقدمة تتضمن العدد الإجمالي للإلكترونات ($N_e$) كمعامل إدخال جديد لتعزيز نمذجة الأنظمة الواجهة في الكيمياء الكهربائية. يسمح هذا النموذج بحساب دالة العمل ويسهل أخذ عينات من مجموعات الكانونية الكبرى (GC) مع أعداد إلكترونات متغيرة. يتم التحقق من دقة نموذج DP-$N_e$ مقابل نتائج نظرية الكثافة الوظيفية (DFT)، مما يظهر أخطاء جذر متوسط المربعات صغيرة (RMSEs) قدرها 0.6 ميلي فولت/ذرة للطاقة و57 ميلي فولت/Å للقوى، مما يشير إلى موثوقيته في وصف سطح الطاقة المحتملة (PES) لخطوات تفاعل تطور الهيدروجين (HER).
تستكشف الأبحاث أيضًا تأثير تأثيرات الكم النووي (NQEs) على طاقات التنشيط لخطوات نقل الإلكترون المرتبط بالبروتون (PCET) في HER. من خلال استخدام طرق مونت كارلو التكامل المساري (PIMC) ضمن إطار GC-HMC، نجد أن إدماج NQEs يقلل بشكل كبير من طاقات التنشيط لخطوات فولمر وهيوروفيسكي بحوالي 0.13 إلكترون فولت و0.09 إلكترون فولت، على التوالي. يشير هذا الانخفاض إلى أن المعالجات الكلاسيكية التقليدية قد تقلل من معدلات التفاعل بمقدار 50-100 مرة عند درجة حرارة الغرفة. تشير نتائجنا إلى أن الانتقال من مسار فولمر-تافل إلى مسار فولمر-هيوروفيسكي يحدث عند جهد مخفض أقل عندما تؤخذ NQEs في الاعتبار، مما يبرز أهمية النفق الكمومي في تعزيز معدلات التفاعل في العمليات التحفيزية الكهربائية. بشكل عام، يؤكد هذا العمل على ضرورة تضمين NQEs في النماذج الحسابية لتحقيق فهم أكثر دقة للآليات الكهروكيميائية.
DOI: https://doi.org/10.1038/s41467-025-58871-7
PMID: https://pubmed.ncbi.nlm.nih.gov/40234425
Publication Date: 2025-04-15
Author(s): Menglin Sun et al.
Primary Topic: Electrocatalysts for Energy Conversion
Overview
Proton-coupled electron transfer (PCET) is a fundamental process in electrocatalysis, particularly in energy conversion applications such as hydrogen evolution reactions (HER). This study addresses the challenges of incorporating nuclear quantum effects (NQEs) into theoretical models of PCET under grand canonical conditions. The authors present a unified computational framework that utilizes extensive grand canonical sampling and a machine learning force field tailored for electrochemical environments. Their findings reveal that NQEs significantly influence PCET simulations, demonstrating that the wave-like nature of protons allows them to tunnel through classical energy barriers, resulting in a notable reduction in activation energy compared to classical approaches.
The implications of this research extend beyond HER, suggesting that insights into NQEs could enhance our understanding of various PCET reactions across different energy conversion processes. The ability to accurately model proton transfer dynamics at the electrode/solvent interface is crucial for advancing electrochemical catalysis and developing innovative green energy technologies. This work highlights the importance of atomic-scale simulations in elucidating the kinetics and thermodynamics of PCET, thereby contributing to the broader field of electrochemistry and chemical physics.
Methods
The section outlines the methodology of the Grand Canonical Hybrid Monte Carlo (GC-HMC) algorithm developed for simulating proton-coupled electron transfer (PCET) processes in complex electrode/solvent environments. The GC ensemble partition function is defined, incorporating variables such as the total electron number ($N_e$), potential energy ($E(R, N_e)$), and electrochemical potential ($\mu_e$). The algorithm employs a constrained Hamiltonian Monte Carlo (HMC) approach to enhance sampling efficiency, allowing for multi-particle displacements while maintaining a constant potential condition through the Metropolis algorithm. This method is particularly advantageous for open systems where particle number can fluctuate, making it more suitable than traditional molecular dynamics (MD) simulations.
The constrained HMC method is further detailed, emphasizing its three-step process: applying a coordinate transformation to monitor the reaction coordinate (RC), integrating the RC within the probability density, and sampling the degrees of freedom using Metropolis and HMC schemes. Additionally, the section discusses the incorporation of quantum effects through the Path Integral Monte Carlo (PIMC) method, which treats the quantum nature of protons in PCET steps. The quantum GC ensemble partition function is presented, highlighting its similarity to the classical formulation, and the acceptance ratio for trial moves in the quantum case is derived. Overall, the methodology aims to provide a comprehensive framework for studying the thermodynamic and kinetic properties of PCET processes while accommodating both classical and quantum mechanical effects.
Results
In this section, the authors present statistical results from their grand canonical (GC) sampling, focusing on the potential and charge variations during proton-coupled electron transfer (PCET) steps. They establish a relationship between the reduction voltage \( U \) relative to the standard hydrogen electrode and the work function \( W(R, N_e) \) of the system, which is influenced by the electrochemical potential parameter \( \mu_e \). The findings indicate that the ensemble average of the work function \( \langle W(R, N_e) \rangle_{GC} \) aligns with \( -\mu_e \), validating their GC sampling approach. The fluctuations in the work function are consistent with previous theoretical models, and the total electron number \( N_e \) exhibits a normal distribution, confirming adherence to the GC ensemble distribution.
The results further reveal that as the reaction progresses from the initial state (IS) to the final state (FS) in both the Volmer and Heyrovsky mechanisms, there is a gradual increase in the total electron number, which is necessary for facilitating reduction reactions. This trend is supported by detailed calculations that show a direct correlation between more negative applied potentials and increased electron numbers, highlighting the system’s electronic response to varying reducing conditions. The authors emphasize the importance of accurately simulating electrochemical open systems, noting that traditional fixed charge or potential schemes may not adequately capture the dynamics of PCET processes.
Discussion
In this study, we present the DP-$N_e$ model, an advanced machine learning potential that incorporates the total number of electrons ($N_e$) as a new input parameter to enhance the modeling of interfacial systems in electrochemistry. This model allows for the calculation of the work function and facilitates the sampling of grand canonical (GC) ensembles with variable electron numbers. The accuracy of the DP-$N_e$ model is validated against density functional theory (DFT) results, demonstrating small root mean square errors (RMSEs) of 0.6 meV/atom for energies and 57 meV/Å for forces, indicating its reliability in describing the potential energy surface (PES) of hydrogen evolution reaction (HER) steps.
The research further explores the impact of nuclear quantum effects (NQEs) on the activation energies of proton-coupled electron transfer (PCET) steps in HER. By employing path integral Monte Carlo (PIMC) methods within the GC-HMC framework, we find that incorporating NQEs significantly lowers the activation energies for the Volmer and Heyrovsky steps by approximately 0.13 eV and 0.09 eV, respectively. This reduction suggests that traditional classical treatments may underestimate reaction rates by 50-100 times at room temperature. Our findings indicate that the transition from the Volmer-Tafel to the Volmer-Heyrovsky pathway occurs at a less reducing potential when NQEs are considered, highlighting the importance of quantum tunneling in enhancing reaction rates in electrocatalytic processes. Overall, this work emphasizes the necessity of including NQEs in computational models to achieve a more accurate understanding of electrochemical mechanisms.