マッカラーの公式

マッカラーの公式（マッカラーのこうしき、仏: formule de MacCullagh）は、球状でない物体が、物体から十分に離れてはいるが、物体の径と比べれば極端には離れていない点で作る万有引力の場を記述するのに重要な公式である。万有引力ポテンシャルの基本的な定義が基になっているが、ポテンシャルの一般的なラプラス展開をするのではなく、2次近似までで切ることで得られる。測地学の理論や地球物理学で広く用いられる。

万有引力ポテンシャルの一般表示

球面座標系で、ポテンシャルが作られる点 $P$ （位置 $\mathbf {r}$ ）、ポテンシャルを作る点 $P'$ （位置 $\mathbf {r'}$ ）の動径成分をそれぞれ $r$ 、 $r'$ とする。座標の原点 $O$ は当面は任意としておく。ベクトル $\mathbf {r}$ と $\mathbf {r'}$ が成す角を $\psi$ 、 $P'$ での体積要素を $\mathrm {d} M(\mathbf {r'} )$ 、物体Bの全質量を $M$ と書く。このとき、物体Bによって作られる万有引力ポテンシャルは一般的に

V(P)=G\int _{M}{\frac {1}{\|\mathbf {r'} -\mathbf {r} \|}}\,\mathrm {d} M(\mathbf {r'} )

と書ける^[1]。ここで

\|\mathbf {r'} -\mathbf {r} \|=r{\sqrt {1-2{\frac {r'}{r}}\cos \psi +\left({\frac {r'}{r}}\right)^{2}}}

で、 $G$ はニュートンの万有引力定数である。

$-1\leq x\leq +1$ のとき $(1+x)^{-1/2}=1-x/2+3x^{2}/8+o(x^{3})$ だから^[2]、 $x=-2(r'/r)\cos \psi +(r'/r)^{2}$ とおくと

{\frac {1}{\sqrt {1-2{\frac {r'}{r}}\cos \psi +\left({\frac {r'}{r}}\right)^{2}}}}=1+{\frac {r'}{r}}\cos \psi -{\frac {1}{2}}\left({\frac {r'}{r}}\right)^{2}+{\frac {3}{2}}\left({\frac {r'}{r}}\right)^{2}\cos ^{2}\psi +o\left[\left({\frac {r'}{r}}\right)^{3}\right]

よって、

