天体力学

天体力学天文学辞典}}は、万有引力の法則に従う天体の運動を古典力学に基づいて扱う学問である。ニュートン力学から成立した物理学の一分野でありCollins、また位置天文学と並び古典天文学の一角を占める。

惑星の公転運動は主に太陽の重力によって支配されているものの、他の惑星などが及ぼす重力が摂動として無視できない影響を及ぼすため、天体力学ではそのような摂動を解析的に取り扱う摂動論が発達した。その最も単純かつ非自明な問題が三体問題である。月の運動は暦の編纂や航海術への応用という実用的な目的のためにとりわけ精確な予測が求められる一方で、惑星の運動に比べ摂動が大きく影響するため、太陰運動論は何世代にも渡って改良されてきた。また天王星の観測データの異常から海王星の存在を予言しその位置を予測したことでも知られる。

天体力学は軌道共鳴、太陽系の安定性、自転軸の歳差と章動、惑星の平衡形状、自転と公転の同期といった問題をも扱う。20世紀には人工衛星・宇宙探査機の軌道設計および軌道制御を扱う軌道力学が派生し、また天体力学の適用対象も太陽系から惑星形成、ブラックホール、そして球状星団および銀河などへと拡大した。

1. ケプラー運動
---------

中心天体からの重力を受ける天体の運動はケプラー運動と呼ばれる天文学辞典-ケプラー運動}}。ケプラー運動では、天体の位置 はニュートンの運動方程式 を満足する木下Murray。 は重力定数と中心天体の質量と問題の天体の質量の和の積であるgroup="注釈" 。なお天体力学では伝統的に質量の単位として太陽質量 が、重力定数 の代わりにその平方根として定義されるガウス引力定数 が採用されるBeutlerPlummerBrouwer。この単位系では、問題の惑星の質量を とすると が成立するPlummerBrouwer。また時刻の単位としては日が、距離の単位としては天文単位が使われるBeutler。

1.1. ケプラーの法則
------------

ケプラーの法則天文学辞典-ケプラーの法則}}は惑星の軌道の最も基本的な性質を述べたものである福島。

第1法則: 惑星は太陽をひとつの焦点とする楕円軌道を描く。
第2法則: 太陽と惑星を結ぶ線分が単位時間に掃く面積は一定である。
第3法則: 惑星の公転周期の二乗は軌道長半径の三乗に比例する。 第1法則が主張する楕円軌道の形状は長半径  、離心率  によって特定される。中心天体との距離が最も小さくなる軌道上の点を近点  と呼ぶが天文学辞典-近点}}、特に太陽のまわりを運動する天体の場合は近日点 天文学辞典-近日点}}、地球のまわりを運動する天体の場合は近地点 天文学辞典-近地点}} などと呼ぶ。中心天体との距離が最も大きくなる軌道上の点が遠点  である天文学辞典-遠点}}。中心天体と問題の天体の距離 は、中心天体と近点を結ぶ線分と動径がなす角 を用いて と表示される。角 は真近点角または真近点離角  と呼ばれる天文学辞典-離心近点角}}木下。なお を半直弦  と呼ぶ福島Murray。

第2法則は角運動量の保存を意味する福島。第3法則に対応して、長半径 は平均角速度を表す平均運動  と次の関係にあるPlummerMurray。

ケプラー運動には楕円軌道の他に放物線軌道、双曲線軌道が存在する福島。これらはいずれも円錐曲線である福島。

1.2. 軌道要素
---------

天体の軌道およびその上の位置を特定するために用いられるパラメータを軌道要素  と呼ぶ天文学辞典-軌道要素}}。上述の楕円軌道の形を特定するために用いられる長半径 と離心率 は軌道要素のひとつである。さらに、軌道面内における楕円軌道の向きを特定するために近点引数 天文学辞典-近日点引数}}福島 が、軌道面を特定するために軌道傾斜角 天文学辞典-軌道傾斜角}} と昇交点黄経 天文学辞典-昇交点}} が用いられる。まず軌道傾斜角 は天体の軌道面が基準面となす角として定義される天文学辞典-軌道傾斜角}}。天体の軌道上の点で軌道面と基準面の双方に乗る点が昇交点であり、昇交点が黄道面内の基準方向となす角が昇交点黄経 である天文学辞典-昇交点}}。最後に近点引数 は昇交点と近点がなす角である天文学辞典-近日点引数}}。近点引数 の代わりに により定義される近点黄経  を採用してもよい天文学辞典-近日点引数}}。

楕円軌道上の天体の位置を表す角度として真近点角 以外に離心近点角 天文学辞典-離心近点角}} 、平均近点角  、平均黄経  がある。離心近点角 は を満足し福島、真近点角 と という関係にある木下。平均近点角 は近点通過時刻を として により定義され、離心近点角 とケプラー方程式 によって結ばれる木下。平均黄経 は により定義される。これらの角 , , , は時間的に変化する量であるが、近点通過時刻 または元期  での平均黄経 を与えればどの角も現在時刻 から計算できるため、例えば近点通過時刻 を軌道要素として用いれば十分である木下。軌道要素の組 はケプラーの軌道要素または軌道6要素木下と呼ばれ、これによって天体の運動状態を完全に特定できる福島。

