deepmodeling / abacus-develop

An electronic structure package based on either plane wave basis or numerical atomic orbitals.
http://abacus.ustc.edu.cn
GNU Lesser General Public License v3.0
169 stars 129 forks source link

Feature: require an automatic workflow of electron-number adjustment for a constant work function #3579

Open sunml99 opened 8 months ago

sunml99 commented 8 months ago

Background

Realistic electrochemical reactions often occur under a constant applied potential. It is therefore crucial to achieve a constant potential condition for DFT-based electrochemcial surface reaction simulations. DFT calculations in ABACUS are currently based on a constant number of electrons instead of a constant potential (i.e., work function).

Describe the solution you'd like

Since we are already able to adjust the total number of electrons in ABACUS slab models to approximate an effective electrical double-layer condition, we would like to request the ABACUS team to construct an automatic workflow of electron-number adjustment for converging the slab model's work function to a certain value. We would also prefer this constant-potential workflow/module to combine with CI-NEB calculations, which can help people do transition-state optimization under a constant potential condition.

Task list only for developers

Notice Possible Changes of Behavior (Reminder only for developers)

No response

Notice any changes of core modules (Reminder only for developers)

No response

Notice Possible Changes of Core Modules (Reminder only for developers)

No response

Additional Context

No response

Task list for Issue attackers (only for developers)

WHUweiqingzhou commented 8 months ago

@sunml99, Thanks for your idea! Indeed, this function is quite important for surface application. We will consider a practical schedule about it.

WHUweiqingzhou commented 8 months ago

Maybe we can implement this method in ABACUS referred by J. Chem. Theory Comput. 2023, 19, 5168−5175.