трудности «плывущих» констант

[Котофеич]

Mike4july1972 писал(а):

Интеграл Фейнмана используется в случае непрерывного поля (теория поля), а в случае конечного числа частиц используется диаграммная техника Фейнмана?

Ну Вы Михаил и дали копоти. Диаграмная техника енто теория возмущений.
Она получается путем разложения бессмысленного фейнмановского интеграла, в
осмысоленный ряд теории возмущений. Почитайте книжку Фейнмана
Квантовая механика и "интегралы" по траекториям.

[Аурелиано Буэндиа]

Mike4july1972 писал(а):

Интеграл Фейнмана используется в случае непрерывного поля (теория поля), а в случае конечного числа частиц используется диаграммная техника Фейнмана?

Интреграл Фейнмана это только способ квантования. Кантовать можно бесконечномерную динамическую систему (поле) и конечномерную (классические частицы). А диаграммная техника это теория возмущения. Полностью согласен с Котофеичем.

[Mike4july1972]

Алексей, Харьков:
Майк, я сейчас занимаюсь изучением фейнмановской диаграммной техники. Эта техника по существу является мощным методом терии многочастичных систем. Мне диаграммная техника нужна для применения в теории сверхпроводимости, которой я занимаюсь.
Смысл диаграммной техники в том, чтобы по теории возмущении записать ряд (с применением графического диаграммного метода) для оператора эволюции системы (пропагатор, функция Грина) по степеням оператора взаимодействия для интересующих нас частиц, и просуммировать этот ряд. За счет взамодействия частиц в многочастичной системе, получается что вместо обычной массы и потенциала для свободной частицы, частица преобретает эффективную массу и эффективный потенциал (н., обычный кулоновский потенциал становится экранированным потенциалом Юкавы).
Т.о., мы вычисляем параметры "квазичастиц" и можем анализировать процессы в твердых телах на языке "квазичастиц", что значительно проще, чем учитывать все частицы и взаимодействия в исходном виде.
Хотя... кажется в вашей задаче о расчете уровней энергии атомов принцип применения диаграммной техники другой... По поводу диаграммной техники в многочастичных системах есть замечательная книга:
A.M. Zagoskin, Quantum theory of many-body systems. Если нужно, я могу вам ее прислать.
Во всяком случае, я могу помочь вам с идеей самого мтода, он оказывается изящным и красивым.
_________
Г.Б. Аракелян, "Числа и величины в современной физике":
Таким образом, анализ размерностей в применении к постоянным c, h, e, G, G_F, m, t_H позволяет среди множества всех вариантов выделить чисто формальным путем тройки c, h, e и G, m , e, из которых только и могут быть составлены безразмерные комбинации alpha и alpha_G(e, m). При этом, фундаментальной величиной безоговорочно следует считать пока лишь постоянную alpha, поскольку истинное значение alpha_G остается неопределенным из-за отсутствия строгих критериев определения фундаментальной массы m_G. В любой системе измерения, построенной на физических постоянных, необходимо принять во внимание отношение e^2/mc ~ 1/137, как некую теоретическую и эмпирическую данность, подсказанную самой природой и вместе с тем выделенную и формально.

[Котофеич]

Ну вот видите Михаил Алексей из Харькова тоже самое говорит.

[Алексей_]

Алексей из Харькова это я Мне очень понравилась книга A.M. Zagoskin, Quantum theory of many-body systems, по применению методов квантовой теории поля в физике конденсировааного состояния, скачанную с вашего сайта. Там последовательно и доступно излагается теория, а также
и приложения (правда встречаются опечатки, надо читать внимательно). Кстати этот А. Загоскин
когда то работал в нашем ФТИНТе. Сейчас я лишен возможности скачивать книги, о чем очень сожалею.

Котофеич писал(а):

Ну вот видите Михаил Алексей из Харькова тоже самое говорит.

[Mike4july1972]

У кого какие есть серьезные публикации и выступления на конференциях по физике? Дайте, пожалуйста, ссылки, чтобы я мог прочитать о Ваших достижениях. Я хоть и дурак и псих, но у меня пару десятков публикаций из них 4 в Physical Review.

[Mike4july1972]