具体的な太陽系惑星の軌道要素の値は#太陽系惑星の軌道要素節および#月の軌道要素節を参照。

1.3. 軌道決定
---------

ある瞬間における天体の座標 および速度 が与えられたならば、その天体の軌道要素は一意に定まりそれを計算することができる福島。しかし実際には1回の観測で得られるのは2つの角度だけであり、天体の軌道要素を決定するためには最低3回の観測を行う必要がある。観測データから軌道要素を決定する方法論は軌道決定として知られている。

2. 摂動論
------

惑星の公転軌道は第一に太陽の重力によって支配されており、0次近似としては太陽-惑星の二体問題とみなすことができる。この近似では惑星の軌道要素は一定であり、時間変化しない。しかし実際には惑星の軌道は他の惑星の摂動  によって変化する天文学辞典-摂動}}福島。そこである瞬間の惑星の軌道について、その瞬間に運動状態が一致するような仮想的なケプラー軌道を考え、その軌道要素を惑星のその時刻の接触軌道要素 天文学辞典-接触軌道要素}} と呼ぶ福島木下。接触軌道要素は他の惑星の摂動によって時間変化するため、それを計算することができれば惑星の軌道が求まることになる。このような摂動手法が定数変化法  である木下Plummer。

2.1. 摂動関数とラグランジュの惑星方程式
----------------------

摂動として働く力が重力などの保存力である場合、天体の運動方程式は摂動関数木下谷川または擾乱関数  として知られる関数 を用いて と書くことができる木下Brouwer。例えば太陽系惑星の場合、 番目の惑星の太陽を中心とする座標  での位置 は、運動方程式 を満足するBeutlerMurray。ここに は惑星 の質量であり、摂動関数の第1項を直接項 、第2項を間接項  と呼ぶBeutlerMurray。

摂動関数 による接触軌道要素 の時間変化はラグランジュの惑星方程式 によって記述されるBoccalettiBeutler。ここに はラグランジュ括弧である。接触軌道要素として を取るとき、ラグランジュの惑星方程式は次のように書き下される木下MurrayBoccaletti。 摂動関数 が与えられたならば、それを摂動展開し惑星方程式を逐次的に解くことにより軌道要素の時間変化が計算できる。

2.2. ガウスの方法
-----------

カール・フリードリヒ・ガウスによる方法は摂動関数ではなく天体に働く力を陽に扱うものでありBeutler、非保存力を扱うことができるBeutler木下。この場合、運動方程式を と書くときBeutler、 を軌道要素として摂動方程式は により与えられるBeutler。摂動 の成分としては以下の2通りの与え方があるBeutler木下。

動径成分 、軌道面内の の法線成分 、軌道面の法線成分 。
軌道の接成分 、軌道面内の の法線成分 、軌道面の法線成分 。前者の立場では、軌道要素 に関するガウスの摂動方程式は次により与えられるBeutler。ここに は半直弦である。

2.3. 永年摂動
---------

離心率 や軌道傾斜角 が小さいときの摂動関数 の展開は literal expansion として知られるMurrayBrouwer。これは摂動関数を角度座標の三角関数の和に分解するものであり、具体的な計算方法がラグランジュ、ラプラス、ルヴェリエ、ニューカムら多くの人の手によって研究されてきたMardling2013}}MurrayBrouwer。例えば中心天体のまわりを公転する2天体について考えるとき、その一方の摂動関数は という形に展開される木下Murray。軌道要素の時間変化は、周期摂動とそれより長い時間スケールでの時間変化を引き起こすに分解できるが、太陽系天体では周期摂動より永年摂動の方が重要である天文学辞典-摂動}}。そのため摂動関数から周期摂動を落としたものをラグランジュの惑星方程式と用いることにより永年摂動の計算が可能となるBeutler木下。例えば近点黄経には時間に比例して増大する項が存在し近点移動が生じる木下。一方で、離心率と軌道傾斜角には永年項が存在せず木下非常に長い周期で時間変化する木下Murray。

2.4. 正準変数
---------

天体力学のいくつかの問題にはケプラーの軌道要素ではなく正準共役量を基本変数として用いるハミルトン力学が適しているMurray。例えばドロネー変数木下 は により定義され、, , が正準共役な組となっているMurray。このときハミルトニアンは により与えられるBrouwerBoccaletti。なおこれらの変数はケプラー問題の作用・角変数と関係しているBoccaletti。

正準形式の摂動論は摂動後のハミルトニアンから角変数を消去するような正準変換を構築することによって実現される木下。このような正準変換を施すと、変換後の作用変数が時間変化しなくなり、問題を自明に解くことができる木下。このような変換は摂動の任意の次数まで続けることができるものの、この摂動級数は収束せず、級数を途中で打ち切る必要がある木下。

3. 応用
-----

