【讨论】太阳系是稳定的吗？

并不那么简单

牛顿第一个认识到了行星轨道的物理本质。他的万有引力定律——描述了物体间的吸引力和它们的质量以及距离之间的关系——漂亮地解释了太阳系行星的运动秩序。利用牛顿定律，根据行星目前的位置和速度以及它所受到的所有引力可以预测这颗行星未来的运动轨迹。太阳占据了太阳系中99.8%的质量，因此做为一个非常好的近似，每颗行星的轨道可以被描述成一个个彼此独立并且焦点位于太阳的椭圆。

如果真能如此操控一切，每一颗行星的轨道都会永远保持不变。但是，行星之间还存在微小的相互吸引。结果是，这些微小的效应并不会一直保持“低调”。只要时间足够长，它们会以复杂并且意料之外的方式累积起来进而产生压倒性的效果。

即便是在牛顿那个年代，对行星运动的观测就已经达到了非常高的水平。开普勒证明，行星在天空中的奇特轨迹可以用三维空间中一个简单的椭圆来解释。但直到三代人之后，牛顿时代的天文学家才以足够的精度确定了行星的运动轨迹，以此揭示出了它们和完美椭圆之间的偏差。

牛顿知道行星之间的吸引会改变它们的轨道，同时他也特别渴望能解释木星和土星轨道的一个显著特性。在整个16和17世纪，天文学家发现，木星正在缓慢地螺旋向内运动，而土星则正在逐渐向外。如果这一趋势继续进行几万年的话，整个太阳系就会陷入危机。虽然付出了巨大的努力，但牛顿无法用他的理论来解释正在发生的这一现象。其中的数学实在是太吓人了。在一封信中，他含蓄地承认了失败：“…如果我没有弄错的话，同时考虑所有对（行星）运动的影响并且使用简单计算即可得的精确定律来定义这些运动超出了人类思维所能及的范围。”

牛顿在解释木星和土星“迷途”轨迹上的失败为18世纪最杰出的数学家提供了巨大的动力。1776年拉普拉斯攻克了这一难题。他证明，木星和土星的轨道会在一个平均值附近出现周期为数千年的振荡，这是由于木星绕太阳公转五圈的时间正好和土星公转两圈的时间相等造成的。这一“共振”使得它们施加在彼此身上的扰动可以连续累加达数百年。一系列看似可以忽略的影响，只要时机恰当也能产生非常实实在在的长期效果。

在建立了他的理论之后，拉普拉斯可以有效地来反演太阳系，进而预报古代行星在天空中的位置。他的结果和2千年前巴比伦人的观测惊人地吻合，根据他自己的说法“公元前228年3月1日4时23分（巴黎时间），土星位于恒星室女γ下方两个手指处”。这一在纵横数千年里确定土星位置的成功使得他对于自己理论的正确性拥有极端的自信，并毫无疑问地导致了拉普拉斯决定论——如果你知道宇宙中所有粒子现在的位置和速度，那么就能精确地知道它们的整个未来。



[图片说明]：完美可预知性的终结。左：拉普拉斯（1749-1827），非凡的数学天才，解决了木星和土星轨道的缓慢变化之谜。这一成功帮助导致了“拉普拉斯决定论”——如果你知道宇宙中每一个粒子目前的状态，就能彻底地预言它们的整个未来。中：勒威耶（1811-1877），他所在时代的杰出数学天才，第一个发现了完美轨道决定论存在问题。右：庞加莱（1854-1912）证明根本无法预言行星在遥远未来的位置。即使是它们目前轨道中微观的不确定性也最终会通过蝴蝶效应迅速演变成压倒性的不可预知性。

可预知的梦想

在拉普拉斯的公式下，太阳系是完全稳定的。他的数学框架不允许行星的周期存在任何长期演化，而且拉普拉斯认为他已经证明了行星的轨道是“牢不可破”的。

然而，到19世纪中太阳系这部精密的机器开始出现了裂缝。天文学家逐渐开始明白，拉普拉斯理论中的基本近似方法使得它无法适用于数千年以上的时间跨度。勒威耶（他根据海王星对天王星的摄动精确地预报了海王星的存在）警告，太阳系可能并不像拉普拉斯所说的那么精准。几十年内，这个问题变得日益紧迫。一个奖金颇丰的国际竞赛也应运而生，第一个证明行星轨道稳定性的人将获得由瑞典国王奥斯卡所颁发的金牌和2,500克朗的奖金。

让人沮丧的是，由庞加莱提交并最终胜出的证明却是这个问题不可解。庞加莱证明，即使一个仅仅包含了一颗恒星和两颗行星的简单系统也是“不可积”的，即找不到一个公式能永远告诉你这两颗行星的位置——进一步也无法对它们轨道的长期稳定性给一个“是”或者“否”的回答。

