Molecular simulation of the hydration of ethene to ethanol using ab initio potentials and free energy curves

Molecular dynamics simulations of aqueous solutions at infinite dilution of the reaction of water with ethene: H2O + CH2CH2 --> CH3CH2OH were performed using Lennard-Jones 12-6-1 potentials to describe the solute-solvent interactions, and TIP3P to describe the water-water interactions. The Morokuma decomposition scheme of ab initio interaction energies at the SCF level and the dispersion component at the MP2 level were used to reproduce the molecular parameters of the solute-water interaction potentials. The results show that the functions that use the EX-PL-DIS-ES interaction model to describe the solvation of the reactant, transition state, and product systems lead to good values of the reaction (Delta G) and acceptable values of the activation (Delta G#) free energy as compared with those from using AMBER-derived parameters, using the available theoretical and experimental data as referents.