[论文解读] Exact Solutions of the Isothermal Lane--Emden Equation with Rotation and Implications for the Formation of Planets and Satellites

Dimitris M. Christodoulou, Demosthenes Kazanas|ArXiv.org|Jun 21, 2007

Thermodynamic properties of mixtures参考文献 2被引用 7

一句话总结

本文推导了在有无差速旋转的情况下，圆柱几何下等温 Lane–Emden 方程的精确解析解，揭示了原行星盘中的径向密度剖面由于边界条件而产生振荡结构。这些振荡形成了密度峰，与太阳系中行星的半长轴惊人地吻合（平均相对误差 ≈ 4%），为行星轨道间距提供了一个数学上严谨的解释，无需依赖提托斯-波得定律。

ABSTRACT

We have derived exact solutions of the isothermal Lane--Emden equation with and without rotation in a cylindrical geometry. The corresponding hydrostatic equilibria are most relevant to the dynamics of the protosolar nebula before and during the stages of planet and satellite formation. The nonrotating solution for the mass density is analytic, nonsingular, monotonically decreasing with radius, and it satisfies easily the usual physical boundary conditions at the center. When differential rotation is added to the Lane--Emden equation, an entire class of exact solutions for the mass density appears. We have determined all of these solutions analytically as well. Within this class, solutions that are power laws or combinations of power laws are not capable of satisfying the associated boundary--value problem, but they are nonetheless of profound importance because they constitute "baselines" to which the actual solutions approach when the central boundary conditions are imposed. Numerical integrations that enforce such physical boundary conditions show that the actual radial equilibrium density profiles emerge from the center close to the nonrotating solution, but once they cross below the corresponding baselines, they cease to be monotonic. The actual solutions are forced to oscillate permanently about the baseline solutions without ever settling onto them because the central boundary conditions strictly prohibit such settling, even in the asymptotic regime of large radii. Based on our results, we expect that quasistatically--evolving protoplanetary disks should develop oscillatory density profiles in their midplanes during their isothermal phase. The peaks in these profiles correspond to local potential minima and their locations are ideal sites for the formation of protoplanets ...

研究动机与目标

• 推导适用于原恒星和原行星盘的圆柱几何下无旋转与有旋转的等温 Lane–Emden 方程的精确解析解。
• 通过识别基于流体静力平衡和边界条件的物理机制，解决太阳系中行星轨道间距这一长期存在的难题。
• 证明由边值问题产生的振荡密度剖面自然地在与观测到的行星轨道半径相匹配的位置形成原行星形成位点。
• 将模型与实证数据进行对比，特别是与 Weidenschilling (1977) 提出的 R−2.5 径向密度剖面进行验证，并检验其对太阳系及系外行星系统的预测能力。
• 提供一种非任意的、基于物理的解释，说明行星系统的有序排列，取代如提托斯-波得定律这类人为经验法则，转而基于在物理约束下的非线性微分方程的解。

提出的方法

• 推导圆柱坐标系下无旋转等温 Lane–Emden 方程的精确解析解，得到一个非奇异且单调递减的密度剖面。
• 将 Lane–Emden 方程扩展以包含差速旋转，导出一类具有幂律基线的新精确解。
• 识别出不满足中心边界条件但可作为渐近参考的基线解（幂律或幂律组合）。
• 通过数值积分强制施加物理中心边界条件，显示实际解围绕基线解振荡，但不会收敛于基线解。
• 将振荡解应用于太阳星云，采用径向密度剖面 ∝ R−2.5（面密度 ∝ R−1.5），与 Weidenschilling (1977) 的结果一致。
• 将模型预测的密度峰与太阳系中行星及矮行星的观测半长轴进行比较，以评估其准确性。

实验结果

研究问题

• RQ1 能否在圆柱几何下，对包含差速旋转的等温 Lane–Emden 方程推导出精确解析解？
• RQ2 在中心处施加物理边界条件时，旋转自引力盘中的径向密度剖面行为如何受到影响？
• RQ3 由边值问题产生的振荡密度剖面是否自然地在与太阳系行星半长轴相匹配的位置形成？
• RQ4 太阳系中行星轨道间距的有序性是否源于流体静力平衡与旋转，而非统计巧合或经验定律？
• RQ5 该模型能否解释其他系外行星系统的轨道构型，例如 55 Cancri 和 HD 37124，其相邻轨道的半长轴比值分别为约 3 和 2，与太阳系相似？

主要发现

• 旋转等温 Lane–Emden 方程存在精确解，并已通过解析方法推导，揭示了一类在物理边界条件下具有振荡行为的新解。
• 由于严格的中心边界条件，原行星盘中的实际径向密度剖面被强制永久地围绕基线幂律解振荡，而不会收敛于基线解。
• 该模型预测的密度峰与太阳系中所有主要行星及矮行星的半长轴匹配，平均相对误差约为 4%。
• 在 0.8 AU 至 11.3 AU 的范围内，平均密度剖面的径向依赖关系必须遵循 R−2.5（面密度 ∝ R−1.5），这一剖面由 Weidenschiller (1977) 经验确定，才能实现该精度。
• 该模型预测在 0.1 AU 至 0.3 AU 之间存在一个凝聚物浓度的内间隙，与观测结果一致：原序星在此区域表现出低近红外发射。
• 该模型的预测与两个系外行星系统（55 Cancri 和 HD 37124）的轨道间距高度一致，其连续半长轴的比值分别约为 3 和 2，与太阳系相似。