Functional Analytic Derivation and CP2K Implementation of the SCCS Model Based on the Solvent-Aware Interface

Ziwei Chai,Sandra Luber

2024-07-30

Abstract:In the self-consistent continuum solvation (SCCS) approach ($\textit{J. Chem. Phys.}$ 136, 064102 (2012)), the analytical expressions of the local solute-solvent interface functions determine the interface function and dielectric function values at a given real space position based solely on the electron density at that position, completely disregarding the surrounding electron density distribution. Therefore, the low electron density areas inside the solute will be identified by the algorithm as regions where implicit solvent exists, resulting in the emergence of non-physical implicit solvent regions within the solute and even potentially leading to the divergence catastrophe of Kohn-Sham SCF calculations. We present a new and efficient SCCS implementation based on the solvent-aware interface ($\textit{J. Chem. Theory Comput.}$ 15, 3, 1996-2009 (2019)) which addresses this issue by utilizing a solute-solvent interface function based on convolution of electron density in the CP2K software package, which is based on the mixed Gaussian and plane waves (GPW) approach. Starting with the foundational formulas of SCCS, we have rigorously and meticulously derived the contributions of the newly defined electrostatic energy to the Kohn-Sham potential and the analytical forces. This comprehensive derivation utilizes the updated versions of the solute-solvent interface function and the dielectric function, tailored to align with the specifics of the GPW implementation. Our implementation has been tested to successfully eliminate non-physical implicit solvent regions within the solute and achieve good SCF convergence, as demonstrated by test results for both bulk and surface models, namely liquid $H_2O$, titanium dioxide, and platinum.

Chemical Physics

What problem does this paper attempt to address?

The paper primarily aims to address an issue in the self-consistent continuum solvation (SCCS) method, specifically the formation of non-physical implicit solvent regions within the solute. In detail, the solute-solvent interface function in the SCCS method determines the dielectric function value based solely on single-point electron density. This can lead to low electron density regions within the solute being incorrectly identified as implicit solvent regions, causing computational issues such as the divergence of Kohn-Sham self-consistent field (SCF) iterations. To solve this problem, the authors propose a new SCCS implementation based on the solvent-aware interface method. This approach introduces a convolution form of the solute-solvent interface function by utilizing the mixed Gaussian and plane wave (GPW) method in the CP2K software package, considering the surrounding electron density distribution. In this way, non-physical implicit solvent regions within the solute can be avoided, ensuring good convergence of the Kohn-Sham SCF calculations. Additionally, the paper provides a detailed theoretical derivation process, including modifications to the generalized Poisson equation in the SCCS model, the definition of electrostatic energy and its contribution to the Kohn-Sham potential, as well as the contribution to analytical forces. These theoretical improvements enable the SCCS method to more accurately simulate solvent effects and enhance computational efficiency. Finally, the authors demonstrate the effectiveness and convergence of the method by testing different systems (such as liquid water, titanium dioxide surfaces, and platinum surfaces). The results show that the new SCCS implementation successfully eliminates non-physical implicit solvent regions within the solute and achieves good SCF convergence.

Functional Analytic Derivation and CP2K Implementation of the SCCS Model Based on the Solvent-Aware Interface

Self-Consistent Solution Of Kohn-Sham Equation By Real-Space Finite-Difference Method

Periodic Continuum Solvation Model Integrated with First-Principles Calculations for Solid Surfaces

Implicit solvation model for density-functional study of nanocrystal surfaces and reaction pathways

Solvent-aware Interfaces in Continuum Solvation

Robust ΔSCF calculations with direct energy functional minimization methods and STEP for molecules and materials

Modeling solvation effects in real-space and real-time within Density Functional Approaches

Spicing up continuum solvation models with SaLSA: the spherically-averaged liquid susceptibility ansatz

A Seamless Grid-Based Interface for Mean-Field QM/MM Coupled with Efficient Solvation Free Energy Calculations

First-principles molecular dynamics simulations at solid-liquid interfaces with a continuum solvent

Grand Canonical Ensemble Approaches in CP2K for Modeling Electrochemistry at Constant Electrode Potentials

Weighted-density functionals for cavity formation and dispersion energies in continuum solvation models

PCMSolver: An open-source library for solvation modeling

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

Analytic First-Order Derivatives of CASPT2 Combined with the Polarizable Continuum Model

Solvation Free Energies in Subsystem Density Functional Theory

Consistent analytical second derivatives of the Kohn–Sham DFT energy in the framework of the Conductor-like screening Model (COSMO) through Gaussian charge distributions

Adaptive Regularized Self-Consistent Field Iteration with Exact Hessian for Electronic Structure Calculation.

Recent developments and applications of reference interaction site model self-consistent field with constrained spatial electron density (RISM-SCF-cSED): A hybrid model of quantum chemistry and integral equation theory of molecular liquids

Implicit self-consistent electrolyte model in plane-wave density-functional theory

Effect of the Solute Cavity on the Solvation Energy and its Derivatives within the Framework of the Gaussian Charge Scheme