Gaussian-Based Periodic Grand Canonical Density Functional Theory with Implicit Solvation for Computational Electrochemistry
arXiv:2508.15705 · doi:10.1021/acs.jctc.5c01403
Abstract
We present a numerical method for grand canonical density functional theory (DFT) tailored to solid-state systems, employing Gaussian-type orbitals as the primary basis. Our approach directly minimizes the grand canonical free energy using the density matrix as the sole variational parameter, while self-consistently updating the electron number between self-consistent field iterations. To enable realistic electrochemical modeling, we integrate this approach with implicit solvation models. Our solvation scheme introduces less than 50% overhead relative to gas-phase calculations. Compared to existing plane wave-based implementations, our method shows improved robustness in grand canonical simulations. We validate the approach by modeling corrosion at silver surfaces, finding excellent agreement with previous studies. Our method is implemented in the quantum chemistry software Q-Chem. This work lays the groundwork for future wavefunction-based simulations beyond DFT under electrochemical operando conditions.
References in corpus (23)
- Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
- CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
- Revised self-consistent continuum solvation in electronic-structure calculations
- Grand canonical electronic density-functional theory: algorithms and applications to electrochemistry
- JDFTx: software for joint density-functional theory
- Joint Density-Functional Theory of the Electrode-Electrolyte Interface: Application to Fixed Electrode Potentials, Interfacial Capacitances, and Potentials of Zero Charge
- The charge-asymmetric nonlocally-determined local-electric (CANDLE) solvation model
- The importance of nonlinear fluid response in joint density-functional theory studies of battery systems
- A unified electrostatic and cavitation model for first-principles molecular dynamics in solution
- A generalized Poisson and Poisson-Boltzmann solver for electrostatic environments
- Applying the Coupled-Cluster Ansatz to Solids and Surfaces in the Thermodynamic Limit
- Evaluating continuum solvation models for the electrode-electrolyte interface: challenges and strategies for improvement
- Function-Space Based Solution Scheme for the Size-Modified Poisson-Boltzmann Equation in Full-Potential DFT
- Interaction Picture Density Matrix Quantum Monte Carlo
- A Phaseless Auxiliary-Field Quantum Monte Carlo Perspective on the Uniform Electron Gas at Finite Temperatures: Issues, Observations, and Benchmark Study
- Toward Accurate Adsorption Energetics on Clay Surfaces
- The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients
- Adsorption and Vibrational Spectroscopy of CO on the Surface of MgO from Periodic Local Coupled-Cluster Theory
- Finite Temperature Auxiliary Field Quantum Monte Carlo in the Canonical Ensemble
- Weighted-density functionals for cavity formation and dispersion energies in continuum solvation models
- Finite-temperature many-body perturbation theory in the grand canonical ensemble
- Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems
- Even Faster Exact Exchange for Solids via Tensor Hypercontraction