A script to compute electrostatic and hydrogen bond contributions to the lattice energy

This script uses Forcite together with the Dreiding forcefield to compute the electrostatic and hydrogen bond contributions to the lattice energy. This can be a useful analysis step after a polymorph prediction sequence. The script adds two new columns to a study table named "Electrostatic energy" and "Hydrogen bond energy" which are calculated as

E_latt = E_crystal/(#molecules per cell) - E_molecule

Geometry optimizations on crystal and molecular structure are performed to obtain forcefield consistent energies.

Minimum version: Materials Studio 4.2 (should work on 4.1, not tested though)
Modules: MS Forcite

Best regards, Max Petersen