Расчет энергии водородной связи в Gaussian

[BaHb_XeJIbCuHr]

Доброго времени суток. Необходимо с помощью Gaussian расчитать энергию водородной связи в димере, образованном из тимола и 1,8-цинеола. У кого был опыт таких расчетов в этой проге прошу поделиться. Пробовал сначала построить обе молекулы в HyperChem, сблизить акцепторный и донорный фрагменты на эффективное расстояние, после чего оптимизировать геометрию в MM+ и сохранить файл в формате, который умеет распознавать Gaussian, а именно .mol. Дальше в самом Gaussian разделял атомные группы и проводил расчет по ряду методик, но значения выходили большие, от 1800 до 2200 кДж/моль.

[spa]

Проведите оптимизацию геометрии отдельно молекулы тимола, затем отдельно молекулы 1,8-цинеола, затем проведите оптимизацию геометрии рядом расположенных молекул при этом они как-то могут сблизиться и не обязательно нужным вам образом. Вычтите из энергии комплекса сумму энергий исходных молекул. Получите что-то похожее на энергию водородной связи.