2-электронные мол. интегралы гауссовых функций

[Суровый]

Подскажите источник, где можно посмотреть подказки по взятию данных интеграллов. Или готовые решения.

[Shorku]

Подскажите источник, где можно посмотреть подказки по взятию данных интеграллов. Или готовые решения.

Гуглить libint за авторством Ed Valeev. Также можно погуглить его лекции по теме. Если кратко, то это либо относительно несложно, но убийственно долго, либо рекурсивные алгоритмы и ... проще готовым воспользоваться. А Вам зачем?

[ Post made via Android ]

[Суровый]

Гуглить libint за авторством Ed Valeev.

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

Также можно погуглить его лекции по теме.

Не получается найти.

Если кратко, то это либо относительно несложно, но убийственно долго, либо рекурсивные алгоритмы и ... проще готовым воспользоваться

Вот-вот. Я месяц не мог вять интегралы перекрывания, пока случайно не наткнулся на подсказку о специальном преобразовании .

А Вам зачем?

Делаю собственную программу с Х-Ф.

[Shorku]

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

А надо обязательно разбираться? Исходники прекрасно компилируются. Для хф этого вполне достаточно. Разбираться же в работе всего этого- еще месяц убъете, а то и больше. Лекции завтра прицеплю. Но, простите за любопытство, зачем Вам писать еще одну хф программу?

[ Post made via Android ]

[Суровый]

Чтобы реально разобраться как работает ХФ.
Расчет Х-Ф, если получится, добавлю в уже сушествующую программу, которая делает пока только анализ.

[Shorku]

Лекции

[Jeffry]

Гауссовских функций много - декартовские, радиальные, Гаусса-Эрмита и пр.
В базисе каких функций надо брать интегралы?

[Shorku]

Гауссовских функций много - декартовские, радиальные, Гаусса-Эрмита и пр.
В базисе каких функций надо брать интегралы?

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

[Jeffry]

декартовские, радиальные, Гаусса-Эрмита

Первые и последние берутся аналитически (первые - в паре статей Taketa, O-ohata, Huzinaga, последние - в статьях Maksić Z.B. и др.), вторые - разложениями в ряды по первым. Коэффициенты разложений надо искать. Где-то видел, но о-очень давно.
С первыми могу подсказать, как получить пару тождеств, используемых при их выводе. В оригинальной статье они просто даны без пояснений.

[VTur]

Journal of Mathematical Chemistry
http://link.springer.com/journal/10910
в помощь

[Суровый]

Лекции

Спасибо!

Гауссовских функций много - декартовские, радиальные, Гаусса-Эрмита и пр.В базисе каких функций надо брать интегралы?

Декартовских.

Всем спасибо за помощь. Сейчас буду разбираться. Может быть потом задам наводящие вопросы.

[АлександрI]

Охх, почитал и вспомнил свою очень давнишнюю мечту: взять какую-то простую-простую молекулу симметрии С1, выписать и взять все 1- и 2-электронные интегралы и пр., о чём говорит топикстартер.
Сделать всё, вплоть до межатомных расстояний по Малликену и Ко.
А потом сделать тот же расчёт SP того же HF, посмотреть листинг аутпута и (с удивлением) сказать: а ведь программа, которую я уже тысячу раз запускал, правильно считает...

[Shorku]

Мне, простите, для общего образования, а что такое межатомные расстояния по Малликену?

[ Post made via Android ]

[surius]

и (с удивлением) сказать: а ведь программа, которую я уже тысячу раз запускал, правильно считает...

пардон за оффтопик
А я бы, с еще большим удивлением сказал: ох, .пи-пи-пи. .пи-пи-пи. я нигде не ошибся.

[Гесс]

и (с удивлением) сказать: а ведь программа, которую я уже тысячу раз запускал, правильно считает...

пардон за оффтопик
А я бы, с еще большим удивлением сказал: ох, .пи-пи-пи. .пи-пи-пи. я нигде не ошибся.

"Садись Иванов, двойка. Ошибка в третьей строке и далее все неверно"

[kbob]

Поскольку 2-e интегралы (ERI) обычно не вычисляются по одиночке, а чаще для наборов примитивных GTO составляющих оболочку, существует несколько оптимизированных алгоритмов вычисления ERI, учитывающих этот факт.
https://rsc.anu.edu.au/~pgill/papers/061COLD.pdf
Схема Obara-Saika реализованная в Эдуардом Валеевым в библиотеке https://github.com/evaleev/libint, является наиболее распространенной, но не единственной.

Кстати, а вы с какой целью интересуетесь 2-e интегралами?
Делаете собственую програмку? На каком языке программирования? Delphy?

[АлександрI]

Мне, простите, для общего образования, а что такое межатомные расстояния по Малликену?

Когда дяди шутят, дети сидят смирно не задают вопросы

[kbob]

Чтобы реально разобраться как работает ХФ.
Расчет Х-Ф, если получится, добавлю в уже сушествующую программу, которая делает пока только анализ.

Расчет RHF тривиален и достаточно подробно описан во многих источниках.
http://sirius.chem.vt.edu/wiki/doku.php ... g:project3

и реализован:
https://github.com/rpmuller/pyquante2/t ... uante2/scf

https://github.com/evaleev/libint/tree/ ... rtree-fock

[АлександрI]

Расчет RHF тривиален и достаточно подробно описан во многих источниках

Step #1: Nuclear Repulsion Energy
........................................
Step #2: One-Electron Integrals
Read the AO-basis overlap,
....................................
(Это я так привёл цитату из первой ссылки)

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

Я ещё помню, в студенческое время я был на практике в Институте, не буду говорить в каком, но по физхимии. Чтобы отвязаться от моих глупых вопросов, мне ткнули книгу Босса. Открываю книгу, на первой странице - интеграл. Я опять за своё: как его взять, что туда поставить, как будет выглядеть результат и т.п. Но, опять,все закатывали глаза вверх и говорили: это знает только Он...
Увы... может уже тогда у меня появилась изложенная выше м е ч т а ...

[Shorku]

Когда дяди шутят, дети сидят смирно не задают вопросы

Простите, что? Вероятно, общаться в таком тоне Вы можете со своими студентами, но Вы не находите, что здесь это уже через край?

[ Post made via Android ]