On the equivalence of grand-canonical and capacitor-based models in first-principle electrochemistry

Georg Kastlunger,Sudarshan Vijay,Andrew Peterson,Shubham Sharma,Xi Chen

DOI: https://doi.org/10.26434/chemrxiv-2023-1t89n

2023-12-11

Abstract:First principles-based computational and theoretical methods are constantly evolving trying to overcome the many obstacles towards a comprehensive understanding of electrochemical processes on an atomistic level. One of the major challenges has been the determination of reaction energetics under a constant applied potential. Here, a theoretical framework was proposed applying standard electronic structure methods and extrapolating to the infinite-cell size limit where reactions do not alter the potential. More recently, grand-canonical modifications to electronic structure methods which hold the potential constant by varying the number of electrons in a finite simulation cell have gotten increasingly popular. In this perspective, we show that these two schemes are thermodynamically equivalent. Further, we link these methods to ones based on capacitive models of the interface, in the limit that the capacitance of the charging components (whether continuum or atomistic) are equal and invariant along the reaction pathway. We further benchmark the three approaches with an example of alkali cation adsorption on Pt(111) proving that all three approaches converge in the cases of Li, Na and K. For Cs, however, strong deviation from the ideal conditions leads to a spread in the respective results. We discuss the latter by highlighting the cases of broken equivalence and assumptions among the approaches.

Chemistry

What problem does this paper attempt to address?

This paper attempts to solve the problem of how to accurately calculate reaction thermodynamics under the condition of constant potential in the electrochemical process. Specifically, the authors explored the thermodynamic equivalence between two main methods - the standard electronic structure method (by derivation to the infinite cell - size limit) and the generalized canonical ensemble method (by changing the number of electrons in the finite simulation cell to keep the potential constant). In addition, they also related these methods to the method based on the capacitor model, especially when the capacitance of the charging component (whether in the continuum or at the atomic level) is equal and unchanged along the reaction path.

### Main research questions:
1. **Thermodynamic equivalence**: Prove that the standard electronic structure method (derived to the infinite cell - size limit) and the generalized canonical ensemble method are thermodynamically equivalent.
2. **Connection of capacitor model**: Relate the above two methods to the method based on the capacitor model and discuss their applicable conditions.
3. **Method verification**: Verify the consistency of these three methods through the adsorption experiments of lithium, sodium, potassium and cesium on the Pt(111) surface, and discuss the special case of cesium.

### Key points:
- **Generalized canonical ensemble method**: Allows the number of electrons in the simulation to change continuously and keeps the potential constant by applying an external charge.
- **Standard electronic structure method**: By derivation to the infinite cell - size limit, the reaction does not change the potential.
- **Capacitor model**: Assume that the electrode - electrolyte interface can be regarded as a parallel - plate capacitor, and correct the free energy in different states by adjusting the capacitance.

### Experimental verification:
- **Lithium, sodium, potassium**: In the adsorption experiments of these three ions, the results of the three methods are consistent, indicating that these methods are equivalent under ideal conditions.
- **Cesium**: Due to the special properties of cesium (such as depolarization), the results deviate, indicating that in some cases, the equivalence of these methods may not hold.

### Conclusions:
- Through theoretical derivation and experimental verification, the paper proves that under certain conditions, the standard electronic structure method, the generalized canonical ensemble method and the capacitor model method are equivalent in calculating the thermodynamics of electrochemical reactions.
- For special cases such as cesium, further research on the underlying physical mechanisms is required to understand why the deviation occurs.

Through these studies, the paper provides a theoretical basis and practical guidance for the simulation of electrochemical processes, which is helpful for more accurately understanding and predicting the thermodynamic behavior of electrochemical reactions.