1 引言

织法理论[1]认为，我们宇宙的时空几何可由一个三维四价张量网络即立方金刚石（3C）网描述，其沿堆垛轴作正交投影到二维平面可得三价（蜂窝）网络。网络的基本单元 a_P = 2√(ln8)·l_P（[1] 推论 3），是一个二进制的 3 比特（8 种状态）的量子纠缠关系。在整个时空，这样的纠缠关系生生灭灭。因荷守恒的约束，某个单元纠缠关系状态的改变，会瞬时扰动整个网络，导致其他单元的状态同步改变，这就是量子纠缠现象。在时空的局部区域，基本单元之间因受规范力约束，可形成结构特定的、处于锁死或受束缚状态的纠缠关系；物质由此演生[2,3]。本文的核心论点是：引力的本质不是一种力，而是时空、物质相互之间 a_P 场密度的相对梯度。本文从时空几何的底层结构出发，推导织法理论下的引力方程与统一度规（密度档、剪切档、拖曳档），给出点源解与强场行为，并以太阳系天体表面的重力加速度、中子星表面的红移与爱因斯坦方程的结果作对照，从而检验方程的准确性。

---

2 引力的本质

广义相对论认为，引力不是一种力，而是时空的弯曲[4]。物质和能量使周围的时空几何发生畸变，物体在弯曲的时空里沿测地线运动，看上去就像彼此吸引。光同样沿测地线传播，因此光线经过大质量天体附近时会发生偏折；又因为弯曲时空中各处的时钟走得快慢不同，光在经过引力场的那段路上会变慢，到达时出现延迟，频率也随之降低。这套几何描述在太阳系尺度的检验中与观测高度吻合[5]，引力就是时空的弯曲也因此成为今天关于引力的标准图像。但广义相对论没有回答：时空为何会弯曲，引力的来源到底是什么，引力波为何以光速传播？在普朗克尺度，广义相对论更与量子力学产生冲突[6]。

我们认为，宇宙的本体是量子纠缠关系的叠加，而非物质。a_P 是一个基本的量子纠缠关系单元，它由二进制的 3 比特表达，是 8 种状态的叠加；故时空几何和物质均通过 a_P 的编织方式表达。时空并不是物质的背景，物质也并非独立于时空。物质是时空局部的 a_P 网络因荷守恒和规范力约束处于锁死或受束缚状态形成[2,3]，真空是处于荷守恒约束的 a_P 关系自由的量子涨落状态的时空。物质的演生过程，会锁定或束缚四周时空 a_P 单元的量子状态。物质要存在于时空中，则必须与时空建立稳定的纠缠关系，而非脱离于时空。引力是时空与物质的固有属性，而非物质单独的属性。天体之所以存在于时空某个相对规律的轨道，是因为其受整个时空几何引力的牵引而达到的与时空，其他天体之间形成的引力平衡状态。

引力来源于 a_P 单元的计数；但在不同能量、动量场景下读出的度规不同，故须按场景各自构建度规。物质锁死或束缚的 a_P 数是可加的整数 N，物质的量就是这个计数：质量 m = N/a_P，单位质量 m₁ = 1/a_P。物质处的 a_P 场密度用 ρ 表示，引力即 a_P 场密度的相对梯度。物质为了与周围时空建立稳定的纠缠关系，其自身锁死或束缚的 a_P 单元同时会决定周围时空 a_P 的量子状态，形成 a_P 场，其强度由其质量决定。物质之间时空的 a_P 场由各自的质量叠加决定。苹果之所以会落地，是因为其脱离于树枝后，必须与周围时空建立稳定纠缠关系，但其周围时空的 a_P 单元更多的被地球束缚，其沿着测地线在空中每个位置都无法建立稳定纠缠关系，只有在落到地球表面后，才与周围时空建立了稳定的纠缠关系。在真空中，因为苹果四周的引力均匀，因此任何位置其都可以与时空建立稳定的纠缠关系，故表现为悬浮在真空中。

设一个鸡蛋与一块石头，体积相同、外形一致，则两者体积内的 a_P 数量完全相同；那么，石头为什么比鸡蛋重？差别在于两者体积内处于锁死状态的 a_P 数量不同。每个 a_P 单元都携带一份 ln8 的熵；要让一个 a_P 进入锁死状态，就必须付出同等的一份能量——色场正是这份能量（质量）的来源。在氢原子的原子核内，胶子将夸克束缚在一起，形成一个六元环（由 6 条 a_P 边闭合而成），每条边的量子状态都处于锁死的稳定状态[3]。原子核内 a_P 的量子状态决定核外的量子涨落状态，它把电子的激发与网络的其余部分关联起来，使电子不断被熄灭，又在别处重新点亮[3]。因此，处于锁死状态的 a_P 数量越多，也意味着 a_P 的可用态数越少，其对电子的束缚也就越强。原子钟数的就是这样一条被荷守恒约束的激发模式——铯-133 的超精细跃迁，数满 9 192 631 770 个周期为 1 秒。而核锁死的 a_P 越多，可用态越少，走完一个周期所需的坐标时间越长，钟就越慢，即 dτ/dt = 1/ρ。

通过织法公式（[1] 推论 4），引力常数与光速由同一个 a_P 锁定：G = a_P²/(4c·ln8)。G 与 c 是同一块织物的两个读数。几何模型上 a_P = 2√(ln8) 就是 2.884 个 lP，真实物理里它是一个量子纠缠关系。c 为光速，其原始含义是纯几何量——熵–面积耦合的饱和强度（[1] 推论 4）。当 c 趋近热力学饱和极限（c∞ = 1，[1] 定理 3）时，该几何量即等同于光速，即信息在时空几何的任意方向都能以单位时间走完一个 a_P 的上限传播。引力也是通过 a_P 传播，故引力波的速度等于光速。大质量天体附近的时间会变慢，是因为光掠过时，大质量天体周围时空因弯曲导致 a_P 场密度大，光相对于真空，需要跨越更多的 a_P 单元，这就是空间导致的时间变慢。

3 引力方程

3.1 引力方程组

单位制（本节起全部推导）：自然单位制——ħ = c = 1，长度以普朗克长度 l_P 计，时间以 l_P/c 计，质量以 m₁ = ħ/(a_P c) 计；网的基本单元是 a_P = 2√(ln8)·l_P（量值 2.884，[1] 推论 3）。§3 起各节的公式一律按此单位制书写，式中不再夹带换算常数；需要国际单位时，见 §3.5 的换算表。

