A New Phase-Field Approach to Variational Implicit Solvation of Charged Molecules with the Coulomb-Field Approximation
For researchers in computational chemistry and biophysics, this work provides a more accurate and stable phase-field method for simulating solvation of charged molecules, addressing unphysical features in prior models.
This paper introduces a new phase-field model for solvation of charged molecules that improves accuracy by keeping phase-field values closer to solute and solvent regions, and demonstrates improved localization of forces and robustness of the hyperbolic tangent profile for larger interfacial widths compared to previous models.
We construct a new phase-field model for the solvation of charged molecules with a variational implicit solvent. Our phase-field free-energy functional includes the surface energy, solute-solvent van der Waals dispersion energy, and electrostatic interaction energy that is described by the Coulomb-field approximation, all coupled together self-consistently through a phase field. By introducing a new phase-field term in the description of the solute-solvent van der Waals and electrostatic interactions, we can keep the phase-field values closer to those describing the solute and solvent regions, respectively, making it more accurate in the free-energy estimate. We first prove that our phase-field functionals $Γ$-converge to the corresponding sharp-interface limit. We then develop and implement an efficient and stable numerical method to solve the resulting gradient-flow equation to obtain equilibrium conformations and their associated free energies of the underlying charged molecular system. Our numerical method combines a linear splitting scheme, spectral discretization, and exponential time differencing Runge-Kutta approximations. Applications to the solvation of single ions and a two-plate system demonstrate that our new phase-field implementation improves the previous ones by achieving the localization of the system forces near the solute-solvent interface and maintaining more robustly the desirable hyperbolic tangent profile for even larger interfacial width. This work provides a scheme to resolve the possible unphysical feature of negative values in the phase-field function found in the previous phase-field modeling (cf. H. Sun, et al. J. Chem. Phys., 2015) of charged molecules with the Poisson--Boltzmann equation for the electrostatic interaction.