============================================================ クリフォード代数とマクスウェル方程式 ============================================================ 2020.03.29追記 うべさんがブログを作りました。ご興味ある方は こちらから_ どうぞ！ 最近、Twitterで知り合った高専生のうべゆうとさんという方に 真空中のマクスウェル方程式の面白い表現を教えていただいたので、 それについて述べてみようと思います。 クリフォード代数というものを使います。 実に興味深いです。結論から言うと、マクスウェル方程式は たった一本の方程式として表現でき、 \dfrac{b_u}{2} \left( b_u A - A d_u \right) + J_e = 0 \tag{##} となります、ここで、 $A$ と $J_e$ は4元ポテンシャルと4元電流密度 A &= \left( \phi, A_1, A_2, A_3 \right) \\ J_e &= (\rho, j_1, j_2,j_3) \tag{##} であり、 $b_u,d_u$ は微分演算子の一種で、微分形式で使う外微分 $d$ と同じ記号なので、 この記事では僕がうべさんの"u"をお借りして $b_u,d_u$ と書くことにします。 これをうべさんはそれぞれ右純然微分作用素、左純然微分作用素と名付けました。 クリフォード代数の基本 ========================= クリフォード代数はその基底の二乗の値に $-1$ を許容する代数です。 複素数の拡張と考えられます。今回我々が使うのは、 四つの基底 $e_0, e_1, e_2, e_3$ に対して、 e_0^2 &= -1 \\ e_1^2 &= e_2^2 = e_3^2 = 1 \\ e_i e_j &= - e_j e_i \tag{##} を満たすものです。クリフォード代数では $e_i$ の二乗は実数にグレードダウンできますが、 最後の式から、 $e_0 e_2, e_2 e_3, e_0 e_3 e_1$ 等は消えずにその形で残ります。 また、この代数は非可換です。 そして、クリフォード代数の商を例えば、次の様に考えます。実数 $a,b,c,d$ に対して、 \dfrac{1}{ a e_0 + b e_1 }(c e_0 + d e_1) &= \dfrac{a e_0 + b e_1}{ (a e_0 + b e_1)^2 }(c e_0 + d e_1) \\ &= \dfrac{ace_0e_0 +ad e_0e_1 +bc e_1 e_0 + bd e_1 e_1}{a^2e_0e_0 +ab e_0e_1 +ab e_1 e_0 + b^2 e_1 e_1} \\ &= \dfrac{1}{-a^2+b^2}\left( (-ac+bd)+(ad - bc)e_0e_1 \right) \tag{##} となります。逆数を右から掛けるか左から掛けるかで結果が変わることにご注意ください。 すると、この代数は逆元を $(i = 0,1,2,3)$ に対して $\dfrac{1}{e_i}e_i = 1$ で定義する、つまり、 \dfrac{1}{e_i} = \begin{cases} -e_0 \ \ \ \ (i=0) \\ e_i \ \ \ \ (i=1,2,3) \end{cases} \tag{##} とすれば良いです。 うべさんの発見 ====================== 僕がうべさんから教えていただいた時には、発見論的な話で聞いたのですが、 紆余曲折は本質を見づらくすると思うので、結論に直接持っていく話をしようと思います。 偏微分微分演算子の記号 $\partial$ を右に作用するか左に作用するかで $p,q$ とします。 つまり、 A \dfrac{p}{px_i} B &= A \dfrac{\partial B}{\partial x_i} \\ A \dfrac{q}{x_iq} B &= \dfrac{\partial A}{\partial x_i} B \tag{##} とします。余談ですが、右と左で対称的な作用をすることから、 左右で対称な記号の対として $p,q,b,d$ を 選んだそうです。 $p$ は、partial differential の $p$ だったり、 彼の尊敬するDiracから $d$ を選んだともおっしゃっていました。 こういう遊び心は良いですね。 そして、いよいよこの記事の冒頭で書いた記号 $b_u,d_u$ を定義します。 その定義は b_u = \sum_{i=0}^3 \dfrac{1}{e_i} \dfrac{p}{px_i} \\ d_u = \sum_{i=0}^3 \dfrac{q}{x_iq} \dfrac{1}{e_i} \tag{##} とします。 マクスウェル方程式のクリフォード代数での記述 ===================================================== 自然単位系、つまり $\varepsilon_0 = \mu_0 = \hbar = 1$ とした(今回は $\hbar$ は使いませんが) 真空中のマクスウェル方程式を書くと &\nabla \cdot \bm{E} = \rho \\ &\nabla \times \bm{B} - \dfrac{\partial \bm{E}}{\partial t} = \bm{j} \\ &\nabla \cdot \bm{B} = 0 \\ &\nabla \times \bm{E} + \dfrac{\partial \bm{B}}{\partial t} = 0 \tag{##} となります。ポテンシャルを $A = (\phi, A_1, A_2, A_3) = (\phi, \bm{A})$ とすれば、 \bm{E} &= - \nabla \phi -\dfrac{\partial \bm{A}}{\partial t} \\ \bm{B} &= \nabla \times \bm{A} \tag{##} となります。ここで、 $b_u A$ を計算してみます。すると、 b_u A &= \left( \sum_{i=0}^3 \dfrac{1}{e_i} \dfrac{p}{px_i} \right) (\phi e_0 + A_1 e_1 + A_2 e_2 + A_3 e_3) \\ &= \left( - \dfrac{p}{px_0}e_0 + \dfrac{p}{px_1}e_1 + \dfrac{p}{px_2}e_2 + \dfrac{p}{px_3}e_3 \right) (\phi e_0 + A_1 e_1 + A_2 e_2 + A_3 e_3) \\ &= \left( \dfrac{\partial \phi}{\partial t} - \dfrac{\partial A_1}{\partial t} e_0 e_1 - \dfrac{\partial A_2}{\partial t} e_0 e_2 - \dfrac{\partial A_3}{\partial t} e_0 e_3 \right. \\ &- \dfrac{\partial \phi}{\partial x} e_0 e_1 + \dfrac{\partial A_1}{\partial x} + \dfrac{\partial A_2}{\partial x} e_1 e_2 - \dfrac{\partial A_3}{\partial x} e_3 e_1 \\ &- \dfrac{\partial \phi}{\partial y} e_0 e_2 - \dfrac{\partial A_1}{\partial y} e_1 e_2 + \dfrac{\partial A_2}{\partial y} + \dfrac{\partial A_3}{\partial y} e_2 e_3 \\ &- \left. \dfrac{\partial \phi}{\partial z} e_0 e_3 + \dfrac{\partial A_1}{\partial z} e_3 e_1 - \dfrac{\partial A_2}{\partial z} e_2 e_3 + \dfrac{\partial A_3}{\partial z} \right) \\ &= \left( \dfrac{\partial \phi}{\partial t} + \dfrac{\partial A_1}{\partial x} + \dfrac{\partial A_2}{\partial y} + \dfrac{\partial A_3}{\partial z} \right)+ E_1 e_0 e_1 + E_2 e_0 e_2 + E_3 e_0 e_3 + B_1 e_2 e_3 + B_2 e_3 e_1 + B_3 e_1 e_2 \tag{##} 最後の行に行くには、式 $(9)$ を使いました。 同様に $A d_u$ を計算すると、 A d_u &= (\phi e_0 + A_1 e_1 + A_2 e_2 + A_3 e_3) \left( \sum_{i=0}^3 \dfrac{q}{x_iq} \dfrac{1}{e_i} \right)\\ &= \left( \dfrac{\partial \phi}{\partial t} + \dfrac{\partial A_1}{\partial x} + \dfrac{\partial A_2}{\partial y} + \dfrac{\partial A_3}{\partial z} \right) -(E_1 e_0 e_1 + E_2 e_0 e_2 + E_3 e_0 e_3 + B_1 e_2 e_3 + B_2 e_3 e_1 + B_3 e_1 e_2) \tag{##} となります。マクスウェル方程式の登場人物が出てきましたね。 $\bm{E},\bm{B}$ です。さらにもう一歩進んで、 \dfrac{1}{2} (b_u A - A d_u) &= E_1 e_0 e_1 + E_2 e_0 e_2 + E_3 e_0 e_3 + B_1 e_2 e_3 + B_2 e_3 e_1 + B_3 e_1 e_2 ( \equiv V ) \tag{##} とします。( $\equiv$ は $V$ をこれで定義するという意味です。) ここで、左からさらに $b_u$ を掛けると、 b_u V &= \left( \sum_{i=0}^3 \dfrac{1}{e_i} \dfrac{p}{px_i} \right) \left( E_1 e_0 e_1 + E_2 e_0 e_2 + E_3 e_0 e_3 + B_1 e_2 e_3 + B_2 e_3 e_1 + B_3 e_1 e_2 \right) \\ &= \dfrac{\partial E_1}{\partial t} e_1 + \dfrac{\partial E_2}{\partial t} e_2 + \dfrac{\partial E_3}{\partial t} e_3 - \dfrac{\partial B_1}{\partial t} e_0 e_2 e_3 - \dfrac{\partial B_2}{\partial t} e_0 e_3 e_1 - \dfrac{\partial B_3}{\partial t} e_0 e_1 e_2 \\ &- \dfrac{\partial E_1}{\partial x} e_0 - \dfrac{\partial E_2}{\partial x} e_0 e_1 e_2 + \dfrac{\partial E_3}{\partial x} e_0 e_3 e_1 + \dfrac{\partial B_1}{\partial x} e_1 e_2 e_3 - \dfrac{\partial B_2}{\partial x} e_3 + \dfrac{\partial B_3}{\partial x} e_2 \\ &+ \dfrac{\partial E_1}{\partial y} e_0 e_1 e_2 - \dfrac{\partial E_2}{\partial y} e_0 - \dfrac{\partial E_3}{\partial y} e_0 e_2 e_3 + \dfrac{\partial B_1}{\partial y} e_3 + \dfrac{\partial B_2}{\partial y} e_1 e_2 e_3 - \dfrac{\partial B_3}{\partial y} e_1 \\ &- \dfrac{\partial E_1}{\partial z} e_0 e_3 e_1 + \dfrac{\partial E_2}{\partial z} e_0 e_2 e_3 - \dfrac{\partial E_3}{\partial z} e_0 - \dfrac{\partial B_1}{\partial z} e_2 + \dfrac{\partial B_2}{\partial z} e_1 + \dfrac{\partial B_3}{\partial z} e_1 e_2 e_3 \\ &= - \left( \dfrac{\partial E_1}{\partial x} + \dfrac{\partial E_2}{\partial y} + \dfrac{\partial E_3}{\partial z} \right) e_0 \\ &+ \left( \dfrac{\partial E_1}{\partial t} - \dfrac{\partial B_3}{\partial y} + \dfrac{\partial B_2}{\partial z} \right) e_1 \\ &+ \left( \dfrac{\partial E_2}{\partial t} - \dfrac{\partial B_1}{\partial z} + \dfrac{\partial B_3}{\partial x} \right) e_2 \\ &+ \left( \dfrac{\partial E_3}{\partial t} - \dfrac{\partial B_2}{\partial x} + \dfrac{\partial B_1}{\partial y} \right) e_3 \\ &+ \left( \dfrac{\partial B_1}{\partial x} + \dfrac{\partial B_2}{\partial y} + \dfrac{\partial B_3}{\partial z} \right) e_1 e_2 e_3 \\ &+ \left( - \dfrac{\partial B_1}{\partial t} - \dfrac{\partial E_3}{\partial y} + \dfrac{\partial E_2}{\partial z} \right) e_0 e_2 e_3 \\ &+ \left( - \dfrac{\partial B_2}{\partial t} - \dfrac{\partial E_1}{\partial z} + \dfrac{\partial E_3}{\partial x} \right) e_0 e_3 e_1 \\ &+ \left( - \dfrac{\partial B_3}{\partial t} - \dfrac{\partial E_2}{\partial x} + \dfrac{\partial E_1}{\partial y} \right) e_0 e_1 e_2 \tag{##} これを b_u V &= - \mathrm{div} \bm{E} e_0 \\ &+ \left( \dfrac{\partial \bm{E}}{\partial t} - \nabla \times \bm{B} \right) e_i \\ &+ \mathrm{div} \bm{B} e_1 e_2 e_3 \\ &+ \left( -\dfrac{\partial \bm{B}}{\partial t} - \nabla \times \bm{E} \right) e_0 e_i e_j \tag{##} と書いても誤解は生じないでしょう。 $e_0 e_i e_j$ は上から順に $e_0 e_2 e_3, e_0 e_3 e_1, e_0 e_1 e_2$ を表します。 すると、これはマクスウェル方程式の左辺になっています。つまり、後は電流密度を $J_e$ として、 J_e = \rho e_0 + j_1 e_1 + j_2 e_2 + j_3 e_3 \tag{##} とすれば、式 $(8)$ の事を b_u V + J_e = 0 \tag{##} と書けるのです。 また、もし磁気モノポールが存在するなら、 J_m = \rho_m e_1 e_2 e_3 + j_{m1} e_0 e_2 e_3 + j_{m2} e_0 e_3 e_1 + j_{m3} e_0 e_1 e_2 \tag{##} によって、マクスウェル方程式を、 b_u V = - J_e + J_m \tag{##} と書けます。右辺の符号が気になるなら、 V d_u = J_e + J_m \tag{##} としても良いです。 最後にうべさんの上げた成果について、述べておきます。 この記事は、うべさんのレビューを基にクロメルが作ったレビューであり、 新規性はどこにあるかと言うと、実をいうと古くから知られていたことを表現しなおしたことの域を得ません。 また、左右の純然微分作用素に関しては、ウェッジ積を用いた同様のものは以前から存在しているらしいです。 しかし、分母に $e_i$ を持ってきたことは、うべさんが初めてではないかと思われます。 もし、この辺りの事情に詳しい方がいらっしゃったら、 クロメル_ までメールでご一報ください。 うべさんは現在、現役の高専生ということで（若い！）、将来がとても楽しみです。 うべさん、面白い話をありがとうございました。今日はここまで、お疲れさまでした。 .. _こちらから: https://ubeyuto.hatenablog.com/ .. _クロメル: http://hooktail.org/wiki/index.php?%A5%AF%A5%ED%A5%E1%A5%EB @@author:クロメル@@ @@accept:2019-06-09@@ @@category:電磁気学@@ @@id:cliffordMaxwell@@