Привожу начало компьютерной программы, реализующей метод наложения конфигураций (автор – Владимир Дзюба из Австралии)
--
Computer code:
Created: Dec 03, 03; 13:34 (probably)
Printed out on August 28, 2004 in Sydney, Australia.
Method which is implemented by this computer code: Configuration interaction for n electrons (CIN).
Name of the file: cin1.f
Implicit doubleprecision(a-h, o-z)
INCLUDE “cin.h”
Common /ne/ ne /nd/nd /mj/mje /nc/nc, ncr /nv/nv
/ /Nlev/Nlev /ER0/ ER0 fk/fk(9) /nso/nlist0
/ /ENERGY/EN0, ICUT, ESIGMA(9) /vsc/vsc(9) /DE/ICUTQ, DE
C common /work/work(IWK)
Dimension nlmax1(5, 5), nlmax2(5, 5), nlmax(9)
Dimension ev(NT1)
Real *8, allocatable :: work(:)
REAL *8 er0, ev, ee
CHARACTER integral1 *20, integral2 *20, integrale *20, LET(6)
Open(unit=1, file=’cin1.dat’,STATUS=’OLD’, form=’formatted’)
READ(1, 12) integral1
READ(1, 12) integral2
READ(1, 12) integrale
12 FORMAT(A20)
PRINT *, ‘File for one-electron integrals-’, integral1
PRINT *, ‘File for two-electron integrals-’, integral2
READ(1, *) ER0
READ(1, *) EN0, ICUT, ICUTQ
READ(1, *) (vsc(I), I=1,9)
READ(1, *)(fk(I), I=1,9)
READ(1, *)(ESIGMA(I), I=1,9)
READ(1, *) amj
READ(1, *) Nlev
READ(1, *) iconf1, lmax1
Do j=1, iconf1
Read(1, *) (nlmax1(i,j), i=1, lmax1+1)
End do
Call readint1(integral1, integral2, integrale, lmax, nlmax)
Call makegaunt(2*lmax+1)
MJe = 2.*ABS(amj)+0.01
Call dinit(nlist0)
Call donrl(nlist0)
Call rnrc(ne, lp)
Call config(ne, lp, iconf1, lmax1, nlmax1)
C call outcnr(ne)
Call mackerel(ne)
C call outcrel(ne)
Call jterm(1)
C call outdets(ne)
Call formj(1)
CALL FORMH
C IWK=2*NXX*LIM+LIM*LIM+(NUME+11)*LIM+2*NUME
IWK=2*NV*LIM+LIM*LIM+(NUME+11)*LIM+2*NUME
C limi=nlev+2
C IWK=2*NV*limi+ limi**2+(NUME+11)* limi+2*NUME+NBLOCK
Allocate(work(IWK))
If Allocate(work(IWK)) then
Print *, ‘work – ok’
Else ‘allocation of work failed’
End if
Print *, ‘work:’, size(work)
Call DODAVE(NV, NLEV, EV, WORK(1), iwk)
Call OUTP(NV, NLEV, EV, WORK(1))
STOP
END
--
Here:
NV – полное количество базисных функций.
DETS - указывает пределы списка (какие номера детерминантов соответствуют каким базисным функциям).
DODAVE - реализует диагонализацию матрицы методом Давидсона.
MATRIJ – вычисляет многоэлектронный матричный элемент для двух детерминантов.
Gdet(N, IDET1): N- номер детерминанта;
IDET1 - список одноэлектронных состояний, соответствующих данному детерминанту.
RSPQ1 – сравнивает детерминанты на предмет того, каким количеством функций они отличаются.
I1, I2, J1, J2 – номера одноэлектронных состояний.

[Аурелиано Буэндиа]

Mike4july1972 писал(а):

Привожу начало компьютерной программы, реализующей метод наложения конфигураций (автор – Владимир Дзюба из Австралии)
--
Computer code:

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

[Mike4july1972]

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

[Аурелиано Буэндиа]

Mike4july1972 писал(а):

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

Как по-вашему любая ли информация представляет интерес? Какой интерес может представлять этот огрызок программы?
Предлагаю Вам (если Вы действительно отзывчивый человек) выложить эту программу целиком с подробными мануалом где-нибудь в интернете и дать здесь ссылку.

[LynxGAV]

Михаил, у Вас очень странный метод написания докторской. Я бы и с кандидатской так не рискнула. Чисто по-человески, складывается впечатление, что Вы своим руководителям делаете рекламу. Как можно с пятью-то дипломами спрашивать что-то? Тут рассказывать, рассказывать и рассказывать..

[Аурелиано Буэндиа]

LynxGAV писал(а):

Чисто по-человески, складывается впечатление, что Вы своим руководителям делаете рекламу. Как можно с пятью-то дипломами спрашивать что-то? Тут рассказывать, рассказывать и рассказывать..

Все проще: он не руководителям делает рекламу, а себе.
Программу так и не выложил

