QM расчет потенциала для МД

[КВГ]

Здравствуйте!
1) Не подскажет ли кто-нибудь как рассчитать квантовохимически (есть PC Gamess) потенциал взаимодействия между ионами (напр., Ta(5+), Cs+ или F-) и поверхностью переходного металла (для конкретности - Ta, Mo, W)? Этот потенциал нужен для последующего использования в какой-либо программе молекулярной динамики (пока не знаю, в какой именно). Поверхность металла моделируем кластером подходящих размеров (2-3 атомарных слоя).
Можно ли как-то учесть кулоновское взаимодействие с зарядом-изображением (который в данном случае допустимо рассматривать как зеркальное изображение заряда "первичного" иона, расположенное в вакууме симметрично относительно плоской поверхности металла)?

2) В методическом отношении имеются ли принципиальные отличия от первого случая при расчете потенциала взаимодействия двух ионов между собой (т.е. никаких поверхностей нет)?

3) Работают ли МД программы с табулированным потенциалом или надо непременно аппроксимировать его аналитич. ф-цией? Какие free программы предпочтительнее использовать, ежели среди оных имеются работающие с табулированными потенциалами?
Спасибо!

[Darth Vasya]

GROMACS может работать с табулированными потенциалами. CHARMM, Насколько знаю, тоже.

[КВГ]

GROMACS может работать с табулированными потенциалами.

Спасибо! Ссылка отличная - этот Gromacs еще и под винду оказывается можно заточить. Осталось разжиться таблицами потенциалов...