проблема с RO-DFT в Gaussian

[Vinianel]

Здравствуйте, уважаемые коллеги.
Столкнулась с проблемой: не сходится scf при попытке обсчета катион-радикала с помощью RO-DFT в Gaussian. С чем может быть связана ошибка, и какие есть возможности для ее решения?
Геометрия комплекса оптимизирована с помощью U-DFT.
Из опробованного: увеличение количества SCF циклов, использование Use=L506, понижение симметрии.
Log файл: https://www.dropbox.com/s/cuwefyz5fgrzz ... A.log?dl=0

[Гесс]

А почему важно именно RO? Для удешевления?
Ну есть вот совет http://www.ccl.net/chemistry/resources/ ... index.html
Использование L506 мне не совсем понятно.
Из собственных идей - использование SCF=xqc (а раз известно что нормальная часть всеравно бесполезна - то можно сразу qc (будет очень долго, возможно vshift из совета выше - лучше.
И возможно IOp(5/6) уменьшить до 6 или еще чегото чтоб хоть как то свести.

[Vinianel]

Гесс, спасибо. Хотела пересчитать с RO-DFT частично в образовательных целях, чтобы понять, каков диапазон получаемых энергий (считаю энергию диссоциации связи Со-С), и как такой расчет отличается от U-DFT. Многие советуют использовать RO-DFT вместо U так как в RODFT нет проблемы spin contamination.
Использование L506 взято из обсуждений на форумах, так как SCF=qc не подходит для open-shell систем (в т.ч указано и в приведенной вами сссылке).
Из того, что советуют по ссылке - different (smaller) basis set - пробовала, проблема остается.
RO-HF - тоже пробовала, безрезультатно.
Другой метод - тут проблематично, там советуют гибридные функционалы, а они не подходят для моей системы.
Попробую vshift, и поиграться с IOp. Вернусь, если не получится.

[Гесс]

Это, добавьте P после #

#P ROB3LYP...

Это даст вам более развернутый аутпут и мы хоть увидим как выглядит эта сходимость.

qc работает для U- систем, это чисто для RO- нет но этого я не знал ибо в RO- не играю обычно.

Если проблема от металла - можно попробовать дать лучший базис на нем.
def2-tzvp c def2svp на лиганде или скажем
lanl2dz с cc-pvdz
Ну а касатльно функционалов - а почему гибридные не идут? Ну и из чистого GGA еще много чего есть. PBEPBE, G97D

[Vinianel]

Еще раз спасибо. Добавила.
По поводу гибридных функционалов - многочисленный бенчмаркинг по энергии связи Сo-C в кобаламинах показал, что гибриды недооценивают энергию связи и плохо соотносятся с экспериментальными данными, в то время как чистые GGA коррелируют хорошо. Я взяла PB86 потому что почти все DFT расчеты по кобаламинам за последние 10 лет ведутся с его помощью. Гримме, правда, использовал TPSS.
У меня как-то не возникало мысли, что проблема может быть в функционале, тем более что чистый HF тоже не сошелся. Я попробую систематично пройтись по вашим предложениям, посмотрю, что сработает.

[Vinianel]

Эм. Похоже я уперлась в стену.
#P - крутая вещь, о ней не знала раньше, теперь везде прописываю ее, стало хоть немного понятно что происходит с расчетом.
Вначале SCF сходимость выглядела как-то так:
https://drive.google.com/file/d/0BwOpBb ... sp=sharing

vshift проблему не решил, лучший базис на металле тоже. Поменяла функционал на гибридный - тоже самое. С мертвой точки сдвинуло следующее - использование гибридного функционала + добавление псевдопотенциала на кобальт. Тогда картинка выглядит намного лучше, + начала считаться геометрия:
https://drive.google.com/file/d/0BwOpBb ... sp=sharing

Однако расчет все равно дохнет после нескольких шагов оптимизации.
"Не беда,- подумала я,- перезапущу оптимизацию с предпоследнего шага!" Как бы ни так. Та же картина, разве что теперь расчет дохнет значительно раньше.
https://drive.google.com/file/d/0BwOpBb ... sp=sharing

Лог файл тут: https://www.dropbox.com/s/sh6r5xwoa9w4v ... 6.log?dl=0

[Гесс]

scf=xqc
Должно помочь.

[Vinianel]

Не могу. " No SCF=QC or DM with ROHF"

[Гесс]

а, черт, забыл, обсуждали же...

[eduard1001]

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

[Vinianel]

eduard1001, вот хороший мануал по Гауссиану для новичков:
http://blogs.cimav.edu.mx/daniel.glossm ... ethods.pdf