- 1 名前：デフォルトの名無しさん [2006/01/24(火) 09:48:23 ]
- このスレッドは、他のスレッドでは書き込めない超低レベル、
もしくは質問者自身何が何だが分からない質問を勇気を持って書き込むスレッドです。
FORTRAN使いが優しくコメントを返しますが、
お礼はFORTRANの布教と初心者の救済をお願いします。
- 401 名前：デフォルトの名無しさん [2006/06/07(水) 15:16:53 ]
- Evca(Ω,x)を求めるプログラムを作りたいのですが教えてください。
Evca(Ω,x)=Ei(Ω,x)+E0(Ω,x)+E1(Ω,x)+E2(Ω,x)
Ei(Ω,x)=-1.79186*Z**(5/3)/r
E0(Ω,x)=Z(2.21/r**2-0.916/r+Ecorr)
Ecorr=-(0.115-0.031logr)
E1(Ω,x)=2πZe**2((1-X)Z**Cs*(Rm**Cs)**2(1+2/3*u**Cs)+xZ**K(Rm**k)**2(1+2/3u**K))
E2(Ω,x)=1/2∑(4πe**2/ΩG**2)**-1(V(G))**2*1/(1-f(G))*(1/ε(G)-1)
∑はG/=0
Gは逆格子ベクトルG(n1,n2,n3)
V(q)=(1-x)V**Cs+xV**k(q)
Vi(q)=-4πZe**2/Ωq**2((1+u)cos(qRm**i)-usin(qRm**i)/qRm) (i=Cs,K)
u=qRmcos(qRm)/sin(qRm)-qRmcos(qRm)
わからない値は後から自分で入れるので教えてください
状態方程式の計算をするらしいです
- 402 名前：デフォルトの名無しさん [2006/06/07(水) 15:24:39 ]
- E2(Ω,x)=1/2∑(4πε**2/ΩG**2)**-1(V(G))**2(1/(1-f(G)))(1/ε(G)-1)
∑はG/=0です
Ｇは逆格子ベクトルです。G(n1,n2,n3)
これを求めるプログラムはどうしたらいいか教えてください
わからない値は自分でいれます。お願いします。
|
 |
