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

[BaHb_XeJIbCuHr]

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