non-bonded cutoffs
I wanted to add to the non-bond cutoff debate. Even though one typically
uses a non-bonded cutoff that is too small, e.g. 8 angstoms, in many
macromolecular simulations, it is likely better to use ewald summation than
no cutoff in crystal simulations..(what does infinite cutoff really mean
when one has periodic boundary conditions..it is really undefined until you
define your unit cell, which you can make larger and larger)...even in
solution simulations without periodic boundary conditions, with many
formal charges, I would argue that infinite cutoffs might lead to
artifacts because one can't run the simulations long enough to get appropriate
averaging of all these formal charges to give realistic ion atmosphere
effects...the cure may be worse than the disease...i do think longer than
8 angstrom cutoffs are appropriate and recommend dual cutoffs(c.f. Berendsen,
van Gunsteren,Darden,Charifson) which can consider effects out to, e.g. 15
to 20 angstroms at little computational overhead...but in summary, there are
many things more important in reducing simulation artefacts than switching
to an infinite non-bonded cutoff--even if we could, i don't think we should.