庞加莱的工作领先了他所在的时代差不多100年，他第一个暗示了混沌理论以及现在被广为人知的“蝴蝶效应”——系统中一个微小的扰动会产生压倒性的效果。（这个词来自大气建模中的发现，即便是一只巴西的蝴蝶在风中搧动翅膀也会在几个月之后改变全球的天气。）庞加莱的工作证明，太阳系的确定轨道无法外推至无穷远，同时任何有关行星命运的结论都必须用概率来表达。



[图片说明]：最坏的情况。如果水星在接下去的60亿年里落入了和木星的共振，它的轨道最终会被拉长并穿越金星的轨道——随后水星可能会被抛射到任何地方，其中也包括了和地球间可能发生的碰撞。但水星在太阳死亡前失控的概率仅有大约1%。

强力解

在当代，有三件事情使得这一曾经晦涩难懂的天体力学领域得以获得新生。空间飞行对行星的位置和空间探测器轨道的精度要求极高。超级计算机通过在大量小时间跨度上反复计算行星的位置和速度进而对未来行星的运动进行高精度强力模拟。第三，太阳系外行星（包括具有共振和混沌迹象的多行星系统）的发现重新唤起了对太阳系命运的兴趣。

20世纪80年代，许多天文学家使用超级计算机研究了太阳系动力学中的微妙之处。这些计算的成果之一就是直接彰显了太阳系确实是混沌的，更重要的是能进行有效轨道预报的时间跨度只有几百万年——远远小于太阳系的年龄。举个例子，即使现在行星的位置能被确定到误差小于一个原子的直径，但你仍然无法准确地预报1亿年之后它们的位置。我们没有办法确切知道公元100,000,000年1月1日元旦这一天是出现在冬季还是夏季以及地球是否还在绕太阳转动。


1995年巴黎天文台勒威耶的“后来人”Jacques Laskar发表了一个有趣的数值实验的结果。在计算机代数的帮助下，他发展出了一种可以在时间上向前推进行星运动的近似方法。实现这一方法的计算机程序包含了超过15万项，它们描述了行星间微妙的引力相互作用，使得他可以在计算机中以200年为一步向前推动太阳系的演化。由此他也成为了第一个在超过十亿年的时间跨度上研究太阳系多种可能演化过程的人。

那么是否有结果表明太阳系中的行星可能会打破它们由来已久的轨道呢？

在向前演化了5亿年之后，Laskar的程序没有显示灾难将要发生的明显趋势。木星、土星、天王星和海王星这些巨行星没有显示出任何不稳定的迹象。（我们现在认为它们在1千万亿年的时间跨度上都是稳定的。）但类地行星则呈现成了微妙的混沌和轨道的漂移，其中尤其是水星最终有可能会失控。

[图片说明]：害群之马。看似无辜的水星却是掌控这太阳系是否失控的唯一关键。

Laskar因此专门进行了一项实验来研究水星的运动情况。为此，他检查了他5亿年的计算，从中挑选出了水星椭圆轨道偏心率达到最大的那一刻。然后，他用这一时刻的轨道构形做为另外四个几乎完全相同的太阳系模型的初始条件，而在这四个模型中仅有地球轨道存在完全无法可观测到的微小差异。当他在时间上向前演化这四个模型时，一开始它们彼此齐头并进，但在几百万年之后蝴蝶效应开始显现，不同的模型出现了不可避免的差异。

在这四个模型都完成了5亿年的演化之后，Laskar会再一次检查行星的轨道，并从中挑选出水星偏心率达到最大的那一刻。接着他会以此做为另外新的四个数值模拟的初始条件，再一次重复上面的步骤。

在重复了十几次之后，就会产生大量可能的结果。Laskar发现，在其中一些结果里水星的轨道会被拉长并且出现和金星轨道相交的危险情况。轨道相交通常会导致灾难：碰撞、会导致潮汐力撕碎一颗或者两颗行星的密近交会或者是把一颗行星完全抛射出太阳系。

对于地球而言，它的结果包括了从经历更频繁的小行星碰撞到直接被摧毁不等。

今天的图景

Laskar1994年的结果第一次有力地证明，即便没有外界的影响，在太阳做为一颗燃烧氢的恒星所剩下的60亿年时间里也会出现轨道混沌。不过，一些重要的问题仍有待回答。重复选择看似最不稳定的轨道到底有多大的“帮助”？如果进行一次从头到尾的模拟，太阳系能维持多久？是什么机制使得水星的轨道失稳？最后，如果在计算中考虑了广义相对论的微妙作用，结果是不是会有不同？

