[论文解读] Comment on "The Mass-Radius relation for Neutron Stars in $f(R)$ gravity" by S. Capozziello, M. De Laurentis, R. Farinelli and S. Odintsov

Stoytcho S. Yazadjiev, Daniela D. Doneva|arXiv (Cornell University)|Dec 17, 2015

Geophysics and Gravity Measurements被引用 7

一句话总结

此评论挑战了Capozziello等人（2015年）在$f(R)$引力中报告的中子星质量-半径关系，通过约旦框架场方程的数值求解，证明其结果与早期经过充分验证的工作不一致。作者表明，他们在约旦框架下的结果与Yazadjiev等人（2014年）的爱因斯坦框架结果完全一致，从而证明差异并非源于框架选择，而是源于原论文方程和数值实现中的错误。

ABSTRACT

In the present comment we discuss the qualitative and the quantitative differences between the neutron stars models in $f(R)$ gravity presented in the recent paper \cite{Capozziello2015} (Capozziello et al. (2015)) and in our paper \cite{Yazadjiev2014} (Yazadjiev et al. (2014)). The authors of \cite{Capozziello2015} claim that these differences come from the fact that they are using the physical Jordan frame while we are working in the Einstein frame. The results in \cite{Yazadjiev2014} are very well tested by two numerical codes using different formulations of the field equations that give us the confidence that they are correct. Thus we decided to perform an additional test. Namely, we derive the relevant set of field equations in the Jordan frame and solve them numerically in order to disprove the claim in \cite{Capozziello2015} that the difference comes from the use of different frame. The obtained results are in perfect agreement with our previous work \cite{Yazadjiev2014} and therefore the numerical results in \cite{Capozziello2015} seem to be incorrect. At the end we point out several possible sources of error in \cite{Capozziello2015} that can lead to the observed large deviations.

研究动机与目标

• 解决Capozziello等人（2015年）报告的$f(R)$引力中中子星模型与Yazadjiev等人（2014年）早期结果之间的显著差异。
• 检验该差异是否源于使用物理约旦框架与爱因斯坦框架的差异。
• 通过推导并求解约旦框架中的完整约化场方程组，验证爱因斯坦框架与约旦框架结果的一致性。
• 识别并纠正原论文在场方程表述中的具体数值与分析错误。
• 通过使用不同框架与代码结构的独立数值验证，重新确认早期结果的正确性。

提出的方法

• 推导静态、球对称中子星在$f(R)$引力约旦框架下的维度约化场方程，特别针对$R^2$引力。
• 使用专用静态代码对推导出的约旦框架方程进行数值求解，确保与物理爱因斯坦框架结果一致。
• 通过交叉验证，确认爱因斯坦框架（来自Yazadjiev等人，2014年）与新求解的约旦框架方程结果在数值误差范围内完全一致。
• 在约旦框架中使用物理解释的物态方程，以避免物质耦合中的框架依赖性模糊。
• 在两个框架中比较不同$f(R)$参数$a$（或$|\alpha|$）值下的质量-半径关系及中心密度依赖性。
• 识别Capozziello等人（2015年）中的具体错误，包括标量曲率的错误无量纲化处理，以及一个无法恢复广义相对论极限的有缺陷的约化$R$方程。

实验结果

研究问题

• RQ1 为何Capozziello等人（2015年）报告的中子星质量-半径关系与Yazadjiev等人（2014年）的结果存在显著差异，尽管两者均使用$f(R)$引力？
• RQ2 两份研究之间的差异是源于使用了不同框架（爱因斯坦框架与约旦框架）所致，还是源于数值或分析错误？
• RQ3 在约旦框架中求解场方程时，其结果是否与爱因斯坦框架中的结果保持一致？
• RQ4 场方程表述或其数值实现中的哪些具体错误可能导致Capozziello等人（2015年）中观察到的显著偏差？
• RQ5 Capozziello等人（2015年）中标量曲率的行为是否物理解释，特别是在$\alpha$为负值且中心密度较高时？

主要发现

• 通过求解约旦框架中完整约化场方程组所获得的数值结果，与Yazadjiev等人（2014年）在爱因斯坦框架中的早期结果完全一致，证实了框架等价性。
• 声称差异源于框架选择（约旦框架与爱因斯坦框架）的说法被推翻，因为两个框架在数值精度范围内产生完全相同的质量-半径关系。
• Capozziello等人（2015年）的质量-半径曲线表现出定性错误的趋势：随着$a$增大，中子星最大质量显著下降，而Yazadjiev等人（2014年）观察到的仅为微小下降（<1%）。
• Capozziello等人（2015年）报告称，在$\alpha$为负值时，星体内部的标量曲率呈正值，这违反了椭圆方程的最大值原理，表明其行为不物理解释。
• 在Capozziello等人（2015年）的方程（14）中发现一个关键错误：由于符号错误，当$\alpha \to 0$时无法恢复正确的广义相对论极限。
• Capozziello等人（2015年）中无量纲标量曲率被错误地定义为$R / r_0^2$，而正确形式应为$R r_0^2$，这导致无量纲方程与数值解中出现系统性误差。