3.1. 太陰運動論
----------

月の満ち欠けは太陰暦および太陰太陽暦の基礎であり、月の運動は古くから記録されてきた。月の軌道は等速円運動ではなく、そこからのずれが存在する暦-月}}。月の軌道が楕円軌道であることによる不等が中心差  であるが暦-月}}、これ以外に例えば太陽の摂動によって次のような不等が存在する。

出差 : 遠地点または近地点が太陽の向きにあるとき、相対的に強い摂動を受ける効果暦-月}}。
二均差 : 1朔望月の間に太陽の摂動によって地球の重力が実効的に変化する効果暦-月}}。
年差 : 地球の離心率のために一年の間に太陽の摂動の強さが変化する効果暦-月}}。これらの不等を説明し、精度よく月の運動を予測することは太陰運動論または月運動論  として古くから調べられてきた。これには純粋な天文学上の興味に加えて、航海術への応用という実用的な目的があったWepster。月の理論は最も一般には他の惑星の摂動や地球や月が球形でないことの効果を考慮する必要があるが、アーネスト・ウィリアム・ブラウンは太陽、地球、月の三体を質点として扱う場合論を太陰運動論の main problem と呼んだBrouwer。月の運動は惑星の運動に比べて顕著に大きな摂動を受けており太陰運動論-コトバンク}}、主な摂動の原因である太陽と月の距離がほとんど変化しないものの太陽が地球と月に及ぼす引力の差異によって主要な摂動が生じるという点で惑星の問題とは大きく異なっているWilson。19世紀末から20世紀初頭にかけて完成したヒル-ブラウンの理論は最も精緻な月の運動論であると評価されているWilson。

またエドモンド・ハレーによって指摘された、古代から続く月食の記録を比較すると月の平均運動が徐々に増大しているように見えるという永年加速の問題があるWepster。ラプラス、アダムズを含む数世代にわたる長い論争を経てKushner、潮汐摩擦によって地球の自転が減速し時刻の定義自体が変化している効果を考慮することによって永年加速の問題は解決されたWepster。

3.2. 軌道共鳴
---------

同一の中心天体のまわりの2つの公転軌道について、その平均運動が簡単な整数比にあるとき尽数関係  にあるという天文学辞典-尽数関係}}。このような軌道は安定化または不安定化し、平均運動共鳴と呼ばれる天文学辞典-平均運動共鳴}}。より正確には、2つの軌道 , が平均運動共鳴にあるとは、, を整数として が成立することを言う谷川Murray。例えば小惑星帯のカークウッドの空隙と呼ばれる小惑星の数が少ない領域は木星と平均運動共鳴にあり不安定化したものだと考えられている天文学辞典-カークウッドの間隙}}。逆に太陽系外縁部には共鳴外縁天体と呼ばれる海王星と平均運動共鳴にある天体群が存在することが知られており、その代表的なものが2:3の平均運動共鳴にある冥王星である天文学辞典-共鳴外縁天体}}。また2つの1:2平均運動共鳴が同時に成立するときラプラス共鳴と呼び、太陽系では木星系のイオ-エウロパ-ガニメデが唯一の例である。

一方、平均運動共鳴とは異なり、永年摂動による近点移動の振動数が摂動天体の固有振動数と尽数関係にあるときは永年共鳴  として知られている谷川天文学辞典-永年共鳴。これは軌道周期に比べ非常に長い時間スケールでの軌道の不安定化を導く天文学辞典-永年共鳴。

3.3. 太陽系の安定性
------------

太陽系惑星の軌道が長期的に安定して保たれるかという太陽系の安定性の問題はアイザック・ニュートン以来研究されてきた伊藤。ニュートンは太陽系は不安定であると考えていたLaskar。

ラグランジュらによる摂動論の研究を経てLaskar、ラプラスは1776年に永年摂動の1次の範囲では惑星の軌道長半径は時間変化せず安定であることを示したLaskar。シメオン・ドニ・ポアソンはラプラスの結果を拡張し、1808年に2次摂動の範囲でも軌道長半径は永年不変量であることを示したLaskar。しかしユルバン・ルヴェリエは1840年から41年にかけて、長期間の軌道進化では高次の摂動が重要であり、摂動の低次の項だけに基づくラプラスらによる安定性の証明は信頼できないと指摘したLaskar。アンリ・ポアンカレはルヴェリエの問題提起を受けて、1880年代に惑星系の軌道は解析的な解の表示が存在しないこと、そして問題の摂動級数は一般に発散することを証明したLaskar。1960年代のコルモゴロフらによるKAM理論は近可積分系の大部分の軌道は摂動が十分に小さければトーラス上の準周期解となることを示しており、太陽系の安定性をこの路線で証明する研究が行われたLaskar。

