Loading [MathJax]/jax/output/CommonHTML/jax.js

2022年11月24日木曜日

軸対称電荷分布(2)

軸対称電荷分布(1)からの続き

前回のz軸対称で,xy平面にも対称な電荷分布のつくるクーロンポテンシャルは,
V(r,s)=e4πε0n=0Pn(s)10R(t)0ρ(r,t)Pn(t){1r>(r<r>)n}4πr2drdt であった(ただし,s=cosθt=cosθ)。

最初に球対称の場合で,上記の式を確認する。このときR(t)=Rρ(r,t)=ρ(r)=H(Rr)ρ0である。H(x) ヘヴィサイドの階段関数4π3R3ρ0=Zeである。
このとき n=0 の項だけが残り,V(r)=e4πε0R0ρ01r>4πr2dr

(1) r>R の場合:
V(r)=e4πε0R0ρ01r4πr2dr=e4πε04πρ01rR0r2dr=Ze24πε01r

(2) r<R の場合:
V(r)=e4πε04πρ0{r01rr2dr+Rr1rr2dr}
=e4πε04πρ0{1rr33+[r22]Rr}=Ze24πε0{r2R3+3(R2r2)2R3}
=Ze24πε01R{32r22R2}

これで,球対称一様電荷分布によるクーロンポテンシャルが再現できた。

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 件のコメント: