Не оптимизируется геометрия

[AlexKsen]

Добрый день, уважаемые форумчане. Подскажите, почему может не оптимизироваться геометрия следующих молекул, out:https://cloud.mail.ru/public/CgM4/RzpppK9j6, https://cloud.mail.ru/public/FmYs/FafkJ2c4y
Причем, с заместителями до C10H21 все отлично соптимизировалось. Не могу понять в чем проблема

[Гесс]

Смотрим в конец файла:

***** FAILURE TO LOCATE STATIONARY POINT, SCF HAS NOT CONVERGED

Смотрим в участок по SCF (проще всего поиском найти слово SCF ибо в этом аутпуте море слабонужной инфы которую задалбываешься проматывать)

99 98 0 -1860.0039610308 -0.0000001531 0.000210349 0.000001589 2790264206415527137142
100 99 0 -1860.0039609945 0.0000000363 0.000920264 0.000001315 2822323171315515699107

SCF IS UNCONVERGED, TOO MANY ITERATIONS

Если построить график изменения энергии от шага (или внимательнее посмотреть на E CHANGE колонку) то видно что начиная примерно с 24ого шага энергию колбасит примерно в интервале 0.4 калорий (не ккал). С другой стороны это не осциляции (незнаю насколько они характерны для гамесса). Вариантов много, имхо наиболее логичный сначала попробовать просто дать больше шагов на scf. Если 400 шагов не решат проблему надо менять алгоритм сведения. Какие алгоритмы есть у гамесса - я не в курсе.

Отвлеченный вопрос - а выбор конформера чем то мотивирован?

[AlexKsen]

Конформера с точки зрения расположения заместителя или длины углеродной цепи?

[Гесс]

С точки зрения изгиба длинных цепочек в пространстве. Для чистых алканов помнится цепочки до 17 атомов предпочтительны в относительно линейной форме а дальше складываются пополам, но тут все может быть по другому. Имхо тут конформационный поиск невреден.