本节建立的引力理论写成下面一组封闭的方程，分三部分：场方程（定场）、源规则（物质到源）、读数规则（场到几何与观测量），另加常数定义与张量扇区（引力波）。式中 u ≡ ln ρ 是场变量；ρ(x) 是 a_P 场密度（= 该处单位体积内锁死和受束缚状态的 a_P 数量 / 真空 a_P 场密度，真空处 ρ = 1，相对真空，锁死只增不减，故 ρ ≥ 1）；n_lock(x) 是单位体积内锁死和受束缚状态的 a_P 数量（即源的计数）；N_lock 是某处（一个源）锁死和受束缚状态的总 a_P 数（N_lock = Σ N_p，见 §4.2）；ρ_m 是质量密度；g_μν 是几何层的度规读数，δ_ij 是三维平直空间度规。

（1）场方程（定场）

∇²(ln ρ) = −(4π/a_P)·n_lock                                 (3.1)

（2）源规则（物质到源）

ρ_m = n_lock/a_P                                            (3.2)

每一条锁死或束缚状态的关系贡献一份质量 m₁（自然单位下 m₁ = 1/a_P），故 ρ_m = n_lock/a_P（源规则的微观依据见 §4.2）。代入式 (3.1)，得几何层的泊松形式

∇²(ln ρ) = −(4πG/c²)·ρ_m                                    (3.3)

（3）读数规则（场到几何与观测量）

ds² = −(1/ρ²) dt² + (ρ² δ_ij + h_ij^TT)(dx^i + N^i dt)(dx^j + N^j dt)        (3.4)

a = c²∇(ln ρ),      1 + z = ρ,      dτ/dt = 1/ρ             (3.5)

式 (3.4) 是统一度规：密度档给出 lapse 1/ρ 与空间标度 ρ²δ_ij，剪切档给出 h_ij^TT（引力波），拖曳档给出 N^i（转动/平动）；静态（N^i = 0、h_ij^TT = 0）即退化为 g₀₀ = −1/ρ²、g_ij = ρ²δ_ij。式 (3.5) 依次给出加速度、红移与时钟速率；各档的推导见 §6。

（4）常数与点源

G = a_P²/(4c·ln8)                                           (3.6)

ρ = e^{μ/r},      μ = N_lock/a_P,      g(r) = GM/r²         (3.7)

G 为引力常数，其量值由 a_P 定出（[1] 推论 4），自然单位下 G = 1。式 (3.7) 是式 (3.1) 到 (3.3) 的点源解（μ 为源的幅度，自然单位下即源的质量 M；弱场、强场同一表达式；国际单位写法见 §3.5）。

（5）张量扇区（引力波）

(∂_t² − ∇²)h_ij^TT = 0                                      (3.8)

真空中无源传播、两个偏振、波速 c；源的耦合常数即 G（推导见 §5.1）。

用法：先数源——由物质分布写出 n_lock，或用式 (3.2) 由 ρ_m 换算；再解场——解泊松方程 (3.1) 得 u = ln ρ；后读数——用式 (3.4)、(3.5) 得度规、加速度与红移。这组方程严格线性、互不耦合：给定源只需解一个泊松方程，再按读数规则读出几何与观测量，无需联立迭代。以下各节逐条导出这组方程。

3.2 a_P 场密度 ρ

面积律 S = (L−1)·ln8（[1] 定理 1）把“几何量”与“计数”接通（纠缠熵按面积标度这一普适行为的既有背景见 [7–9]）：面积 A 对应的不是连续的长度平方，而是 L 个 a_P——每个 a_P 携带一份面积 a_P²（[1] 推论 3）：

A = L·a_P²                                                  (3.9)

于是定义 a_P 场密度：

ρ(x) ≡ （单位体积内锁死和受束缚状态的 a_P 数量）/（真空 a_P 场密度）(3.10)

每个 a_P 携带 ln8 的熵，是 8 种量子状态的叠加。a_P 处于锁死状态，即其量子状态完全不自由；锁死状态直接决定了束缚状态，即其量子状态自由度有限。真空里没有锁死的关系——关系全部自由、可自由重排——因此真空场密度取 1。相对真空，锁死只增不减，故 ρ ≥ 1。锁死状态的 a_P 数量越大，场密度越大，引力越大。

引力的本体（本文的核心命题）：

引力不是一种力，而是 a_P 场密度的相对梯度。

加速度取这个比值的相对梯度。场变量 u ≡ ln ρ 无量纲（式 (3.13)），故 ∇u 的量纲是长度的倒数，要配成加速度只能乘 c²；这也就是几何层的读法——u = −Φ/c²、a = −∇Φ（§3.3（4））：

a = c²∇(ln ρ) = c²(∇ρ)/ρ                                    (3.11)

弱场下它退化为牛顿引力（§3.4）。

3.3 场方程

网形：三维网络取立方金刚石（3C）网——每个节点由四条键连到四个近邻，四条键的方向指向正四面体的四个顶点，两两夹角 arccos(−1/3) = 109.4712°；它由两套面心立方（fcc）子晶格沿体对角线错开四分之一互穿而成，键长取 a_P；其晶体学标准结构见 [10]。沿堆垛轴作正交投影，即得三价（蜂窝）网——面积律用的就是这张平面三价网。三维网之所以取立方金刚石（3C）网，是两个条件一起定出来的：基本单元有 4 条关系，4 条键必须等价，键长同为 a_P；4 条键的方向必须对称，否则网不是各向同性的，而键角的唯一对称解是正四面体配位，即 arccos(−1/3) = 109.4712°。满足这两条、且整网宏观各向同性的网就是立方金刚石（3C）网。它的最短闭合回路是六元环：6 个节点、6 条键（最短回路长度 = 6，且不含四元环，已逐条数值核对）。六元环的 6 个节点不严格共面（沿体对角线看呈椅式起伏）。时空几何的网形讨论，详见《织法公式：时空几何的涌现机制》[1]。

图 1　三维四价张量网络、即立方金刚石（3C）网的局部细节。

基本单元 a_P = 2√(ln8)·l_P（[1] 推论 3）即

a_P² = 4 ln8                                                (3.12)

（1）场变量。一个节点就是 2³ = 8 个态（3 比特），相互独立的键的贡献相乘 ⇒ ln8 对键可加。源是锁死和束缚状态的纠缠关系数 N——整数计数、可加。可加源的多极叠加要求场量可加，局部标量中唯一可加的是

u ≡ ln ρ                                                    (3.13)

