Релятивистские базисы и порядок связи

[eFo]

Добрый день. Только начал осваиваться в квантовой химии, и возник такой вопрос.. Корректно ли использовать релятивистские базисы (типа L2) для расчёта порядка связей? Просто так получилось, что, например, при переходе от нерелятивистского базиса к релятивистскому длины связей уменьшаются, плотность в критических точках соответствующих связей увеличивается, а порядок связи существенно уменьшается (бывает даже в 2 раза), что, как мне кажется, нелогично.. Или всё норм? Спасибо за любой совет, что почитать \ ответ )

[Лечащий Врач]

Уточните, как именно считаете порядок связи, а также какие именно используете метод и базис.

[amge]

Раз релятивистский базис типа L2, то, скорее всего, это Природа. Природа считает порядки связей по Малликену. Поэтому я бы этим порядкам связей не доверял (по крайней мере, малликеновские же заряды при достаточно хорошем базисе ничего разумного из себя не представляют).

[Лечащий Врач]

Интересно тогда, каким именно образом в Природе реализована теория Бейдера. Потому что если есть возмжность получить wfn-файл, то дальше можно воспользоваться делокализационными индексами для расчета порядков связей. А если только внутренний анализ есть, то увы.

[sanya1024]

Откуда в Природе Бейдер? нет его там, и не было... Есть возможность получить орбитали в формате MOLDEN. Если из MOLDEN-а после этого можно что-то похожее на wfn вытащить, то можно работать дальше. Еще можно тот же молденовский файл считать ChemCraft-ом, но там из орбиталей можно получить только cube файлы (орбиталей по отдельности, разных их комбинаций, а также полной электронной и спиновой плотности). Не знаю, может ли это как-то пригодиться...

eFo, когда делали релятивистский расчет, не забыли использовать файл basis4.in и four=1?

[EvgeniX]

Есть, есть Бейдер в ПРИРОДе. Что-то вроде print=+density.
QCC Front-End умеет визуализировать критические точки.

Надо ещё разобраться, как лаплассиан анализировать.

[eFo]

sanya1024, да, и базис такой юзал, и four=1 ставил. Как сказал EvgeniX, критические точки можно находить с помощью print=+density и смотреть в QCC.
Т.е. дело всё-таки в способе расчёта порядка связей?

[sanya1024]

Интересно век живи -- век учись )) А где QCC берут? (у меня когда-то была давнишняя версия, да давно протухла...)

[amge]

EvgeniX , eFo, спасибо за информацию! Я тоже не знал, что в Природе есть Бейдер. В инструкции написано только про print=+bonds (малликеновские порядки связей). Надо будет поразбираться.

[EvgeniX]

Интересно век живи -- век учись )) А где QCC берут? (у меня когда-то была давнишняя версия, да давно протухла...)

Старую и используем

EvgeniX , eFo, спасибо за информацию! Я тоже не знал, что в Природе есть Бейдер. В инструкции написано только про print=+bonds (малликеновские порядки связей). Надо будет поразбираться.

Бейдер+DFT+металлокомплексы дают не совсем правду.

[amge]

Бейдер+DFT+металлокомплексы дают не совсем правду.

Только Природовская реализация или вообще?

[EvgeniX]

Я только в природовской реализации смотрел, но поверхностно.

[eFo]

Ладно, спасибо за обсуждение, будем думать... sanya1024, у меня QCC v2.09, не знаю насколько стар... Просто я слышал, что создатель этого QCC слишком рано умер..

[sanya1024]

Ладно, спасибо за обсуждение, будем думать... sanya1024, у меня QCC v2.09, не знаю насколько стар...

Я уж теперь и не помню, какая была версия. Лицензия была привязана к железу, и при смене компа, ес-но, протухла. На noname была ссылка на депозит, но оттуда файл удалили. Может, где еще лежит?

Просто я слышал, что создатель этого QCC слишком рано умер..

Не знала

[eFo]

Ладно, спасибо за обсуждение, будем думать... sanya1024, у меня QCC v2.09, не знаю насколько стар...

Я уж теперь и не помню, какая была версия. Лицензия была привязана к железу, и при смене компа, ес-но, протухла. На noname была ссылка на депозит, но оттуда файл удалили. Может, где еще лежит?

Просто я слышал, что создатель этого QCC слишком рано умер..

Не знала

Я могу послать, если надо... С keygen'ом

[sanya1024]

Спасибо.

[EvgeniX]

При первом запуске вводим любой номер.

[sanya1024]

Спасибо и Вам, EvgeniX.

[EvgeniX]

Не за что. Жаль, что автор не сможет доделать программу, но хотя бы она в свободном доступе будет.
(Беседин Д.В. и другими интересными проектами для iXBT.com занимался)

[amge]

Коллеги, научите, пожалуйста, как извлекать из природовского расчета print=+density какую-нибудь информацию, кроме координат критических точек (порядки связей, заряды, etc.)