一方で、1950年頃からは電子計算機による太陽系の長時間高精度シミュレーションが行われるようになった。初期のものとしては1951年の W. J. Eckert らによる5惑星シミュレーションがある伊藤。Laskar は1989年の論文でシミュレーションの結果リャプノフ時間500万年で不安定化すると主張した谷川。しかしリャプノフの意味での不安定性にもかかわらず、伊藤孝士と谷川清隆は&pm;40億年のシミュレーションでは惑星軌道は安定に存在し続けたと報告している谷川。太陽系の安定性に関する一般的な理論は2009年現在未だ存在しない。

4. 自転と潮汐
--------

4.1. 自転: 歳差と章動
--------------

多くの天体は公転に加えて自転しており、自転運動はオイラーの運動方程式によって記述されるBeutler。測地学では地球の自転を地球に対して固定された座標系で議論することが多いものの、天文学分野では慣性系を用いて議論することが好まれるBeutler。惑星の自転はある軸まわりの回転として表現でき、その軸を 、自転角速度を とするとき自転は角速度ベクトル により記述されるBeutler福島。角速度ベクトルは自転角運動量 と という関係にある。ここに は慣性モーメントテンソル である福島。しばしば座標系として慣性主軸を取り、そのとき慣性モーメントテンソルは主慣性モーメント , , を固有値とする対角行列となるMurray。

地球の自転軸は月と太陽および他の惑星による摂動を受け、複雑に変化する福島。このうち長周期での軸の移動を歳差 福島、より短周期での振動を章動 福島 と呼ぶ。歳差の周期は約2万6000年であり、春分点の移動をもたらす福島。章動のうちもっとも振幅の大きな成分は周期18.6年であり、月の昇交点がこの周期で移動していることによる福島。歳差および章動は木下宙によって1977年に精密な理論が構築された伊藤。

4.2. 潮汐力
--------

潮汐力 天文学辞典-潮汐力}}は重力の非一様性のために生じる非一様な重力の作用でありMurray、月および太陽による潮汐力は海の潮汐の原因として知られている天文学辞典-潮汐}}。潮汐力はまた天体の潮汐変形、潮汐トルク、潮汐加熱天文学辞典-潮汐加熱}}といった現象を引き起こす天文学辞典-潮汐}}。例えば地球の表面における月による潮汐力は、ポテンシャル により与えられるMurray。ここに は地球の半径、, は地球と月の質量、 は地球と月の距離、 は地球の表面重力、 は月の公転面を基準に計った地球上の点の緯度、 はルジャンドル多項式である。

潮汐による海水の移動が生じる摩擦は地球の自転を減速させる天文学辞典-潮汐摩擦}}。この結果、全角運動量の保存により月は地球から遠ざかる天文学辞典-潮汐摩擦}}。

4.3. 惑星の平衡形状
------------

惑星は厳密には球形ではなく、自転による変形および潮汐力による潮汐変形天文学辞典-潮汐}}を被るMurray。このような変形は軸対称であり、近似的に中心軸から計った角度 の関数として という形に表現できるMurray。また潮汐変形の程度はラブ数によって定量化されるMurray。

主慣性モーメント , , を持つ天体がその外部につくる重力ポテンシャル の表式 はマッカラーの公式と呼ばれる福島Murray。ここに は天体の重心とポテンシャルの評価点を結ぶ軸まわりの慣性モーメントであり福島、評価点の座標を とするとき により与えられる福島Murray。

4.4. 自転と公転の同期
-------------

月は常に同じ面を地球に向けているが、これは月の自転周期と公転周期が同期しているためである。これは地球の重力による月の潮汐変形が原因でありMurray、潮汐ロックと呼ばれる天文学辞典-潮汐ロック}}。

5. その他のトピック
-----------

5.1. 三体問題
---------

重力相互作用する3天体の運動を求める問題は三体問題として知られる。第三体の質量が他の二体に比べて極めて小さく、二体に及ぼす重力が無視できるとき制限三体問題と呼び、特に二体が円運動するときを円制限三体問題と呼ぶ木下。この問題は多くの人の手によって調べられてきておりMurray、三体問題は求積法により解くことはできないものの木下、特殊解のひとつであるラグランジュ点はよく知られている。

5.2. 環
------

土星や天王星に存在する環は衛星と相互に重力を及ぼし合うMurray。環の構造や安定性、羊飼い衛星天文学辞典-羊飼い衛星}}といった問題が取り扱われるTiscarenoMurray。

5.3. 彗星と太陽系小天体の軌道
-----------------

彗星は大きな離心率を持ち、特に極端なものはサングレーザーと呼ばれる。しばしば彗星は木星との近接散乱により大きな摂動を受けるが、これは円制限三体問題とみなすことができ、ティスランの判定式によって彗星の同一性を判定できる福島。また彗星が大きな離心率を獲得する機構として古在メカニズムが提案されている。

小惑星などの太陽系小天体の軌道はカオスを示すことでも注目されるMurray。小惑星帯の小惑星の多くは小惑星-木星系の、または小惑星-木星-土星系の平均運動共鳴に由来するカオス軌道を持つTsiganis2007}}。これは軌道要素のカオス拡散といった効果を生じるTsiganis2007}}。

また宇宙塵などの小天体の場合、輻射圧などの重力以外の摂動が軌道進化において重要である場合がある。