其中 ρ 是 a_P 场密度（§3.2）；关系被锁死一份，可用态数就摊薄一份，故 ρ ≥ 1。

（2）四腿展开（网 → 连续）。四面体四腿的单位矢满足

Σ_a nᵃ = 0,      Σ_a nᵃnᵇ = (4/3)δᵃᵇ                       (3.14)

把四条腿写成 (1,1,1)、(1,−1,−1)、(−1,1,−1)、(−1,−1,1) 再归一化到单位长：对角元 (1+1+1+1)/3 = 4/3，非对角因正负相消为零。沿各腿把 u 展开到 O(a_P²)，线性项由 Σ_a nᵃ = 0 整体消去：

u(x + a_P nᵃ) − u(x) ≈ a_P (nᵃ·∇)u + (a_P²/2)(nᵃ·∇)²u
Σ_a (nᵃ·∇)²u = Σ_a nᵃ_α nᵃ_β ∂_α∂_β u = (4/3)∇²u
⇒ Σ_j (u_j − u_i) = (a_P²/2)·(4/3)∇²u = (2/3)a_P²∇²u        (3.15)

（式 (3.15) 的展开在三维四价网上；面积律在平面三价网上研究；两层只在涌现度规处对接（§6）。）

（3）二次能量。弹性应变取相对变化。此处出现的密度按体积计数：ρ_体 ≡ 每单位坐标体积内的 a_P 条数比（即 a_P 场密度，真空取 1；与本文其余处的线计数 ρ 的关系为 ρ_体 = ρ³，见 §6）。于是单位关系摊到的体积 ∝ 1/ρ_体 ⇒ 坐标格距比 b/b₀ = ρ_体^{−1/3}，故相对应变 δb/b = −(1/3)ln ρ_体。能量在应变上二次，以 u = ln ρ 写出仍是二次型，其刚度系数就是式 (3.15) 的系数：

W = (κ_reg/2)∫(∇u)²,      κ_reg = (2/3)a_P²                 (3.16)

（归一约定：节点上的梯度算符由四腿展开 (3.15) 给出，即 Σ_j(u_j − u_i) = (2/3)a_P²∇²u；取该算符系数为刚度，κ_reg = (2/3)a_P²（数值 5.545177444）。）格点能量按节点求和 W = (κ_reg/2)Σ_i(∇u)²·v₁ 与其连续极限 (κ_reg/2)∫(∇u)² 出自同一关系计数（v₁ 见式 (3.22)）。

（4）源项与场方程。源可加 ⇒ 源项只能线性：W_src = −∫σ·u。变分 δW/δu = 0 给出

∇²u = −σ/κ_reg                                              (3.17)

源强系数由引力常数 G 定出。几何层的泊松读法是 ∇²Φ = 4πGρ_m，场变量与它的关系是 u = −Φ/c²（等价地 a = c²∇u = −∇Φ，式 (3.11)）；代入得

∇²u = −(4πG/c²)·ρ_m                                        (3.18)

再用源规则 ρ_m = n_lock/a_P（§4.2）把源写成本层的计数：

∇²u = −(4πG/(c²a_P))·n_lock                                (3.19)

G/c² 是纯数：由 [1] 推论 4 的 G = a_P²/(4c·ln8) 与式 (3.12)，

G/c² = a_P²/(4c³·ln8) = l_P²/c³ = 1      （自然单位）        (3.20)

于是 4πG/(c²a_P) = 4π/a_P——源强系数就是 G（换成质量密度源写即式 (3.28)）。又 σ = σ₁·n_lock，其中 σ₁ = (4π/a_P)·κ_reg = (8π/3)a_P = 24.161392390530，故

∇²u = −(4π/a_P)·n_lock                                     (3.21)

（5）节点体积与格点形式。金刚石晶胞的棱 a = 4a_P/√3、每个晶胞 8 个节点，故每个节点对应的体积

v₁ = a³/8 = 8a_P³/(3√3)                                    (3.22)

节点 i 上的关系数 N_i = n_lock·v₁；把式 (3.21) 的 ∇²u 代入式 (3.15)：

Σ_j (u_j − u_i) = −q·N_i                                   (3.23)
q = (8π/3)a_P/v₁ = π√3/(4 ln8)      （q = 0.654189836987）

q 是把连续系数落成格点计数（每条关系、每份体积）的换算因子。反过来把 N_i = n_lock·v₁ 代入式 (3.23)、并用式 (3.12) 与式 (3.22)：

∇²u = −(4q/√3)·a_P·n_lock = −(4π/a_P)·n_lock                (3.24)

即回到式 (3.21) 的系数（也就是 §3.1 的场方程 (3.1)）。但这是代数恒等式，不是独立检验：把 q 与 a_P² = 4 ln8 代入，(4q/√3)·a_P = (4/√3)·(π√3/(4 ln8))·a_P = π·a_P/ln8 = 2π/√(ln8) = 4π/a_P（数值 4.35719012284867）；它只说明格点形式与连续形式出自同一关系计数。PPN 参数由同一个 ρ 给出：时钟读数给 g₀₀ = −1/ρ²、几何给 g_ij = ρ²δ_ij，展开 ρ = 1 + U 即得 γ = β = 1（§6、§7.1）。

3.4 点源解

点源：一个格点上 N_lock 条不完全自由的关系，n_lock = N_lock·δ³(x)。用 ∇²(1/r) = −4πδ³(x) 对 ∇²u = −(4π/a_P)·N_lock·δ³(x) 直接积分：

ln ρ = N_lock/(a_P · r)      ⟹      ρ(r) = e^{μ/r},   μ = N_lock/a_P    (3.25)

换成观测语言：由源规则（§4.2）源幅度 μ = N_lock/a_P，而 N_lock/a_P 就是源的质量 M；自然单位下 G = 1（式 (3.20)），于是由式 (3.11) 读出加速度

g(r) = GM/r²                                                (3.26)

弱场、强场同一表达式；国际单位下 μ = GM/c²（式 (3.29)）。

3.5 国际单位制下的换算

本节之前的方程（§3.1–§3.4）一律写自然单位制（ħ = c = 1，长度以 l_P 计、时间以 l_P/c 计、质量以 m₁ = ħ/(a_P c) 计），式中不含换算常数。以下按量纲把 c 与 ħ 的幂补回，写成国际单位制——长度、时间、质量分别以 l_P、l_P/c、ħ/(a_P c) 折算。

