Cleaned up alchemy.py a bit, but running into an issue with OpenMM's numerical integration of long-range dispersion correction for TIP3P with reaction field and switch.