{\frac {r}{\|\mathbf {r'} -\mathbf {r} \|}}=1+{\frac {r'}{r}}\cos \psi +\left({\frac {r'}{r}}\right)^{2}-{\frac {3}{2}}\left({\frac {r'}{r}}\right)^{2}\sin ^{2}\psi +\dots

が得られる^[3] 。これより $P$ での万有引力ポテンシャルは次のようになる。

V(P)={\frac {G}{r}}\int _{M}\mathrm {d} M+{\frac {G}{r^{2}}}\int _{M}r'\cos \psi \mathrm {d} M+{\frac {G}{r^{3}}}\int _{M}r'^{2}\mathrm {d} M-{\frac {3G}{2r^{3}}}\int _{M}r'^{2}\sin ^{2}\psi \mathrm {d} M+\dots

マッカラーの定理

このように万有引力ポテンシャルの2次までの近似 $V(P)$ は、 $V_{0}(P)$ , $V_{1}(P)$ , $V_{2}(P)$ , $V_{3}(P)$ の4つの項に分解できることがわかる。つまり

V(P)=V_{0}(P)+V_{1}(P)+V_{2}(P)+V_{3}(P).

第1項

V_{0}=Gr^{-1}\int _{M}\,\mathrm {d} M=GMr^{-1}

は球対称な質量分布のポテンシャルに対応する。これは物体の全質量が1点に集められている場合と同じである。減少の度合 1/r が他の項と比べ緩やかであるため、非常に遠距離ではこの項が支配的となる。これは惑星のケプラー軌道を確立するのに用いられた、古典的なニュートンポテンシャルである。単極子（monopôle）、または単極子的質量分布（distribution de masse monopolaire）と呼ばれる。

第2項

V_{1}=Gr^{-2}\int _{M}r'\cos \psi \,\mathrm {d} M

は双極子（dipôle）、または双極子的質量分布（distribution de masse dipolaire）に対応する。座標の原点を物体の質量中心に一致させると、この項は消える^[4]。実際、右の図を参照すると次のように書くことができる。

V_{1}=Gr^{-2}\int _{M}r'\cos \psi \,\mathrm {d} M=Gr^{-2}\int _{M}x\,\mathrm {d} M=GMr^{-2}x_{0}

ここで $x_{0}$ は物体の質量中心の Ox-軸成分。

第3項と第4項、 $V_{2}$ , $V_{3}$ は四極子（quadrupôle）、または四極子的質量分布（distribution de masse quadrupolaire）である。 $V_{2}$ の計算を続けると

{\begin{aligned}V_{2}&=Gr^{-3}\int _{M}r'^{2}\,\mathrm {d} M=Gr^{-3}\int _{M}(x^{2}+y^{2}+z^{2})\mathrm {d} M\\&={\frac {1}{2}}Gr^{-3}\int _{M}[(y^{2}+z^{2})+(x^{2}+z^{2})+(x^{2}+y^{2})]\,\mathrm {d} M={\frac {1}{2}}Gr^{-3}(A+B+C)\end{aligned}}

となる。ここで $A$ , $B$ , $C$ はそれぞれOx, Oy, Oz-軸のまわりの慣性モーメントを表す。慣性モーメントの平均 ${\bar {I}}$ を

{\bar {I}}={\frac {1}{3}}(A+B+C)

,

とおくと、

V_{2}={\frac {3}{2}}Gr^{-3}{\bar {I}}

となる。 $V_{3}$ に対しては、

V_{3}=-{\frac {3}{2}}Gr^{-3}\int _{M}r'^{2}\sin ^{2}\psi \,\mathrm {d} M=-{\frac {3}{2}}Gr^{-3}\int _{M}r'^{2}\cos ^{2}\alpha \,\mathrm {d} M=-{\frac {3}{2}}Gr^{-3}I

ここで $I$ は ${\overline {OP}}$ -軸のまわりの慣性モーメントである。これより、次のマッカラーの公式（マッカラーの定理とも）が得られる。

V(P)=GMr^{-1}-{\frac {3}{2}}Gr^{-3}(I-{\bar {I}})

この公式は、地球や他の惑星等の球状に近い物体に適用するのに非常に便利である。物体に何らかの対称性があれば、近似は十分遠距離でも有効である。明らかに、2次までの近似であることから有効な適用範囲には限界があり、外部重力場をさらに正確に表すためにはより高次の調和関数が必要になる。

対称性のある物体に対するマッカラーの公式

質量中心を通る軸 ${\overline {OP}}$ のまわりの慣性モーメントは、慣性テンソルを表す 3×3行列を対角化することで得られる主慣性モーメント $A$ , $B$ , $C$ を使って

I=n_{x}^{2}A+n_{y}^{2}B+n_{z}^{2}C

と書くことができる。ここで $n_{x}$ , $n_{y}$ , $n_{z}$ は ${\overline {OP}}$ （単位ベクトル）の慣性主軸 Ox, Oy, Oz への方向余弦を表す。

物体が Oz-軸を中心に回転する（Oz について軸対称である）と仮定する。このとき $A=B$ である。 $\phi$ を ${\overline {OP}}$ が Oxy-平面となす角とする^[5]。このとき $n_{z}=\sin \phi$ であり、 $n_{x}^{2}+n_{y}^{2}+n_{z}^{2}=1$ だから方向余弦は

\sin ^{2}\phi =1-n_{x}^{2}-n_{y}^{2}

である。 $A=B$ としているので、マッカラーの公式は

V(P)=GMr^{-1}-{\frac {1}{2}}Gr^{-3}(3I-2A-C)

となる。慣性モーメント I は

I=A(n_{x}^{2}+n_{y}^{2})+C\sin ^{2}\phi =A(1-\sin ^{2}\phi )+C\sin ^{2}\phi

となるから、回転体に対するマッカラーの公式は最終的に次の形になる ^[6]。

V(P)=GMr^{-1}+{\frac {1}{2}}Gr^{-3}(C-A)(1-3\sin ^{2}\phi )

地球の力学的形状係数

軸対称性があるとき、万有引力ポテンシャルの多重極展開を2次まで行ったものは常に次の一般式で書ける。

{\begin{aligned}V(P)&=GMr^{-1}-J_{2}GMa^{2}r^{-3}P_{2}(\cos \theta )\\&=GMr^{-1}-{\frac {1}{2}}J_{2}GMa^{2}r^{-3}(1-3\sin ^{2}\phi )\end{aligned}}

ここで

P_{2}(\cos \theta )={\frac {3}{2}}\sin ^{2}\phi -{\frac {1}{2}},\quad (\phi ={\frac {\pi }{2}}-\theta )

は2次のルジャンドル多項式である。地球の場合、2次の帯調和重力係数 $J_{2}$ はしばしば地球の力学的形状係数（facteur de forme géodynamique）と呼ばれる。この式をマッカラーの公式と等しいと置くと、

J_{2}={\frac {C-A}{Ma^{2}}}

この最後の関係式もマッカラーの公式と呼ばれることがある。

注釈

^ ここでは $V$ を仕事として定義しているが、物理学（の他分野）では、万有引力（重力）ポテンシャルは $\Phi =-V$ を満たす相互作用エネルギー $\Phi$ と定義されるのが慣習である。
^ $\textstyle {\frac {1}{\sqrt {1-2{\frac {r'}{r}}\cos \psi +\left({\frac {r'}{r}}\right)^{2}}}}=1-{\frac {1}{2}}\left[-2{\frac {r'}{r}}\cos \psi +\left({\frac {r'}{r}}\right)^{2}\right]+{\frac {3}{8}}\left[-2{\frac {r'}{r}}\cos \psi +\left({\frac {r'}{r}}\right)^{2}\right]^{2}+\dots$ $\textstyle =1+{\frac {r'}{r}}\cos \psi -{\frac {1}{2}}\left({\frac {r'}{r}}\right)^{2}+{\frac {3}{8}}\left[4\left({\frac {r'}{r}}\right)^{2}\cos ^{2}\psi -4\left({\frac {r'}{r}}\right)^{3}\cos \psi +\left({\frac {r'}{r}}\right)^{4}\right]+\dots$ $\textstyle =1+{\frac {r'}{r}}\cos \psi -{\frac {1}{2}}\left({\frac {r'}{r}}\right)^{2}+{\frac {3}{2}}\left({\frac {r'}{r}}\right)^{2}\cos ^{2}\psi +o\left[\left({\frac {r'}{r}}\right)^{3}\right]$
^ $\textstyle {\frac {r}{\|\mathbf {r'} -\mathbf {r} \|}}=1+{\frac {r'}{r}}\cos \psi +{\frac {1}{2}}\left({\frac {r'}{r}}\right)^{2}(3\cos ^{2}\psi -1)+\dots =1+{\frac {r'}{r}}\cos \psi +{\frac {1}{2}}\left({\frac {r'}{r}}\right)^{2}(3-3\sin ^{2}\psi -1)+\dots$ $\textstyle =1+{\frac {r'}{r}}\cos \psi +\left({\frac {r'}{r}}\right)^{2}-{\frac {3}{2}}\left({\frac {r'}{r}}\right)^{2}\sin ^{2}\psi +\dots$ .
^ このように、（近似を続けつつ）万有引力ポテンシャルの双極子項が消えるように座標軸の取り方を決めることで、地球の質量中心を正確に決定できる。
^ 地球物理学の文脈では、 $\phi$ は一般に地点の緯度（もしくは緯度にほぼ一致する角）を表す。
^ $V(P)=GMr^{-1}-{\frac {1}{2}}Gr^{-3}(3I-2A-C)$ $=GMr^{-1}-{\frac {1}{2}}Gr^{-3}(3A-3A\sin ^{2}\phi +3C\sin ^{2}\phi -2A-C)$ $=GMr^{-1}-{\frac {1}{2}}Gr^{-3}[A-C-3(A-C)\sin ^{2}\phi ]$ $=GMr^{-1}+{\frac {1}{2}}Gr^{-3}(C-A)(1-3\sin ^{2}\phi )$