Prediction of clathrate structure type and guest position by molecular mechanics.

The clathrate hydrates occur in various types in which the number, size, and shape of the various cages differ. Usually the clathrate type of a specific guest is predicted by the size and shape of the molecular guest. We have developed a methodology to determine the clathrate type employing molecular mechanics with the MMFF force field employing a strategy to calculate the energy of formation of the clathrate from the sum of the guest/cage energies. The clathrate type with the most negative (most stable) energy of formation would be the type predicted (we mainly focused on type I, type II, or bromine type). This strategy allows for a calculation to predict the clathrate type for any cage guest in a few minutes on a laptop computer. It proved successful in predicting the clathrate structure for 46 out of 47 guest molecules. The molecular mechanics calculations also provide a prediction of the guest position within the cage and clathrate structure. These predictions are generally consistent with the X-ray and neutron diffraction studies. By supplementing the diffraction study with molecular mechanics, we gain a more detailed insight regarding the details of the structure. We have also compared MM calculations to studies of the multiple occupancy of the cages. Finally, we present a density functional calculation that demonstrates that the inside of the clathrates cages have a relatively uniform and low electrostatic potential in comparison with the outside oxygen and hydrogen atoms. This implies that van der Waals forces will usually be dominant in the guest-cage interactions.