ベクトルラプラシアン

POINT

  • 作成中...(未完成,未検算)
  • ベクトルラプラシアンの計算方法と注意点.

ベクトルラプラシアンとは

ベクトルラプラシアン
ベクトル場(ベクトル値関数)$\boldsymbol{A}(\boldsymbol{x})=\begin{pmatrix}A_x(\boldsymbol{x}) \\ A_y(\boldsymbol{x}) \\A_z(\boldsymbol{x})\end{pmatrix}$に対して,デカルト座標系のベクトルラプラシアンを
\begin{aligned}
\boldsymbol{\nabla}^2 \boldsymbol{A}(\boldsymbol{x})
=
\begin{pmatrix}
\boldsymbol{\nabla}^2 A_x(\boldsymbol{x}) \\
\boldsymbol{\nabla}^2 A_y(\boldsymbol{x}) \\
\boldsymbol{\nabla}^2 A_z(\boldsymbol{x})
\end{pmatrix}
\end{aligned}
と定義する.但し,右辺のラプラシアンはスカラーに作用する通常のラプラシアンである.

$\mathrm{grad}$, $\mathrm{div}$, $\mathrm{rot}$や通常のラプラシアンと同じように,他の座標系では上を座標変換して得られた関係式でベクトルラプラシアンが定義されます.

ベクトルラプラシアンは,ベクトル解析の関係式

\begin{aligned}
\boldsymbol{\nabla}^2 \boldsymbol{A}
=\mathrm{grad\,}\mathrm{div\,}\boldsymbol{A}-\mathrm{rot\,}\mathrm{rot\,}\boldsymbol{A}
\end{aligned}
を利用して計算することもできます.この関係式は,次の記事で導出しています($\mathrm{rot\,}\mathrm{rot\,}\boldsymbol{A}$$=\mathrm{grad\,}\mathrm{div\,}\boldsymbol{A}-\boldsymbol{\Delta\,A}$):ベクトル解析の公式 - Notes_JP

あるいは,この式をベクトルラプラシアンの定義とすることもできます.

一般式

以下,$-g=|\mathrm{det\,}(g_{\mu\nu})|$と表すことにする.
ベクトルラプラシアン
共変微分(参考:共変微分の計算法 - Notes_JP)を用いてベクトルラプラシアンを計算すると
\begin{aligned}
&\boldsymbol{\nabla}_\alpha \boldsymbol{\nabla}^\alpha \boldsymbol{A} \\
&= \Biggl[ \frac{1}{\sqrt{-g}} \partial_\alpha (\sqrt{-g} g^{\alpha \beta} \partial_{\beta} A^\mu) \\
&\qquad + \frac{1}{\sqrt{-g}} \partial_\alpha (\sqrt{-g} g^{\alpha \beta} \Gamma^{\mu}_{\:\lambda\beta} A^{\lambda}) \\
&\qquad +g^{\alpha \beta} \Gamma^{\mu}_{\:\nu\alpha}(\partial_\beta A^{\nu} + \Gamma^{\nu}_{\:\lambda \beta} A^{\lambda} )
\Biggr] \boldsymbol{e}_\mu \\
&=\Biggl[\frac{1}{\sqrt{-g}} \partial_\alpha (\sqrt{-g} g^{\alpha \beta} A^{\mu}_{\: ;\beta}) \\
&\qquad +\Gamma^{\mu}_{\:\nu\alpha} g^{\alpha \beta} A^{\nu}_{\: ;\beta}
\Biggr]\boldsymbol{e}_\mu
\end{aligned}
となる.

【証明】
共変微分の計算法 - Notes_JPで計算したように,

