The Effective Molecular Method (EMO) Method

When a molecule is introduced into an aqueous solvent, the water molecules arrange themselves about the solute in a "solvent shell". This solvent shell is not static and the solvent molecules move around, and even enter and exit the solvent shell over time.

The changing environment impacts on the electronic structure of the solute and solvent. However, the Kohn-Sham orbitals obtained for a solute in liquid water are highly delocalised over all the molecules in solution (see the figure on the right).

Until now it has been impossible to isolate information pertinent to a specific solute. We have developed a new method for extracting effective solute MOs, as well as for analyzing local individual solute solvent interactions via a beta parameter.

This takes us a large way forward in the study of solvation effects, and gives us ability to distinguish and investigate the reactivity of solutes in solution. We are now able to examine the molecular orbitals (MOs) of molecules in solution, and to examine how they change over time.

In this graph the light blue line represents the change in energy of the HOMO (highest occupied orbital) of one water molecule over time. This water molecule represents a "solute", and the remaining water molecules the "solvent". The "solvent" HOMOs form a band of energies which lie between the two sets of blue diamond marked lines.

Using Car-Parrinello molecular dynamics coupled with the new method for resolving MOs we are now able to pose and in part answer more penetrating questions, for example; Where does the HOMO of a solution lie, on the solvent or the solute? How do the molecular orbitals of a molecule change on solvation? How does the reactivity of a molecule change in solution (as compared to the gas phase)?
For example we are able to "regenerate" isolated molecular orbitals from the extended delocalised KS orbitals. Above and to the right is shown the molecular orbital(MO) energy diagram for a single isolated molecule of water. Adjacent and directly above are shown the new "effective MOs" which we have determined from the same solution that gave the KS orbitals pictured at the top.