棒エレクトレットの分極の回転

 電場  \mathbf{E}電束密度(電気変位) \mathbf{D} の違いを説明するときに棒エレクトレット(棒電石)がしばしば登場する.棒エレクトレットは棒磁石の「電気版」である.理想的な棒エレクトレットは円柱状の物質で,軸方向に一様な分極をもつ.自由に動ける電荷が存在しない場合,分極  \mathbf{P} しか存在しないので電束密度の発散は0になる.


 \begin{align} \boldsymbol{\nabla}\cdot\mathbf{D}=0 \end{align}


これからただちに  \mathbf{D}=0 としてしまうと,エレクトレット外部では  \mathbf{E}=0 となってしまい矛盾する.電束密度の発散が0でも,回転が0であるとは限らない.


 \begin{align} \boldsymbol{\nabla}\times\mathbf{D}=\epsilon_0  \boldsymbol{\nabla}\times\mathbf{E} +  \boldsymbol{\nabla}\times\mathbf{P} =  \boldsymbol{\nabla}\times\mathbf{P}   \end{align}


つまり分極の回転が0でなければ電束密度の回転も0にならない.しかし,棒エレクトレットは一様な分極をもつので一見すると  \boldsymbol{\nabla}\times\mathbf{P} は0になるように思える.グリフィス「電磁気学I」ではこの問題提起だけをして,解答は読者に任せる,と書いてある.


 エレクトレット内部では分極は一様なので明らかに回転は0である.0でない回転をもつとしたらエレクトレットの表面部分しかありえない.エレクトレットを半径  a,長さ  L の円柱とし,分極の方向を  z 軸方向として分極の大きさを  P_0 と書くと,任意の位置の分極は


 \begin{align} \mathbf{P} = P_0 \, \theta(z) \theta(L-z) \theta(a-\rho) \hat{\mathbf{z}} \end{align}


と書ける.ここで円柱座標  (\rho, \phi, z) を使った. \mathbf{P} z 成分しかもたず, \phi によらないことから,


 \begin{align} \boldsymbol{\nabla}\times \mathbf{P} = -\frac{\partial P_z}{\partial \rho} \hat{\boldsymbol{\phi}} = P_0 \, \theta(z) \theta(L-z) \delta(a-\rho) 
 \hat{\boldsymbol{\phi}}  \end{align}


となる.つまりエレクトレットの側面部で0でない(大きさ無限大の)回転をもつことになる.