5.4. 重力ポテンシャルの高次成分
------------------

厳密には天体は球形ではなく、それに対応して天体の重力ポテンシャルには単極子項への補正が存在する。これは特に地球を周回する人工衛星の軌道に最も大きな摂動として寄与するため、軌道力学では重力ポテンシャルの補正を考慮する必要がある。軸対称な天体の場合には、重力ポテンシャル は、 を天体の質量、 を天体の半径、 を質量分布に関する定数として、ルジャンドル多項式 を用いて

と書ける木下。

5.5. 一般相対論
----------

強重力場のもとでは一般相対性理論によるニュートン重力からの補正が必要となる。これは水星の近日点移動の要因のひとつとして有名であるMisner。例えばシュワルツシルト時空におけるハミルトン–ヤコビ方程式は と書けるMisner。一般相対論効果はブラックホールなどのコンパクト天体で顕著でありMisner、銀河中心の恒星の運動は超大質量ブラックホールの一般相対論効果を強く受ける。また連星パルサーを代表とするコンパクト星連星では重力波放出により軌道が収縮する。

5.6. 惑星形成
---------

惑星形成理論は微惑星の集積として惑星が形成される過程を議論するものであり、微惑星の合体成長過程は天体力学と関係している伊藤。

5.7. 恒星系力学
----------

恒星系力学は多数の重力相互作用する恒星からなる系を取り扱う理論であり、球状星団や銀河の力学的な性質の基礎となる。この理論は天体力学と顕著な繋がりがある他、統計力学およびプラズマ物理学とも関係しているBinneyTremaine}}。

6. 歴史
-----

6.1. ケプラーの法則
------------

1576年から1601年にかけて、ティコ・ブラーエ  はデンマーク、次いでチェコのプラハにおいて太陽と惑星を観測し、望遠鏡がない当時としては最高精度の誤差1 - 2分角でその位置をしたBeutler。ヨハネス・ケプラー  はブラーエの観測結果をもとにケプラーの法則に到達し、1609年の『新天文学』 、1619年の『宇宙の調和』  においてこれらの法則を公刊したBeutler。

6.2. ニュートンの『自然哲学の数学的諸原理』
------------------------

エドモンド・ハレー  の勧めもありTimberlake、1687年にアイザック・ニュートン  は『自然哲学の数学的諸原理』を出版しTimberlake、ニュートン力学および天体力学の基礎を築いた。なおニュートンがプリンキピアを書き上げるにあたって、ロバート・フック Timberlake やジョン・フラムスティード Timberlake ら同時代の研究者の業績に大きく影響を受けている。

まず第1巻でニュートンは質量  および運動量  を定義し、力  について論じているTimberlake。続いて運動の法則を定式化しTimberlake、中心力場のもとでは面積速度が一定であることTimberlake、円錐曲線を描いて運動する物体には距離の二乗に反比例する中心力が作用していることTimberlake、その場合に楕円軌道を描く物体の周期は楕円の長半径の1.5乗に比例することTimberlakeを示した。

さらにニュートンは互いに引力を及ぼす二体問題についても論じ、その重心まわりの運動に帰着できることを示し、逆二乗則の場合には重心まわりの軌道は円錐曲線となることを主張したTimberlake。また、ニュートンはその理論を月の運動に適用し三体問題の一般解を求めようとしたものの見出すことができず、プリンキピアでは近似解についてのみ記述しているTimberlake。

プリンキピアの第2巻は空気抵抗などの抵抗力のもとでの物体の運動を扱っているTimberlake。//The System of the World// と題された第3巻は前2巻とは異なり自然哲学を扱ったもので、ニュートンはそれまでの巻で展開した数学理論を天界の物体の運動に適用したTimberlake。木星の衛星、土星の衛星、そして惑星がいずれもケプラーの法則を満たすことから、天体間には逆二乗則の引力が働いていること、そして地球-月間に働くこの引力は地球上の物体が地球の中心に向かって落下しようとする力と同じものであると論じているTimberlake。そしてこのことからすべての物体間に重力が作用することを主張したTimberlake。さらに第3巻では自転する球体は扁平な形に変形することTimberlake、潮汐が月の引力によるものであることTimberlake、月の運動Timberlake、月と太陽の重力による地球の歳差の計算Timberlake、彗星の軌道Timberlakeといった内容が扱われている。

1693年にハレーは古代バビロニアおよび中世アラブ界の月食の記録を当時の記録と比較し、月の永年加速を指摘したWepster。1749年に :en:Richard Dunthorne は永年加速の大きさを1平方世紀あたり10秒と求めたWepster。

6.3. 解析力学
---------

