古典力学 において、重力ポテンシャル と は、空間の 各点に、保存重力場 における固定された基準点から物体をその点に移動させるために必要な単位質量あたりの仕事(伝達されるエネルギー)を関連付けるスカラーポテンシャルである。これ は 、質量 が電荷 の役割を果たす電位 に類似している 。ポテンシャルがゼロとなる基準点は、慣例的にあらゆる質量から無限に離れているため、 有限 距離において負のポテンシャルとなる。これらの類似性は、両方の関連場が 保存力を 持つことと相関している。
数学的には、重力ポテンシャルはニュートンポテンシャルとも呼ばれ、 ポテンシャル理論 の研究において基礎となる。また、均一に帯電または分極した楕円体によって生成される静電場および静磁場を解く際にも用いられる。[ 1 ]
位置エネルギー ある場所における重力ポテンシャル(V )は、その場所における単位質量あたりの 重力ポテンシャルエネルギー (U )です。
V = あなた メートル 、 {\displaystyle V={\frac {U}{m}},}
ここで、m は物体の質量です。位置エネルギーは、重力場が物体を無限遠から空間内の所定の位置に移動させる仕事に等しく(大きさは負ですが)、物体の質量が1キログラムの場合、その物体に割り当てられる位置エネルギーは重力位置エネルギーに等しくなります。したがって、位置エネルギーは、重力場が単位質量を無限遠から移動させる仕事の負の値として解釈できます。
状況によっては、位置からほぼ独立した場を仮定することで方程式を簡略化できます。例えば、地球の表面に近い領域では、重力加速度 g は一定とみなすことができます。この場合、ある高さから別の高さへの位置エネルギーの差は、ほぼ正確には、高さの差に比例します。 Δ あなた ≈ メートル グラム Δ h 。 {\displaystyle \Delta U\approx mg\Delta h.}
質量Mの 質点 から距離x にある重力ポテンシャル Vは 、単位質量を無限遠からその点まで運ぶために外部のエージェントが行う必要がある仕事Wとして定義できます。 [ 2 ] [ 3 ] [ 4 ] [ 5 ]
V ( × ) = W メートル = 1 メートル ∫ ∞ × F ( × ′ ) ⋅ d × ′ = 1 メートル ∫ ∞ × G メートル M × ′ 2 d × ′ = − G M × 、 {\displaystyle V(\mathbf {x} )={\frac {W}{m}}={\frac {1}{m}}\int _{\infty }^{x}\mathbf {F} \left(\mathbf {x} '\right)\cdot d\mathbf {x} '={\frac {1}{m}}\int _{\infty }^{x}{\frac {GmM}{x'^{2}}}dx'=-{\frac {GM}{x}},} ここで、G は重力定数 、F は重力の力です。積GMは 標準的な重力パラメータ であり、 G またはMを 個別に用いるよりも高い精度で知られていることがよくあります。ポテンシャルは質量あたりのエネルギーを単位 とし、例えばMKS 単位系ではJ/kgとなります。慣例により、定義されている場所では常に負の値をとり、xが 無限大に近づくにつれて0に近づきます。
重力場 、ひいては質量の大きい物体の周囲の空間における小物体の加速度は、重力ポテンシャルの 負の勾配 である。したがって、負の勾配の負の値は、質量の大きい物体に向かう正の加速度を生み出す。ポテンシャルには角度成分がないため、その勾配は次のように表される。 ここで、 x は質量点から小物体に向かう長さx のベクトルであり、は質量点から小物体に向かう単位ベクトル である。したがって、加速度の大きさは反二乗則 に従う。 1つの = − G M × 3 × = − G M × 2 × ^ 、 {\displaystyle \mathbf {a} =-{\frac {GM}{x^{3}}}\mathbf {x} =-{\frac {GM}{x^{2}}}{\hat {\mathbf {x} }},} × ^ {\displaystyle {\hat {\mathbf {x} }}} ‖ 1つの ‖ = G M × 2 。 {\displaystyle \|\mathbf {a} \|={\frac {GM}{x^{2}}}.}
質量分布 に関連するポテンシャルは、質点のポテンシャルの重ね合わせである。質量分布が有限個の質点の集合であり、質点が点x 1 , ..., x n に位置し、質量がm 1 , ..., m n である場合、点x における分布のポテンシャルは V ( × ) = ∑ 私 = 1 n − G メートル 私 ‖ × − × 私 ‖ 。 {\displaystyle V(\mathbf {x} )=\sum _{i=1}^{n}-{\frac {Gm_{i}}{\|\mathbf {x} -\mathbf {x} _{i}\|}}.}
点x とr 。r は分布質量(灰色)に含まれ、微分質量dm ( r ) は点r に位置します。 質量分布が3次元ユークリッド空間 R 3上の質量 測度 dm として与えられている場合、ポテンシャルは− G /| r | とdm の畳み込み です。多くの場合、これは積分に等しくなります。 ここで| x − r | は点x とr間の 距離 です。rにおける分布の密度を表す関数ρ ( r ) が あり、dm ( r ) = ρ ( r ) dv ( r ) ( dv ( r ) はユークリッド体積要素 )である場合、重力ポテンシャルは体積積分です。 V ( × ) = − ∫ R 3 G ‖ × − r ‖ d メートル ( r ) 、 {\displaystyle V(\mathbf {x} )=-\int _{\mathbb {R} ^{3}}{\frac {G}{\|\mathbf {x} -\mathbf {r} \|}}\,dm(\mathbf {r} ),} V ( × ) = − ∫ R 3 G ‖ × − r ‖ ρ ( r ) d v ( r ) 。 {\displaystyle V(\mathbf {x} )=-\int _{\mathbb {R} ^{3}}{\frac {G}{\|\mathbf {x} -\mathbf {r} \|}}\,\rho (\mathbf {r} )dv(\mathbf {r} ).}
V が連続質量分布 ρ ( r )から導かれるポテンシャル関数である場合、ρ は ラプラス演算子 Δ を用いて復元できます。 これは、ρ が連続で、有界集合の外側ではゼロである場合には常に点ごとに成り立ちます。一般 に、ラプラス演算子を分布の意味でとれば、質量測度 dm も 同様の方法で復元できます。結果として、重力ポテンシャルはポアソン方程式 を満たします。3変数ラプラス方程式のグリーン関数 とニュートンポテンシャル も参照してください。 ρ ( × ) = 1 4 π G Δ V ( × ) 。 {\displaystyle \rho (\mathbf {x} )={\frac {1}{4\pi G}}\Delta V(\mathbf {x} ).}
対称楕円体や退化楕円体を含むあらゆる楕円体形状について、この積分は既知の超越関数を用いて表すことができます。[ 6 ] これらには、3つの半軸が等しい球面、2つの半軸が等しい扁平楕円体(楕円体参照 )と長楕円体、1つの半軸が無限大の退化楕円体(楕円柱と円筒)、そして2つの半軸が無限大の無限シートが含まれます。これらの形状はすべて、重力ポテンシャル積分(定数G 、ただし定数電荷密度 𝜌は除く)の電磁気学への応用において広く用いられています。
球面対称性 球対称の質量分布は、分布の外側にいる観測者にとっては、質量の全てが中心に集中しているかのように振る舞い、殻定理 により実質的に質点 として振る舞います。地球表面では、加速度はいわゆる標準重力 g によって与えられ、約9.8 m/s 2ですが、この値は緯度と高度によってわずかに変化します。地球は 扁平回転楕円体 であるため、加速度の大きさは赤道よりも極でわずかに大きくなります。
球対称の質量分布の下では、球座標系でポアソン方程式 を解くことができる。半径R 、密度ρ、質量m の均一球体内部の重力g は中心からの距離r に比例して変化し、球体内部の重力ポテンシャルは[ 7 ] [ 8 ] となり、 これは球体外部のポテンシャル関数と微分的に結びつく(上図参照)。 V ( r ) = 2 3 π G ρ [ r 2 − 3 R 2 ] = G メートル 2 R 3 [ r 2 − 3 R 2 ] 、 r ≤ R 、 {\displaystyle V(r)={\frac {2}{3}}\pi G\rho \left[r^{2}-3R^{2}\right]={\frac {Gm}{2R^{3}}}\left[r^{2}-3R^{2}\right],\qquad r\leq R,}
一般相対性理論 一般相対論 では、重力ポテンシャルは計量テンソル に置き換えられます。重力場が弱く、源が光速に比べて非常にゆっくりと移動している場合、一般相対論はニュートン力学の重力に帰着し、計量テンソルは重力ポテンシャルによって展開されます。[ 9 ]
多極展開 点x における電位は次のように与えられる。 V ( × ) = − ∫ R 3 G | × − r | d メートル ( r ) 。 {\displaystyle V(\mathbf {x} )=-\int _{\mathbb {R} ^{3}}{\frac {G}{|\mathbf {x} -\mathbf {r} |}}\ dm(\mathbf {r} ).}
質量中心をベクトルx とr の原点とし、ポテンシャルが計算される点をベクトルx の先頭に置いた質量分布 (灰色) の図。ポテンシャルは、一連のルジャンドル多項式 で展開できます。点x とr を、 質量中心 に対する位置ベクトル として表します。積分の分母は平方根として表され、次のようになります。 ここで、最後の積分では、r = | r | であり、θは x とr の間の角度です。 V ( × ) = − ∫ R 3 G | × | 2 − 2 × ⋅ r + | r | 2 d メートル ( r ) = − 1 | × | ∫ R 3 G 1 − 2 r | × | コス θ + ( r | × | ) 2 d メートル ( r ) {\displaystyle {\begin{aligned}V(\mathbf {x} )&=-\int _{\mathbb {R} ^{3}}{\frac {G}{\sqrt {|\mathbf {x} |^{2}-2\mathbf {x} \cdot \mathbf {r} +|\mathbf {r} |^{2}}}}\,dm(\mathbf {r} )\\&=-{\frac {1}{|\mathbf {x} |}}\int _{\mathbb {R} ^{3}}{\frac {G}{\sqrt {1-2{\frac {r}{|\mathbf {x} |}}\cos \theta +\left({\frac {r}{|\mathbf {x} |}}\right)^{2}}}}\,dm(\mathbf {r} )\end{aligned}}}
( § 数学的形式 を参照。) 積分関数は、係数を明示的に計算することにより、Z = r /| x | のテイラー級数として展開できます。同じ結果を得るためのより簡単な方法は、一般化 二項定理 を使用することです。[ 10 ] 結果として得られる級数は、ルジャンドル多項式の 生成関数です。これは 、 | X | ≤ 1 かつ| Z | < 1 に有効です。係数P n は、次数n のルジャンドル多項式です。したがって、積分関数のテイラー係数は、X = cos θのルジャンドル多項式で与えられます。したがって、システムのすべての質量要素に対して r < | x | となるような位置x について収束する級数にポテンシャルを展開できます(つまり、システムを囲む、質量中心を中心とする球の外側)。 積分は、質量中心のx 方向の成分です。ベクトルx は 質量中心から放射されるため、この成分は消えます。したがって、積分を和の符号の下に持ってくると、 ( 1 − 2 X Z + Z 2 ) − 1 2 = ∑ n = 0 ∞ Z n P n ( X ) {\displaystyle \left(1-2XZ+Z^{2}\right)^{-{\frac {1}{2}}}\ =\sum _{n=0}^{\infty }Z^{n}P_{n}(X)} V ( × ) = − G | × | ∫ ∑ n = 0 ∞ ( r | × | ) n P n ( コス θ ) d メートル ( r ) = − G | × | ∫ ( 1 + ( r | × | ) コス θ + ( r | × | ) 2 3 コス 2 θ − 1 2 + ⋯ ) d メートル ( r ) {\displaystyle {\begin{aligned}V(\mathbf {x} )&=-{\frac {G}{|\mathbf {x} |}}\int \sum _{n=0}^{\infty }\left({\frac {r}{|\mathbf {x} |}}\right)^{n}P_{n}(\cos \theta )\,dm(\mathbf {r} )\\&=-{\frac {G}{|\mathbf {x} |}}\int \left(1+\left({\frac {r}{|\mathbf {x} |}}\right)\cos \theta +\left({\frac {r}{|\mathbf {x} |}}\right)^{2}{\frac {3\cos ^{2}\theta -1}{2}}+\cdots \right)\,dm(\mathbf {r} )\end{aligned}}} ∫ r cos ( θ ) d m {\textstyle \int r\cos(\theta )\,dm} V ( x ) = − G M | x | − G | x | ∫ ( r | x | ) 2 3 cos 2 θ − 1 2 d m ( r ) + ⋯ {\displaystyle V(\mathbf {x} )=-{\frac {GM}{|\mathbf {x} |}}-{\frac {G}{|\mathbf {x} |}}\int \left({\frac {r}{|\mathbf {x} |}}\right)^{2}{\frac {3\cos ^{2}\theta -1}{2}}dm(\mathbf {r} )+\cdots }
これは、質量中心からの距離が同じ場合を比較すると、物体が伸長すると、球状の質量による電位と比較して、伸長方向の電位は低くなり、垂直方向の電位は高くなることを示しています。(表面からの距離が同じ場合を比較すると、 その逆になります。)
単位と数値 重力ポテンシャルのSI単位は平方メートル毎秒(m 2 /s 2 )、またはジュール毎キログラム(J/kg)です。地球 、 太陽、 そして天の川 銀河の 質量に対する 、様々 な地点における重力ポテンシャルの絶対値は 、以下の表に示されています。つまり、地球表面にある物体が地球の重力場から「脱出」するには60 MJ/kg 、太陽の重力場から「脱出」するにはさらに900 MJ/kg、天の川銀河の重力場から「脱出」するには130 GJ/kg以上の重力ポテンシャルが必要です。このポテンシャルは 脱出速度 の2乗の半分です。
位置 に関して 地球 太陽 天の川 地球の表面 60 MJ/kg 900 MJ/kg ≥ 130 GJ/kg レオ 57 MJ/kg 900 MJ/kg ≥ 130 GJ/kg ボイジャー1号 (地球から170億km)23 J/kg 8 MJ/kg ≥ 130 GJ/kg 地球から0.1光年 0.4 J/kg 140 kJ/kg ≥ 130 GJ/kg
これらの場所での重力 を比較してください。
参照
注記 ^ Solivérez, CE (2016).偏極楕円体の静電気学と静磁気学:偏極テンソル法 (第1版). 無料科学情報. ISBN 978-987-28304-0-3 。 ^ マリオン, JB; ソーントン, ST (1995). 粒子とシステムの古典力学 (第4版). ハーコート・ブレース・アンド・カンパニー. p . 192. ISBN 0-03-097302-3 。^ Arfken, George B.; Weber, Hans J. (2005). Mathematical Methods For Physicists International Student Edition (6th ed.). Academic Press . p. 72. ISBN 978-0-08-047069-6 。^ サング, デイビッド; ジョーンズ, グラハム; チャダ, グリンダ; ウッドサイド, リチャード; スターク, ウィル; ギル, エイダン (2014). ケンブリッジ国際AS・Aレベル物理学コースブック (イラスト入り). ケンブリッジ大学出版局 . p. 276. ISBN 978-1-107-69769-0 。^ マンカスター、ロジャー (1993). Aレベル物理学 (イラスト入り). ネルソン・ソーンズ . p. 106. ISBN 978-0-7487-1584-8 。^ マクミラン、WD (1958). 『潜在能力の理論 』ドーバー・プレス. ^ ローリー、ウィリアム・ローリー(2011年) 『地球物理学方程式の学生ガイド 』ケンブリッジ大学出版局、p.69、 ISBN 978-1-139-49924-8 。68ページの抜粋 ^ サンチェス=ラベガ、アグスティン (2011). 『惑星大気入門 (イラスト入り)』CRC Press. p. 19. ISBN 978-1-4200-6735-4 。19ページの抜粋 ^ Grøn, Øyvind; Hervik, Sigbjorn (2007)、 アインシュタインの一般相対性理論:宇宙論への現代的応用 、Springer Science & Business Media、p. 201、 ISBN 978-0-387-69200-5 ^ Wylie, CR Jr. (1960). Advanced Engineering Mathematics (第2版). New York: McGraw-Hill . p. 454 [定理2、セクション10.8].
参考文献 ウラジミロフ, VS (1971), 『数理物理学の方程式』 , ロシア語からオードリー・リトルウッド訳. アラン・ジェフリー編. 『純粋・応用数学』第3巻, ニューヨーク: マルセル・デッカー社, MR 0268497 。Wang, WX (1988). 「回転楕円体座標系における同次回転楕円体のポテンシャル.I. 外点における」J. Phys. A: Math. Gen. 21 ( 22): 4245– 4250. Bibcode : 1988JPhA...21.4245W . doi : 10.1088/0305-4470/21/22/026 . Milon, T. (1990). 「楕円体座標における同次楕円体のポテンシャルに関する一考察」J. Phys. A: Math. Gen. 23 ( 4): 581– 584. doi : 10.1088/0305-4470/23/4/027 . ラストール、ピーター (1991). 『ポストプリンキピア:物理学者と天文学者のための重力』 ワールド・サイエンティ フィック 7頁以降ISBN 981-02-0778-6 。 Conway, John T. (2000). 「不均質回転楕円体族の重力ポテンシャルの厳密解」 . Mon. Not. R. Astron. Soc . 316 (3): 555– 558. Bibcode : 2000MNRAS.316..555C . doi : 10.1046/j.1365-8711.2000.03524.x . Cohl, HS; Tohline, JE; Rau, ARP (2000). 「トロイダル関数を用いた重力ポテンシャル決定法の発展」. Astron. Nachr . 321 (5/6): 363– 372. Bibcode : 2000AN....321..363C . doi : 10.1002/1521-3994(200012)321:5/6<363::AID-ASNA363>3.0.CO;2-X . ソーントン、スティーブン・T.; マリオン、ジェリー・B. (2003) 『粒子とシステムの古典力学』 (第5版)、ブルックス・コール、ISBN 978-0-534-40896-1 。Zhu, Lupeia (1988). 「重力と地球の密度構造」 . 地球大気科学科. EAS-437 地球ダイナミクス . セントルイス大学. カリフォルニア工科大学. 2011年7月26日時点のオリジナルよりアーカイブ. 2009年3月25日 閲覧 . Charles D. Ghilani (2006年11月28日). 「地球の重力場」 . ペンシルベニア州立大学測量工学プログラム. 2011年7月18日時点のオリジナルよりアーカイブ。 2009年3月25日 閲覧 。 福島俊夫 (2014). 「重力場の長球状調和展開」 .天体物理学ジャーナル. 147 ( 6): 152. Bibcode : 2014AJ....147..152F . doi : 10.1088/0004-6256/147/6/152 .