A unified electrostatic and cavitation model for first-principles molecular dynamics in solution

TL;DR

This paper introduces a combined electrostatic and cavitation model for first-principles molecular dynamics in solution, enabling efficient and accurate simulations of solvated systems with applications to molecular dimerization.

Contribution

It presents a novel unified model integrating electrostatic and cavitation effects for first-principles MD, improving efficiency and accuracy in solvated system simulations.

Findings

Cavitation energy matches complex algorithms with fewer parameters.

Model achieves accuracy comparable to established quantum-chemistry methods.

Application to tetracyanoethylene dimers reveals structural and dynamical insights.

Abstract

The electrostatic continuum solvent model developed by Fattebert and Gygi is combined with a first-principles formulation of the cavitation energy based on a natural quantum-mechanical definition for the surface of a solute. Despite its simplicity, the cavitation contribution calculated by this approach is found to be in remarkable agreement with that obtained by more complex algorithms relying on a large set of parameters. Our model allows for very efficient Car-Parrinello simulations of finite or extended systems in solution, and demonstrates a level of accuracy as good as that of established quantum-chemistry continuum solvent methods. We apply this approach to the study of tetracyanoethylene dimers in dichloromethane, providing valuable structural and dynamical insights on the dimerization phenomenon.