A molecular-dynamics study of a model S(N)1 dissociation reaction at the water liquid/vapor interface

The thermodynamics and dynamics of a model S(N)1 reaction: t-BuCl --> t-Bu+ + Cl- is studied at the water liquid/vapor interface using molecular-dynamics computer simulations. The empirical valence bond approach is used to couple two diabatic states, covalent and ionic, in the electronically adiabatic limit. Umbrella sampling calculations are used to calculate the potential of mean force along the reaction coordinate (defined as the t-Bu to Cl distance) in bulk water and in several locations at the interface. We find a significant increase of the dissociation barrier height and of the reaction free energy at the interface relative to the bulk. This is shown to be due to the reduced polarity of the interface. Reactive flux correlation function calculations show significant deviation of the rate constant from the transition-state theory: The transmission coefficients range from 0.49 in the bulk to 0.05 above the Gibbs surface. The low transmission coefficient at the interface despite the lower friction is shown to be due to slow vibrational relaxation.