这些问题最近都得到了回答。多个彼此独立的小组进行了数值模拟，发现了太阳系的唯一致命弱点——使得水星失控的机制。正是由于木星的引力作用，水星一直处于危险之中。

行星间的引力会造成一系列轨道相对于完美椭圆的偏离效果。在这些扰动中最显著的就是轨道进动。当一个轨道发生进动的时候，其椭圆的长轴会改变指向，因此行星的近日点也会沿顺时针或者逆时针方向缓慢但稳定地移动。水星轨道目前的进动速度为每年0.16°，而木星的则为每年0.23°。但数值模拟显示，在非常长的时间下，引力作用会使得水星的进动速率大幅提升。最具有戏剧性的是，如果水星的进动速率接近木星的话，就会引发长期共振，这会迫使水星轨道的进动和木星的同步。


[图片说明]：灾难之路。水星具有一条偏心率适中的轨道。木星的轨道（不在此图之内）也有偏心率，但较之更小。这两条椭圆轨道的长轴分别被称为“拱线”，它们会缓慢地沿着逆时针方向（蓝色箭头）进动。如果两者的速率达到同步，那么如图所示水星的轨道会被拉得越来越长，直到它和金星的轨道相交。此时，水星和金星间的密近交会会把它们朝任何一个方向甩出去，由此也把混沌传播遍整个内太阳系。

长期共振的形成还会带来其他的麻烦。在数百万年的时间里，木星会逐渐地提取水星轨道的角动量。这一效应在大质量的木星身上几乎看到不到效果，但对于小得多的水星而言它的轨道偏心率会增大到能引发灾难的程度。

现在的计算机已经发展到了能够进行数千次太阳系模拟的程度，其中每一个都不采用近似方法而且还包含了广义相对论以及诸如小行星谷神星和月球等天体的影响。2009年夏，Laskar和他的合作者Mickael Gatineau公布了迄今所进行的最庞大的研究结果，他们共进行了2,501组太阳系模拟。

[图片说明]：建模水星。上图：在绝大多数水星的长期模拟中，一切都波澜不惊。这里显示的是水星轨道偏心率的典型变化，从现在开始到未来的20亿年中没有剧烈的起伏。（偏心率为0意味着以太阳为中心的圆轨道。偏心率为1则是一条以太阳为端点的、形似直线的椭圆轨道。目前水星的偏心率为0.21。）
下图：然而在少数情况下，水星会失控。在这一模拟中，其作者会从中挑选出水星偏心率较大的时刻，并以此做为“分支点”开始四组新的模拟。当这些新的模拟中水星又出现大偏心率时，再一次将其做为分支点，循环往复。在这些模拟中内太阳系偶然会出现失控的结果。

这些模拟包含了大量的细节。例如，它们显示，我们很幸运爱因斯坦是对的。著名的广义相对论效应会为水星的进动速率增加额外的每年0.43个角秒。这使得水星的轨道很难与木星发生危险的长期共振。如果没有这一广义相对论效应，水星轨道失稳的概率为10%。但有了广义相对论的作用，在太阳死亡前太阳系失控的可能性仅有大约1%。

虽然很小，但1%的概率仍然是不可忽略的。Laskar和Gastineau计算出了几种太阳系未来可能的轨道，其中地球的情况都不妙。在其中尤其剧烈的一种情况下，地球会和火星发生毁灭性的直接碰撞。在另一种情况下，火星会在距离地球表面仅几百千米的地方掠过。虽然没有撞上，但这也绝非是一件好事。在地球和火星逐渐靠近的过程中潮汐拉伸和挤压作用会把地球加热到完全能熔化地壳和地幔的程度。地球的海洋会蒸发成一个包围着全球性岩浆海的水蒸汽大气。

[图片说明]：密近交会。行星间的密近交会和碰撞一样糟糕。如图所示的计算机模拟显示，如果火星从距离地球足够近的地方掠过，潮汐力会使得较小的行星严重变形、自转加快，并且还会有大量的岩幔被拽出，形成小行星大小的岩石和残块流——其中一些会被地球俘获，引发灾难性的后果。在这里地球被不现实地显示为一个球形。在真实的情况下它也会出现暂时的变形。事实上，地球和火星间极端的密近交会耗散地球内部大量的潮汐能，使得整颗行星熔化。即使地球和火星根本没碰上，地球也会成为一个岩浆星球。

我个人更愿意看着那99%而不是剩下的1%。现在对于太阳系稳定性这一数百年之久的问题我们已经有了一个确定的概率解。在接下去的5千万年里任何行星都不会失控。在太阳面对其60亿年之后的终点时，行星则确实会展现出它们不同寻常的轨迹。

稳定是一个相对的概念，没有绝对的稳定

太阳系还是比较稳定的。
对于宇宙来说人类活动的时间太短了，不足以判断稳定与否。