ニュートンのプリンキピアは当時考案されたばかりの微分法および積分法の使用を避け幾何学的な考察に基づくものであり極めて難解なものであった。プリンキピアの出版後18世紀初頭にかけてピエール・ヴァリニョン 、ヨハン・ベルヌーイ 、Jakob Hermann  らはプリンキピアの内容をゴットフリート・ライプニッツ  らによる微積分学の言葉を用いて理解するようになったCaparrini。1730年頃からはダニエル・ベルヌーイ 、レオンハルト・オイラー 、アレクシス・クレロー 、ジャン・ル・ロン・ダランベール らによって保存則やポテンシャルの概念などが導入され、1760年頃までには現在の力学に近い形にまで整備されたCaparrini。ダランベールは1743年に //Traité de dynamique// を出版したGrattan-Guinness。オイラーは1749年にニュートンの運動方程式を初めて現在知られている形で書き下しているBeutlerBogolyubov。ジョゼフ＝ルイ・ラグランジュ  は1750年代から統一的な原理に基づく力学の再構築に取り組み、現在解析力学として知られる体系を1788年の著書 //Mécanique analytique// にまとめ上げたCaparriniGrattan-Guinness。

6.4. 二体問題と三体問題
--------------

上述のように、アイザック・ニュートンはプリンキピアにおいて惑星軌道が円錐曲線であるならば逆二乗則に従う中心力が作用していることを示したものの、逆に逆二乗則の重力を受けて運動する物体の軌道がどのようなものかという問題に対しては十分な回答を著述しなかった。この問題は1710年の Jakob Hermann の研究、そしてそれに続くヨハン・ベルヌーイの研究によって解決された。

1730年代にピエール・ルイ・モーペルテュイ  率いる観測隊は地球が赤道付近で膨らんでいる扁球であることを証明したTerrall。これにより地球の形状に関するジャック・カッシーニ  の測量Terrallが棄却され、それと対立していたニュートンの理論の正しさが明らかになったBodenmannTerrall。この観測に参加していたアレクシス・クレローは地球の形状に関する1743年の著書 //Théorie de la figure de la terre// を出版した後に天体力学の研究を始め、1747年11月にパリで三体問題に関する口頭発表を行い、月の近地点移動を説明するためには万有引力の法則に逆三乗則に従う付加項が必要であると主張したBodenmann。この主張は激しい拒否反応を引き起こし、短距離側ではなく遠距離側で万有引力の法則を修正する必要があると考えていたレオンハルト・オイラーとの間で論戦となったBodenmannWilson。ダランベールもこの問題に興味を示し、独自のアイデアで研究に参入したBodenmann。1714年に英国が定めた経度法の懸賞金に繋がる可能性からBodenmann月の近地点移動はこの三者による研究競争となったものの、1749年にクレローは当初の主張を撤回し当時は無視されていた太陽による高次摂動を考慮することによって月の近地点移動を説明できることを示しBodenmann、この成果によって帝国サンクトペテルブルク科学アカデミーの賞を1750年に獲得したBodenmannBogolyubov。その後クレローはハレー彗星の軌道の摂動計算などの研究を行っている。

1748年にパリの科学アカデミーは木星と土星の相互摂動に関するコンテストを開催し、レオンハルト・オイラーが優勝したWilson。彼は木星と土星の運動のケプラー軌道からの逸脱を完全に説明することはできなかったもののWilson、その後の天体力学の研究において極めて重要な役割を果たす三角級数の方法を導入したWilson。またオイラーの研究には観測データからのパラメータ推定に関する先駆的な業績が含まれているWilson。

トビアス・マイヤー  はオイラーの木星と土星の理論を発展させ太陽-地球-月系に応用することによりWilson、月の天文表を作成し1753年に出版したWilson。その正確さは1760年までにジェームズ・ブラッドリー  の観測によって裏付けられ、1767年に創刊された航海年鑑の基礎となったWilson。 レオンハルト・オイラーは三体問題を求積するために運動の積分を探し求めたものの、必要な数の積分を得ることはできなかったBogolyubov。そこで三体が同一直線に乗る配位の特殊解に目を向け、1766年に三体問題に関する論文 //Considerationes de motu corporum coelestium// の中で制限三体問題の平衡点であるラグランジュ点のうち直線解と呼ばれる L1, L2 を発見したGeigesBogolyubov。ラグランジュは1772年にすべての平衡点、特に正三角形解を発見した。ラグランジュはまた一般三体問題の18本の方程式を7本の方程式に帰着できることを示している。

円制限三体問題におけるヤコビ積分は1836年にカール・グスタフ・ヤコブ・ヤコビ  によって導入されたGeiges。

6.5. 摂動論の開発
-----------

摂動論の基本的な道具立てはジョゼフ＝ルイ・ラグランジュによって整備されWilson、ピエール＝シモン・ラプラス  によって発展した。接触軌道要素はレオンハルト・オイラーによって厳密に定義されたBogolyubov。ラグランジュは月の秤動に関する研究を行い、1764年にフランス科学アカデミーの賞を獲得した。またラグランジュは1779年に摂動関数を導入したBeutler。