加速度        a = c²∇(ln ρ) = c²(∇ρ)/ρ                         (3.27)
场方程        ∇²(ln ρ) = −(4πG/c²)·ρ_m                         (3.28)
点源          μ = GM/c²,   ρ = e^{μ/r},   g(r) = GM/r²        (3.29)
引力常数      G = a_P²c³/(4ħ·ln8) = l_P²c³/ħ                    (3.30)

各式依次是式 (3.11)、(3.21)、(3.25)–(3.26) 的国际单位写法；质量单位 m₁ = ħ/(a_P c)。[1] 推论 4 的自然单位写法为 G = a_P²/(4c·ln8)。

数值例（国际单位）：g = GM/R² 代入太阳得 ≈ 274.2 m/s²，代入地球得 ≈ 9.82 m/s²，与观测一致（自洽性验证，§7）。

---

4 源：a_P 量子状态不完全自由的关系数

4.1 物质

方程 ∇²(ln ρ) = −(4π/a_P)·n_lock 的源只有一个：n_lock。它是什么？

a_P 是时空几何的基本单元，每个 a_P 携带 ln8 的熵。因此 a_P 的量子纠缠程度有完全锁死、不完全自由和完全自由三种；而锁死状态直接决定了不完全自由的量子状态，换句话说，锁死状态的 a_P 会束缚其周围 a_P 的量子状态，使其处于不完全自由的状态。因此 n_lock 指的是单位体积内锁死状态和因此受束缚状态的 a_P 的计数。真空中，a_P 的量子状态完全自由。物质内部，a_P 则是处于锁死和不完全自由的状态，物质的质量主要由锁死状态的 a_P 贡献。 n_lock 不是能量、不是动量、不是任何外部输入——它是方程唯一的源。

4.2 源规则

n_lock 与质量之间只差唯一尺度 a_P：

ρ_m = n_lock / a_P                                          (4.1)

N_lock = Σ N_p,   N_p = m/m₁
（m 为质量，m₁ 为质量单位；与质量出自同一关系计数；N_p 为单个粒子的质量份数）

等价读法（与 a_P 场密度读法出自同一关系计数，[1] 推论 3）：质量 = m₁ × 质量份数；质量份数的比值由《质量比值的织法起源》[11] 的内部相空间体积读数给出（质量 ∝ 内部相空间体积），其绝对取值列为开放；粒子谱本身（谁在谱上、带什么荷与色）见《物质谱的织法结构》[3]。

这带来一个可核验的推论：两个粒子的质量比 = 它们的源幅度比 = 它们的质量份数之比。以质子/电子为例：

m_p/m_e = 1836.152673440（观测）                              (4.2)

即：引力方程的源，读出的正是观测质量谱；m_p/m_e 不进引力场方程本身，它只进“源”。

4.3 唯一尺度

全篇只出现一个基本量：a_P（以 l_P 计时为 2.884，见 §3.2）。它同时是

• 基本单元 a_P，即一个量子纠缠关系（3 比特、共 8 种状态）（[1] 推论 3）；
• n_lock ↔ 质量的换算率（§4.2）；
• 唯一质量单位 m₁ = 1/a_P。

---

5 张量扇区与标量扇区

网的弹性动力学在张量扇区与标量扇区上给出的结论完全相反——张量扇区管引力波，标量扇区管“力”。

5.1 张量扇区：(∂_t² − ∇²)h_ij^TT = 0

网的能量规则只有一条：能量 = −ln 态数（态数是纯数，形变下按 Jacobian 变换，§5.2）。把它在金刚石网上展开，得到两类贡献——键长项（态数只依赖键长）与体积项（态数 ∝ 形变 Jacobian det F）。键长项单独作用时 C11 = C12、C44 = 0，剪切出现零模 ⇒ 横支速度为零、无引力波；体积项的弹性张量比例是 (C11−C12) : C12 : C44 = 2 : 0 : 1。织法没有第二条计数规则，键长项不是独立项，它就是体积项的形变分量 ⇒ 弹性张量取纯体积型：

C12 = 0,      C11 = 2C44       （即 β/α = 1/2 那一点）        (5.1)

动力学矩阵。两类贡献写成标准的两项（价力场模型的键长–键角展开 [12]）——键长项（键长平方的偏差）与键角项（键角平方的偏差）；每个节点 4 条键、6 个键角，b′₁、b′₂ 为同一节点上两条形变后的键矢（b′ = b + Δu）：

W = (α/2) Σ_{键} (b′²/a_P² − 1)² + (β/2) Σ_{角} (b′₁·b′₂/a_P² + 1/3)²        (5.2)

（α、β 为键长项、键角项的刚度；织法只有一条计数规则 ⇒ 取体积型那一点 β/α = 1/2。）把式 (5.2) 展开到位移的二次项：键长项给出 (2α/a_P²)·Σ_键 (n·Δu)²，键角项给出 (β/2a_P²)·Σ_角 (n₁·Δu₂ + n₂·Δu₁)²（n 为单位键方向、Δu 为键或角两端的位移差；键角项的指标是交叉配对）。两者都是位移差的二次型，按 D(q) = Σ_R Φ(R)e^{iq·R}（Φ 为二次型的力常数矩阵，R 取到次近邻）写成 6×6 动力学矩阵并对角化，三个声学本征值为

v² = (2, 2, 4)      （[100]、[110]、[111] 三个方向完全一致）   (5.3)

二重简并应对两个横向位移，纵支对应 ∇·u。v² 的无量纲数值取决于键刚度与长度单位的归一（(2,2,4) 只是其中一个归一），物理内容是比值：横支二重简并、v_L²/v_T² = 2 ⇒ c_L/c_T = √2，且与方向无关（各向同性）。绝对标度由时标约定固定——把“一次 tick 走完一个 a_P”取作时间单位，横向（物理）支的波速就定为 c（[1] 定理 3 的热力学饱和极限），纵支为 √2 c。

两类贡献与弹性常数的对应：C11、C12、C44 是立方对称弹性的三个独立常数（Voigt 记法：C11 = ∂²w/∂ε₁²、C12 = ∂²w/∂ε₁∂ε₂、C44 = ∂²w/∂ε₄²，w 为应变能密度，ε₄ = 2ε_yz）。沿 [100] 方向分解，纯纵支 v_L² = C11/ρ_m、纯横支 v_T² = C44/ρ_m；代入式 (5.3) 的 (v_L², v_T²) = (4, 2) 得 C11 = 2C44，代入各向同性条件 C11 − C12 = 2C44 得 C12 = 0。各向异性因子（Zener 因子 [13]）A = 2C44/(C11−C12) = 1 ⇒ 宏观各向同性；ν = 0（无泊松比）、c_L/cT = √2 精确。同一能量按均匀小应变的二阶展开、并对胞内 A、B 两套子晶格的相对位移取极小（w = −H{εw}H_{ww}^{−1}ε），给出弹性常数之比 C11 : C12 : C44 = 2 : 0 : 1（纯体积项 (1/2)tr ε² 单独作用给出同一比例：C11 = 1、C12 = 0、C44 = 1/2），与它们的取值自洽。

