[论文解读] Construction of Highly Accurate Models of Rotating Neutron Stars - Comparison of Three Different Numerical Schemes

Tetsuo Nozawa, Nikolaos Stergioulas+2 more|CERN Bulletin|Apr 21, 1998

Geophysics and Gravity Measurements被引用 4

一句话总结

本文比较了三种数值代码在广义相对论框架下构建快速旋转中子星高精度模型的表现。在使用相同的方程态、插值方法和积分区域的条件下，这些代码对大多数模型的计算结果一致性优于0.1%，验证了其在多方指数 N > 0.5 及无相变的现实方程态下的可靠性，确立了此类模型的标准物理参数集。

ABSTRACT

We conduct a direct comparison of three different representative numerical codes for constructing models of rapidly rotating neutron stars in general relativity. Our aim is to evaluate the accuracy of the codes and to investigate how the accuracy is affected by the choice of interpolation, domain of integration and equation of state. In all three codes, the same physical parameters, equations of state and interpolation method are used. We construct 25 selected models for polytropic equations of state and 22 models with realistic neutron star matter equations of state. The three codes agree well with each other (typical agreement is better than 0.1% to 0.01%) for most models, except for the extreme assumption of uniform density stars. We conclude that the codes can be used for the construction of highly accurate initial data configurations for polytropes of index N > 0.5 (which typically correspond to realistic neutron stars), when the domain of integration includes all space and for realistic equations with no phase transitions. With the exception of the uniform density case, the obtained values of physical parameters for the models considered in this paper can be regarded as ``standard'' and we display them in detail for all models.

研究动机与目标

• 评估三种不同数值代码在广义相对论框架下构建快速旋转中子星模型的精度。
• 研究插值方法、积分区域和方程态对数值精度的影响。
• 确定代码在一系列物理参数下是否产生一致的结果。
• 为数值相对论模拟中的初始数据，建立一组高精度中子星模型的基准集。
• 识别代码在极端情况（如均匀密度星）下失效的条件。

提出的方法

• 本研究采用三种不同的数值代码，每种代码使用不同的数值格式求解轴对称、定态旋转中子星的爱因斯坦场方程。
• 所有代码使用相同的物理参数、方程态（包括多方和现实方程态）以及插值技术，以确保一致性。
• 积分区域被扩展至整个空间，以最小化边界效应并提高精度。
• 代码通过迭代松弛法同时求解相对论流体动力学方程和引力场方程。
• 针对25种多方方程态（N > 0.5）和22种现实方程态构建了模型，覆盖了不同的中心密度和旋转速率。
• 通过比较所有三种代码的关键物理参数（质量、半径、旋转频率和紧凑度）来评估收敛性和精度。

实验结果

研究问题

• RQ1 当构建快速旋转中子星模型时，三种不同数值代码在精度上如何比较？
• RQ2 插值方法和积分区域的选择在多大程度上影响模型的数值精度？
• RQ3 与现实方程态相比，代码在多方指数 N > 0.5 的多方方程态下表现如何？
• RQ4 在极端情况（如均匀密度星）下，代码在何种条件下出现分歧？
• RQ5 这些代码生成的模型是否可被视为标准且可靠的初始数据，适用于数值相对论模拟？

主要发现

• 对于使用多方指数 N > 0.5 和现实方程态构建的大多数模型，三种代码的一致性优于0.1%。
• 对于无相变的现实方程态模型，所有代码的一致性进一步提升至优于0.01%。
• 唯一显著的不一致出现在均匀密度星的极端情况，此时差异超过1%。
• 当积分区域包含整个空间且方程态平滑无相变时，代码产生一致且高度精确的结果。
• 所有模型的物理参数（质量、半径、旋转频率、紧凑度）均被证实可靠，可作为未来模拟的标准参数。
• 本研究建立了一个包含47个模型（25个多方，22个现实）的基准数据集，具有高精度，适用于数值相对论中的初始数据。