[論文レビュー] 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.

研究の動機と目的

• 電極の電位変化に対する電気的に起因する界面の幾何的および電子的応答を厳密に分析すること。
• カノニカル（一定電荷）およびグランドカノニカル（一定電位）の集合が、界面の停留点を記述する際の違いを明確にすること。
• カノニカル集合における局所的最小値が、グランドカノニカル集合においても最小値でなくなる条件を調査すること。
• 幾何的効果が電気的に起因する界面の総静電容量にどのように寄与するかを定量化すること。
• カノニカルシミュレーションにおける負の静電容量の表面的パラドックスを解消し、それが系のサイズ効果に起因することを示すこと。

提案手法

• 電極の電位 $U$ および原子座標 $\vec{r}$ の関数として、グランドカノニカルのポテンシャルエネルギーサーフェス（gcPES）を導出する。ここで過剰電荷は $q = -\partial \mathcal{E}/\partial U$ で与えられる。
• シュワルツの定理を適用し、混合2階微分を等置することで、電荷の配置勾配 $\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$ が、カノニカルの局所的最小値がグランドカノニカル集合でも最小値のままであるかどうかを決定する。
• 無限大のセルサイズ極限において、カノニカルとグランドカノニカルのヘッセ行列の差が消えるため、特性の切り替えや負の静電容量は、一定電荷シミュレーションにおける小規模系のアーティファクトであると示唆される。
• 電位変化に対する幾何的応答は、システムの静電容量への追加的寄与として同一にマッピング可能であり、構造的リラクゼーションが界面静電的性質に与える役割を強調している。
• 分析により、グランドカノニカル集合では局所的最小値において常に正の静電容量が得られるが、カノニカル集合では特性の切り替えにより負の静電容量が得られることがあることが明らかになった。これは理論的モデリングにおいて重要な差異である。