[논문 리뷰] 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명|arXiv (Cornell University)|2023. 07. 19.

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$를 포함한다.
• Schwarz의 정리를 적용하여 혼합 두 번째 도함수를 동치로 연결함으로써, 전하의 구성 기울기 $\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$은 캐논니컬 국소 최솟값이 그랜드 캐논니컬 앙상블에서도 국소 최솟값으로 유지되는지를 결정한다.
• 무한한 세포 크기 근사에서 캐논니컬과 그랜드 캐논니컬 헤시안 간의 차이가 사라지며, 이는 정적 점 특성 전환과 음의 용량이 고정된 전하 시뮬레이션에서 작은 시스템에서 기인한 결과임을 시사한다.
• 전위 변화에 대한 기하 반응은 동일하게 시스템의 용량에 추가 기여로 매핑될 수 있으며, 이는 구조적 완화가 인터페이스 전기적 성질에 미치는 역할을 강조한다.
• 분석 결과, 그랜드 캐논니컬 앙상블에서는 국소 최솟값에서 항상 양의 용량을 얻지만, 캐논니컬 앙상블은 특성 전환로 인해 음의 용량을 유도할 수 있으며, 이는 이론적 모델링에서의 핵심적 차이점이다.