连续极限（λ = C12 = 0）给出

∂_t²u_i = c_T²(∇²u_i + ∂_i∇·u)                              (5.4)

网位移场在几何层的读数是应变，即 h_ij := ∂_i u_j + ∂_j u_i（应变的二倍；桥在 §6）；代入式 (5.4) 得 ∂_t²h_ij = c_T²(∇²h_ij + ∂_i∂_j h)，取横无迹部分（∇·u = 0 即横支）：

(∂_t² − ∇²)h_ij^TT = 0                                      (5.5)

即：真空中的无源传播、两个偏振、波速 = c_T = c（一次 tick 走完一个 a_P 就是光速，[1] 定理 3 的热力学饱和极限）。纵向支恰好是纯规范（h_ij = ∂_i u_j + ∂_j u_i 就是坐标变换）⇒ 物理谱 = 2 个横偏振，与引力的观测自由度一致。

带源时，张量扇区的方程是

□ h_ij^TT = −κ·τ_ij^TT,      κ = 4π/a_P                      (5.6)

这里的 κ 与标量扇区式 (3.1) 的系数是同一个（也就是几何层的 G，§3.3）：源 τ_ij = n_lock·⟨n_i n_j⟩ 是锁死关系的方向张量（⟨·⟩ 表示对该处锁死键的方向取平均，n 为单位键方向），其迹 tr τ = n_lock 即标量（牛顿）扇区的源、无迹的横向部分给出引力波的源。张量扇区与标量扇区在关系层用同一套计数、同一个耦合，不需要额外的换算因子；纵向支为纯规范，物理谱只有 2 个横偏振。

为什么必须有体积项：若弹性只依赖键长（中心力），则 β = 0，剪切出现零模 ⇒ 无引力波——这已被引力波观测排除[14]。而“态数只随局域体积（Jacobian）变”这一条自然给出各向同性：即使把键长项当作独立项按任意权重 w 混入（C11 = 8+w、C12 = 8、C44 = w/2），A = 2C44/(C11−C12) 也严格等于 1，各向同性与权重无关；只有把非中心力项换成 Keating 键角刚度（C12 = 8 − 16β ≠ 0）才需要把 β/α 调到 1/2。

5.2 标量扇区：全阶严格线性

标量（“牛顿”）扇区不能用局域体积泛函。局域体积型全阶泛函（张量扇区用）由统计力学的读数规则给出：自由能 = −ln 态数（[1] 的熵–面积律就是态数的对数），而态数是局域体积的函数、形变下按 Jacobian 变换 ⇒ 态数 ∝ det F，故

W = Σ_i [ −ln det F_i + tr(F_i − I) ]                       (5.7)

Fi 由节点 i 的 4 条键矢构造，−ln det F 就是 −ln 态数；它的展开首项为 −tr(F−I)，加 tr(F−I) 正为抵消该项，使 W 在无应变处为零、最低阶为二次。它只含一阶导 F，故 Euler–Lagrange 方程 ∂_i(∂f/∂F\{ij}) = 0 退化为“∂f/∂F 为常数”的代数条件，结合无穷远处无应变（F → I）即得 F ≡ I，也就是 u′ ≡ 0（u′ = du/dr；数值残差在机器精度内为零）⇒ 局域体积泛函的真空必给 u = 常数 ⇒ 不产生静态牛顿场。所以它只管引力波，不管力。

标量扇区来自梯度泛函；而 u = ln ρ 是计数的对数、计数是乘性的 ⇒ 梯度泛函天然是二次型 ⇒ 方程严格线性。因此

∇²(ln ρ) = −(4π/a_P)·n_lock     （全阶严格成立）                    (5.8)
ρ = e^{μ/r}                     （全阶精确解，不是外插）

这一点正是织法与 GR 的真正分野：GR 把非线性放在方程里（泊松只是弱场线性化近似）；我们把非线性放在读数里（方程本身严格线性，非线性全部搬进了 ρ = e^{μ/r} 这一层）。

---

6 涌现度规：关系层到几何层的桥

关系层的场量是 a_P 场密度 ρ（场变量 u = ln ρ）与关系数的计数，几何层的语言是度规。两层只在涌现度规处对接一次——把 ρ 与网的形变读成时钟、长度与拖曳。读数只有两条动作：数 a_P（长度）与数 tick（时间）；a_P 是常数（键长 a_P 恒定、网不可压），两条读数都与坐标无关：

• 长度（线计数）：单位坐标长度内有 ρ/a_P 条关系，每条关系的本征长度恒为 a_P ⇒ 本征长度 dl = ρ·dx（等价地：单位坐标体积内有 ρ³/a_P³ 条、本征体积 = ρ³d³x）；
• 时间（态数摊薄）：一次 tick 走一个 a_P（时标约定，[1] 定理 3 的热力学饱和极限）；锁死占掉一份态 ⇒ 可用态数摊薄 ⇒ dτ/dt = 1/ρ。

两条相乘 ⇒ dτ·dl = dt·dx ⇒ (−g₀₀)·g_ij = 1 ⇒ 本征光速 = c 自动成立（§6.5）。据此，读数按网的局域形变分档给出：体积（迹）部分 → 时钟与空间标度（§6.1、§6.2），无迹（剪切）部分 → 度规的空间无迹部分（§6.3），流（拖曳）部分 → 度规的时间–空间交叉项（§6.4）。

6.1 时钟读数

一次 tick 走完一个 a_P（时标约定，[1] 定理 3 的热力学饱和极限）⇒ 局部钟的周期由可用态数决定：一个节点 8 个态，关系被锁死一份、可用态数就摊薄一份 ⇒ 单位坐标时间内的固有时（局部钟速率）∝ 1/ρ（ρ ≥ 1：锁死越多、钟越慢）。时间分量由固有时定义（dτ² = −g₀₀ dt²），故钟速率的读数就是

dτ/dt = 1/ρ      ⇒      g₀₀ = −(dτ/dt)² = −1/ρ²             (6.1)

