orca ccsd diis

[VTur]

СС - это не одноконфигурационный расчет. Базис aug-cc-pVTZ адекватен методу. Дело в чем-то другом.
Я бы выкинул все синглеты и учет триплетов по теории возмущений. Все равно все состояния, за исключением основного, - фиктивные, и нужны только для увеличения гибкости функции.

[kbob]

СС - это не одноконфигурационный расчет. Базис aug-cc-pVTZ адекватен методу. Дело в чем-то другом.
Я бы выкинул все синглеты и учет триплетов по теории возмущений. Все равно все состояния, за исключением основного, - фиктивные, и нужны только для увеличения гибкости функции.

Я имел ввиду что singlereference methods могут быть неприменимы в конкретном случае, поскольку SCF расчет тоже плохо сходится.

[VTur]

А какая разница? Ну, можно написать UCCD, возможно станет лучше, но не уверен. Все равно, то что есть - это надстройка над Х-Ф.

[kbob]

А какая разница? Ну, можно написать UCCD, возможно станет лучше, но не уверен. Все равно, то что есть - это надстройка над Х-Ф.

Возможно, но поскольку формула соединения и геометрия не извесна - гадать не буду.

[Гесс]

геометрически это куб из 8 атомов.

[kbob]

геометрически это куб из 8 атомов.

Кубан, без водородов?

[Гесс]

типа да. С поправкой на тип атомов. В литературе ужо есть углерод и B4N4, а также водород и литий, но они посчитаны как не минимумы. Экспериментально есть красный кислород http://www.membrana.ru/particle/10584, но он не очень кубический.

[VTur]

А это, случаем не радикал? Есть что-нибудь про магнитные свойства?

[ Post made via Android ]

[Гесс]

Триплет выше, по ДФТ и МП2.
Нет ничего по экспериментальным свойствам. Такого кластера вообще быть недолжно.

[VTur]

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

[ Post made via Android ]

[Гесс]

Какие нестыковки? У мну нет нестыковок между расчетными данными, у меня есть расчет который не сходится без видимых на то причин.
У мну есть расчет в CCSD(T)/cc-pvdz в MOLPRO, который сошелся ни на что не жалуясь. Просто он содержит некоторые странные фичи, которых я непонимаю, и по которым я не могу получить консультацию. Если здесь есть люди разбирающиеся или хотя бы работавшие в молпро - я создам отдельную тему. Форум молпро можно не предлагать, я там пару раз постил вопросы, они обычно остаются без ответа.

[VTur]

А нельзя еще в чем-нибудь посчитать?

[ Post made via Android ]

[Гесс]

Программ как бы хватает. Гауссианом вроде как долго, это моя любимая прога, но типа нерекомендованная, я ее оставлю напоследок. Счас согласно теме kbobа попробую поставить и запустить PSI4. Писать письмо создателям CFOUR у меня сегодня вдохновения нет.

[Гесс]

Возможно жизнь наладится и в орке. Я поднял критерий сходимости до 5e-5 а число итераций до 2000, это не помогло, 4 или 5 процессов так и не сошлись, однако задача на этом не упала а пошла дальше, уведомив о том что
WARNING: It seems that the Numerical Gradient ISN'T COMPLETE!
Please check your results VERY carefully!
For possible error information check files CCSD.last*.dpl*
Дальше был второй шаг оптимизации в котором не сошлось 3 процесса, опять предупреждение и уход на третий шаг.
Третий шаг закончился нормально и перешел на 4ый, счас идет четвертый и пока не так и плохо идет. На третьем шаге ни один из критериев Energy change, RMS gradient, MAX gradient, RMS step, MAX step не сошелся, но посмотрим как будет дальше.

PSI4 - когда я впервые поставлю прогу без приключений - я напьюсь. Она компилируется, но десяток тестов падает, в мануале сказано "если упал хоть один тест жалуйтесь разработчикам". пока не хочу продолжать пока не починил тесты. Я правильно понял что форума или чего то такого у них нет и единственный механизм получения помощи - письма мелдким почерком этим людям: http://www.psicode.org/developers.php ?

[sanya1024]

Я в PSI4 разочаровалась после того, как простой обычный бензол с симметрией сошелся в СС быстро, а без симметрии -- повесил 8-ядерную ноду.

А киньте, плз, сюда геометрию -- может, в Природе получится?

[Гесс]

Начиная с третьего шаги пошли все веселее и веселее - на 11 шаге градиенты сходятся уже за менее 40 итераций (ну правда критерий сходимости у них по прежнему 5e-5). За ночь прошло 6 шагов, в то время как за первые двое суток - всего 4. Молекулу конкретно разфигачило в разные стороны, увеличив длины связей к четвертому шагу в полтора раза против нормальных, но счас она пошла сползаться обратно, разумеется никакой симметрии там уже нет, но вцелом исходная геометрия еще с горем поплам угадывается.

[Enchanter]

CCSD(T)/aug-cc-pVDZ : (S) -303.4431711; (T) -303.4038478, без ZPE пока. Чего вы мучаетесь с Molpro и Orca если есть G09?

[Гесс]

CCSD(T)/aug-cc-pVDZ : (S) -303.4431711; (T) -303.4038478, без ZPE пока. Чего вы мучаетесь с Molpro и Orca если есть G09?

Когда я столкнулся с первой проблемой в молпро один умный человек отправил меня в турбомоль и орку и сказал что нормальные люди не делают CCSD оптимизации в гауссиане. Наверное мне стоит уточнить у этого человека что он имел в виду.

[sanya1024]

А что нам скажет мультиплетность 9 в том же варианте расчета?

[Гесс]

Цичас узнаем.