Electrostatic component of solvation: comparison of SCRF continuum models

Area of Science:

• Computational Chemistry
• Theoretical Chemistry
• Quantum Mechanics

Background:

• Accurate solvation free energy calculations are crucial for understanding chemical processes in solution.
• Quantum mechanical self-consistent reaction field (SCRF) methods are widely used to model solvent effects.
• Different SCRF models exist, employing various approximations for the solvent reaction field.

Purpose of the Study:

• To systematically compare electrostatic contributions to solvation free energy from three distinct SCRF methods.
• To evaluate the liquid-phase dipole moments as indicators of solute response to the solvent reaction field.
• To assess the impact of solute cavity definitions and charge density representations on SCRF results.

Main Methods:

• Comparison of the generalized Born model (SM5.42R), the Nancy multipolar expansion model, and the MST polarizable continuum model.
• All calculations performed at the Hartree-Fock (HF) level with the 6-31G(d) basis set.
• Analysis of 18 molecules with polar groups across three dielectric permittivity values.

Main Results:

• Significant differences were observed in the calculated electrostatic solvation free energies among the three SCRF methods.
• The liquid-phase dipole moments showed variations, indicating different sensitivities of the solute to the solvent reaction field.
• The choice of solute cavity and charge density representation influenced the results for all tested methods.

Conclusions:

• No single SCRF model universally reproduces electrostatic solvation energies and dipole moments.
• The study underscores the importance of carefully selecting SCRF methods and their parameters for accurate solvation predictions.
• Further development is needed to refine SCRF models for better quantitative agreement with experimental or higher-level theoretical data.