（同一个 ρ 也给出力的读数：a = c²∇(ln ρ)，式 (3.11)、(3.27)——时钟与力出自同一个 ρ。）

6.2 长度读数

几何层的面积由关系数读出（[1] 推论 3：A = 关系数 × a_P²，式 (3.9)）⇒ 关系数的比 ρ 是空间度量的标量输入（无迹的剪切读数见 §6.3）。空间部分取共形平直、共形因子只依赖 ρ，即 g_ij = ρ^{2n}δ_ij；各向同性坐标（GR 弱场的标准坐标[5]）要求它与 g₀₀ 的一次项大小相同：

g_ij = (1 + 2nU + …)δ_ij,   g₀₀ = −(1 − 2U + …)   ⇒   n = 1
⇒      g_ij = ρ² δ_ij                                        (6.2)

（几何层面 g_ij 的整体标度本来就是坐标自由，各向同性坐标把它定死；空间部分不携带独立方程。）

6.3 剪切读数（张量半边）

网的局域形变按不可约分解，两半各归一处。迹（体积）部分改变单位体积内的关系数，也就是改变 ρ——这一半已由式 (6.1)、(6.2) 承担。无迹（剪切）部分一阶不改体积（det(I + h) = 1 + tr h + O(h²)、tr h = 0 ⇒ 不改 ρ），它只能通过方向读取不一致显形：同一段坐标间隔，沿不同方向数到的关系数不同 ⇒ 长度读数方向依赖。把“沿 n̂ 的长度读数”与各向同性值相比，偏差就是无迹形变在 n̂ 上的投影：

g_ij = ρ²(δ_ij + h_ij^TT)                                   (6.3)
h_ij^TT = [∂_i u_j + ∂_j u_i]^TT,  δ^{ij}h_ij^TT = 0,  ∂^i h_ij^TT = 0

u_i 是 §5.1 的网位移场，h_ij := ∂_i u_j + ∂_j u_i 是它的应变读数（应变的二倍）。三个性质各有出处：无迹 = 不碰体积（体积已归 ρ，两档互不污染）；横向 = 纵向支为纯规范（§5.1），物理谱只剩 2 个横偏振；对称 = 度规与键张量 Σ n_i n_j 都对称。探测器读的就是这个读数之比：ΔL/L = ½·h_ij^TT n̂^i n̂^j（n̂ 为探测器臂方向）——它是可直接测量的量，不是内部记号。

6.4 拖曳读数（流档）

物质要与周围时空保持稳定的纠缠（§2），它一动，那层关系就跟着换手 ⇒ 网在运流。设网的运流速度为 w：静态时 Δt = ρ²Δx（本征光速 1/ρ²，即 §6.1、§6.2 的 A·B = 1）；网以 w 运流时，光在坐标时间 Δt 内相对网的位移是 Δx − wΔt，故 Δt = ρ²(Δx − wΔt)，解出坐标光速 dx/dt = w ± 1/ρ²。与 ADM 形式逐项对齐（g_{0i} = (ρ²δ_ij + h_ij^TT)N^j，弱场即 ρ²N_i），得

N^i = −w^i                                                  (6.4)

即：拖曳就是网的运流速度，系数 1 来自读数规则本身——这一步不借任何坐标变换。（平动的那一份是坐标性的，转动的那一份不可被同时消去。）

运流由谁给：锁死要素把它的两通道读数份额一起带走。同一个 U 给出两个共轭读数——空间通道 B = e^{2U}（长度怎么读）、时间通道 A = e^{−2U}（tick 怎么读），A·B = 1；要素以速度 v 运动时，两通道的份额被一起带走：

w_i = (B_i − A_i) v_i = (e^{2U_i} − e^{−2U_i}) v_i          (6.5)

份额就是源自己的关系数（可数 ⇒ 互不重叠、直接相加，w = Σ_i w_i）。两通道方向必然相反，各向同性的部分不产生移位 ⇒ 移位只吃两者的失配 B − A：弱场下 B − A = ρ² − ρ^{−2} ≈ 4U，故

N^i = −(ρ² − ρ^{−2}) v^i   ⟶  弱场 N^i ≈ −4U v^i           (6.6)

那个 4 是读数本身的弱场系数（A·B = 1 的后果），不是额外输入。网不可压（a_P 常数 ⇒ 键长 a_P 恒定），稳恒运流满足 ∇·w = 0；源进这一档的只有角动量 J，故最低阶解只能是偶极：

N_i(x) = −2 (J × r̂)_i / r²,     量值 2J/r²                 (6.7)

数值核：均匀转动球的远场 N_y·r² 收敛到 −2J（网格 160 给 −3.350439，偏差 −0.018%；2J = 3.351032）、∇·N = 0，且 ½∇×N = −(J/r³)[3(Ĵ·r̂)r̂ − Ĵ] 与 Lense–Thirring 张量结构逐项一致（拖曳与转动同向）。读数换算中的格点离散项无偏、不累积，物理残留为 O((a_P/λ)²)（λ = 1 μm 时 3.2×10⁻⁵⁷）⇒ 任何波长都测不到。陀螺进动的观测对照见 §7.4。

6.5 统一度规

把各档合起来，几何层的度规读数就是式 (3.4) 的统一形式：

ds² = −(1/ρ²) dt² + (ρ²δ_ij + h_ij^TT)(dx^i + N^i dt)(dx^j + N^j dt)        (3.4)

表 1　统一度规的三个槽及其读出

槽	由谁读出	管什么
lapse 1/ρ 与空间标度 ρ²δ_ij	密度 ρ（τ_00）	物质、能量、各向同性压强
h_ij^TT	剪切（τ_ij 的无迹横部分）	引力波（2 偏振）
N^i	流（τ_0i）	转动、平动

各档由同一个 ρ 与同一条计数规则给出，不引入新的可调量；N^i = 0、h_ij^TT = 0 即退化为静态式 g₀₀ = −1/ρ²、g_ij = ρ²δ_ij：ρ 越大（锁死越多），钟越慢，而空间间隔的读数越大。一阶展开 g₀₀ = −1 + 2U、g_ij = (1 + 2U)δ 与 GR 的各向同性弱场形式同形 ⇒ γ = 1（二阶给 β = 1，见 §7.1）。两层的对接只有这一处：关系层的 ρ 与网的形变在此一次换成度规；§3.2 的面积律与 §5.1 的弹性常数都在各自层内使用，不跨层比系数。

7 与观测的对照（自洽性验证）

7.1 PPN 与太阳系

