Растворитель (диполь) и MP2 в гауссиане, где взять Molcas

[Зара]

Добрый день!

Вопрос простой, но нигде не могу найти ответ:
Почему при расчете в гауссиане методом MP2 с использованием модели Онзагера если задавать радиус полости и диӕлектрическую константу, расчет не проходит
( #P MP2/6-31+G(d) SCRF=(Dipole, A0=4.10, Dielectric=2.023))
, тогда как если опустить A0 и Dielectric
( #P MP2/6-31+G(d) SCRF=(Dipole) )
- расчет идет на ура

В других базисах (напр, b3lyp) такой картины не наблюдается

ошибка, выдаваемая гауссианом: Only single-point PSCF with reaction field.

И еще вопрос:
Кто-нибудь знает, как и где можно скачать сорсы (или уже скомпиленную версию) Molcas под Линух

Заранее спасибо

[Nord]

Для расчета дипольных моментов нужно знать матрицу плотности в первом порядке ТВ. Соответственно, если ее вычисление не запрограммировано (а оно, например, требуется и для расчета градиентов по смещениям ядер), то рассчитать дипольные моменты нельзя. Хотя это можно сделать, например, конечно-разностным способом дифференцируя энергию системы во внешнем электричнском поле по его напряженности.

Что касается Molcas, то самый простой способ "достать" его -- купить

[Зара]

Nord, спасибо, что откликнулись, но я только начинаю считать, поӕтому такой вопрос:
Как узнать (или указать ключ в гауссиане для расчета) матрицу плотности в первом порядке ТВ? И почему в методе функционала плотности, ни в хартри фоке такого не происходит?

спасибо

[helicase]

это несколько странно, если явно опустить А0 и е, то они, видимо, задаются по умолчанию, что-то типа е=78, А0=1. так?

[Зара]

Я тоже так поняла, что они берутся по умолчанию. Но мне бы хотелось задать a0 и e в явном виде.
Nord пишет , что для расчета дипольных моментов нужно знать матрицу плотности в первом порядке ТВ. Можно как-нибудь задать ключ, чтобы гауссиан рассчитал ӕту матрицу?

[helicase]

Этим я хотел сказать вот что: например, в 98 гауссиане были проблемы в PCM с полостью, точнее с параметром разбиения TSNUM (расчет не удавалось завершить). Это было исправлено в 03 версии. Возможно, и здесь нечто подобное c Онзагером -попробуйте поварьировать А0.

[Зара]

Да, обязательно попробую.
helicase, думаете ӕто просто глюк?

[helicase]

да, не исключаю этого

[Зара]

да, попробую. helicase, спасибо за участие!

Все же, подскажите, плиз!!!, где взять molcas

[Nord]

Molcas - платная программа. Вот нажалуюсь на Вас Русу (я его, надеюсь, в пятницу увижу)

[Yu/2]

Все же, подскажите, плиз!!!, где взять molcas

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

[Зара]

По поводу molcas приношу свои извинения. Кое-где в интернете она числится как условно-бесплатная, в одном хранилище хим-совта почти дали скачать. Честно говоря, не знаю, кто такой Рус, но, в таком случае - то же ему мои извинения.
Что касается сложности - меня попросили разобраться что ӕто такое и как с ӕтим работать. Поскольку у нас в лаборатории никто никогда ее в глаза не видел - поручили ӕто мне. Yu/2, спасибо, что предупредили

[Вик]

По поводу molcas приношу свои извинения. Кое-где в интернете она числится как условно-бесплатная, в одном хранилище хим-совта почти дали скачать. Честно говоря, не знаю, кто такой Рус, но, в таком случае - то же ему мои извинения.
Что касается сложности - меня попросили разобраться что ӕто такое и как с ӕтим работать. Поскольку у нас в лаборатории никто никогда ее в глаза не видел - поручили ӕто мне. Yu/2, спасибо, что предупредили

Сейчас текущая версия Molcas 7 . Релиз был сравнительно недавно и в ближайшие 3-4 года нового не будет (это цикл разработки molcas),imho. Стоит он денег и условия лицензирования следующие:
вы покупаете лицензию на 1 год или на 3 года. В течении этого времени Molcas будет работать. Цена 3-х годичной лицензии - 2200 euro.

Денег своих стоит. Преимущества:

1) Удобный API;
2) Стабильный API (переход 6->7 менее болезнен нежели 5->6 )
Как следствие - приятно программировать под Molcas и имплементировать post SCF методы

3) Гибкая работа с базисами;
4) Оптимизирован под ANO базисы;
5) релиз 7-го был связан завершением работ по CD (Cholesky Decompositon) . Если нужны SP расчеты - отличный метод.

6) Отличная поддержка от Molcas's team ...

Это мое личное мнение

Home page: www.teokem.lu.se/molcas/

2Yu/2 в 7-ке уже более менее автоматизирована работа различных модулей. Но и в 6-ом уже не страшно и не сложно с ним работать -)

[Зара]

Вик, спасибо, что все так подробно рассказали. Буду надеяться, что когда-нибудь институт решится купить молкас.
Теперь, по крайней мере, благодаря участникам форума имею кое-какое представление о программе.

Всем спасибо !!!