Abstract: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 $\Gamma$-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.

A NUMERICAL METHOD FOR THE PHASE-FIELD MODEL OF THE VESICLE DYNAMICS

A Hybrid Numerical Method For The Phase-Field Model Of Fluid Vesicles In Three-Dimensional Space

Numerical approximations of the Navier–Stokes equation coupled with volume-conserved multi-phase-field vesicles system: Fully-decoupled, linear, unconditionally energy stable and second-order time-accurate numerical scheme

Totally Decoupled Implicit–explicit Linear Scheme with Corrected Energy Dissipation Law for the Phase-Field Fluid Vesicle Model

Adaptive Finite Element Method for a Phase Field Bending Elasticity Model of Vesicle Membrane Deformations

Decoupled Energy Stable Schemes for Phase-Field Vesicle Membrane Model.

An Energy Stable $C^0$ Finite Element Scheme for A Phase-Field Model of Vesicle Motion and Deformation

A Phase Field Method for Simulating Morphological Evolution of Vesicles in Electric Fields

Numerical Simulation of Cell Migration and Deformation with Flow Using a Coupled Lattice Boltzmann and Phase-Field Method

A New Phase-Field Approach to Variational Implicit Solvation of Charged Molecules with the Coulomb-Field Approximation

An Efficient Scheme for A Phase Field Model for the Moving Contact Line Problem with Variable Density and Viscosity

A phase-field model of elastic and viscoelastic surfaces in fluids

Discontinuous Finite Volume Element Method for a Coupled Navier-Stokes-Cahn-Hilliard Phase Field Model

Topology‐preserving phase‐field modeling of elastic bending energy

A Phase-Field Fluid Modeling and Computation with Interfacial Profile Correction Term.

Phase Field Modeling and Numerical Algorithm for Two-Phase Dielectric Fluid Flows

Modelling and Simulations of Multi-component Lipid Membranes and Open Membranes via Diffusive Interface Approaches

A Phase Field Method for the Numerical Simulation of Rigid Particulate in Two-Phase Flows

A fictitious domain method with a hybrid cell model for simulating motion of cells in fluid flow

High-volume fraction simulations of two-dimensional vesicle suspensions

A hybrid phase field method for fluid-structure interactions in viscous fluids