一阶（γ）：度规 (6.1)、(6.2) 的一次项就是各向同性弱场形式——ρ = 1 + U（U 为小量）⇒ g₀₀ = −1 + 2U、g_ij = (1 + 2U)δ ⇒ γ = 1 自动成立。力取相对梯度 ∇(ln ρ)、时钟取 1/ρ，两者同出这一个 ρ。

二阶（β）：把源规则 n_lock = a_P·ρ_m 代入式 (3.21)，即得 ∇²(ln ρ) = −4π·ρ_m（u = ln ρ），点源解为 ρ = e^{μ/r}（式 (3.25)）。代入式 (6.1)、(6.2)，按 U = μ/r 展开到 O(U²)：

g₀₀ = −ρ^{−2} = −e^{−2U} = −(1 − 2U + 2U² − …)              (7.1)

PPN 的标准形式是 g₀₀ = −(1 − 2U + 2βU² + …)，逐项对比即得

β = 1,      γ = 1                                           (7.2)

β = 1 不是调出来的：它来自“计数的对数是 ln ρ”这唯一结构（§3.3）——非线性全部在读数 ρ = e^{μ/r} 里，方程本身严格线性（§5.2）。太阳系全部一阶、二阶检验通过（γ、β 分别由卡西尼号测距与月距激光测距等实验定出，见 [5,15,16]），差别退到 O(U³)。

数值（太阳系）：一阶、二阶 PPN 参数与 GR 逐项相同（γ = β = 1）；差异出现在 u³ 阶，同一面半径下 −g₀₀ 的相对差 ≈ u³/3——太阳 ≈ 3×10⁻¹⁸、白矮 ≈ 6×10⁻¹²、S2 近星点 ≈ 1.5×10⁻¹¹——远低于当前测量误差（红移的相对差见 §7.2 的表）。

7.2 强场：必须统一到同一面半径

比较强场必须统一到同一个【面半径 R = areal】——因为 g_ij = ρ²δ_ij 是空间共形平直，ρ = e^{μ/r} 中的 r 是 isotropic 半径；而面积由 g_θθ = ρ²r² 读出，4πR² = 4πρ²r²，故

R = ρ·r = r·e^{μ/r}                                         (7.3)

喉部极小由 dR/dr = e^{μ/r}(1 − μ/r) = 0 定出，在 r = μ 处：

R_min = eμ = 2.718282 μ                                     (7.4)

记 u = μ/r（R = r·e^u = μ·e^u/u）。在同一面半径下解

R/μ = e^u/u                                                 (7.5)

红移由时钟读数读出：1 + z = 1/√(−g₀₀) = ρ ⇒ z_织法 = e^u − 1；GR 在同一面半径下 z_GR = (1 − 2U)^{−1/2} − 1，其中 U = μ/R = u·e^{−u}：

z_织法 = e^u − 1,      z_GR = (1 − 2u·e^{−u})^{−1/2} − 1       (7.6)

比对的口径：两边取同一个质量 M（由 r → ∞ 处 ρ = 1 + μ/r 的 1/r 系数定出，两种理论给出同一个 μ = GM/c²）与同一个面半径 R；给定 μ、R，式 (7.5) 唯一解出 u，再代入式 (7.6)。逐天体代入（太阳取 IAU 标称值[17]；白矮取天狼星 B 的动力学质量与半径 0.008098 R☉[18]；S2 取 GRAVITY 测得的 120 AU 近星点距离[19]与 Sgr A* 质量[20]；中子星两例为 NICER 的 M–R 结果[21,22]）：

表 2　同一面半径下各天体表面的引力红移对比

对象	M（M☉）	R	u	Δz/z
太阳	1	6.957×10⁸ m	2.1×10⁻⁶	< 10⁻¹¹
天狼星 B（白矮）	1.018	5.634×10⁶ m	2.67×10⁻⁴	+1.2×10⁻⁸
S2 近星点	4.30×10⁶	1.795×10¹³ m	3.54×10⁻⁴	+2.1×10⁻⁸
J0030+0451	1.34	12.710 km	0.1879	+0.78%
J0740+6620	2.072	12.390 km	0.3507	+3.54%

差别从 u³ 阶开始：同一面半径下 −g₀₀ 之差 = −(1/3)u³ + O(u⁴)，故 −g₀₀ 的相对差 ≈ u³/3、红移的相对差 ≈ u²/6（u ≪ 1 的渐近）⇒ 强场“温和”：中子星表面也只差百分之几，与面半径的观测误差同量级（8%–10%） ⇒ 现有中子星观测尚不能分辨。

一切强场比较必须以同一面半径为准：把面半径当作坐标半径比，会得到假偏差。

---

7.3 引力波

张量扇区给出 (∂_t² − ∇²)h_ij^TT = 0（§5.1）：两个偏振、无源传播、波速 = c[14]——与 GR 的引力波自由度一致，且偏振数 2 是网动力学的输出（纵向支为纯规范，自动离场）。

7.4 转动拖曳：陀螺进动（GP-B）

拖曳档给出的 g₀ᵢ（式 (6.6)、(6.7)）就是 Lense–Thirring 势。拖曳是物理量（不可被坐标变换消去），也是本理论与 GR 在转动这一档的可对照之处。陀螺在极轨上的平均拖曳进动（极轨对 Lense–Thirring 张量取平均，⟨cos²θ⟩ = ½，数值核 0.500003）为

Ω_LT = ½ · GJ/(c²r³)                                         (7.7)

取 GP-B 轨道 r = 7018 km、地球 J_⊕ = 5.861×10³³ kg·m²/s：GJ/(c²r³) = 81.96 mas/yr ⇒ 轨道平均进动幅值 Ω_LT ≈ 41 mas/yr。GP-B 引的 GR 预测 39.2 mas/yr[23] 是陀螺自转轴的漂移率，即上述幅值经导星几何投影后的值：自转轴指向导星 IM Pegasi（赤纬 16.841°），投影因子 cos δ = 0.957，41 × 0.957 ≈ 39.2 mas/yr ⇒ 同一口径下本理论与 GR 同值，两者都是同一 Lense–Thirring 场的输出；与 GP-B 实测 37.2 ± 7.2 mas/yr 相差 0.28σ，落在实测之内。

7.5 a_P / l_P 的地位

理论内唯一确定的长度比是

a_P / l_P = 2√(ln8)      （≈ 2.884）                    (7.8)

