|
|
Вопрос к математикам, физикам и химикам
11.10.2003 23:13
|
|
|
Задача(из области квантовой химии) :
Квантовохимическая формулировка:
Рассчитать(численно) интералл перекрывания орбиталей двух электронов различных атомов атомов( базис STO)
Математическая формулировка:
Есть две функции изначально заданные в сферических координатах и ортогональные друг другу при условии если точки соответствующие R=0 совпадают F1(R,teta,phi) и F2(R,teta,phi). Теперь проводим транслюцию одной из функций в пространстве на вектор (x,y,z). Необходимо вычислить скалярное произведение этих функций, т.е. найти(численно) итеграл произведения функции комплексно сопряженно к первой на вторую по всему пространству.
Конечно, можно перейти к координатам x,y,z и рассчитывать уже в них разбивая пространство на маленькие кубики и пересчитывая для каждой такой ячейки из (x,y,z) в (R1,teta1,phi1) и (R2,teta2,phi2) - получится весьма громоздко.
Существует ли какой нибудь более элегантный способ рассчета?
|
|
|
KOHTPA
|
|
Carpal Tunnel
|
|
|
|
|
|
|
Рег.: 22.01.2003
|
|
Сообщений: 33647
|
|
|
|
Рейтинг: 2374
|
|
|
Прочитай, что такое STO.
Это вполне определенного вида функции.
---
"Vyroba umelych lidi, slecno, je tovarni tajemstvi."
Karel Capek
|
|
|
|
Re: Вопрос к математикам, физикам и химикам
[re: KOHTPA]
12.10.2003 01:16
|
|
|
Да действительно ...
F(n,l,m;R, teta, phi) = const * r ^ f(n) *exp (alpha * r)*Y(l,m;teta,phi), где
Y(l,m; teta,phi) = exp(i*m*phi) * (1-cos(teta)^2)^(abs(m)/2) * d^(l+abs(m)) (cos(teta)^2-1)^l / d(cos(teta)) d - дифференцирование
|
|
|
who
|
|
old hand
|
|
|
|
|
|
|
Рег.: 03.12.2002
|
|
Сообщений: 841
|
|
Из: В1021пр
|
|
Рейтинг: 0
|
|
|
Виталий, с появлением в форуме. Самый естественный ответ - заглянуть в исходники gamess'а. По поводу аналитического взятия подобных интегралов, кажется, что-то есть в Абаренкове-Братцеве-Тулубе.
|
|
|
|
Re: Вопрос к математикам, физикам и химикам
[re: who]
12.10.2003 16:04
|
|
|
Спасибо Борислав.
А у тебя есть эта книжка?
Если нет, то где она может быть? (Химфак, ГЗ-7эт)
Можно ли у тебя взять Степановскую книжку?
Исходники гамесса ... да я представляю сколько там строк на фортране ...
... у меня есть исходники mpqc ... он на cpp там я смотрел ...  - в общем хрен разберешься ...
... названия переменных не отражают их смысла или по крайней мере он там хитро завуалирован ...
... полно неочевидных конструкций ... отсуствие комментариев.
|
|
|
|
Re: Вопрос к математикам, физикам и химикам
[re: who]
12.10.2003 18:16
|
|
|
Можно ли в этом случае перейти от комплексных Y(l,m,teta,phi) к действительным - т.е. поступить так как мы это делаем для водородоподобного атома?
|
|
|
who
|
|
old hand
|
|
|
|
|
|
|
Рег.: 03.12.2002
|
|
Сообщений: 841
|
|
Из: В1021пр
|
|
Рейтинг: 0
|
|
Re: Вопрос к математикам, физикам и химикам
[re: Anonymous]
12.10.2003 19:04
|
|
|
>Можно ли в этом случае перейти от комплексных Y(l,m,teta,phi) к действительным
Разумеется.
E[F(m, \ldots)]=E[F(-m,\ldots)]
Поэтому
G_1(m,\ldots)=(1/2)[F(m,\ldots)+F(-m,\ldots)]=cos(m\phi)*g_1
G_2(m,\ldots)=(i/2)[F(m,\ldots)-F(-m,\ldots)]=-sin(m\phi)*g_1
- вещественные функции из Span{F(m,\ldots),F(-m,\ldots)} c исходной E.
Книга Абаренкова et al. есть на 4-м этаже.
|
|
|
|
Re: Вопрос к математикам, физикам и химикам
[re: who]
12.10.2003 19:56
|
|
|
Хорошо.
Вот например если мы возьмем какой нибудь кластер из атомов B.
Если я правильно понял в STO каждый электрон представляется одной функцией.
Так вот рассмотрим его электрон n=2 l=1 m=-1 ему соответствует определенная функция - и всей системе определенное значение энергии.
А если мы присвоим ему не Y(1, -1) а линейную комбинацию (Y(1, -1) + Y(1, 1)) / sqrt(2) то ведь значение энергии ... как мне кажется получится другое ... насколько правомерно так поступать?
|
|
|
KOHTPA
|
|
Carpal Tunnel
|
|
|
|
|
|
|
Рег.: 22.01.2003
|
|
Сообщений: 33647
|
|
|
|
Рейтинг: 2374
|
|
Re: Вопрос к математикам, физикам и химикам
[re: Anonymous]
12.10.2003 20:20
|
|
|
Почему???
Энергия остается той же.
Ты просто "смотришь" с другой стороны, поэтому у тебя "оси повернуты",
а энергия-то остается.
---
"Vyroba umelych lidi, slecno, je tovarni tajemstvi."
Karel Capek
|
|
|
who
|
|
old hand
|
|
|
|
|
|
|
Рег.: 03.12.2002
|
|
Сообщений: 841
|
|
Из: В1021пр
|
|
Рейтинг: 0
|
|
Re: Вопрос к математикам, физикам и химикам
[re: KOHTPA]
12.10.2003 21:44
|
|
|
Оператор момента входит только в одноэлектронные компоненты гамильтониана
H=\sum_a h_1(a) + \sum_{a,b} h_2(a,b).
[Смысл предыдущей реплики]:
h_1[g(m)]=h_1[g(-m)]
|
|
|
who
|
|
old hand
|
|
|
|
|
|
|
Рег.: 03.12.2002
|
|
Сообщений: 841
|
|
Из: В1021пр
|
|
Рейтинг: 0
|
|
Re: Вопрос к математикам, физикам и химикам
[re: who]
12.10.2003 22:11
|
|
|
Кстати, нашел более общее объяснение, почему можно пользоваться любыми (унитарными) преобразованиями базиса (в частности, F-->G). Правда, только для метода ХФ: существует формулировка метода в терминах матрицы плотности \rho(r,r'), а \rho однозначно опр. _подпространством_ функций.
См. Степанов-Пупышев (Green Book), 269-270.
|
|
|
KOHTPA
|
|
Carpal Tunnel
|
|
|
|
|
|
|
Рег.: 22.01.2003
|
|
Сообщений: 33647
|
|
|
|
Рейтинг: 2374
|
|
Re: Вопрос к математикам, физикам и химикам
[re: who]
12.10.2003 22:27
|
|
|
> См. Степанов-Пупышев (Green Book), 269-270.
; ) Он у меня: а) есть дома; б) есть в лабе. 10 раз. : )
Но "Greeb Book" мне понравилось.
Надо завести тоже какие-нибудь обозначения: "Heavy Book," "Light Book."
---
"Vyroba umelych lidi, slecno, je tovarni tajemstvi."
Karel Capek
|
|
|
top
|
|
|
|
(10.0.0.73)
|
|
|
|
|
|
|
|
|
|
|
Рейтинг: 3392
|
|
Re: Вопрос к математикам, физикам и химикам
[re: KOHTPA]
14.10.2003 00:23
|
|
|
сразу видно что вы химики 
|
|
|
KOHTPA
|
|
Carpal Tunnel
|
|
|
|
|
|
|
Рег.: 22.01.2003
|
|
Сообщений: 33647
|
|
|
|
Рейтинг: 2374
|
|
Re: Вопрос к математикам, физикам и химикам
[re: top]
14.10.2003 16:05
|
|
|
И что из этого?
Мы все трое очень разные химики.
---
...Я работаю антинаучным аферистом...
|
|
|
top
|
|
|
|
(10.0.0.73)
|
|
|
|
|
|
|
|
|
|
|
Рейтинг: 3392
|
|
Re: Вопрос к математикам, физикам и химикам
[re: KOHTPA]
15.10.2003 14:40
|
|
|
мб вы и разные но вас объединяет то , что вы химики !
Вам замечание (+). Флуд в тематическом разделе. Хватит флудить.
Редактировал Morly (15.10.2003 14:47)
|
|
Список
Плоский
Дерево
