Three-Dimensional RISM Integral Equation Theory for Polarizable Solute Models

Area of Science:

• Computational Chemistry
• Physical Chemistry
• Molecular Modeling

Background:

• Accurate modeling of solvation phenomena is crucial for understanding chemical processes.
• Classical fixed-charge force fields have limitations in describing solute polarizability.
• Polarizable force fields offer improved electrostatic descriptions but integrating them with implicit solvent models is challenging.

Purpose of the Study:

• To develop and validate a novel implicit solvent model for the polarizable AMOEBA force field.
• To adapt the 3D Reference Interaction Site Model (3D RISM) integral equation theory for use with polarizable force fields.
• To assess the performance of the new model in calculating excess chemical potentials and solvation energies.

Main Methods:

• Integration of the polarizable AMOEBA force field with the 3D RISM integral equation theory using an embedding cluster (EC-RISM) strategy.
• Development of the conceptual, physical, and algorithmic basis for the EC-RISM approach.
• Application of the EC-RISM method to benchmark systems for validation.

Main Results:

• The EC-RISM method successfully couples the polarizable AMOEBA force field with an implicit solvent model.
• Calculated free energies in solution show reasonable agreement between the AMOEBA force field and quantum-chemical reference calculations.
• The model allows for the distinct evaluation of energetic and solvation contributions.

Conclusions:

• The developed EC-RISM approach provides a viable method for modeling solvation with polarizable force fields.
• The AMOEBA force field, when parameterized appropriately, demonstrates good performance in reproducing chemical potentials.
• This work suggests potential avenues for systematic improvement of polarizable force fields through better parametrization.