它是 [1] 推论 3 的定义关系（a_P 与 l_P 互为定义的读数），本身不携带独立的可检验预言；织法成立与否，靠本文方程与强场读数（§7.2）这类物理量的输出说话。

8 结论

1. 引力在关系层的场方程为 ∇²(ln ρ) = −(4π/a_P)·n_lock：引力不是一种力，而是 a_P 场密度的相对梯度；其中 G 与 c 的涌现公式由 [1] 推论 4 给出，本文直接引用。
2. 源规则：ρ_m = n_lock/a_P，m = m₁ × 质量份数（m₁ = ħ/(a_P c)）；源幅度比即观测质量比（m_p/m_e）。
3. 张量扇区 (∂_t² − ∇²)h_ij^TT = 0（2 个偏振、波速 = c，各向同性由单条计数规则保证）；标量扇区严格线性，ρ = e^{μ/r} 全阶精确。
4. 度规读数统一为式 (3.4) 的一个式子：密度档给出 lapse 与空间标度，剪切档 h_ij^TT 给出引力波（§5.1、§6.3），拖曳档 N^i 由网的运流给出、远场为 Lense–Thirring（§6.4）。
5. 与观测的对照是自洽性验证：PPN 中 γ = β = 1 自动通过；陀螺拖曳进动与 GP-B 相符（§7.4）；强场偏差落在中子星面半径观测误差之内；引力波自由度与 GR 一致。

---

参考文献

[1] Xinkuang Liu. The Weaving Formula: Emergence Mechanism of Spacetime Geometry[Z]. Zenodo. DOI: 10.5281/zenodo.22899821.
[2] Xinkuang Liu. The Weaving Origin of the Gauge Forces[Z]. Zenodo. DOI: 10.5281/zenodo.22710656.
[3] Xinkuang Liu. The Weaving Structure of the Matter Spectrum[Z]. Zenodo. DOI: 10.5281/zenodo.22928553.
[4] Albert Einstein. Die Grundlage der allgemeinen Relativitätstheorie[J]. Annalen der Physik, 1916, 354(7): 769–822. DOI: 10.1002/andp.19163540702.
[5] Clifford M. Will. Theory and Experiment in Gravitational Physics[M]. 2nd ed. Cambridge: Cambridge University Press, 2018. DOI: 10.1017/9781316338612.
[6] Steven J. Carlip. Quantum gravity: a progress report[J]. Reports on Progress in Physics, 2001, 64(8): 885–942. DOI: 10.1088/0034-4885/64/8/301.
[7] Luca Bombelli, Rabinder K. Koul, Joohan Lee, Rafael D. Sorkin. Quantum source of entropy for black holes[J]. Physical Review D, 1986, 34(2): 373–383. DOI: 10.1103/PhysRevD.34.373.
[8] Mark Srednicki. Entropy and area[J]. Physical Review Letters, 1993, 71(5): 666–669. DOI: 10.1103/PhysRevLett.71.666.
[9] Shinsei Ryu, Tadashi Takayanagi. Holographic derivation of entanglement entropy from AdS/CFT[J]. Physical Review Letters, 2006, 96(18): 181602. DOI: 10.1103/PhysRevLett.96.181602.
[10] Charles Kittel. Introduction to Solid State Physics[M]. 8th ed. Hoboken: Wiley, 2005.
[11] Xinkuang Liu. The Weaving Origin of the Proton–Electron Mass Ratio: 1836[Z]. Zenodo. DOI: 10.5281/zenodo.22951070.
[12] P. N. Keating. Effect of invariance requirements on the elastic strain energy of crystals with application to the diamond structure[J]. Physical Review, 1966, 145(2): 637–645. DOI: 10.1103/PhysRev.145.637.
[13] Clarence Zener. Elasticity and Anelasticity of Metals[M]. Chicago: University of Chicago Press, 1948.
[14] B. P. Abbott, et al. (LIGO Scientific Collaboration and Virgo Collaboration). Gravitational waves and gamma-rays from a binary neutron star merger: GW170817 and GRB 170817A[J]. The Astrophysical Journal Letters, 2017, 848(2): L13. DOI: 10.3847/2041-8213/aa920c.
[15] B. Bertotti, Luciano Iess, Paolo Tortora. A test of general relativity using radio links with the Cassini spacecraft[J]. Nature, 2003, 425(6956): 374–376. DOI: 10.1038/nature01997.
[16] James G. Williams, Slava G. Turyshev, Dale H. Boggs. Progress in lunar laser ranging tests of relativistic gravity[J]. Physical Review Letters, 2004, 93(26): 261101. DOI: 10.1103/PhysRevLett.93.261101.
[17] Andrej Prša, Petr Harmanec, Guillermo Torres, et al. Nominal values for selected solar and planetary quantities: IAU 2015 Resolution B3[J]. The Astronomical Journal, 2016, 152(2): 41. DOI: 10.3847/0004-6256/152/2/41.
[18] Howard E. Bond, Gail H. Schaefer, Ronald L. Gilliland, et al. The Sirius system and its astrophysical puzzles: Hubble Space Telescope and ground-based astrometry[J]. The Astrophysical Journal, 2017, 840(2): 70. DOI: 10.3847/1538-4357/aa6af8.
[19] GRAVITY Collaboration. Detection of the gravitational redshift in the orbit of the star S2 near the Galactic centre massive black hole[J]. Astronomy & Astrophysics, 2018, 615: L15. DOI: 10.1051/0004-6361/201833718.
[20] GRAVITY Collaboration. Mass distribution in the Galactic Center based on interferometric astrometry of multiple stellar orbits[J]. Astronomy & Astrophysics, 2022, 657: L12. DOI: 10.1051/0004-6361/202142465.
[21] Thomas E. Riley, Anna L. Watts, Slavko Bogdanov, et al. A NICER view of PSR J0030+0451: Millisecond pulsar parameter estimation[J]. The Astrophysical Journal Letters, 2019, 887(1): L21. DOI: 10.3847/2041-8213/ab481c.
[22] Thomas E. Riley, Anna L. Watts, Paul S. Ray, et al. A NICER view of the massive pulsar PSR J0740+6620 informed by radio timing and XMM-Newton spectroscopy[J]. The Astrophysical Journal Letters, 2021, 918(2): L27. DOI: 10.3847/2041-8213/ac0a81.
[23] C. W. F. Everitt, et al. (Gravity Probe B Collaboration). Gravity Probe B: final results of a space experiment to test general relativity[J]. Physical Review Letters, 2011, 106(22): 221101. DOI: 10.1103/PhysRevLett.106.221101.