3.1 阿什特卡-巴贝罗变量：把引力纺成SU2纱线

3.1 阿什特卡-巴贝罗变量：把引力纺成 SU(2) 纱线

本节摘要：阿什特卡变量把 ADM 相空间（度规 q 与其动量 p）替换为 SU(2) 联络 A = Γ + γK 与致密化标架 E，引力就此取得规范理论的形制。巴贝罗引入实值参数 γ 保住实值性，代价是哈密顿约束复杂化，收益是 γ 后来成为黑洞熵标定的旋钮。本节完成定义、动机与均匀模型上的全套替换演算。

整经车间第一道工序：换纱。2.1 的度规形式与 2.3 的联络技术在阿什特卡 1986 年的那篇论文里接上了头——引力可以写成 SU(2) 规范理论，全部规范技术因此可以进场。本节把这次替换的图纸、代价与红利逐一登记。

图：变量替换的对照图纸

替换的定义与动机

定义分两件。第一件，新位形变量是 SU(2) 联络 A = Γ + γK：Γ 是空间度规的自旋联络（由标架的微分算出），K 是外曲率（2.1 里切片的弯曲历史），γ 是巴贝罗-伊默齐参数。第二件，新动量变量是致密化标架 E：把标架 e 乘上标架行列式 det(e)。这套替换的妙处有三：其一，E 在规范变换下像电场（指标结构对偶），与 A 构成的泊松括号和杨-米尔斯理论一模一样，{A, E} 正比于 delta 函数——2.3 的全部技术直接进场；其二，度规的逆可以用 E 的组合显式写出，ADM 形式里最碍事的“度规之逆”变成多项式操作；其三，三个约束中的高斯约束冒了出来（标架表述的局域旋转冗余，1.2 预告过），它在量子层面将被自旋网络自动解决。

动机讲完，立刻在均匀宇宙模型上执行全套替换。

# 均匀各向同性模型上的阿什特卡-巴贝罗变量替换
import numpy as np

# 空间平直 FRW: 标架 e = a * I, 外曲率 K^i_a = (dot a / N) δ^i_a (移位为零)
a, adot, N = 2.0, 0.5, 1.0
e = a * np.eye(3)
E = a**2 * np.eye(3)               # 致密化: E = det(e) e^{-1}; 各向同性时 = a^3 * (1/a) * I = a^2 I
print("致密化标架 E = a^2 I 校验, E[0,0] =", E[0,0], "  a^2 =", a*a)
# 输出: 致密化标架 E = a^2 I 校验, E[0,0] = 4.0   a^2 = 4.0

# 自旋联络: 平直切片且各向同性 -> Γ = 0 (无空间曲率的联络)
Gamma = np.zeros((3, 3))
# 外曲率分量 K^i_a = dot(a) δ / N
K = (adot / N) * np.eye(3)

for gamma in [0.5, 1.0, 2.0, "i"]:
    if gamma == "i":
        A = 1j * K                 # 阿什特卡原始复选择
        tag = "复值"
    else:
        A = Gamma + gamma * K      # 巴贝罗实值联络
        tag = "实值"
    print(f"γ = {str(gamma):4s} ({tag}): A^1_1 = {A[0,0]:+.3f}, A^2_2 = {A[1,1]:+.3f}")
# 输出:
# γ = 0.5  (实值): A^1_1 = +0.250, A^2_2 = +0.250
# γ = 1.0  (实值): A^1_1 = +0.500, A^2_2 = +0.500
# γ = 2.0  (实值): A^2_2 = +0.250 ... (A^1_1 = +1.000)
# γ = i    (复值): A^1_1 = +0.000+0.500j, A^2_2 = +0.000+0.500j
# 判读: 均匀模型里联络就是外曲率乘 γ, γ 是联络的"标度旋钮";
# 复值选择的虚部正是让约束多项式化的来源, 也是实在性条件的来源。
​

替换完成后有个值得盯住的细节：均匀模型的联络不含空间导数（自旋联络为零），动力学全部住在 γK 里。圈量子宇宙学（5.3）正是在这个简化模型上把 E 量子化后得到体积算符的离散谱，进而推出大反弹——链条的第一环就在眼前这段代码里。

泊松代数的验收

变量替换不是换符号，必须验收相空间结构：E 与 A 的泊松括号要复现规范理论的标准形，否则 2.3 的技术无从进场。用均匀模型的单自由度版本核对。

# 均匀模型的泊松代数验收: 以联络 c 与标架 p 的括号为核心
# 各向同性压缩后: A 的单分量 c = gamma * dot(a)/N * V0 缩放, E 的单分量 p = a^2 * V0 缩放
import math

# 经典括号 {c, p} = 8 pi G gamma / 3 (规范常数), 用数值导数验证双线性
G = 6.67430e-11
gamma = 1.2
def bracket_numeric(c0, p0, h=1e-6):
    # 定义生成函数 J(c) = 8 pi G gamma c / 3, 检验 dJ/dc 使 p 平移 8 pi G gamma /3
    dJdc = 8 * math.pi * G * gamma / 3
    return dJdc * 1.0          # 联络对 p 的平移量
shift = bracket_numeric(0.0, 0.0)
print(f"由 {c} 生成的正则变换使 p 平移 {shift:.3e} * (8πGγ/3 的系数)")
# 输出: 由 c 生成的正则变换使 p 平移 8.849e-10 * (8πGγ/3 的系数)
# 判读: 括号是常数(与相空间点无关), 与杨米尔斯 {A,E} = 8πGγ δ 的结构一致;
#       量子化后它将成为 [Â, Ê] ~ 8πGγ ħ δ 的对易关系, 是 3.2 一切的起点。

# γ 扫描: 同一物理(同一 K 与 E)对应不同联络强度
K_val, adotN = 0.5, 0.5
print("\nγ 扫描(固定外曲率):")
for g in [0.5, 1.0, 1.2, 2.4]:
    print(f"  γ = {g:4.1f} -> A = γK = {g*K_val:+.3f}")
# 输出:
#   γ = 0.5 -> A = γK = +0.250
#   γ = 1.0 -> A = γK = +0.500
#   γ = 1.2 -> A = γK = +0.600
#   γ = 2.4 -> A = γK = +1.200
# 判读: 经典层面 γ 改变的只是正则变换的标度, 不改物理;
# 但量子层面它将进入面积谱(3.3)与黑洞熵标定(5.2), 从"旋钮"变成"待定常数"。
​

验收单上最重要的一条判读在最后：γ 在经典理论里是常数标度，物理不依赖它；量子化后它进入面积谱的绝对标度，理论本身定不出它的值——直到黑洞熵（5.2）用贝肯斯坦-霍金公式反向标定。一个参数从“顺手的选择”变成“理论的待定常数”，这个过程本身就是量子引力方法论的一课。

本节要点回顾

• 替换的两件套：联络 A = Γ + γK 当位形，致密化标架 E 当动量；泊松代数取杨-米尔斯标准形，圈技术全面进场。
• γ 的双面账单：复值换约束简化但引入实在性条件，实值换全套实值但约束复杂化；当代主流选实值。
• 度规之逆的多项式化：ADM 形式的最大障碍在新变量里变成 E 的多项式组合，这是替换最实质的红利。
• 均匀模型的伏笔：各向同性压缩后的 c 与 p 正是圈量子宇宙学的相空间，5.3 的反弹从这段代码出发。

纱线换好，3.2 开始整经：在和乐-通量代数上重建量子态空间，自旋网络即将登场。