Abstract:The polarizable continuum model (PCM) is a computationally efficient way to incorporate dielectric boundary conditions into electronic structure calculations, via a boundary-element reformulation of Poisson's equation. This transformation is only rigorously valid for an isotropic dielectric medium. To simulate anisotropic solvation, as encountered at an interface or when parts of a system are solvent-exposed while other parts are assumed to be in a nonpolar environment, ad hoc modifications to the PCM formalism have been suggested, in which a dielectric constant is assigned separately to each atomic sphere that contributes to the solute cavity. The accuracy of this "heterogeneous" PCM (HetPCM) method is tested here for the first time, by comparison to results from a generalized Poisson equation solver. The latter is a more expensive and cumbersome approach to incorporate arbitrary dielectric boundary conditions, but one that corresponds to a well-defined scalar permittivity function, epsilon(r). We examine simple model systems for which a model function epsilon(r) can be constructed in a manner that maps reasonably well onto a dielectric constant for each atomic sphere, using a solvent-exposed dielectric constant epsilon_solv = 78 and a variety of smaller values to represent a hydrophobic environment. For nonpolar dielectric constants epsilon_nonp <= 2, differences between the HetPCM and Poisson solvation energies are large compared to the effect of anisotropy on the solvation energy. For epsilon_nonp = 4 and epsilon_nonp = 10, however, HetPCM and anisotropic Poisson solvation energies agree to within 2 kcal/mol in most cases. As a realistic use case, we apply the HetPCM method to predict solvation energies and pKa values for a blue copper protein. The HetPCM method affords pKa values that are more in line with experimental results as compared to either gas-phase calculations or homogeneous (isotropic) PCM results.

Continuous dielectric permittivity II: An Iterative Method for Calculating the Polar Component of the Molecular Solvation Gibbs Energy Under a Smooth Change in the Dielectric Permittivity of a Solution

Production of Charm, Charmonium, and Beauty in 800$\,$GeV Proton-Nucleon Collisions

Polarisable soft solvent models with applications in dissipative particle dynamics

A Molecular Contact Theory for Simulating Polarization: Application to Dielectric Constant Prediction.

Dielectric anisotropy in polar solvents under external fields

On the analogy between the restricted primitive model and capacitor circuits. Part II: A generalized Gibbs-Duhem consistent extension of the Pitzer-Debye-Hückel term with corrections for low and variable relative permittivity

Kawamata-Viehweg vanishing and quint-canonical map of a Complex threefold

A new method to calculate broadband dielectric spectra of solvents from molecular dynamics simulations demonstrated with polarizable force fields

Some challenges of diffused interfaces in implicit-solvent models

Dielectric continuum methods for quantum chemistry

Entropic Origin of Ionic Interactions in Polar Solvents

A fast-Fourier transform method to solve continuum-electrostatics problems with truncated electrostatic interactions: Algorithm and application to ionic solvation and ion–ion interaction

Dielectric Boundary Force in Molecular Solvation with the Poisson-Boltzmann Free Energy: A Shape Derivative Approach.

Incorporating Dipolar Solvents with Variable Density in Poisson-Boltzmann Electrostatics

Multistate, polarizable QM/MM embedding scheme based on the direct reaction field method: Solvatochromic shifts, analytical gradients and optimization of conical intersections in solution

Accounting for electronic polarization in nonpolarizable force fields

Ion Transport in Dipolar Medium I: A Local Dielectric Poisson–Nernst–Planck/Poisson–Boltzmann Model

Equation of Motion for the Solvent Polarization Apparent Charges in the Polarizable Continuum Model: Application to Time-Dependent CI

On the theory of electric double layer with explicit account of a polarizable co-solvent

Polarizable Continuum Models and Green's Function GW Formalism: On the Dynamics of the Solvent Electrons

Testing a heterogeneous polarizable continuum model against exact Poisson boundary conditions