ERROR in dH calculation. Not working. Blowing up.
Some dH blow up, probably b/c of numerical error in (sig/rij). Implemented a cutoff for (sig/rij)