Проблема с файлом контрольной точки в Gaussian 09w

[Андрей Васин]

Добрый день уважаемые коллеги, начал работать с программой Gaussian 09w A-02 SMP ver. 7.0 (для Windows) и при расчёте кристаллической решётки апатит германата столкнулся с тем, что, насколько я понял, программа не может получить данные из файла контрольной точки (с расширением .chk). Если не использовать команды типа "Guess=Read" и "Geom=Allcheck", то расчёт запускается нормально. Файлы input (в виде текстового документа) и .log прилагаются. Заранее благодарю всех отписавшихся.

[VIPer]

попробуйте убрать пробел из имени файла и пути к rwf и chk вообще, либо всё в кавычки поставьте, там где пробел.

[Гесс]

Как называется чекпоинт файл который вы до этогосоздали в C:\g09w\ ?
Переименуйте так чтоб не было пробелов и укажите это имя при считывании в инпуте.

Из общего - расчет периодических систем в гауссиане идея вообще плохая, а систем такого размера - однозначно плохая. Подумайте о кодах нацеленных на периодические расчеты.
Кроме того жуткий зоопарк ключей в инпуте opt=(calcfc,cartesian,maxcycle=200,restart) freq=anharmonic b3lyp/genecp guess=read maxdisk=50GB geom=(connectivity,allcheck) 10f 6d extrabasis integral=grid=600974 iop(4/5=7) scf=(conver=4,xqc,maxcycle=200) вызывает подозрения. Вы либо оооочень хорошо знаете что вы делаете либо совсем наоборот.

[Андрей Васин]

To Гесс. Спасибо, чекпоинт файл назыается также. В общем задача была оптимизировать структуру элементарной ячейки с моделированием ИК-спектра. Restart был добавлен для того,чтобы перезапустить расчёт из контрольной точки в случае если условие сходимости (которое я в этот раз сделал более грубым) не будет выполнено так как при 1-ой попытке расчёта (с более грубой grid) после 11-и дней расчёт не сошёлся. А какие программы вы порекомендуете для расчёта таких систем?

[Гесс]

Как "также"? Гауссиан не любит пробелы в названии или после названия чекпоинт файла. Это должно решаться простой переименовкой файла и редактированием названия в инпуте.
Restart вызывает перезапуск, его ненадо вставлять в инпут который запускается в первый раз. Да и вообще restart редко когда нужен, обычно достаточно считывать данные из чекпоинта.
Расчет мог несойтись по многим разным причинам и сетка лишь одна из них. Вы одновременно меняет scf-алгоритм, критерий сходимости, почему то число scfциклов, сетку интегрирования, делаете чтото через iopы, заказываете очень дорогой ангармонический спектр, расчет в картезиановых координатах вместо внутренних и кучу прочей мелочи. Я не уверен что это заработает как вам хочется даже когда вы решите проблему со сходимостью. Даже если мы скажем что вы будете проходить scf процедуру за день, то оптимизация займет недели, а ангармонические частоты - страшно подумать.
Я пока не занимаюсь периодическими расчетами, но насколько я знаю, к наиболее популярным относятся бесплатный Quantum ESPRESSO, платные VASP и WIEN2k, а также Crystal и Castep условия которых я не помню.
Как минимум большинство этих кодов plane-wave, то есть заточенные к периодическим системам в противовес орбитальным прогам обычно используемых для конечных молекул.
Чисто на всякий случай упомяну что если у вас нет лицензии на Гауссиан то публикация с гауссианом вызовет сложности.

[Андрей Васин]

Также в смысле 63 gapatite.chk. Про пробелы понял, спасибо. Советы учту.

[Гесс]

касательно использованного IOP - Read intermediate SCF results which are on the checkpoint file - может быть не лучшей идеей если scf осцилировал или сходился к чему то нездоровому.

Ну назовите его в обоих местах 63_gapatite.chk, может и заведется.