[论文解读] A Theoretical Investigation of the Grand- and the Canonical Potential Energy Surface: The Interplay between Electronic and Geometric Response at Electrified Interfaces

Simeon D. Beinlich, Georg Kastlunger+2 more|arXiv (Cornell University)|Jul 19, 2023

Electrochemical Analysis and Applications被引用 5

一句话总结

本文通过分析在恒电势（巨参量系综）和恒电荷（正则系综）势能面上的驻点，为理解电极界面处电子响应与几何响应之间的相互作用提供了严格的理论框架。研究发现，正则系综中的局部极小点可能因特征切换而在巨参量系综中变为鞍点，从而导致负电容——这一现象被证明是小体系尺寸的产物，对基于DFT的电催化体系精确模拟具有重要意义。

ABSTRACT

How does an electrochemical interface respond to changes in the electrode potential? How does the response affect the key properties of the system - energetics, excess charge, capacitance? Essential questions key to ab-initio simulations of electrochemical systems, which we address in this work on the basis of a rigorous mathematical evaluation of the interfacial energetics at constant applied potential. By explicitly taking into account the configurational and electronic degrees of freedom we derive important statements about stationary points in the electronically grand canonical ensemble. We analyze their geometric response to changes in electrode potential and show that it can be mapped identically onto an additional contribution to the system's capacitance. We draw similar conclusions for the constant charge ensemble which equally allows to assess the respective stationary points. Our analysis of the relation between the canonical and grand canonical energetics reveals, however, one key difference between both ensembles. While the constant potential ensemble yields in general positive capacitances at local minima, the capacitance of local minima in the constant charge ensemble might become negative. We trace back this feature to the possibility of character switching of stationary points when switching between the ensembles causing the differences in the response to perturbations. Our systematical analysis not only provides a detailed qualitative and quantitative understanding of the interplay between electronic and configurational degrees of freedom and their contributions to the energetics of electrified interfaces but also highlights the similarities and subtle dissimilarities between the canonical and grand canonical description of the electronic degrees of freedom, which is crucial for a better understanding of theoretical calculations with and without potentiostat.

研究动机与目标

• 严格分析电极电势变化下电极界面的几何与电子响应。
• 阐明正则系综（恒电荷）与巨参量系综（恒电势）在描述界面驻点时的差异。
• 研究正则系综中局部极小点在巨参量系综中不再为极小点的条件。
• 量化几何效应对电极界面总电容的贡献。
• 解决正则模拟中负电容现象的表观悖论，并将其与体系尺寸效应关联。

提出的方法

• 推导巨参量系综势能面（gcPES）作为电极电势 $U$ 和原子坐标 $\vec{r}$ 的函数，其中过剩电荷 $q = -\partial \mathcal{E}/\partial U$。
• 应用施瓦茨定理使混合二阶导数相等，从而关联电荷的构型梯度 $\vec{\nabla}q$ 与力的导数 $\partial \vec{\mathcal{F}}/\partial U$。
• 利用gcPES与cPES的海森矩阵分析驻点的曲率与稳定性，特征值决定其为局部极小点或鞍点。
• 引入判据 $H_{n,n}^* > C_{\text{el}}^* \left( \partial U / \partial r_n \right)^2$ 以判断正则系综中的局部极小点是否在巨参量系综中仍为极小点。
• 分析无限晶胞尺寸极限，表明系综间差异消失，从而将负电容问题解释为有限尺寸效应的产物。
• 比较正则系综与巨参量系综的海森矩阵，量化驻点特征差异及其对电容的影响。

实验结果

研究问题

• RQ1 电极界面在电极电势变化下的几何响应如何与电子响应及总电容相关？
• RQ2 在何种条件下，正则系综（恒电荷）势能面上的局部极小点不再为巨参量系综（恒电势）势能面上的局部极小点？
• RQ3 为何在正则模拟中总电容可能为负？这是物理效应还是人为误差？
• RQ4 电子与几何对电容的贡献如何相互作用，以决定界面结构的稳定性？
• RQ5 驻点特征的系综差异在多大程度上依赖于体系尺寸？这对DFT模拟的可靠性有何影响？

主要发现

• 若几何响应超过电子响应，正则系综中的局部极小点可能在巨参量系综中变为鞍点，导致gcPES中出现鞍点。
• 正则系综中出现负总电容，是因为几何曲率 $H_{n,n}^*$ 小于静电曲率 $C_{\text{el}}^* (\partial U / \partial r_n)^2$，表明在电势控制下，正则系综的局部极小点无法保持稳定。
• 判据 $H_{n,n}^* > C_{\text{el}}^* (\partial U / \partial r_n)^2$ 决定了正则系综中的局部极小点是否在巨参量系综中仍为局部极小点。
• 在无限晶胞尺寸极限下，正则系综与巨参量系综的海森矩阵差异消失，表明特征切换与负电容是恒电荷模拟中较小体系的产物。
• 对电势变化的几何响应可等价映射为系统电容的额外贡献，凸显了结构弛豫在界面静电学中的作用。
• 分析表明，巨参量系综在局部极小点处仅产生正电容，而正则系综可能因特征切换产生负电容，这是理论建模中的关键差异。