写在前面
这段时间有点闲，不是因为事不多，主要也是不想干，所以看了看膜拜已久的相对论的知识，特此记录，文中涉及到的公式大多是从别的网站找到的，链接我也会在文章底部留下，感觉国内网站上简单分析相对论的并不多，或者是晦涩难懂的长篇大论，或者上升到了玄学的高度，这也是我写这篇文章的目的，这篇文章适合大多数理工科非物理专业（学过高数和线代），想要知道相对论的推导细节的人。其他人也可以看跳过公式即可。 按理说不应该将此文章放到博客园中，但是博客园对页面定制做的比较好，所以比较适合写各种文章。
一、基本概念
1.正交变换
在线性代数中，正交变换是线性变换的一种。对一个由空间 投射到同一空间 的线性转换，如果转换后的向量长度与转换前的长度相同，则为正交变换。两个向量经过正交变换后内积仍然不变
$$\begin{pmatrix}x_1‘\\\\ x_2‘\end{pmatrix}= \begin{pmatrix} a_{11}&a_{12} \\\\ a_{21}&a_{22} \end{pmatrix} \begin{pmatrix} x_1\\\\x_2 \end{pmatrix} $$ 就是一个变换，其中的一种正交变换为： $$\begin{pmatrix} a_{11}&a_{12} \\\\ a_{21}&a_{22} \end{pmatrix}=\begin{pmatrix} cos\phi & sin\phi \\\\ -sin\phi&cos\phi \end{pmatrix} $$
2.闵可夫斯基空间
设 V是实数域上的四维空间，若η 是一个非退化的对称型矩阵且其正惯性指数等于3，则称(V,η)是一个闵可夫斯基空间。 洛伦兹变换 V上的正交变换成为洛伦兹变换。
3.黎曼空间
黎曼空间 是一种矢量空间，它满足空间中存在在度规张量（一种度量距离的矩阵），使临近两点的距离由正定二次型决定，看起来也不知道这是什么意思，后面我会在用到的时候继续提到，到时候再说，它非常重要。这里先说说黎曼几何和其他的集合有什么区别：
| 空间类型 | 空间曲率 | 平行线 | 三角形三内角之和 | 圆周率 | 例子 |
|---|---|---|---|---|---|
| 黎氏空间 | 正 | 无 | >180° | <π | 球面 |
| 欧式空间 | 零 | 一条 | =180° | =π | 平面 |
| 罗氏空间 | 负 | 无 | >=两条以上 | >π | 伪球面，马鞍面 |
我们早知道在欧式空间中，也就是普通的平面直角坐标系中，过一点只能做一条平行线，而在球面上过一点，一条平行线都做不出来，因为球面上的一条直线线必定是最大圆周，这样必然会相交，而在马鞍面上会做出来看起来类似双曲线一样的两条平行线
3.张量
因为张量的性质是它可以满足一切物理定律必须与坐标系的选择无关的特性，两个张量如果在一个参考系下相等，那么经过洛伦兹变换后，再任意参考系下都相等。 广义相对论的做法是将张量先在一个参考系下表示，然后推导出相应的基本方程，因为在不同参考系下都相等，因此得到的物理规律也适用于一切参考系。 即一切物理定律与坐标系选择无关，且符合爱因斯坦的等效原理(后面会提)
4.能动张量
能量动量张量stress-energy-momentum tensor(SEM): 描述能量与动量的密度和通量，其为一个二阶张量(就是矩阵，后面提到二阶或者二维都是这个意思)，关于SEM后面我会详细介绍
5. 爱因斯坦取和原则(Einstein summation convention)
爱因斯坦取和原则是一种标记的约定，在处理关于坐标的方程式中非常有用，该约定是爱因斯坦对于同一个符号同时出现在上或下标中相当于对这个符号出现的所有可能进行加和例如 $$ x^iy_i = \sum_{i}x_iy_i 而对于x^{\mu\nu} 只给出代表符号，那么实际上等于\sum_{i}\sum_{j}x^{ij} $$
我们称这个符号为dummy index
$$ 连续出现两次相同的符号，不强调他们是上标还是下标,都是满足爱因斯坦取和原则的 $$$$- 只出现一次的称为自由索引 $$$$- 但是注意不会出现三个相同的标记在同一项中 $$$$如x_{ii} \cdot y_{ji} 但是x_{ii}+y_{ji}却是可以的 $$$$- 当同时出现上下标时，也表示他们的取和如 x^i_j = \sum_{i,j}x^i_j $$$$- 一个等式两边的自由索引和虚拟索引个数应该一致 $$
二、狭义相对论
狭义相对论建立在惯性系的基础上，惯性系就是一个参考系，它相对于绝对空间静止或匀速直线运动，牛顿通过木桶思想实验说明了绝对空间的存在(感觉并不科学), 狭义相对论旨在描述高速运动物体的物理特性。
相对论历史:伽利略的相对性原理描述了不同的惯性参考系下观察同一事件的相对速度就是单纯的叠加,其坐标变换公式为: $$ \begin{cases} t‘=t\\\\ x‘=x-vt\\\\ y^{‘}=y\\\\ z‘=z\\\\ \end{cases}$$ 爱因斯坦则使用了洛伦兹变换： $$ \begin{cases} t‘=\frac{t-\frac{v^2}{c^2} x}{\sqrt{1-\frac{v^2}{c^2}}}\\\\ x‘=\frac{x-vt}{\sqrt{1-\frac{v^2}{c^2}}}\\\\ y‘=y\\\\ z‘=z\\\\ \end{cases} $$
洛伦兹变换的推导
假设最简单的线性变换， (t,x,y,z)代表s系，另外的一个代表s‘, 其他情况更复杂的情况都可以从基本情况推出来 $$ \begin{cases} t‘=\alpha t + \beta x \\\\ x‘=\gamma t + \sigma x \\\\ y‘=y\\\\ z‘=z\\\\ \end{cases} $$
$$这里的\alpha,\beta,γ和\sigma 都只依赖于v，如果 x=vt,则 x′ = 0$$
因为相对我s系是vt，就说明我是起点，因为是匀速直线运动 代入第二个式子得
$$x‘ = \sigma (v)⋅(x−vt ) ①$$
$$相对另一物体速度为0时，说明x=x^′,则此时\sigma =1 ，即\sigma (0)=1$$
$$根据相对性原理,s观察s‘的相对速度是-v,代入 x = \sigma (−v)(x^′+vt^′ ) ②$$
$$由2式可得 : \frac{x−\sigma x^′}{\sigma v}= t‘ $$
$$将1式代入2式可得 (\frac{1}{\sigma v}−\frac{\sigma }{v})x+\sigma t=t^{‘} ③式$$
$$将该式与方程组中第一个方程对比，可得 \sigma =\alpha$$
$$再由四维时空中两点的间隔对于两种参照系下相等:$$
$$ t^′2−x^′2= t^2−x^2 $$
代入特殊点 x = vt, x‘ = 0 有：
$$(\frac{x−\sigma x}{\sigma v})^2+0 =−t^2+x^2$$
$$两边同除以t得:\frac{1}{\sigma} =−1+v^2 即\sigma =\frac{1}{\sqrt{1−v^2 }}$$
$$将\sigma 代入③式可得\beta=0$$
于是得到洛伦兹坐标变换公式：
$$ \begin{cases} t‘=t- \frac{\frac{v^2}{c^2} x}{\sqrt{\frac{1-v^2}{c^2}}}\\\\ x‘=\frac{x-vt}{\sqrt{\frac{1-v^2}{c^2}}}\\\\ y‘=y\\\\ z‘=z\\\\ \end{cases} $$
对时间乘因子ic,c是具有速度量纲的常数，ict就是具有长度量纲的一个虚数了，将ict第四维度，规定坐标变换时，变换矩阵必须正交，用来保证任意两点的距离不变，比如洛伦兹变换： $$t^′= \frac{\frac{t− vx}{c^2} }{\sqrt{1−(\frac{v}{c})^2 } }$$
通过洛伦兹变换就可以解释 钟慢效应，尺短效应等高速运动物体物理特性发生的差异，这里不再详细解释。
洛伦兹变换的矩阵形式 $$如果把x,y,z依次记为x_1, x_2 , x_3 ，记ict为x^4,写成矩阵形式就是： $$$$ \begin{pmatrix}x_1‘\\\\ x_2‘\\\\x_3‘\\\\x_4‘ \end{pmatrix}= \begin{pmatrix} \gamma & 0 & 0 & i\beta\gamma \\\\ 0 & 1 & 0 & 0 \\\\ 0 & 0 & 1 & 0 \\\\ -i\beta\gamma & 0 & 0 & \gamma \\\\ \end{pmatrix} \begin{pmatrix}x_1\\\\ x_2\\\\x_3\\\\x_4 \end{pmatrix} $$ $$ 上式中 \beta = \frac{v}{c} , \gamma =\frac{1}{\sqrt{1-(\frac{v}{c})}} $$$$ $$此时 s^2=x_1^2+x_2^2+x_3^2+x_4^2=r^2+(ict)^2 $$ $$注意这里说的x是\Delta x，即物体的距离s^2= r^2−(ct)^2$$ $$s就是四维时空两点中的间隔， 为了后面统一表述这里将s^2 替换成 ds^2 $$ $$他们并没有什么差别，d是距离而不是微分， 将ds^2 表示成矩阵相乘的形式 $$ $$ ds^2 = \begin{pmatrix}-1 & & &\\\\ & 1 & & \\\\ & & 1 &\\\\& & & 1 \end{pmatrix} \begin{pmatrix}c^2t^2\\\\ x_2^2\\\\x_3^2\\\\x_4^2 \end{pmatrix}=\eta⋅(c^2 t^2, x_1^2,x_2^2,x_3^2 ) $$ $$ 其中 \eta_{\alpha\beta}=diag(−1,1,1,1)，这里 \eta_{\alpha\beta} 对角线有三正一负，$$ $$这里的系数矩阵称为\eta,也被称为闵可夫斯基度量(度规)， 那么距离也可以表示成 $$ $$ds^2= \sum_{\mu=0}^3\sum_{\nu=0}^3\eta_{\mu\nu} x^\mu x^\nu$$ $$如果用爱因斯坦取和标记表示 ds^2=\eta_{\mu\nu} x^\mu x^\nu （\mu、\nu为矩阵的行和列的取值）$$ $$在闵可夫斯基空间V中的适合 s^2=0的x称为类光向量或者零矢量，$$ $$ 或称x是类光的(lightlike)(因为此时空间距离x^2=(ict)^2 )，$$ $$V中适合 \eta(x,x)>0的向量 x 成为类空向量或称x是类空的(spacelike )，$$ $$适合\eta(x,x)\lt 0 的向量x称为类时向量或称x是类时的(timelike) , $$$$闵可夫斯基空间中所有零矢量的集合构成了该时间所有事件的光锥$$
三、广义相对论
介绍
狭义相对论重点还有有名的质能方程，推导也不难，我的重点在广义相对论的核心--场方程， 爱因斯坦在改造引力理论为协变形式(就是通过某种变换始终等价)时发现，力学方程(比如F=ma)在不同惯性参考系中经过洛伦兹变换后的形式不一样了，当有了加速度就脱离了惯性系，此外引力和非惯性系有深刻关联，于是爱因斯坦想既然在惯性系上无法满足，我们就把物理定律推广到一切参考系。 而其中最重要的工具就是张量了，前面提到张量的性质是它可以满足一切物理定律必须与坐标系的选择无关的特性。 想一想向量是不是与参考系的选择无关呢。
等效原理
等效原理有不同的版本表述：比较多的说法是 弱等效原理: 即引力质量和惯性质量是相等的： 观测者不能再局部区域分辨加速度所产生的的惯性力和物体所产生的引力， 比如星际穿越中通过旋转的太空舱提供的加速度=g，使得宇航员可以感受不到失重。 $$ 对于引力有 F_G= \frac{ (GM_G m_G)}{r} = m_G g = m_1 a, $$$$其中m_G 是引力质量， m_1 是惯性质量，$$在高中时我们默认两个m是一样的，但引力和惯性力是,牛顿第二定律中描述了没有施力物体也有可能产生加速度（爱因斯坦认为加速度是因为质量过大造成时空弯曲，物体仅仅是在弯曲的时空中沿测地线匀速运动，但是放到正常的欧式空间观察则是加速了），惯性力是由物体的质量来决定他的大小，表征物体惯性的大小，m越大，惯性越大。 但是精确的实验得到两个质量是相等的，所以可以将他们进行等效。后面我们将不再区分这两种质量和力。
爱因斯坦场方程
右侧给出了一张推导场方程的先后顺序，下面我会按照从前到后的顺序进行介绍， 因为场方程实际上是边推边猜出来的，最后会介绍爱因斯坦的思路和最后怎么验证他是对的 
1.闵可夫斯基距离
$$s^2= \sum_{\mu=0}^3\sum_{\nu=0}^3 \eta_{\mu\nu} x^\mu x^\nu $$
2.闵可夫斯基度规
$$\eta_{\mu\nu}$$
3.黎曼|伪黎曼流形
流形是局部具有欧几里得空间性质的空间，流形的定义体现了坐标系是局部的,并限制了这些局部范围内坐标允许的变换， 一个表面是一个二维流形，时空是一个四维流形，经典力学空间和广义相对论所用时空都是流形，$$黎曼流形具有正定的量度，即 ds^2 始终为正$$
4.洛伦兹流形(Lorentzian manifold)
他是度规为(p,q)= (1,N-1)的一个流形，因此闵可夫斯基时空是满足四维时空下(N=4)的一个洛伦兹流形，而洛伦兹流形属于伪黎曼流形的一种，接下来将闵可夫斯基时空的度规张量在伪黎曼流形中表示成伪黎曼度规，
5. 度规张量($g_{μν}$):
我们命名时空时空坐标$$\xi^\alpha 即 \xi^0=CT,\xi^1=Y, \xi^2=Y,\xi^3=Z $$ 我们继续表述在新参考系下坐标进行微分： $$d\xi^\alpha =\frac{\partial \xi^\alpha }{\partial x^0 }dx^0+\frac{\partial \xi^\alpha }{\partial x^1 } dx^1+\frac{\partial \xi^\alpha }{\partial x^2 } dx^2+\frac{\partial \xi^\alpha}{\partial x^3} dx^3 $$$$ d\xi^\alpha = \frac{\partial \xi^\alpha}{\partial x^\mu} dx^\mu,d\xi^\beta= \frac{\partial \xi^\beta}{\partial x^ \nu} dx^\nu $$ $$那么 ds^2=\eta_{\alpha \beta} \frac{\partial \xi^\alpha}{\partial x^\mu} \frac{\partial \xi^\beta}{\partial x^ \nu} dx^\nu dx^\mu =g_{\mu\nu} dx^ \nu dx^\mu $$ $$其中 g_{\mu\nu} = \eta_{\alpha \beta} \frac{\partial \xi^\alpha }{\partial x^\mu } \frac{\partial \xi^\beta}{\partial x^\nu } $$ $$ 可以看到g_{\mu\nu}=g_{\nu\mu}，因为\eta是一个对角矩阵 $$ $$ 我们将g的逆表示成 g^{\mu\alpha} $$ $$ g^{\mu\alpha} g_{\alpha\nu}=(\eta_{\gamma\beta} \frac{\partial \xi^\gamma}{\partial x^\alpha} \frac{\partial \xi^\beta}{\partial x^\mu })^{-1}⋅ (\eta_{\gamma\beta} \frac{\partial \xi^\gamma}{\partial x^\alpha} \frac{\partial \xi^\beta}{\partial x^\nu }) $$$$=\frac{\partial \xi^\beta}{\partial x^\mu }^{-1}⋅\frac{\partial \xi^\gamma}{\partial x^\alpha }^{-1}⋅\eta_{\gamma\beta}^{-1}⋅ \eta_{\gamma\beta}⋅\frac{\partial \xi^\gamma}{\partial x^\alpha } \frac{\partial \xi^\beta}{\partial x^\nu} = \delta_\nu^\mu （Kronecker delta） $$ $$ 上面这个表示和结果非常重要，后面也要用到 g_{\mu\nu} 非常重要就在于 $$$$ 它包含了所有关于时空和引力场的信息，g_{\mu\nu}一共有16个项$$$$ 但是由于对称的原因合并下来只有10个项就是上三角或下三角那10个 $$
6. 测地线方程:
对等效原则的考虑，如果我们要描述一个物体在地球引力场中运动，必须遵循以下步骤： $$ 1.用局部惯性自由落体参考系来描述运动 $$ $$ 2.对局部惯性参考系进行变换，变换到大地参考系 $$ $$ 1.用自由落体参考系来描述运动 $$ $$ 在自由下落的参考系中，我们命名时空坐标为\xi^\alpha ,其中\alpha \in {0,1,2,3} $$ $$ \frac{d^2 \xi^\alpha}{d\tau ^2}=0 , ds^2=cd\tau ^2 $$$$ (\because t‘ = \frac{t − \frac{v^2}{c^2} x}{\sqrt{\frac{1−v^2}{c^2}}} \therefore dt‘ = cdt \Rightarrow ds^2=cd\tau ^2) $$ $$ 这里\tau 值得是在静止参考系中观察者测量的时间，也称为固有时间 $$ $$ 2.在大地参考系中描述运动(向上加速) $$ $$ 令x^\mu 是任意非惯性参考系的各轴 $$ $$ 有： \frac{d\xi^\alpha }{d\tau} =\frac{\partial \xi^\alpha }{\partial x^\mu} \frac{dx^\mu}{d\tau} $$ $$ 这就是爱因斯坦取和标记: 对于上式实际上为 \frac{d\xi^\alpha }{d\tau} =\sum_{i=0}^{3}\frac{\partial \xi^\alpha }{\partial x^i} \frac{dx^i}{d\tau} $$ $$ 将链式法则运用到自由落体方程中有： $$ $$ 0 = \frac{d}{d\tau} (\frac{d\xi^\alpha }{d\tau}) =\frac{d}{d\tau}( \frac{\partial \xi^\alpha }{\partial x^\mu } \frac{dx^\mu}{d\tau}) $$ $$ = \frac{d}{d\tau} (\frac{\partial \xi^\alpha }{\partial x^\mu }) \frac{dx^\mu}{d\tau} +\frac{\partial \xi^\alpha}{\partial x^\mu } \frac{d^2 x^\mu}{d\tau ^2 } $$ $$ = \frac{d}{d\tau} (\frac{\partial \xi^\alpha }{\partial x^\mu }) \frac{dx^u}{dx^\nu } \frac{dx^ \nu}{d\tau} +\frac{\partial \xi^\alpha}{\partial x^\mu } \frac{d^2 x^\mu}{d\tau ^2 } $$ $$ = \frac{dx^\nu}{d\tau} (\frac{\partial \xi^\alpha}{\partial x^\mu \partial x^ \nu }) \frac{dx^\mu}{d\tau} +\frac{\partial \xi^\alpha }{\partial x^\mu } \frac{d^2 x^\mu}{d\tau ^2 } $$ $$ 两边同时乘\frac{\partial x^\beta}{\partial \xi^\alpha} 得： $$ $$ 0 = \frac{\partial x^\beta}{\partial \xi^\alpha } \frac{dx^\nu}{d\tau} (\frac{\partial \xi^\alpha}{\partial x^\mu \partial x^ \nu }) \frac{dx^\mu}{d\tau} +\frac{\partial \xi^\alpha }{\partial x^\mu} \frac{d^2 x^\mu}{d\tau ^2 } \frac{\partial x^\beta}{\partial \xi^\alpha} $$ $$ 令\frac{\partial x^\alpha }{\partial x^\mu }=\delta_\mu^\beta (Kronecker delta 克罗内科函数， $$$$即如果\beta和\mu一致值为1，不一样则值为0 ) 带入可得： $$ $$ 0 = \frac{\partial x^\beta}{\partial \xi^\alpha} \frac{dx^\nu}{d\tau} \frac{\partial \xi^\alpha }{\partial x^\mu \partial x^ \nu } \frac{dx^\mu}{d\tau} +\delta_\mu^\beta \frac{d^2 x^\mu}{d\tau ^2 } 因为只有相等时为1，所以可以用\beta代替\mu ,即 $$$$ 0 = \frac{\partial x^\beta}{\partial \xi^\alpha} \frac{dx^\nu}{d\tau} \frac{\partial \xi^\alpha }{\partial x^\mu \partial x^ \nu } \frac{dx^\mu}{d\tau} +\frac{d^2 x^\beta}{d\tau ^2 } $$ $$ 最终： 0 = \frac{d^2 x^\beta}{d\tau ^2 }+\Gamma_{\mu \nu}^\beta \frac{\partial x^\mu}{\partial\tau} \frac{\partial x^ \nu}{\partial \tau} (测地线方程Geodesic equation) $$ $$ 其中 \Gamma_{\mu \nu}^\beta= \frac{\partial x^\beta}{\partial \xi^\alpha} \frac{\partial \xi^\alpha}{\partial x^\mu \partial x^ \nu } (Christoffel symbol) $$ $$ 上式被称为测地线方程 $$ $$ 通过方程可知: $$ $$ 1. 在加速的参考系下，粒子的运动不再是直线运动了，因为某种惯性力施加在它上面： $$$$f^\beta=−m \Gamma_{\mu \nu}^\beta \frac{\partial x^\mu}{\partial \tau } \frac{\partial x^ \nu}{\partial \tau} $$$$ 上述环境如果实在重力场中进行自由落体，那么惯性力就是重力了， $$$$ 如果我们知道christoffel符号的组成部分，我们就可以确定上述中物体的运动速度了。 $$$$ - 这个等式将直线运动推广到了弯曲的时空中，等会我们会介绍如按照度量张量 $$$$看到从时空度两种计算christoffel符号 $$$$ - 该方程实际上代表10项的4个不同方程\beta从0到3， \mu和 \nu 的组合$$$$会有16种，但其中有些对称 所以只有10个 $$$$- 在同一个测地线方程中，我们将在牛顿极限测地线方程中演示如何在$$$$特定条件下推导张量和重力之间的直接关系。 $$
7.Christoffel 符号：
$$ 即测地线方程中的 \Gamma_{\mu \nu}^\beta= \frac{\partial x^\beta}{\partial \xi^\alpha} \frac{\partial \xi^\alpha }{\partial x^\mu \partial x^\nu } $$
8.Christoffel张量
$$ 即测地线方程中的 \Gamma_{\mu \nu}^\beta= \frac{\partial x^\beta}{\partial \xi^\alpha} \frac{\partial \xi^\alpha }{\partial x^\mu \partial x^\nu } $$ $$ 为了能在爱因斯坦场方程中用张量表示，所以这里的christoffel符号 $$$$需要表示成张量 $$$$我们已经将度量张量(metric tensor)和Christoffel符号定义为： $$$$ g_{\mu \nu}= \eta_{\alpha \beta} \frac{\partial \xi^\alpha }{\partial x^\mu } \frac{\partial \xi^\beta}{\partial x^ \nu } $$$$ \Gamma_{\mu \nu}^\beta = \frac{\partial x^\beta}{\partial \xi^\alpha } \frac{\partial \xi^\alpha }{\partial x^\mu \partial x^ \nu } $$ $$ 用不同形式改写 度规张量g_{\nu \mu},$$$$ g_{\mu \nu} = \eta_{\alpha \beta} \frac{\partial \xi^\alpha }{\partial x^\mu } \frac{\partial \xi^\beta}{\partial x^ \nu} , $$$$ 对称地 ：g_{\nu\mu} = \eta_{\alpha \beta} \frac{\partial \xi^\alpha }{\partial x^\mu} \frac{\partial \xi^\beta}{\partial x^ \nu } $$$$ 将 \alpha 替换成\sigma 可以得到 g_{\nu\mu}= \eta_{\alpha \beta} \frac{\partial \xi^\alpha }{\partial x^\mu } \frac{\partial \xi^\beta}{\partial x^ \nu } $$ $$将 \nu 替换成\alpha ,可以得到： g_{\alpha \mu}= \eta_{\alpha \beta} \frac{\partial \xi^\sigma}{\partial x^\mu} \frac{\partial \xi^\beta}{\partial x^\alpha } $$ $$ 计算 g_{\alpha \mu} 对x^ \nu 的偏导数 \frac{\partial g_{\alpha \mu}}{\partial x^ \nu }=\frac{\partial }{ \partial x^ \nu } (\eta_{\alpha \beta} \frac{\partial \xi^\sigma}{\partial x^\mu } \frac{\partial \xi^\beta}{\partial x^\alpha }) $$ $$其中\eta_{\alpha \beta} 是一个常量提出来可以得到: $$$$ \frac{\partial g_{\alpha \mu}}{\partial x^ \nu } =\eta_{\alpha \beta} \frac{\partial ^2 \xi^\sigma}{\partial x^ \nu \partial x^\mu } \frac{\partial \xi^\beta}{\partial x^\alpha } + \eta_{\sigma\beta} \frac{\partial \xi^\sigma}{\partial x^\mu } \frac{\partial ^2 \xi^\beta}{\partial x^ \nu \partial x^\alpha } $$ $$ 接下来给Christoffel符号 两边同时乘以\frac{\partial \xi^\sigma}{\partial x^\beta } 有： $$$$ \frac{\partial \xi^\sigma}{\partial x^\beta } \Gamma_{\mu \nu}^\beta= \frac{\partial \xi^\sigma}{\partial x^\beta } \frac{\partial x^\beta}{\partial \xi^ \lambda } \frac{\partial ^2 \xi^ \lambda }{\partial x^\mu \partial x^ \nu } =\frac{\partial \xi^\sigma}{\partial \xi^ \lambda } \frac{\partial \xi^ \lambda }{\partial x^\mu \partial x^ \nu } =\delta_ \lambda ^\sigma \frac{\partial ^2 \xi^ \lambda }{ \partial x^\mu \partial x^ \nu } $$ $$ 所以将\sigma被 \lambda 替换有: $$$$\frac{\partial \xi^ \lambda }{\partial x^\beta }v \Gamma_{\mu \nu}^\beta= \frac{\partial ^2 \xi^ \lambda }{ \partial x^\mu \partial x^ \nu } $$ $$ 这样就可替换\frac{\partial g_{\alpha \mu}}{\partial x^ \nu } 中的一些项 $$$$ \frac{\partial g_{\alpha \mu}}{\partial x^ \nu } = \eta_{\alpha \beta} \frac{\partial ^2 \xi^\sigma}{\partial x^ \nu \partial x^\mu } \frac{\partial \xi^\beta}{\partial x^\alpha } + \eta_{\sigma\beta} \frac{\partial \xi^\sigma}{\partial x^\mu } \frac{\partial ^2 \xi^\sigma}{\partial x^ \nu \partial x^\alpha } $$$$ = \eta_{\alpha \beta} \frac{\partial \xi^\beta}{\partial x^\alpha } \frac{\partial \xi^\alpha }{\partial x^\rho } \Gamma_{\mu \nu}^\rho +\eta_{\sigma\beta} \frac{\partial \xi^\sigma}{\partial x^\mu } \frac{\partial \xi^\beta}{\partial x^\rho } \Gamma_{ \nu\alpha} ^\rho =g_{\rho \alpha} \Gamma_{\mu \nu}^\rho +g_{\mu\rho} \Gamma_{\nu\alpha} ^\rho (3) $$ $$ 我们交换\mu 和 \nu 可以得到 $$$$ \frac{\partial g_{\alpha \nu}}{\partial x^\mu } = g_{\rho \alpha} \Gamma_{\nu\mu}^\rho +g_ {\nu\rho} \Gamma_{\mu\alpha} ^\rho （4） $$ $$ 交换\mu 和 \alpha 可以得到 $$$$ \frac{\partial g_{\mu \nu}}{\partial x^\alpha } =g_{\rho \mu} \Gamma_{ \nu\alpha} ^\rho +g_{\nu\rho} \Gamma_{\alpha \mu}^\rho (5) $$ $$ (3)+(4)-(5)可以得到\frac{\partial g_{\alpha \mu}}{\partial x^\nu }+\frac{\partial g_{\alpha \nu}}{\partial x^\mu } − \frac{\partial g_{\mu \nu}}{\partial x^\alpha }= $$$$ g_{\rho \alpha} \Gamma_{\mu \nu}^\rho +g_{\mu\rho} \Gamma_{\nu\alpha} ^\rho +g_{\rho \alpha} \Gamma_{\nu\mu}^\rho +g_{\nu\rho} \Gamma_{\mu\alpha} ^\rho − g_{\rho \mu} \Gamma_{\nu\alpha} ^{\rho} −g_{\nu\rho} \Gamma_{\alpha \mu}^{\rho} =2g_{\rho \alpha} \Gamma_{\mu \nu}^\rho $$ $$ 因此g_{\rho \alpha} \Gamma_{\mu \nu}^\rho = \frac{1}{2} (\frac{\partial g_{\alpha \mu} }{\partial x^ \nu }+\frac{\partial g_{\alpha \nu}}{\partial x^\mu } − \frac{\partial g_{\mu \nu}}{\partial x^\alpha }) $$ $$ 两边同时乘以g^{\beta\alpha} ,得 g^{\beta\alpha} g_{\rho \alpha} \Gamma_{\mu \nu}^\rho = \delta_\rho ^\beta \Gamma_{\mu \nu}^\rho = $$$$ \frac{1}{2} g^{\beta\alpha} (\frac{\partial g_{\alpha \mu}}{\partial x^ \nu }+\frac{\partial g_{\alpha \nu}}{\partial x^\mu } − \frac{\partial g_{\mu \nu}}{\partial x^\alpha }) $$ $$ 将\beta 替换成\rho \Gamma_{\mu \nu}^\rho = \frac{1}{2} g^{\beta\alpha} (\frac{\partial g_{\alpha \mu}}{\partial x^ \nu }+\frac{\partial g_{\alpha \nu}}{\partial x^\mu } − \frac{\partial g_{\mu \nu}}{\partial x^\alpha }) $$ $$ 通常我们将\frac{\partial g_{\alpha \mu}}{\partial x^ \nu } \equiv g_{\alpha \mu, \nu} $$ $$ 所以 \Gamma_{\mu \nu}^\rho = \frac{1}{2} g^{\beta\alpha} (g_{\alpha \mu,v}+ g_{\alpha \nu,\mu}− g_{\mu \nu,\alpha } ) $$
9.Christoffel符号2(连接系数)
$$ 当表述一个矢量场V(x,y,z)时，举个例子 V = (3xy) e_x + (x+ 4y + 3z) e_y + (2y) e_z , $$ $$ 它表示了一个场的运动轨迹，当我们对各分量求导时，我们不用担心基向量e,$$$$因为它就是一个常量，但是如果是在一般的坐标系上，如果各分量不正交，$$$$那么e就不能直接省略了，因为e_x 并不一定仅仅对应x^\alpha ，一般的 $$$$ 我们令V = V^\alpha e_\alpha 注意这里使用的是爱因斯坦取和标记 $$$$ 对于一个向量V = V^\alpha e_\alpha ，其对x^\beta 求偏导 为： $$$$ \frac{\partial V}{\partial x^\beta } =\frac{\partial V^\alpha }{\partial x^\beta } e_\alpha +V^\alpha \frac{\partial e_\alpha }{\partial x^\beta } $$$$ 这样Chirstoffel 符号又可定义成:\Gamma_{\alpha \beta}^\gamma (e_\gamma ) = \frac{\partial e_\alpha ?}{\partial x^\beta } $$$$ 其中\alpha 表示要求到的矢量，\beta表示 坐标变化引起\alpha 变化的矢量， $$$$ \gamma 表示该分量该点导数的指向 $$$$ 这种形式的christoffel符号又称为 连接系数 $$$$ 最后，Christoffel符号具有以下特征：$$$$ -它们是在下部索引对称，即，\Gamma^{\gamma}_{\alpha \beta} =\Gamma{\gamma} _{\beta\alpha} （这是从上述定义清楚）[1] $$
10. 协变导数：
在微分几何里，人们希望推广这个概念到一般微分流形上。首先求导（或求微）的对象从函数推广到向量场（就是向量丛的截面，如切向量场和余切向量场）， 定义域则移到了整个流形上（不再是平坦的空间）， 求导的方向可以是任何切向量的方向。 这样得到的导数就称为协变导数，其微分称为协变微分。
向量的协变微分
$$ 对于一个向量V = V^\alpha e_{\alpha} 其对x^{\beta} 求偏导 为：\frac{ \partial V}{\partial x^\beta }=\frac{\partial V^\alpha }{\partial x^\beta} e_\alpha +V^{\alpha} \frac{\partial e_\alpha }{\partial x^\beta } $$$$ 之前我们得到连接系数为： \Gamma_{\alpha \beta}^\gamma (e_\gamma ) = \frac{\partial e_{\alpha} }{\partial x^\beta } $$$$ 将其带入上式可得： \frac{\partial V}{\partial x^\beta} =\frac{\partial V^\alpha }{\partial x^\beta } e_\alpha +V^\alpha \Gamma_{\alpha \beta}^\gamma e_\gamma $$$$ 等式的右侧有两个虚拟符号\alpha 和\gamma ，交换他们 (为什么能交换) $$$$ \frac{\partial V}{\partial x^\beta}=\frac{\partial V^\alpha }{\partial x^\beta} e_\alpha +V^\gamma \Gamma_{\gamma\beta}^\alpha e_\alpha ， 提出e_\alpha $$$$ \frac{\partial V}{\partial x^\beta}=(\frac{\partial V^\alpha }{ \partial x^\beta } +V^\gamma \Gamma_{\gamma\beta}^\alpha )e_\alpha $$$$ 此表达式表示V在x的\beta方向上的变化率 $$$$ 可以用 \nabla来表示 协变导数值: $$$$\nabla_\beta V^\alpha =\frac{\partial V^\alpha }{ \partial x^\beta } +V^\gamma \Gamma_{\gamma\beta}^\alpha $$$$ \nabla_\beta V^\alpha 的等价符号是 D_\beta V^\alpha 和 V_{;\beta}^\alpha $$$$ 备注: $$$$ - 度规张量的偏导数消失在M的任何点等价于Christoffel符号的消失， $$$$这样我们可以在任何惯性或局部惯性狂阿基背景下重写整个基本等式: $$$$\nabla_\beta V^\alpha =\frac{\partial V^{\alpha} }{\partial x^\beta } $$$$ - Christoffel符号本身不会转换成张量是因为，总可以找到令Christoffel为0 $$$$ 的局部惯性系，但对于张量却不会。 $$$$ 此外还可以推导出协变向量的协变微分。 由于此等式对于每个任意A向量都成立，$$$$因此方括号中的数量必须为null。因此，我们证明了向量B的协变分量的$$$$协变导数的表达式如下： $$$$ \nabla_\alpha B_\mu = \frac{\partial B_\mu}{\partial x^\alpha} - B_\sigma \Gamma^\sigma_{\alpha \mu} $$$$ 其中B为协变向量 $$
11. 黎曼张量
换向差(commutator)：表示一个操作整箱和方向的差异如:[AB] = AB-BA ,如果更改顺序后没有影响，则换向差为0 ，因为在平坦的空间求协变微分操作先后都一样将换向差应用于协变微分的任何非零结果都可以归因与空间的曲率，从而归因于黎曼张量， 黎曼张量实际上衡量了时空的曲率
$$ 我们通过计算协变微分的换向差（commutation）\nabla_c \nabla_b V_\alpha −\nabla_b \nabla_c V_\sigma $$$$来推导黎曼张量 或者用分号表示:V_{\alpha ;bc} −V_{\alpha ;cb} $$$$ 前面推导的协变微分为:\nabla_b V_\alpha =\partial _b V_\alpha −\Gamma_{ab}^d V_d $$$$ 于是： \nabla_c \nabla_b V_\alpha = \partial _c (\nabla_b V_\alpha )−\Gamma_{\alpha c}^e \nabla_b V_e −\Gamma_{bc}^e \nabla_e V_\sigma $$$$ 其中右侧第一项：\partial _c (\nabla_b V_\alpha )=\partial _c \partial _b V_\alpha −\partial _c (\Gamma_{ab}^d V_d ) $$$$ \partial _c (\nabla_b V_\alpha )=\partial _c \partial _b V_\alpha −\partial _c (\Gamma_{ab}^d ) V_d−\partial _c \Gamma_{ab}^d (V_d ) $$$$ 右二\Gamma_{ac}^e \nabla_b V_e=\Gamma_{ac}^e (\partial _b V_e −\Gamma_{be}^d V_d) $$$$ 右三\Gamma_{bc}^e \nabla_e V_\alpha =\Gamma_{bc}^e (\partial _e V_\alpha −\Gamma_{\alpha e}^d V_d ) $$$$ 带入公式可得\nabla_c \nabla_b V_\alpha = $$$$\partial _c \partial _b V_\alpha - \partial _c (\Gamma_{ab}^d ) V_d −\partial _c (V_d ) \Gamma_{ab}^d −\Gamma_{ac}^e (\partial _b V_e−\Gamma_{be}^d V_d )−\Gamma_{bc}^e (\partial _e V_\alpha −\Gamma_{ae}^d V_d ) $$$$ 然后交换bc可得\nabla_b \nabla_c V_\alpha $$$$ = \partial _b \partial _c V_\alpha - \partial _b (\Gamma_{ac}^d ) V_d −\partial _c (V_d ) \Gamma_{ac}^d −\Gamma_{ac}^e (\partial _c V_e−\Gamma_{ce}^d V_d )−\Gamma_{cb}^e (\partial _e V_\alpha −\Gamma_{ae}^d V_d ) $$$$ 计算\nabla_c \nabla_b V_\alpha −\nabla_b \nabla_c V_\alpha $$$$ = - \partial _c (\Gamma_{ab}^d ) V_d −\partial _c (V_d ) \Gamma_{ac}^d −\Gamma_{ac}^e (\partial _c V_e−\Gamma_{ce}^d V_d )−\Gamma_{cb}^e (\partial _e V_\alpha −\Gamma_{ae}^d V_d ) $$$$ 展开： \nabla_c \nabla_b V_\alpha − \nabla_b \nabla_c V_\alpha $$$$ = −\alpha _c (\Gamma_{ab}^d ) V_d −\alpha _c (V_d ) \Gamma_{ab}^d − \Gamma_{ac}^e \partial _b V_e+\Gamma_{ac}^e \Gamma_{be}^d V_d+\partial _b (r_ac^d ) V_d+\partial _b (V_d ) \Gamma_{ac}^d+\Gamma_{ab}^e \partial _c V_e+\Gamma_{ce}^d V_d $$$$ = \partial _b (V_d ) \Gamma_{ac}^d −\partial _c (V_d ) \Gamma_{ab}^d+\Gamma_{ab}^e \partial _c V_e −\Gamma_{ac}^e \partial _b V_e+(\partial _b (\Gamma_{ab}^d )−\partial _c (\Gamma_{ab}^d )+\Gamma_{ac}^e \Gamma_{be}^d −\Gamma_{ab}^e \Gamma_{ce}^d ) V_d $$$$ 前面我们计算出Christoffel的另一种表示\Gamma_{\alpha \beta}^\Gamma e_\gamma = \frac{\partial (e_{\alpha } )}{\partial x^\beta } $$$$ 因此： \Gamma_{ab}^e \partial _c V_e − \Gamma_{ac}^e \partial _b V_e= \partial _c \Gamma_{ab}^e V_e −\partial _b \Gamma_{ac}^e V_e （偏导顺序无关）$$$$ =\partial _c \partial _b V_a−\partial _b \partial _c V_\alpha =0 $$$$ 通过交换b,c 可以得到： \partial _b (V_d ) \Gamma_{ac}^d −\partial _c (V_d ) \Gamma_{ab}^d=0 $$$$ 最后，commutation为： $$$$ \nabla_c \nabla_b V_\alpha − \nabla_b \nabla_c V_\alpha =(\partial _b (\Gamma_{ac}^d ) −\partial _c (\Gamma_{ab}^d )+\Gamma_{ac}^e \Gamma_{be}^b −\Gamma_{ab}^e \Gamma_{ce}^d ) V_d=R_{abc}^d V_d $$$$ 将R_{abc}^d 称为黎曼张量。 $$$$ 简化表示成:R_{abc}^d= \Gamma_{ac,b}^d −\Gamma_{ab,c}^d+\Gamma_{ac}^e \Gamma_{be}^d −\Gamma_{ab}^e \Gamma_{ce}^d $$$$ 因为有4个索引，所以 每个索引有4种取值，所以一种有4^4=256个项 $$
12.张量收缩
坐标系变换
$$对坐标系旋转\theta 可以得到新的坐标系，则向量A在两个坐标系下的坐标分别为： $$$$ A=A^x e_x+A^y e_y $$$$ A=A^′x e_x+a^′y e_y $$$$ 其中A^{′x}=L_1+L_2 $$$$ L_1=A^x cos?(\theta ) $$$$ L_2=A^y cos?(\pi/2−\theta ) $$$$ 则(A^′)^x= A^x cos?(\theta )+A^y cos?(\pi/2−\theta ) $$$$ 同理 A^′y= A^x cos?(\pi/2+\theta )+A^y cos?(\theta ) $$$$ 则A‘就可以表示成矩阵的乘法: $$$$ \begin{pmatrix} A^{‘}_{x} \\\\ A^{‘}_{y} \end{pmatrix} = \begin{pmatrix} cos(\theta) & sin(\theta) \\\\ -sin(\theta) & cos(\theta) \end{pmatrix} \begin{pmatrix}A^{x}\\\\ A^{y} \end{pmatrix} $$
基向量变换
$$在新的基向量上表示旋转后的向量，对于向量逆时针旋转\theta 相当于顺时针将坐标系旋转\theta $$ $$ \begin{pmatrix} A^{‘}_{x} \\\\ A^{‘}_{y} \end{pmatrix} = \begin{pmatrix} cos(\theta) & -sin(\theta) \\\\ sin(\theta) & cos(\theta) \end{pmatrix} \begin{pmatrix}A^{x}\\\\ A^{y} \end{pmatrix} $$$$ dx^{′1}=(\partial x^{′1})/(\partial x^{1} ) dx^1+(\partial x^{′1})/(\partial x^2 ) dx^2 $$$$ dx^{′2}=(\partial x^{′2})/(\partial x^1 ) dx^1+(\partial x^{′2})/(\partial x^2 ) dx^2 $$$$ dx^{‘i}=\frac{\partial x^{‘i}}{\partial x^j} dx^j $$$$ 计算适当的导数我们得到: $$$$ \frac{\partial x^{‘1}}{\partial x^1} = cos\theta \qquad \frac{\partial x^{‘2}}{\partial x^1} = sin \theta $$$$ \frac{\partial x^{‘1}}{\partial x^2} = -rsin\theta \qquad \frac{\partial x^{‘2}}{\partial x^2} = rcos\theta $$$$ 因为dx^{‘1}和d^{‘2}实际上反映了e_1 和e_2 的大小所以: $$$$ e_1 = \frac{\partial x^{‘1}}{\partial x^1} \vec{i} +\frac{\partial x^{‘2}}{\partial x^1} \vec{j} = cos \theta \vec{i}+\ sin\theta vec{j} e_2 = \frac{\partial x^{‘1}}{\partial x^2} \vec{i} +\frac{\partial x^{‘2}}{\partial x^2} \vec{j} = -rsin\theta \vec{i} + rcos\theta \vec{j} $$$$这样我们就得到了分量A^{‘i}的定义,又称为逆变(反变或contravariant)分量 $$$$ A^{‘i} = \frac{\partial x^{‘i}}{\partial x^j} A^j $$
协变分量和逆变分量
$$ 图中 X = e_1 x^1+e_2 x^2 $$$$ 我们构造了e_1 和 e_2 的对偶向量e^1 和 e^2，我们希望用对偶向量表示原向量X $$$$ 我们令 协变分量为 x_1= e_1 x ，注意协变分量是标量 $$$$ 我们期望的结果是: X = e^1 X_1+e^2 X_2 $$$$ 对偶向量是对于e^i满足\lt e^i,e_j \gt =\delta_j^i 的向量e^j $$$$ 举个例子:e^1 是 e_1 的对偶向量，他满足 e^1 \cdot e_1=1 $$$$ 即对偶空间的相同下标的基向量内积为1， 关于对偶空间和对偶基可以看： $$ https://blog.csdn.net/sumx2015/article/details/78909635 解释的很形象了 $$ 对偶基的要求还蕴含了任意非相同标记的基向量和对偶基向量内积为0，就是说他们垂直 $$$$ 下面我们来证明为什么这样的向量可以构造出来新的基底 $$$$ 由图根据平行四边形法则我们可知: $$$$ |x|sin\theta =|e^1 x_1 |cos\theta $$$$ sin\theta =\frac{|e^1 (e_1 x)|}{|x|} cos\theta $$$$ sin\theta =\frac{|e^1 | e_1 |cos?(\pi/2−\theta ) |}{|x|} cos\theta $$$$ 1 = |e^1| |e_1|cos\theta $$$$ 即e^1 e_1=1 , 同理可以证明e^2 e_2=1 $$$$ 反过来 也可以推导出 |x|sin\theta =|e^1 x_1 |cos\theta $$$$ 所以\lt e^i,e_j\gt =\delta_{j}^{i} 是满足协变分量为x_1= \lt e_1 ,x\gt 的原向量X在对偶空间上的表示 $$$$ 尤其是当 e_1=1时，A点就和 X_1 重合了，此时X_1 就是x在 e_1 轴上的投影 $$$$ 对于逆变分量， 如果改变基向量e_1 的大小，为了保证X= e则逆变分量X^1 的大小会反向改变, $$ 比如基向量减小 一半，逆变分量就会增加一倍， 而对于协变分量，如果原空间的基向量模增加一倍， 对偶空间的基向量模减小一倍，因为内积为1的约束， 对应协变分量就增加一倍，所以 协变指的是，随着原空间的基向量的变化而正相关变化 之前我们推导对于逆变张量的分量: $$ A^{‘i} = \frac{\partial x^{‘i}}{\partial x^{‘i}} A^j $$$$ 同理对协变张量的分量效果也是一致的： $$$$ A_{‘i} = \frac{\partial x^{‘i}}{\partial x^{‘i}} A_j $$
张量的缩放
当把上面的逆变张量和协变张量推广到普通情况就相当复杂: $$ 令原张量表示成: S_{j_1 j_2?j_q}^{i_1 i_2?i_p } $$$$ 则经过变换后的矩阵\bar{S}可以表示成： $$$$ \bar{S}_{j_1 j_2?j_q}^{i_1 i_2?i_p}= S_{\beta_1 \beta_2?\beta_q}^{\alpha_1 \alpha_2?\alpha_p}\frac{\partial \bar{x}^{i_1 }}{\partial x^{\alpha_1}} \frac{\partial \bar{x}^{i_2 }}{\partial x^{\alpha_2}}\cdots \frac{\partial \bar{x} ^{i_p }}{\partial x^{\alpha_p } } \frac{\partial \bar{x} ^ {\beta_1}}{\partial x^{j_1}} \frac{\partial \bar{x} ^{\beta_2}}{\partial x^{j_2 }}\cdots \frac{\partial \bar{x} ^{\beta_p }}{\partial x^{j_q }} $$$$ 当我们尝试让某个上标i_k 和某个下标j_b 相等时，就可以在右式中约去， $$$$ 反映到左式中就可以直接约去，我们把这种操作叫做contraction $$$$ 张量压缩，当然我们也可以进行扩张，制造两个相同的下标。 $$$$ 度量张量最有用的功能之一就是在其他张量的协变量和范变量分量之间进行转换 $$$$ 因为对 度规张量进行变换可以得到： $$$$ g_{\mu\nu} = \eta_{\alpha \beta} \frac{\partial \xi ^{\alpha}}{\partial x^\nu} \frac{\partial \xi ^\beta}{\partial x^\nu} $$$$ g^{‘}_{\mu\nu} = \eta_{\alpha \beta} \frac{\partial \xi ^{\alpha}}{\partial x^{‘\nu}} \frac{\partial \xi ^\beta}{\partial x^{‘\nu}} = \eta_{\alpha \beta} \frac{\partial \xi ^ \alpha}{\partial x^\sigma}\frac{\partial x^\sigma}{\partial x^{‘\mu}} \frac{\partial \xi^\beta}{\partial x^\rho} \frac{\partial x^\rho}{\partial x^{‘\nu}} $$$$ g^{‘}_{\mu\nu} =\eta_{\alpha \beta} \frac{\partial \xi ^ \alpha}{\partial x^\sigma} \frac{\partial \xi^\beta}{\partial x^\rho} \frac{\partial x^\sigma}{\partial x^{‘\mu}} \frac{\partial x^\rho}{\partial x^{‘\nu}} $$$$ g^{‘}_{\mu\nu} = g_{\sigma \rho}\frac{\partial x^\sigma}{\partial x^{‘\mu} }\frac{\partial x^\rho}{\partial x^{‘\nu}} $$$$ 之所以能够进行压缩，是因为度规张量恰好 $$$$ 需要对下标求偏导，如上面的压缩一样，对相同的标记进行了约分。 $$$$ 因此，给定张量的逆变分量和原始基向量，并希望确定协变分量， $$$$一种方法可能是确定对偶基向量， 如g_{ij} A^j=A_i ，g^{ij} A_i=A^j , g^{ij} A_i=A^j $$
13.Ricci张量和Ricci标量
$$ 通过上一节12中介绍到的张量收缩: $$$$ 可以通过收缩黎曼张量的第一和第三下标来得到ricci张量 $$$$ R^{\mu}_{\beta\mu\nu} = g^{\mu\alpha}R_{\alpha \beta \mu \nu} = R_{\beta \nu} $$$$ 同时使用黎曼张量和度规的对称性可以得到 $$$$ R_{\beta \nu} = g^{\mu\alpha}R_{\alpha \beta \mu \nu} = g^{\alpha\mu}R_{\mu \nu \alpha \beta} = R^{\alpha}_{\nu\alpha\beta}=R_{\nu\beta} $$$$ 由此可见ricci张量是对称的 $$$$可以依次收缩Ricci张量以获得Ricci标量R R^{\nu}_{\nu} = g^{\nu \beta}R_{\beta \nu} = R $$
14. 能动张量：
$$ 应力能量动量张量(Stress \quad energy \quad momentum \quad tensor,SEM):T^\alpha \beta $$$$ 简称能动张量，用来描述能量和动量在时空中的密度与通量(flux), $$$$其为牛顿物理中应力张量的推广，能动张量在广义相对论中是重力场的源， T^\alpha $$$$ \beta 给出四维动量或4-动量P(\frac{E}{c},p_x,p_y,p_z)之\alpha 分量$$$$(\alpha ,\beta依次代表1,2、、)通过一个坐标为常数x^b 表面的通量， $$$$ 另外要注意的是SEM 是对称的(当自旋张量为0时), 即T^{ab}=T^{ba} $$$$ T^{00} 代表能量密度 ， T^{0i} 代表能量通过x^i 表面的能量密度，$$$$等同于 T^{i0} 第i动量之动量密度，当\alpha \beta $$$$ 不同时表示剪应力，相同时表示正向应力。 $$$$ 由于我们最终的目标是在时空几何图形及其内容之间建立关系，$$$$因此我们首先必须找到用什么来描述，$$$$由于质量，能量和动量都相关(E^2=p^2 c^2+m_0^2 c^4)， $$$$ 证明:p^2=(m_0 c)^2 , $$$$ p^2= \eta_u \nu p^\mu p^ \nu=\frac{E^2}{c^2} $$$$ −p_x^2−p_y^2−p_z^2=\frac{E^2}{c^2} −p^2⇒(m_0 c)^2=\frac{E^2}{c^2 −p^2}\Rightarrow E^2 = p^(2c^2 ) $$$$ 所以似乎假设: 广义相对论的引力场来源应该包括动量，能量和质量。 $$$$另一方面，用泊松方程( 泊松方程：\nabla^2 \phi=f ，就是各维的二阶偏导数之和， ) $$$$ 表示的描述由质量密度引起的重力的等效牛顿方程为： $$$$− \nabla g=−\nabla(−\nabla\phi)=4\pi G\rho $$$$ 那么能量密度是张量还是标量呢? $$$$ 先考虑相对静止的非相互作用例子相对于彼此的提及dxdydz,,也被称作 dust cloud $$$$ 在它自己的参考系，这个dust有自己的能量密度（质量密度） \rho _0=m_0 n_0 ， $$$$在这里m_0 代表了尘粒的质量, n_0 $$$$ 代表了单位提及例子的数量，在洛伦兹变换下，会变成: $$$$ m_0 \Rightarrow m_0 \Gamma , n_0→n_0 \Gamma\Rightarrow \rho _0 \Rightarrow\rho =\rho _0 \Gamma^2 u^\mu=\Gamma(1,v) $$$$ =(\Gamma,v_x \Gamma, v_y \Gamma,v_z \Gamma)^T=(u^t,v_x u^t, v_y u^t,v_z u^t )^T $$$$ 所以设置\Gamma= u^t 使用这种相等性，并且每个粒子的能量等于p^t=mu^t,$$$$ 我们改写S‘中的总能量密度： $$$$ \rho ≡np^t=(n_0 u^t )(mu^t )=(n_0 m) u^t u^t=\rho _0 u^t u^t $$$$ 可以将\rho 解释成对称的二维张量，也就是我们说的能动张量， T^{\mu \nu}=T^ \nu\mu=\rho _0 u^\mu u^ \nu $$$$ 上述可知T^{00} 可表示成能量密度 ，对于4x4的其他项来说： $$$$ T^tx=\rho _0 u^t u^x=(n_0 m) u^t u^x=(n_0 u^t )(mu^t ) v_x=np^t v_x=\frac{(nAv_x dt) p^t}{Adt} $$$$ 其中A_vx dt 代表了 在dt时间内垂直x方向穿过表面A的dust体积， $$$$ 可以表示通过该表面的颗粒总数， $$$$ 则T^tx=T^xt 除了T的第一行和第一列外，其他的部分我们表示成 $$$$ T^kl=\rho_0 u^k u^l=(n_0 m) u^k u^l=(n_0 m) u^t v_k u^l=(n_0 u^t ) v_k (mu^l )= $$$$nv_k (mu^l )=nv_k p^l=mv_k p^l=\frac{(nAv_k dt) p^l}{Adt} $$$$ 其中p^l 是动量，则T^{kl} 表示l动量在k方向的通量， $$
15.牛顿极限
牛顿极限： 是一种数学近似值，用于表示(1)弱引力(2)与光速相比运动缓慢的物体(3)缓慢改变引力场的物理系统，在这种情况下，牛顿万有引力定律可以用于获得准确的值， 在牛顿极限中，时空近似平坦，并且可以在有限距离内使用闵可夫斯基度量，
16.测地线方程2和 度规张量分量g的00分量.：
$$ 牛顿引力由两个方程组成:一个高速我们我们物质如何对引力做出反应，$$$$另一个告诉我们物质如何产生引力。 $$$$ 第一个方程是从牛顿第二定律中得到：$$$$ a = -\nabla\phi 表示加速度为a, \phi为力 （加速度为力的负梯度？？） $$$$ 测地线方程描述的是仅受中立作用的粒子的世界线，因此我们的目标是证明，$$$$ 在牛顿极限的情况下，测地线方程可简化成第一个牛顿重力方程 $$$$ 测地线方程为：$$$$ 0 = \frac{d^2 x^\beta}{d\tau ^2 }+\Gamma_{\mu \nu}^\beta \frac{\partial x^\mu}{\partial \tau} \frac{\partial x^ \nu}{\partial \tau} $$$$ 因为粒子运动缓慢所以时间分量支配了其他分量， $$$$ 因此我们可以近似的得出： 由于: \frac{dx^i}{d\tau } \ll \frac{dt}{d\tau} $$$$ 0 = \frac{d^2 x^\beta}{d\tau ^2 }+\Gamma_{\mu \nu}^\beta \frac{\partial x^\mu}{\partial \tau} \frac{\partial x^ \nu}{\partial \tau} $$$$≈ \frac{d^2 x^\beta}{d\tau ^2 }+\Gamma_{00}^\beta (\frac{cdt}{d\tau} )^2 \frac{d^2 x^i}{d\tau ^2 }+\Gamma_{00}^i (\frac{cdt}{d\tau })^2= 0 公式1 $$$$ 前面提到Christoffel符号： \Gamma_{00}^i=\frac{1}{2} g^{ij} (g_{j0,0}+g_{j0,0}−g_{00,j} ) $$$$ 但是因为引力场很弱所以，所以时间导数g_{j0,0} 为0，$$$$ 因此，可以将christoffel 符号简化成 \Gamma_{00}^i=\frac{1}{2} g^{ij} g_{00,j} $$$$ 如果引力场足够弱，时空仅略微从狭义相对论的闵可夫斯基空间变形，$$$$ 我们可以考虑时空度量为一个小的扰动闵可夫斯基规度\eta_{\mu \nu} $$$$ 则g_{\mu \nu}=\eta_{\mu \nu}+h_{\mu \nu} , |{\mu \nu}|\ll 1 $$$$ g_{00,j}={(\eta_{00}+h_{00} )}_{,j} g_{00,j}= {\eta_{00}}_{,j} +h_{00,j}=0+h_{00,j} $$$$ 由于\eta_{00}=1 g_{00,j}=h_{00,j} $$$$ 在此步骤中，公式1变成: \frac{d^2 x^i}{d\tau ^2 }=− \Gamma_{00}^i \frac{cdt}{d\tau }^2 $$$$ \frac{d^2 x^i}{d\tau ^2 }= \frac{1}{2} g^{ij} h_{00,j} (\frac{cdt}{d\tau} )^2 $$$$ 如果我们定义 g^{ij}=\eta^{ij} −h^{ij}, 我们发现g^\mu\sigma g_{\sigma v}=\delta_ \nu^\mu, $$$$ 于是我们得到： $$$$ \frac{d^2 x^i}{d\tau ^2 }= \frac{1}{2} \eta^{ij} h_{00,j} (\frac{cdt}{d\tau} )^2 (\eta来度量g) $$$$ \eta^{ii}=−1, 其他值为0，所以： $$$$ \frac{d^2 x^i}{d\tau ^2 }= −\frac{1}{2} h_{00,i} (\frac{cdt}{d\tau} )^2 $$$$ 将左侧的导数从\tau 改成t，为此，首先将上式中的i替换成0: $$$$\frac{c^2 d^2 t}{d\tau ^2 }=−\frac{1}{2} h_{00},0 \frac{cdt}{d\tau}^2 $$$$ c^2 \frac{d^2 t}{d\tau ^2 }=0 （因为重力场是常数 h_{00},0=0) $$$$ \frac{d^2 t}{d\tau ^2 }=0, $$$$ 有了这个结果，我们仍然需要对\tau 使用偏导数 $$$$ \frac{d^2 x^i}{d\tau ^2 }=\frac{d}{d\tau} \frac{dx^i}{d\tau} =\frac{d}{d\tau} (\frac{dt}{d\tau} \frac{dx^i}{dt}) = $$$$ \frac{d}{d\tau} (\frac{d}{d\tau} \frac{dx^i}{d\tau} )+\frac{dx^i}{dt} (\frac{d}{d\tau} \frac{dx^i}{dt})=\frac{dt}{d\tau } (\frac{dt}{d\tau} \frac{d}{d\tau } \frac{dx^i}{d\tau})+\frac{dx^i}{dt} (\frac{d}{d\tau} \frac{dt}{d\tau} ) $$$$=(\frac{dt}{d\tau} )^2 (\frac{d^2 x^i}{dt^2} )+\frac{dx^i}{dt} (\frac{d^2 t}{d\tau ^2 }) $$$$ 所以 \frac{d^2 x^i}{d\tau ^2 }=(\frac{dt}{d\tau} )^2 (\frac{d^2 x^i}{dt^2} ) （制造出对t的偏导数） $$$$ 因为(\frac{d^2 t}{d\tau ^2} ) =0 $$$$ 那么 \frac{d^2 x^i}{d\tau ^2 }=(\frac{dt}{d\tau} )^2 (\frac{d^2 x^i}{dt^2} )=−\frac{1}{2} h_{00,i} (\frac{cdt}{d\tau} )^2 $$$$ \frac{d^2 x^i}{dt^2}=−\frac{c^2}{2} h(00,i) $$$$ 牛顿力学为\frac{d^2 x^i}{dt^2} =−\partial \frac{\phi}{\partial x^i} (\phi为牛顿引力势 ，可以把这个公式看成已知) $$$$ 最后，以向量形式表示: \frac{d^2 r}{dt^2} =−grad\phi 其中\phi=\frac{c^2 h_{00}}{2} $$$$ 备注： 将度规标准分量g_{00} 表示为： g_{00} =\eta_{00}+h_{00}=1+\frac{2\phi}{c^2} （就是将用\phi表示h_{00}） $$$$ 我们可以看到左侧的公制张量 与右侧的重力势之间存在直接联系 $$$$ 我们可计算地球上的h_{00} 可以得到其为10^(−9) 这意味着因重力场 $$$$而对于minkowski度量的偏差可以忽略不计。 $$
17.爱因斯坦的猜想与验证：
$$ 爱因斯坦张量微分形式比较，于是推测出 G_{\mu \nu}= \kappa T_{\mu \nu} $$$$ 爱因斯坦想要构造一个等式，左边是时空的弯曲情况，右边是时空中的能量。 $$$$ 可以看到已经有了T_{\mu \nu} 现在需要左边描述时空弯曲的程度， $$$$于是 G_{\mu \nu} 需要 满足 是一个二维张量就像T_{\mu \nu} 一样， $$$$ 应该描述使用曲率，并且相对黎曼张量是线性的，具有零散度，最后再牛顿极限下，$$$$应该减小到4\pi G\rho $$$$ 爱因斯坦的第一次尝试使用了Ricci张量作为解决方案: $$$$ 由 在牛顿极限下的测地线方程：\frac{d^2 r}{dt^2 }=−grad\phi 其中\phi=\frac{c^2 h_{00}}{2} $$$$ 可得到 将公制张量的分量 g_{00} 可以表示成 g_{00}= \eta_{00}+h_{00} =1+\frac{2\phi}{c^2} $$$$ 前面说到g_{\mu \nu} 非常重要就在于 它包含了所有关于时空和引力场的信息 ， $$$$ 所以方程中应该也是类似g_{uv} 这种二阶导数的形式， $$$$ 爱因斯坦先考虑的是使用黎曼张量的收缩形式：Ricci张量： $$$$ 前面得到的黎曼张量的表达式为： $$$$ R_{abc}^d= \Gamma_{ac,b}^d −\Gamma_{ab,c}^d+\Gamma_{ac}^e \Gamma_{be}^d −\Gamma_{ab}^e \Gamma_{ce}^d $$$$ 在静态和弱引力场的情况下，只有 R_{00} 有贡献： $$$$ R_{00}=R_{00}a^\alpha =\Gamma_{0\alpha ,0}^\alpha −\Gamma_{00,a}^\alpha +O(h2 ) = −\Gamma_{00,i}^i $$$$ 前面提到牛顿极限下，$$$$ \Gamma_{00}^i=\frac{1}{2} g^{ij} (g_{j0,0}+g_{j0,0}−g_{00,j} ) （\eta^{ii}=−1 ？ 其他值为0，所以：） $$$$ \Gamma_{00}^i=−\frac{1}{2} g^{ij} g_(00,j) , 因为g^{ij} 受\eta^{ij} 支配，并且 g_{00,j}=h_{00,j} $$$$ 所以\Gamma_{00}^i=−\frac{1}{2} \eta ^{ij} h_{00,j}=\frac{1}{2} \sigma _j^i h_{00,j} $$$$ \Gamma_{00,i}^i=\frac{1}{2} \sigma _j^i h_{00,ij} $$$$ R_{00}= −\Gamma_(00,i)^i =−\frac{1}{2} (\partial _1^2 h_{00}+\partial _2^2 h_{00}+\partial _3^2 h_{00} ) （00分量占主导） $$$$ =−\frac{1}{2} \nabla^2 h_{00}=−\frac{1}{c^2} \nabla^2 \phi $$$$ 爱因斯坦试图将R_{ab} 放在爱因斯坦方程的左侧，使他等于kT, $$$$ 他甚至用此解决了水星进动的问题，但是因为实际两边不等，所以最终爱因斯坦不得不放弃 $$$$ 爱因斯坦的第二次尝试: $$$$令等式左边为 G^{\mu \nu}=R^{\mu \nu} −\frac{1}{2} Rg^{\mu \nu} ， $$$$其中 R 是ricci标量,这种形式的G_{ab} 是对称的，显然描述了时空曲率， $$$$因此只需证明总导数为0就像右侧的能动张量的导数为0 $$$$(表明能量守恒),为此，我们从bianchi \quad identity 开始： $$$$ \nabla_\sigma R_{\alpha \beta}{\mu \nu}+\nabla_ \nu R_{\alpha \beta \sigma \mu}+\nabla_\mu R_{\alpha \beta \nu\sigma} =0 $$$$ 每项都乘上 g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} (因为 g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} 的导数为0 ，所以可以把它们作为常数放到\nabla中) $$$$ \nabla_\sigma g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} R_{\alpha \beta}{\mu \nu}+\nabla_ \nu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} R_{\alpha \beta \sigma \mu}+\nabla_\mu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} R_{\alpha \beta \nu\sigma} =0 $$$$ \nabla_\sigma g^{\gamma\sigma} R+\nabla_ \nu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} R_{\alpha \beta\sigma \mu}+\nabla_\mu g^{\gamma\sigma} g^{\alpha \mu } g^{\beta \nu} R_{\alpha \beta \nu\sigma} =0 $$$$ \nabla_\sigma g^{\gamma\sigma} R+\nabla_ \nu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu }R_{\sigma \mu\alpha \beta}+\nabla_\mu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} R_ {\nu\sigma \alpha \beta}=0 $$$$ \nabla_\sigma g^{\gamma\sigma} R−\nabla_ \nu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} R_{\mu\sigma \alpha \beta}+\nabla_\mu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} R_ {\nu\sigma \alpha \beta}=0 $$$$ Ricci 张量可以表示成 R^{\mu \nu} =g^{\mu\beta} g^{v\sigma} r_{\beta\sigma } 上式改写为： $$$$ \nabla_\sigma g^{\gamma\sigma} R−\nabla_ \nu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu } R_{\mu\sigma \alpha \beta}+\nabla_\mu g^{\gamma\sigma} g^{\alpha \mu} g^{\beta \nu} R_{\nu\sigma \alpha \beta}=0 $$$$ \nabla_\sigma g^{\gamma\sigma} R−\nabla_ \nu g^{\gamma\sigma} g^{\beta \nu} g^{\beta \nu} R_{\sigma \beta } +\nabla_\mu g^{\gamma\sigma} g^{\alpha \mu} R_{\sigma \alpha} =0 $$$$ \nabla_\sigma g^{\gamma\sigma} R −\nabla_ \nu R^{\gamma \nu}−\nabla_\mu R^{\gamma\mu}=0 $$$$ \nabla_\sigma g^{\gamma\sigma} R −\nabla_\sigma R^{\gamma\sigma} −\nabla_\mu R^{\gamma\sigma} =0 $$$$ \nabla_\sigma g^{\gamma\sigma} R −2\nabla_\sigma R^{\gamma\sigma} =0 $$$$ \nabla_\sigma (g^{\gamma\sigma} R −2R^{\gamma\sigma} )=0 $$$$ \nabla_\sigma (\frac{g^{\gamma\sigma} R}{2} −R^{\gamma\sigma} )=0 $$$$ 我们证明了爱因斯坦张量的散度为0 ，所以我们已经找到了爱因斯坦方程的左边部分 $$
18.常量K：
$$ 前面我们已经得到了爱因斯坦方程是G^{\mu \nu}=kT^{\mu \nu}，需要确定常数k $$$$ 谓词需要证明对于静态和弱重力场下(牛顿极限)下，爱因斯坦方程可以简化成牛顿重力定律。 $$$$ 首先使用不同形式的爱因斯坦方程，在计算时能够更加容易一些。 $$$$ R^{\mu \nu}−\frac{1}{2} g^{\mu \nu} R=kT^{\mu \nu} $$$$ 两边同时乘g_{\mu \nu} 可以得到： $$$$ g_{\mu \nu} R^{\mu \nu} −\frac{1}{2} g_{\mu \nu} g^{\mu \nu} R=kg_{\mu \nu} T^{\mu \nu} $$$$ g_{\mu \nu} R^{\mu \nu}=R 并且 g_{\mu \nu} T^{\mu \nu}=T 所以： $$$$ R−\frac{1}{2} Rg_{\mu \nu} g^{\mu \nu}=kT $$$$ 因为g^{\mu \nu} 是g_{\mu \nu} 的逆 ，他们的乘积是I，因为使用contraction $$$$所以 g_{\mu \nu} g^{\mu \nu}=\delta_ \nu^ \nu =1+1+1+1=4 $$$$ 因此 R −\frac{1}{2} R× 4=kT $$$$ R = -kT ， 将它代入到原式中 $$$$ R^{\mu \nu} −\frac{1}{2} g^{\mu \nu}× (−kT )= kT^{\mu \nu} $$$$ R^{\mu \nu}+\frac{1}{2} kg^{\mu \nu} T= kT^{\mu \nu} $$$$ R^{\mu \nu}=k(T^{\mu \nu} −\frac{1}{2} g^{\mu \nu} T) $$$$ 前面我们证明了分量R_{00}≈−\frac{1}{c^2} \nabla^2 \phi $$$$ 但是我们可以指出，当度量从固定时空减少到度量\eta时 ，我们可以通过Ricci张量的定义得到 $$$$ R^{00}=g^0\mu g^0 \nu R_{\mu \nu}≈\eta^0\mu \eta^0 \nu R_{\mu \nu} =(−1)(−1) R_{00} =R_{00} $$$$ R^{00}≈−\frac{1}{c^2} \nabla^2 \phi $$$$ 牛顿限制也意味着 T^{00}=\rho c^2 （因为T = \rho u^u u^ \nu 并且 u^i \ll u^0=c） $$$$ T = g_{\mu \nu} T^{\mu \nu}≈g_{00} T^{00}≈\eta_{00} T^{00}=T^{00} =\rho c^2 $$$$ K再代入R^{\mu \nu}=k(T^{\mu \nu} \frac{1}{2} g^{\mu \nu} T) $$$$ R^{00}=k(T^{00} −\frac{1}{2} \eta^{00} T) $$$$ −\frac{4\pi G\rho }{c^2} =K(\rho c^2 −\frac{1}{2}×1×\rho c^2 ) $$$$ −\frac{4\pi G\rho }{c^2} =\frac{1}{2} K\rho c^2 $$$$ K = -\frac{8\pi G}{c^4} $$$$ 终于我们得到了完整的爱因斯坦方程R_{\mu \nu} −\frac{1}{2} g_{\mu \nu} R=−8\pi GT_{\mu \nu} $$
当完全写出时，EFE（爱因斯坦场方程）由10个耦合的非线性二阶偏微分方程组成， 本来是16个，前面讲到对于对称的进行了合并，实际上解爱因斯坦方程就是在解度规张量g，所以相应的方程数为未知数的个数.文章末尾会有github上一个repo提供了对解的验证，对于正确的解会提供对应的方程中特殊部分的值如黎曼张量等
本文提及的场方程仅仅是最基础版本的，后来爱因斯坦又在方程中加入了一项宇宙常数，因为他认为宇宙是静态的，而原始方程解出来是确实膨胀的，所以加了一项常数抵消，但是哈勃证明了他的错误，但宇宙常数最后还是被留下来，用来表示宇宙中存在的暗物质和暗能量
结语
可以看到关于场方程的推导数学上并不涉及非常难的部分，但是对于非物理、数学专业有太多的概念需要了解， 同时我们也看到，爱因斯坦是通过牛顿极限来得到场方程的形式和常数K的值，是将其特殊化求解的方法，所以牛顿体系和爱因斯坦的相对论是完全兼容的，爱因斯坦描述了在任意参照系下的能量时空之间的关系，包含了其特例牛顿的体系。 爱因斯坦在构造场方程的时候遇到的最大困难在于等式的左边，用来描述时空的曲率，使用了黎曼张量的收缩形式，但使用 R^μν−1/2 g^μν R 感觉一般人也不会想到吧。 爱因斯坦在给别人的信中说“推导狭义相对论花了我好几天呢”，果然厉害的人喜欢装逼， 但是，广义相对论花了爱因斯坦10年的时间中间还走进了死胡同。 在这之中还有许多数学大家的帮忙，所以难怪我们看不懂100年前人的结晶。
参考资料
参考文章
爱因斯坦场方程推导过程的逻辑梳理:
https://wenku.baidu.com/view/5c305d42ce84b9d528ea81c758f5f61fb7362863.html
本文主要来源，强推，有非常完整的描述:
Einsteinrelativelyeasy http://einsteinrelativelyeasy.com/index.php
另一些类似的网站:
还有对伽利略和麦克斯韦的成果介绍:
http://www.einstein-online.info/spotlights/gr.html
http://newt.phys.unsw.edu.au/einsteinlight/
对偶空间和对偶基的理解:
https://blog.csdn.net/sumx2015/article/details/78909635
维基百科
百度百科
参考视频:
场方程推导:
协变与逆变向量:
张量的基本运算(包括收缩):
代码：
spacetimeengine:
github上一个python写的repo提供了对 解的验证，对于正确的解会提供对应的方程中特殊部分的值如黎曼张量等
solution_array = [ metric, coordinate_set, index_config, cosmological_constant ]
仅需要输入解数组,其中包括了度规和，坐标表示，上下标配置,和宇宙常数即可进行验证
https://github.com/spacetimeengineer/spacetimeengine
this is another part
原文地址：https://www.cnblogs.com/captainhook/p/11784432.html