前回のz軸対称で,x−y平面にも対称な電荷分布のつくるクーロンポテンシャルは,
V(r,s)=−e4πε0∞∑n=0Pn(s)∫10∫R(t)0ρ(r′,t)Pn(t){1r>(r<r>)n}4πr′2dr′dt であった(ただし,s=cosθ,t=cosθ′)。
このとき n=0 の項だけが残り,V(r)=−e4πε0∫R0ρ01r>4πr′2dr′
(1) r>R の場合:
V(r)=−e4πε0∫R0ρ01r4πr′2dr′=−e4πε04πρ01r∫R0r′2dr′=−Ze24πε01r
(2) r<R の場合:
V(r)=−e4πε04πρ0{∫r01rr′2dr′+∫Rr1r′r′2dr′}
=−e4πε04πρ0{1rr33+[r′22]Rr}=−Ze24πε0{r2R3+3(R2−r2)2R3}
=−Ze24πε01R{32−r22R2}
これで,球対称一様電荷分布によるクーロンポテンシャルが再現できた。
Mathematicaでは次のように計算する。
v[r_] := -3*NIntegrate[s^2/(r*HeavisideTheta[r-s] + s*HeavisideTheta[s-r]), {s, 0, 1}]Plot[v[r], {r, 0, 5}, PlotRange -> {0, -1.6}]
図:球対称一様電荷分布のクーロンポテンシャル(R=1.0, V(R)=-1.0)
0 件のコメント:
コメントを投稿