[论文解读] Newtonian and relativistic polytropes

Fernanda Machado Araujo, Cecilia Chirenti|arXiv (Cornell University)|Feb 11, 2011

Elasticity and Material Modeling参考文献 7被引用 3

一句话总结

本文通过兰恩-艾梅德方程和TOV方程，对牛顿力学与广义相对论下的多级星进行了数值比较研究，表明只有广义相对论能正确预测中子星的最大质量极限。研究识别出当n=1时，临界致密参数σ_CR ≈ 0.42为径向不稳定的起始点，凸显了对致密天体必须引入相对论修正的必要性。

ABSTRACT

The equation of state inside very compact objects like neutron stars is still largely unkown. Even though a lot progress has been made in recent years to develop the so-called realistic equations of state, a lot of insight can be gained by using polytropic equations of state to integrate the stellar equations of structure. In this work we provide a brief review of the Newtonian and relativistic equations of structure and present some numerical results, which we believe that can be useful for students starting to work on this field. The internal structure of the Newtonian polytropes is obtained by the numerical integration of the Lane-Emden equation, and we used our results which can also be used to study the Chandrasekhar limit for white dwarfs. However, Newtonian physics cannot correctly describe very compact and massive objects. There is an upper mass limit for neutron stars, supported by observations, which is not predicted by the Newtonian equations. Neutron stars are best described under the framework of General Relativity. The introduction of the TOV equations (as well as the relativistic Lane-Emden equations) is, therefore, necessary to correctly identify stable and non-stable models via the mass-radius relation. The analysis of the EoS also becomes relevant to discard models that allow a possible violation of causality (sound speed larger than c).

研究动机与目标

• 为进入天体物理学领域的学生提供牛顿与相对论恒星结构方程的教学性综述。
• 对不同多级指数n的多级状态方程，数值求解兰恩-艾梅德方程与相对论TOV方程。
• 通过临界致密参数σ_CR识别相对论模型中径向不稳定的起始点。
• 证明牛顿物理无法预测中子星的最大质量，而广义相对论可以。
• 通过检查状态方程模型中可能出现的超光速声速（因果性破坏）来评估物理一致性。

提出的方法

• 使用无量纲变量x和y对牛顿兰恩-艾梅德方程进行数值积分，适用于多级指数n。
• 通过托尔曼-奥本海默-沃尔科夫（TOV）方程实现相对论恒星结构，采用无量纲坐标ξ、ν和θ。
• 采用多级状态方程P = Kρ^γ，其中γ = 1 + 1/n，以关联压强与密度。
• 定义无量纲质量˜M = M / (4πρ_c A^3)与致密参数σ = GM/(Rc^2)，用于分析稳定性和因果性。
• 在MATLAB中通过迭代数值方法求解n = 3.0, 2.0, 1.0, 0.5, 0.1, 0.0的相对论兰恩-艾梅德系统。
• 基于dM/dρ_c > 0的稳定性分析，将σ_CR识别为径向不稳定的阈值。

实验结果

研究问题

• RQ1 牛顿与相对论多级星的质量-半径关系有何不同，特别是在预测最大质量极限方面？
• RQ2 在相对论多级星中，径向不稳定性首次出现的致密参数σ为何值？该值如何依赖于多级指数n？
• RQ3 牛顿多级模型能否恢复并理解白矮星的钱德拉塞卡极限？
• RQ4 与牛顿模型相比，相对论模型在多大程度上避免了超光速声速（因果性破坏）？
• RQ5 相对论兰恩-艾梅德方程的数值解与不同n值下的已知文献结果相比如何？

主要发现

• 当n = 1时，临界致密参数σ_CR ≈ 0.42标志着径向不稳定的起始点，该值将稳定（σ < σ_CR）与不稳定（σ > σ_CR）模型分隔开来。
• 对于n = 3.0的多级星，当σ = 1.0时，最大无量纲质量约为˜M ≈ 0.2307，与白矮星模型的已知结果一致。
• 牛顿模型预测质量随半径减小而无限增长，无法再现中子星中观测到的物理质量上限。
• 相对论模型能正确预测每种状态方程的最大质量，其中˜M在σ ≈ 0.4–0.5附近达到峰值，适用于n = 1和n = 0.5的情况。
• 当n = 0.1时，最大无量纲质量˜M ≈ 0.3977出现在σ = 1.0处，表明在高致密状态下存在强烈的相对论效应。
• 分析结果证实，当σ > σ_CR时，模型表现出dM/dρ_c < 0，表明存在径向不稳定性，因此为非物理模型。