ピエール＝シモン・ラプラスは1773年頃から天体力学の研究を始め、天体の運動および地球の形状・海の潮汐に取り組んだGrattan-Guinness。ラプラスは1776年に永年摂動の1次の範囲では惑星の軌道長半径は時間変化しないことを示したLaskar。また1787年に木星および金星の摂動によって地球軌道の離心率が変化することにより月の永年加速が説明できると主張したWepster。1789年のフランス革命に伴う環境の激変もありながらGrattan-Guinness、ラプラスは1796年に //Exposition du système du monde// をGrattan-Guinness、1799年から1827年にかけて5巻からなる『天体力学論』  を出版したGrattan-Guinness。この著作は以下の内容を取り扱っているGrattan-Guinness。

第1巻:  平衡と運動に関する一般論、 重力と物体の運動。
第2巻:  天体の形状、 海洋と大気の運動、 天体の重心まわりの運動。
第3巻:  惑星の運動、 月の運動。
第4巻:  木星、土星、天王星の衛星、 彗星、 世界観について。
第5巻:  地球の自転と形状、 弾性流体の運動、 惑星を覆う流体運動、 歳差と秤動、 惑星と彗星の運動、 衛星の運動ラグランジュは1814年に出版した //Mécanique analytique// の第2版の中で摂動関数およびラグランジュの惑星方程式といった天体力学の基本的な道具立てをまとめ、高次摂動の系統的な計算が可能であることを示したWilsonGrattan-Guinness。

6.6. 軌道決定
---------

ティティウス・ボーデの法則は1766年にヨハン・ティティウス  によって発見され、1772年にヨハン・ボーデ  によって紹介されたことで知られるようになったMurray。これは、太陽系惑星の軌道長半径が簡単な数列 により与えらえるというものであるMurray。1781年のウィリアム・ハーシェル  による天王星の発見が の予測に一致したため、この法則は一層興味を集めるようになった。

1801年1月にジュゼッペ・ピアッツィ  は に対応するMurrayケレスを発見し2月上旬まで観測を続けたものの、見失ったBritanica-Asteroid}}Grattan-Guinness。そこでカール・フリードリヒ・ガウス  は同年9月からケレスの軌道計算に取り組みGrattan-Guinness、11月にケレスの軌道計算に成功したGrattan-Guinness。ガウスはフランツ・フォン・ツァハ  へ計算結果を送り、ツァハとヴィルヘルム・オルバース  はガウスの予測通りの位置にケレスを再発見したGrattan-Guinness。さらに翌年に発見された小惑星パラスの軌道計算にも成功しGrattan-Guinness、ガウスはゲッティンゲン大学の天文台のポストを得たGrattan-Guinness。ガウスはさらに天体力学の研究を進め、その成果を1809年に『天体運行論』  として出版したGrattan-Guinness。軌道決定に関するガウスの方法はラグランジュやラプラスによるものよりコンパクトであるGrattan-Guinness。なお『天体運行論』は最小二乗法に関する解説が含まれていることでも知られるGrattan-Guinness。

6.7. 正準理論
---------

ウィリアム・ローワン・ハミルトン  は自身の光学に関する研究から着想を得てNakane1834年から35年にかけての一連の論文においてハミルトン力学を創始した。1836年に円制限三体問題に新しい運動の積分を発見したNakaneヤコビはこの論文を書き上げた後にハミルトンの論文を読んだと考えられておりNakane、彼は力が時間に依存する場合へとハミルトンの理論を拡張し1837年に現在ハミルトン–ヤコビ方程式として知られる単一の偏微分方程式を書き下したNakane。ハミルトンの方程式はヤコビによって「正準」  と命名された中根。

6.8. さらなる発展
-----------

1821年にアレクシス・ブヴァール  は天王星の天文表を出版したが、その後の観測はブヴァールの計算と食い違った。これは未知の惑星の摂動によるものであると考え、ジョン・クーチ・アダムズ  とユルバン・ルヴェリエ  は独立にこの未知の惑星の軌道を計算し、ルヴェリエの予測をもとにヨハン・ゴットフリート・ガレ  が1846年に海王星を発見したNASA-Neptune}}。

1833年にシメオン・ドニ・ポアソン  は独立変数として真近点角 ではなく時刻 を取ることを提案し、Philippe Gustave le Doulcetはこの方法を発展させたWilson。

ルヴェリエは摂動関数の7次までの literal expansion を遂行し、1855年に出版した。ルヴェリエの計算結果は最も広く用いられてきたものであるMurray。Felix Boquet は1889年にルヴェリエの結果を8次に拡張したMurrayほか、サイモン・ニューカム  らはさらに理論を発展させたMurray。ペーター・ハンゼン  も摂動論に多くの貢献を行ったBrouwer。

1856年にジェームズ・クラーク・マクスウェル  は土星の環が固体であるならば不安定であることを証明し、無数の粒子からできているであろうことを指摘したTiscareno。

1889年にフェリックス・ティスラン  は彗星の同一性に関するティスランの判定式を提案した。ティスランはまた4巻からなる『天体力学概論』中根  を出版した。

6.9. 太陰運動論
----------

