Дай дураку каплкластер...

[amge]

кстати, кто пробовал подсовывать Природе неродные базисы?

Я пробовал. Ничего хорошего. С небольшими базисами типа 6-31G* нормально, а если побольше, да с диффузными функциями - и проблемы со сходимостью SCF.

[Droog_Andrey]

Для DFT пользуюсь def2-xZVP. А диффузные функции для больших систем можно добавлять только наружным атомам.

[Гесс]

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

[kbob]

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

Вы что вообще ищите? Геометрию основного состояния, геометрию переходного состояния, энергию ионизации, энергию активации, еще что.

[Гесс]

На данный момент - геометрию основного состояния, энергию ионизации, сродство к электрону и протону, теплоту образования, энергию напряжения. Далее список будет расширен.

Куплю квартиру - большую, дешево, в хорошем состоянии, возле метро.

[sanya1024]

Я пробовал. Ничего хорошего. С небольшими базисами типа 6-31G* нормально, а если побольше, да с диффузными функциями - и проблемы со сходимостью SCF.

Ага, я помню что вроде у Вас был такой опыт. Спасибо.

[Гесс]

MN12SX/maug-cc-pvQz//MN12SX/cc-pvdz для моего тестового набора наврал в потенциале ионизации и сродстве к электрону не более 20%, в сродстве к протону до 10%.

По NO - IP совпал с точностью до 0.07 eV, сродство к протону до 4 ккал/моль, сродство к электрону налажало - ушло в минуса до -0.92eV, при экспериментальном значении вроде как в +0.03eV.
Вероятно лажает и на димерах воды и гидратированных водных ионах, но пока лучшее из того что я получил.

[Гесс]

У нас в лабе считают (то ли предрассудок, то ли не совсем), что у каждой программы есть свой набор базисов, более "естественный", что ли, для нее. Что алгоритмы, скажем, Гауссиана и Гамесса, больше приспособлены для работы с Поппловскими базисами, а у немцев (Турбомоль и Орка) "родные" базисы -- Альрихсовские.

Угу, согласно мануалу Орки, раздел 8.1 мануала 2.9

• In our own research we almost exclusively use the basis sets of the kralsruhe group. They have been updated to the ‘def2’ set that is more consistent than the older basis sets.

Но "все в этой книге может оказаться ошибкой" (откуда?). В смысле, универсального решения нет, каждый раз надо смотреть на задачу конкретно.

Предположительно Ричард Бах, "Иллюзии" или "Карманный справочник мессии"

В смысле, универсального решения нет, каждый раз надо смотреть на задачу конкретно.

С одной стороны жалко, а с другой стороны если бы это было просто то все бы сделали раньше.

[Гесс]

Обновляя свое предыдущее утверждение от первого июля:
MN12SX оказался катастрофически слабым в геометриях аниона тримера воды, результаты бешено зависят от базиса, не похожи друг на друга, непохожи на литературные минимумы полученные другими DFTфункционалами а так же в MP2, и пр.
К коллекции Droog_Andrey "проблемных" молекул, могу предложить добавить анион NNO, у которого весьма неприятно пляшут длины связей (от метод-базисов), а угол по результатам МП2 и каплкластера отличается примерно на 10 градусов.

[sanya1024]

Предположительно Ричард Бах, "Иллюзии" или "Карманный справочник мессии"

Угадали

Беда еще в том, что как ни выбирай супер-пупер-базис для маленькой системы, как только перейдешь к большой, придется "урезать осетра" -- и к чему тогда были все пляски с бубнами?

[Yu/2]

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

А где эта коллекция?

[Droog_Andrey]

Делокализована по форуму

[Droog_Andrey]

анион NNO

Дык ето... он же не существует.

P.S. В смысле, там кислород почти отваливается: http://dx.doi.org/10.1063/1.1364679

[Гесс]

Я про состояние А1. его достаточно активно считали.
А вообще для молекулы есть скажем сродство к электрону и я не уверен к кому оно относится. http://webbook.nist.gov/cgi/cbook.cgi?I ... Energetics

[Droog_Andrey]

Не осиливаю что значит состояние A1, но догадываюсь, что это что-то изоэлектронное молекуле диоксида азота. По идее, такая метастабильная штука существует. Посчитаю-ка сам...

[Гесс]

Это обозначение из статьи которую вы дали. Оно же используется как минимум в трех последующих.
Я имею в виду молекулу с углом 130-143 градуса.

UPD
Вот тут просто собрано много разных цифр (правда больше по нейтралу)

lowdin-may.pdf

А вот тут есть разные рассчеты геометрии аниона

jp049204i.pdf

[Droog_Andrey]

Короче всё это фигня. Не получается такого аниона. Состояние его, кстати, обозначается ²A'.

А расчёты все можно смело выкидывать на помойку. Потому что сколько ни добавляй диффузных функций, электронное облако тупо разрастается до предела, а валентный угол приближается к 180 градусам. В итоге получаем молекулу NNO, окружённую очень далеко делокализованным, почти оторвавшимся электроном, и всё равно проигрываем по энергии нейтральной NNO.

Так что нет у этой частицы сродства к электрону. Максимум, что можно получить, - вышеупомянутый комплекс N₂·O^—, но это через разрыв связи NO, поэтому уже не то.

[Гесс]

Чето нифига не понял.
Неполучается согласно чего? У меня анион в CCSD/aug-cc-pvtz пришел на 132 градуса, в MP2 c тем же базисом в 142, Согласно приложенным статьям он сходится туда же и в других методах и базисах. Вы ж даже статьи не скачали.
Экспериментальное сродство к электрону как я уже писал есть, причем адиабатическое мелкое а VDE полтора эВ.

[Droog_Andrey]

А толку-то их скачивать. Сравните энергию, которая получилась для аниона, с энергией нейтрального NNO на этом же уровне теории. Пока энергия аниона выше, расчёт не имеет смысла.

Экспериментальное - это дважды вышеупомянутый комплекс.

[Гесс]

Сравните энергию, которая получилась для аниона, с энергией нейтрального NNO на этом же уровне теории.

Это я знаю. Это ничего не доказывает. Невероятное количество расчетов дает энергию аниона димера воды выше нейтрального димера. Что не мешает двум молекулам воды успешно локализовать электрон. Кстати я попробовал это несчастное (H2O)2 соптимизировать вашим U500+ на MP2. За двое суток на 40 процах система прошла меньше 7 SCFциклов.

Пока энергия аниона выше, расчёт не имеет смысла.

ИМХО достаточно спорное утверждение. Этот занятный минимум обсчитывается всеми кому не лень.

Экспериментальное - это дважды вышеупомянутый комплекс.

Как видно выше я неотрицал что это может быть именно так. Но я пока в этом не уверен однозначно.

есть скажем сродство к электрону и я не уверен к кому оно относится.