\begin{aligned}
& \boldsymbol{\nabla}_\alpha \boldsymbol{\nabla}_\beta \boldsymbol{A}
= A^{\mu}_{\:\: ;\beta\,;\alpha} (\boldsymbol{e}_\mu \otimes \boldsymbol{f}^\alpha \otimes \boldsymbol{f}^\beta)\\
&A^{\mu}_{\:\: ;\beta\,;\alpha} \\
&=\partial_\alpha \partial_\beta A^{\mu}
+ ( \partial_\alpha\Gamma^{\mu}_{\:\lambda \beta} + \Gamma^{\mu}_{\:\nu\alpha} \Gamma^{\nu}_{\:\lambda \beta} )A^{\lambda} \\
&\quad + (\Gamma^{\mu}_{\:\lambda \beta} \partial_\alpha A^{\lambda} + \Gamma^{\mu}_{\:\lambda \alpha} \partial_\beta A^{\lambda}) \\
&\quad \textcolor{red}{ - \Gamma^{\nu}_{\:\beta\alpha} (\partial_\nu A^\mu + \Gamma^{\mu}_{\:\lambda\nu} A^\lambda) } \\
&=(\partial_\alpha \partial_\beta A^{\mu} \color{red}{- \Gamma^{\nu}_{\:\beta\alpha} \partial_\nu A^\mu} ) \\
&\quad + ( \partial_\alpha\Gamma^{\mu}_{\:\lambda \beta}
\textcolor{red}{ - \Gamma^{\nu}_{\:\beta\alpha} \Gamma^{\mu}_{\:\lambda\nu} }
+ \Gamma^{\mu}_{\:\nu\alpha} \Gamma^{\nu}_{\:\lambda \beta}
) A^{\lambda} \\
&\quad + (\Gamma^{\mu}_{\:\lambda \beta} \partial_\alpha A^{\lambda} + \Gamma^{\mu}_{\:\lambda \alpha} \partial_\beta A^{\lambda})
\end{aligned}
が成り立つ.ここで,今求めたいものは
\begin{aligned}
\boldsymbol{\nabla}_\alpha \boldsymbol{\nabla}^\alpha \boldsymbol{A}
&= g^{\alpha \beta} A^{\mu}_{\:\: ;\beta\,;\alpha} \boldsymbol{e}_\mu
\end{aligned}
である.各項にわけて計算する.

【第1項】
関係式

\begin{aligned}
g^{\alpha \beta} \Gamma^{\nu}_{\:\beta\alpha}
&=-\frac{1}{\sqrt{-g}}\partial_{\alpha} (\sqrt{-g} g^{\nu\alpha})
\end{aligned}
(例えば場の古典論$\S 86$)を用いれば,
\begin{aligned}
&g^{\alpha \beta} (\partial_\alpha \partial_\beta A^{\mu} - \Gamma^{\nu}_{\:\beta\alpha} \partial_\nu A^\mu ) \\
&=\frac{1}{\sqrt{-g}} \partial_\alpha (\sqrt{-g} \partial^{\alpha} A^\mu)
\end{aligned}
となる.これは,$A^\mu$にラプラシアンを作用させた形になっている.

【第2項】
第1項と同様にして

\begin{aligned}
& g^{\alpha \beta} ( \partial_\alpha\Gamma^{\mu}_{\:\lambda \beta} - \Gamma^{\nu}_{\:\beta\alpha} \Gamma^{\mu}_{\:\lambda\nu} ) \\
&=\frac{1}{\sqrt{-g}} \partial_\alpha (\sqrt{-g} g^{\alpha \beta} \Gamma^{\mu}_{\:\lambda\beta})
\end{aligned}
だから,
\begin{aligned}
& g^{\alpha \beta} ( \partial_\alpha\Gamma^{\mu}_{\:\lambda \beta} - \Gamma^{\nu}_{\:\beta\alpha} \Gamma^{\mu}_{\:\lambda\nu} )A^{\lambda} \\
&=\frac{1}{\sqrt{-g}} \partial_\alpha (\sqrt{-g} g^{\alpha \beta} \Gamma^{\mu}_{\:\lambda\beta} A^{\lambda})
- g^{\alpha \beta} \Gamma^{\mu}_{\:\lambda\beta} \partial_\alpha A^{\lambda}
\end{aligned}
となる.よって,
\begin{aligned}
& g^{\alpha \beta}
( \partial_\alpha\Gamma^{\mu}_{\:\lambda \beta}
- \Gamma^{\nu}_{\:\beta\alpha} \Gamma^{\mu}_{\:\lambda\nu}
+ \Gamma^{\mu}_{\:\nu\alpha} \Gamma^{\nu}_{\:\lambda \beta}
) A^{\lambda} \\
&= \frac{1}{\sqrt{-g}} \partial_\alpha (\sqrt{-g} g^{\alpha \beta} \Gamma^{\mu}_{\:\lambda\beta} A^{\lambda})
- g^{\alpha \beta} \Gamma^{\mu}_{\:\lambda\beta} \partial_\alpha A^{\lambda} \\
&\qquad +g^{\alpha \beta} \Gamma^{\mu}_{\:\nu\alpha} \Gamma^{\nu}_{\:\lambda \beta} A^{\lambda}
\end{aligned}

【第3項】

\begin{aligned}
&g^{\alpha \beta} (\Gamma^{\mu}_{\:\lambda \beta} \partial_\alpha A^{\lambda} + \Gamma^{\mu}_{\:\lambda \alpha} \partial_\beta A^{\lambda}) \\
&=2g^{\alpha \beta} \Gamma^{\mu}_{\:\lambda \beta} \partial_\alpha A^{\lambda}
\end{aligned}
となる.

以上で示された.//

円筒座標

