Просьба перепроверить результат NwChem

[kbob]

Добрый день.
Запускаю расчет RI-MP2 энергии для молекулы кубана (C8H8) в различных базисах (геометрия равновесная).
Для cc-pCVTZ, cc-pCVQZ получаю нормальные значения, но для базиса cc-pCV5Z значение MP2-энергии не реально большое.

Входные параметры
[cut]title "cubane cc-pV5Z MP2-RI energy"

geometry
symmetry Oh
H 1.41063149 1.41063149 1.41063149
C 0.78371973 0.78371973 0.78371973
end

basis
H library cc-pV5Z
C library cc-pCV5Z
end

basis "ri-mp2 basis"
H library cc-pV5Z-ri
C library cc-pV5Z-ri
end

task rimp2[/cut]

Выходные значения энергий
[cut]-------------------------------------------
SCF energy -307.507255162945
Correlation energy -4369.270156917300
Total MP2 energy -4676.777412080244
-------------------------------------------[/cut]

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

[endlesslake]

nwchem branch = Development
nwchem revision = 23558

SCF energy -307.507255155134
Correlation energy -487.606592219138
Total MP2 energy -795.113847374272

больше на ошибку похоже

[VTur]

А зачем сжатые базисы использовали?

[endlesslake]

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

тот же инпут, с использованием сферических гармоник

-------------------------------------------
SCF energy -307.507287526028
Correlation energy -1.815754093360
Total MP2 energy -309.323041619388
-------------------------------------------

[kbob]

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

тот же инпут, с использованием сферических гармоник

-------------------------------------------
SCF energy -307.507287526028
Correlation energy -1.815754093360
Total MP2 energy -309.323041619388
-------------------------------------------

Покажите ваш input-файл если не трудно. Можно в топик на форуме NwChem отписать.
http://www.nwchem-sw.org/index.php/Spec ... o_big.html

[endlesslake]

Код: Выделить всё

title "cubane cc-pV5Z MP2-RI energy"

geometry
symmetry Oh
H 1.41063149 1.41063149 1.41063149
C 0.78371973 0.78371973 0.78371973
end

basis spherical
H library cc-pV5Z
C library cc-pCV5Z
end

basis "ri-mp2 basis"
H library cc-pV5Z-ri
C library cc-pV5Z-ri
end

scf
direct
end

task rimp2