シャルル＝ウジェーヌ・ドロネー  は1860年および1867年に二巻からなる //La Théorie du mouvement de la lune// を出版し、月の運動について論じたWilson。その中でドロネーは Jacques Binet  が1841年に導入した変数をもとにドロネー変数として知られる正準変数を定義しているWilson。ただしドロネーの理論は級数の収束が遅く十分な精度を得るためには多大な計算を要するという難点があったWilson。

ジョージ・ウィリアム・ヒル  は1870年代からドロネーの理論を発展させたWilson。彼は月の軌道を楕円軌道ではなく三体問題の近似解である卵形の軌道として扱いWilson、またそれまで天体力学ではあまり普及していなかった複素指数関数 を全面的に採用したWilson。アーネスト・ウィリアム・ブラウン  は1896年に //An Introductory Treatise on the Lunar Theory// を出版した後も月の理論についての研究を続け、1919年に月の天文表を完成させたWilson。

6.10. 力学系の理論
------------

19世紀末に三体問題の求積不可能性が Heinrich Bruns  によるブルンスの定理、そしてアンリ・ポアンカレ  によるポアンカレの定理によって明らかになったBarrow-Green。ポアンカレはこの定理および関連する彼の研究成果を1892年から1899年にかけて出版された3巻からなる著書『天体力学の新しい方法』中根  にまとめているBarrow-Green。その後ポアンカレは微分方程式の解を解析的に求めるのではなく、その定性的な性質を明らかにする力学系の理論を創始したBarrow-Green。なお、先行する1880年代にはアレクサンドル・リャプノフ  が力学系の先駆的な研究を行っているBarrow-Green。ポアンカレの力学系の理論はジョージ・デビット・バーコフ  らによって受け継がれ20世紀に大きく発展したBarrow-Green。バーコフは1927年に //Dynamical Systems// を出版しているBarrow-Green。

20世紀に入ってからもエードヴァルド・ヒューゴ・フォン・ツァイペル  は伝統的な摂動論の研究を行い、リンドステッド -ポアンカレの方法を発展させた。またトゥーリオ・レヴィ＝チヴィタ Barrow-Green、カール・スンドマン Barrow-Green らは三体問題の数学的な研究を継続した。

6.11. 一般相対性理論
-------------

アルベルト・アインシュタイン  は1915年に一般相対性理論を完成させた。この理論は強重力場中でニュートン理論への補正項を生じ、アインシュタインはこれによって水星の近日点移動の予測値と観測値の不一致が説明できることを示したMisner。後に5巻からなる //Celestial mechanics// を出版したことで知られるMurray萩原雄祐  は1930年代に一般相対論的天体力学の研究を行った伊藤。アインシュタインは1938年にレオポルト・インフェルト、バーネッシュ・ホフマン とともにポスト・ニュートン展開による補正項を含むN体系の運動方程式であるアインシュタイン・インフェルト・ホフマンの方程式を導出した。

6.12. 準惑星と太陽系小天体
----------------

20世紀には観測技術の進展によって太陽系天体が多く発見され、またその理論も進展した。平山清次  は1918年に小惑星の族の概念を導入した天文学辞典-族}}。クライド・トンボー  は1930年に冥王星を発見した天文学辞典-冥王星}}。Mikhail Lidov  と古在由秀  は1961年から62年に彗星が大きな離心率を獲得する機構を説明し得る古在メカニズムを提案した。ピーター・ゴールドレイク  とScott Tremaine  は1979年に環における羊飼い衛星の存在を理論的に予想した。

6.13. 現代の天体力学
-------------

ソ連のスプートニク計画による1957年に世界初の人工衛星スプートニク1号の打ち上げ以降、宇宙空間における人工物の軌道制御を扱う軌道力学が急速に進展した伊藤。また同時期に電子計算機が実用化されたことにより数値シミュレーションによる軌道計算が可能となった伊藤。一方で理論的研究も続けられ、アンドレイ・コルモゴロフ  らによるKAM理論、堀源一郎  らによるリー変換摂動論の開発などの進展があった伊藤。特にKAM理論は摂動論の有効性を一般的に示すものとみなされ、20世紀天体力学最大の成果と評されている伊藤。また数値計算に関して1980年代から90年代に開発された吉田春夫らによるシンプレクティック数値積分法、そして杉本大一郎らによるGRAPEなどの特筆に値する発展がある伊藤。

7. 付録
-----

7.1. 太陽系惑星の軌道要素
---------------

理科年表2021年版による、太陽系惑星の質量および元期2021年7月5.0日TT＝JD2 459 400.5TTの平均軌道要素。黄道と平均春分点は2001年1月1.5日TT＝JD2 451 545.0TTの値による。



7.2. 月の軌道要素
-----------

国立天文台による月の平均軌道要素。 はユリウス世紀数で、 をユリウス日として により与えられる。

8. 脚注
-----

8.1. 注釈
-------

8.2. 出典
-------

9. 参考文献
-------

9.1. 書籍 
-------------

9.2. 書籍 
-------------

9.3. 論文など
---------

10. 関連項目
--------

天球座標系
暦法
三体問題と多体問題
N体シミュレーション