对广义相对论中计算能量条件的Python代码存疑

计算科学 计算物理学 智者
2021-12-13 02:52:58

在广义相对论中,一种可能的方法来确定时空 [即洛伦兹流形(M,g)在哪里g是一个任意度量张量。] 是一个“合理的物理”[1]一,然后是检查所谓的能量条件[2], 举个例子:

物理时空是同时满足两者的时空:

{Gμν=:Rμν12Rgμν=8πGc4TμνTμνuμuν0

在哪里Tμνuμuν0对于所有类时向量uα被称为弱能量条件(WEC)

现在,如果你有一个度量和能量-动量张量,你可以计算几乎所有 GR 中的基本量,包括 WEC。知道这个事实,然后考虑这个 SageManifolds 的 Python 代码:

https://nbviewer.jupyter.org/github/sagemanifolds/SageManifolds/blob/master/Worksheets/v1.3/SM_Friedmann_equations.ipynb

如何根据此代码计算 WEC?(我使用这个例子是因为作者使用了一个完美的流体能量-动量张量,另外,他定义了 4-(类时间)向量和 4-(类时间)对偶向量(covectors or one-forms)

[1]WALD.R.,广义相对论芝加哥大学出版社。第 211-219 页。芝加哥,1984

[2]LOBO.FSN,物理学基础理论:虫洞、曲速驱动和能量条件施普林格。第 193-196.vol 186 页。瑞士。2017 年。

1个回答

给定你必须在一个正交基中计算它的分量,然后找到与类时空特征向量对应 Tμνλ0,λ1,λ2,λ3

λ00,λ0+λi0,(i=1..3)

在您的情况下,使用明显的标准正交基,变为对角线因此条件为 这也在 Wald 中进行了解释和维基百科Tμν(ρ,p,p,p)

ρ0,ρ+p0

我对 Sage 不是很熟悉(我今天刚下载)。但我对xAct非常熟悉,它是 Mathematica的强大微分几何和张量操作包套件,我是其中的贡献者之一(xprint包)。

您可能会在AskSage获得更多关于圣人的帮助