Спиновая плотность

[Neitrino]

Подскажите какие keywords позволяют получать числовые значения спиновой плотности в Gaussian03W.

[alxyppv]

Спиновую плотность где? Если в ядрах, то это выдается автоматом (и в единицах плотности, и пересчитанные в константы СТВ) - если в молекуле меньше 100 атомов. Если в молекуле больше 100 атомов, то надо добавить "IOp(6/82=1)". Ну и конечно не забываем ставить #p в начале командной строки.

[Neitrino]

Спиновая плотность электронов для радикалов.

[alxyppv]

Еще раз повторю вопрос. В каких точках Вас интересует значение спиновой плотности? Или Вас интересует не значение спиновой плотности, а спиновая "заселенность" - т.е. интеграл спиновой плотности по всему атому?

[VTur]

В какой модели?

[Neitrino]

Спиновая "заселенность" - т.е. интеграл спиновой плотности по всему атому.

[VTur]

Попробуйте AIM(tight,BO). Но в ГАУССИАНЕ 03 это, то работает, то не работает. Работает в ГАУССИАНЕ 98.

[Neitrino]

Спасибо, попробую. Не могли бы объяснить почему при расчетах РМ3 методом полная энергия положительна. Например, для ацетилена 0,0807 a.u.

[amge]

Спасибо, попробую. Не могли бы объяснить почему при расчетах РМ3 методом полная энергия положительна. Например, для ацетилена 0,0807 a.u.

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

[Neitrino]

А как тогда сравнить полную энергия полуэмперического метода с неэмперическими? От чего нужно отнимать энергию полученную полуэмперическим методом, чтобы получить отрицательные значения?

[amge]

А как тогда сравнить полную энергия полуэмперического метода с неэмперическими?

А это и не нужно (бессмысленно) делать. Сравнивать можно только энергии соединений, посчитанных одним и тем же методом, а в случае неэмпирики - только изомеров, посчитанных одним и тем же методом в одинаковом базисе.

[pokoleg]

а в случае неэмпирики - только изомеров, посчитанных одним и тем же методом в одинаковом базисе.

Золотые слова!

[Neitrino]

Значения полной энергии расчитанные полуэмперикой такими как АМ1 и РМ3 для одной и той же молекулы могут различаться. Полная энергия рассчитанная РМ3 методом более точная или близка к экспериментальным значениям. Несомненно что неэмперика дает более точные значения. Но как быть со спиновой плотностью? А еще не поскажете как устранить ошибку связанную с l9999?

[VTur]

Значения полной энергии расчитанные полуэмперикой такими как АМ1 и РМ3 для одной и той же молекулы могут различаться.

Несомненно, это же разные схемы параметризации.

Полная энергия рассчитанная РМ3 методом более точная или близка к экспериментальным значениям.

Не очень понятно, о каких экспериментальных значениях идёт речь.

Несомненно что неэмперика дает более точные значения. Но как быть со спиновой плотностью? А еще не поскажете как устранить ошибку связанную с l9999?

Можно немного поподробнее?

[Neitrino]

Error termination request processed by link 9999.
Error termination via Lnk1e in C:\G03W\l9999.exe at Thu May 21 10:46:16 2009.
Job cpu time: 0 days 0 hours 5 minutes 25.0 seconds.
File lengths (MBytes): RWF= 11 Int= 0 D2E= 0 Chk= 4 Scr= 1

[VTur]

Это не то, ошибка где-то выше по файлу. Этот оверлей просто завершает выполнение и архивирует данные.

[Neitrino]

Не подскажете в какой проге метода функционала плотности можно посчитать ширину запрещенной зоны?

[Darth Vasya]

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

[Neitrino]

Т. е. можно и PBEPBE, PBE1PBE, B3LYP, LSDA и другие. Я правильно понял?

[Neitrino]

У меня вот такая проблема:
Expectation value S**2: 0.7500
Original(UHF): 5.27
Дублетное состояние. Как можно добиться чтобы original был тождественен E. v.