Solution reaction space Hamiltonian based on an electrostatic potential representation of solvent dynamics.

Quantum chemical solvation models usually rely on the equilibrium solvation condition and is thus not immediately applicable to the study of nonequilibrium solvation dynamics, particularly those associated with chemical reactions. Here we address this problem by considering an effective Hamiltonian for solution-phase reactions based on an electrostatic potential (ESP) representation of solvent dynamics. In this approach a general ESP field of solvent is employed as collective solvent coordinate, and an effective Hamiltonian is constructed by treating both solute geometry and solvent ESP as dynamical variables. A harmonic bath is then attached onto the ESP variables in order to account for the stochastic nature of solvent dynamics. As an illustration we apply the above method to the proton transfer of a substituted phenol-amine complex in a polar solvent. The effective Hamiltonian is constructed by means of the reference interaction site model self-consistent field method (i.e., a type of quantum chemical solvation model), and a mixed quantum/classical simulation is performed in the space of solute geometry and solvent ESP. The results suggest that important dynamical features of proton transfer in solution can be captured by the present approach, including spontaneous fluctuations of solvent ESP that drives the proton from reactant to product potential wells.