\begin{aligned}
(\boldsymbol{\nabla}^2 \boldsymbol{A})_r
&=\boldsymbol{\nabla}^2 A_r
-\frac{1}{r^2}A_r
-\frac{2}{r^2}\frac{\partial A_\theta}{\partial \theta} \\
(\boldsymbol{\nabla}^2 \boldsymbol{A})_\theta
&=\boldsymbol{\nabla}^2 A_\theta
-\frac{1}{r^2}A_\theta
+\frac{2}{r^2}\frac{\partial A_\theta}{\partial \theta} \\
(\boldsymbol{\nabla}^2 \boldsymbol{A})_z
&=\boldsymbol{\nabla}^2 A_z
\end{aligned}

【証明】
円筒座標系のクリストッフェル記号でゼロでないものは

\begin{aligned}
\Gamma^{r}_{\:\theta\theta} =-r,\quad
\Gamma^{\theta}_{\:\theta r} = \Gamma^{\theta}_{\:r \theta} =\frac{1}{r}
\end{aligned}
だけである.(参考:クリストッフェル記号 - Notes_JP).

また,$A^r, A^\theta, A^z$でテンソル解析での成分を表し,$A_r, A_\theta, A_z$で(共変ベクトルの成分ではなく)ベクトル解析の成分を表すとすると

\begin{aligned}
A^r &= A_r/\sqrt{g_{rr}} = A_r \\
A^{\theta} &= A_{\theta}/\sqrt{g_{\theta\theta}} = A_{\theta}/r \\
A^{z} &= A_{z}/\sqrt{g_{zz}} = A_{z}
\end{aligned}
である(参考:「テンソル記法」から「ベクトル解析の記法」への変換方法 - Notes_JP).

【$r$成分】

\begin{aligned}
& (\boldsymbol{\nabla}^2 \boldsymbol{A})_r = (\boldsymbol{\nabla}^2 \boldsymbol{A})^r \\
&= \frac{1}{\sqrt{-g}} \partial_\alpha (\sqrt{-g} g^{\alpha \beta} \partial_{\beta} A^r) \\
&\qquad + \underline{\frac{1}{\sqrt{-g}} \partial_\theta (\sqrt{-g} g^{\theta\theta} \Gamma^{r}_{\:\theta\theta} A^{\theta})} \\
&\qquad +g^{\theta\theta} \Gamma^{r}_{\:\theta\theta}
(\underline{\partial_\theta A^{\theta}} + \Gamma^{\theta}_{\:r \theta} A^{r} ) \\
&=\boldsymbol{\nabla}^2 A_r \\
&\qquad + 2\cdot \frac{1}{r^2} (-r) \partial_\theta \biggl(\frac{ A_\theta}{r}\biggr) \\
&\qquad +\frac{1}{r^2} (-r) \frac{1}{r} A_{r} \\
&=\boldsymbol{\nabla}^2 A_r
-\frac{1}{r^2} A_{r}
-\frac{2}{r^2} \partial_\theta A_\theta
\end{aligned}
となる.途中で,$\sqrt{-g}$, $g^{\theta\theta}$, $\Gamma^{r}_{\:\theta\theta}$が$\theta$依存しないことから,下線部が等しくなることに注意.

【$\theta$成分】

【$z$成分】
$\Gamma^{z}_{\:\mu\nu}=0\:(\forall \mu,\nu)$より$(\boldsymbol{\nabla}^2 \boldsymbol{A})_z=\boldsymbol{\nabla}^2 A_z$.

以上で,導出が完了した.//

極座標

\begin{aligned}
(\boldsymbol{\nabla}^2 \boldsymbol{A})_r
&=\boldsymbol{\nabla}^2 A_r
-\frac{2}{r^2}A_r
- \frac{2}{r^2 \sin\theta}\frac{\partial}{\partial\theta}(A_\theta \sin\theta)
- \frac{2}{r^2 \sin\theta}\frac{\partial A_\varphi}{\partial \varphi} \\
(\boldsymbol{\nabla}^2 \boldsymbol{A})_\theta
&=\boldsymbol{\nabla}^2 A_\theta
-\frac{1}{r^2 \sin^2\theta} A_\theta
+\frac{2}{r^2} \frac{\partial A_r}{\partial\theta}
-\frac{2\cos\theta}{r^2 \sin^2\theta}\frac{\partial A_\varphi}{\partial \varphi} \\
(\boldsymbol{\nabla}^2 \boldsymbol{A})_\varphi
&=\boldsymbol{\nabla}^2 A_\varphi
-\frac{1}{r^2 \sin^2\theta} A_\varphi
+\frac{2}{r^2 \sin\theta}\frac{\partial A_r}{\partial \varphi}
+\frac{2\cos\theta}{r^2 \sin^2\theta}\frac{\partial A_\theta}{\partial \varphi}
\end{aligned}

【証明】
//

参考文献/記事