[Mike4july1972]

Почему сразу надо подозревать в какой-либо нечестности?
Мне слишком мало осталось жить, чтобы я стремился специально поступать нечестно.
Я просто пытаюсь приблизиться к пониманию истины.
Я не стремлюсь делать никому рекламы.
Мой уровень низкий, я с этим не спорю.
В Австралии нет кандидатской диссертации, есть только докторская.
Программу по частям меня просили не выкладывать. Я вынужден набирать текст программы с распечатки, поэтому получается медленно.
Раз Вы настаиваете, привожу еще один фрагмент программы.

Code (additional(2-d part)):
Subroutine readint1(integral1, integral2, integrale, lmax1, nlmax)
Implicit doubleprecision(a-h,o-z)
Include “cin.h”
Common/some/jmin, jmax, nhint /NSO/NSO
Common/nnn/ee(NVM), nnn(NVM), kk(NVM), ll(NVM), jj(NVM), nlist
Common/rint1/rint1(NRI1), rsig(NRI1), rdsig(NRI1), iint1(NRI1)
Common/qqq/qq(NIMAX2), nad(NIMAX2), kat(NKAT2), ngint, katomax
Common/qpar/nvl(NVM), nvlist
Common/qqqe/qe(NIMAX1), qqe(NIMAX1), nade(NIMAX1),
. kate(NKAT1), ngintq, katomq /fk/fk(9), fkk(9)
Common/val/nlistq, nnv(NVM), kapv(NVM) /ns/ns /de/icutq
Dimension nlmax(9)
Real qq, qe, qqe
CHARACTER integral1 *20, integral2 *20, integrale *20
OPEN(UNIT=13, FILE= integral1, STATUS =’old’,
, FORM=’UNFORMATTED’)
Do i=1,9
Nmax(i)=0
End do
Read(13) jmin, jmax, jmax, nlist
Jmin=jmin
Jminq= jmin
NSO= jmin
Ns=nlist
If(nlist.gt.NVM) then
Write(6, 106) nlist, NVM
106 format(‘Too many basis states: nlist=’,i4, ‘while NVM=’, i4)
Stop
End if
Read(13)(nnn(I), kk(I), LL(I), JJ(I), EE(I), I=1, nlist)
Lmax1=0
Do i=1, nlist
Write(6,6)i, nnn(i), kk(i), ll(i), jj(i), ee(i)
C Write(10,6)i, nnn(i), kk(i), ll(i), jj(i), ee(i)
6 format(i4,’)n=,i2, ‘kap=’,i2,’l=’,i2, ‘j=’,i2,’/2en’,
, f12.6)
If(i.eq.nso) write(6,7)
C If(i.eq.nso) write(10,7)
7 format(‘______________’)
If(kk(i).lt.0.and.i.gt.jmin) then
Nlmax(ll(i)+1)= Nlmax(ll(i)+1)+1
Lmax1=max0(lmax1,ll(i))
End if
End do
Print 5, lmax1, (nlmax(l), l=1, lmax1+1)
5 format(/’Number of basis functions for l=0 to’,i3,’:’/
/ 10i4)
Read(13) NHINT, NLIST11
If(nhint.gt.NRI1) then
Write(6, 107) nhint, NRT1
107 format(‘Too many integrals: nhint=’, i6, ‘’,i6)
Stop
End if
Read(13)(RINT1(I), I=1,NHINT)
Read(13)(IINT1(I), I=1,NHINT)
Read(13)(RSIG(I), I=1,NHINT)
Read(13)(RDSIG(I), I=1,NHINT)
Close(13)
*****************************
OPEN(UNIT=13, FILE=integral2, STATUS=’old’,
, FORM=’UNFORMATTED’)
Read(13) jmin, nlist, ngint, katom
If(ngint.gt.NIMAX2.or.katom.gt.NKAT2) then
Print 110, ngint, NIMAX2, katom, NKAT2
110 format(‘Too much date to read:’,/’ngint=’,i7,
, ‘NIMAX2=’,i7, ‘katom=’,i6, ‘NKAT2=’,i6)
Stop
End if
Read(13)
Print *, ngint, NIMAX2, katom, NKAT2
Read(13)(kat(i), i=1,katom)
Read(13)(nad(i), i=1,ngint)
Read(13)(qq(i), i=1, ngint)
Close(13)
*****************************
If(icut.eq.0) then
Nvl=0
Nvlist=0
Fkk=1.
Else
OPEN(UNIT=13, FILE=integrale, STATUS=’old’,
, FORM=’UNFORMATTED’)
Read(13)nlistq,ngintq,katomq
Read(13)(nnv(i), kapv(i), i=1, nlistq)
Read(13)(qqe(i), i=1, ngintq)
Read(13)(nade(i), i=1, ngintq)
Read(13)(kate(i), i=1, katomq)
Read(13)(fkk(i), i=1, 9)
Close(13)
Nvlist=nlistq
C print *, ‘nlistq=’, nlistq
C print 117,(nnv(i), kapv(i), i=1, nlistq)
C 117 format(‘nnv=’,i4,’kapv=’,i4)
Do in=jmin+1, nlist
Do iv=1, nlistq
If(kk(in).eq.kapv(iv).and.nnn(in).eq.nnv(iv))nvl(iv)=in
End do
End do
Print 17, (fkk(i), i=1,9)
17 format(‘fk:’,5f10.5)
Return
901 stop
end
в файле cin.h видимо заданы максимальные границы изменения ряда параметров компьютерной программы, таких например, как максимальный размер матрицы (NXX) для диагонализации методом Давидсона.

[Mike4july1972]

Я пытаюсь восстановить компьютерные программы, реализующие метод наложения конфигураций, много-частичную теорию возмущений, метод Хартри-Фока.
Подпрограмма FORMJ формирует конфигурации для данного момента J?
(part 3 of the code:)
Subroutine RSPQ1(N, IDET1, IDET2, IS, NF, I1, I2)
DIMENSION IDET1(128), IDET2(128)
c-----------------------------
c write(*, 5)(idet1(i), i=1,n)
c write(*, 5)(idet2(i), i=1,n)
5 format(15x, 10I3)
IS=1
Ni=0
Nj=0
I2=0
J2=0
Nf=3
I=1
J=1
L0=0
200 l1=idet1(i)
L2= idet2(j)
If(l1.LT.l0) Goto 1100
L0=l1
If(l1.NE.l2) Goto 210
I=i+1
J=j+1
If(i.LE.n.and.j.LE.n) Goto 200
Goto 230
210 if(l1.GT.l2) Goto 220
Ni=ni+1
If(ni.GT.2) return
I1=i2
I2=i
I=i+1
If(i.LE.n) Goto 200
Goto 230
220 nj=nj+1
If(nj.GT.2) return
J1=j2
J2=j
J=j+1
If(j.LE.n) Goto 200
230 nf=ni
If(i.LE.j) Goto 240
Do 10 k=j,n
J1=j2
J2=k
10 continue
Goto 250
240 if(i.EQ.j) Goto 250
Nf=nj
Do 20 k=I,n
I1=i2
I2=k
20 continue
250 if(nf.EQ.0) return
If(j2.EQ.i2) Goto 260
K=iabs(j2-i2)
If(k.NE.2*(k/2)) is=-is
K2=idet2(j2)
260 if(nf.EQ.1) Goto 1000
If(j1.EQ.i1) Goto 1000
K=iabs(j1-i1)
If(k.NE.2*(k/2)) is=-is
K=idet2(j1)
C
C ????? error ???!!!
Idet2(j1)= Idet2(j1)
C ????? error ???!!!
C
Idet2(i1)=k
1000 if(j2.ne.i2) idet2(i2)=k2
Return
1100 write(*, 15) l0, l1
15 format(1x, ‘Rspq: wrong order of positions’,
! I3, ‘and’, I3, ‘in one of configurations’)
Stop
end
Здесь в набранном мной тексте программы возможна путаница буквы l и цифры 1, а так же буквы о и цифры 0.

[Аурелиано Буэндиа]

Mike4july1972 писал(а):

Я не стремлюсь делать никому рекламы.
Мой уровень низкий, я с этим не спорю.
В Австралии нет кандидатской диссертации, есть только докторская.

Если это Philosophy Doctor то это и есть кандидатская (читайте здесь)! Уточните у Фламбума. Если у Вас уже есть звание кандидата, то лучше поискать postdoc чем стремиться получить PhD. Какой у вас уровень это еще не известно. Может и не слабый. Подозреваю что Вы просто прикидываетесь.

Mike4july1972 писал(а):

Программу по частям меня просили не выкладывать. Я вынужден набирать текст программы с распечатки, поэтому получается медленно.
Раз Вы настаиваете, привожу еще один фрагмент программы.

Еще раз говорю: прекратите заваливать форум бессмысленными фрагментами кода. И это не просьба.

Вот посмотрите:
http://web.am.qub.ac.uk/DARC/ (Код написанный в Белфасте)
Возможно это вам пригодится