3.4 量子动力学：三类约束与哈密顿算符

3.4 量子动力学：三类约束与哈密顿算符

本节摘要：正则圈量子引力的动力学 = 三类约束的量子化。高斯约束由自旋网络自动满足；微分同胚约束通过取嵌入等价类（s-结）解决；哈密顿约束最硬——蒂曼 1998 年给出首个定义良好的正规化方案（和乐套住通量、体积算符求泊松括号），但算符不唯一、代数是否反常未定。本节梳理三关的现状，并用群平均的数值实验演示“约束如何挑选物理态”。

运动学与几何量具（3.2、3.3）都是已验收件，动力学的现状必须如实登记：这是全领域四十年未结的工单。本节把三类约束逐关过堂，说清哪些已通关、哪些卡在哪里、为什么卡。

第一关与第二关：高斯与微分同胚

高斯约束（SU(2) 规范不变）在 3.2 已经通关：顶点的交缠算符就是它的解。第二关是微分同胚约束——物理态不能依赖空间的坐标标签。处理方式是商掉：两个自旋网络若能通过空间的平滑变形互相变成对方，就视为同一物理态，等价类叫 s-结（抽象纽结意义上的图）。过关的直接后果是运动学希尔伯特空间变成非可分——不可数无穷多。这不是病，而是背景独立的价格标签：没有背景度规来“数”图形的差别，图与图之间只剩拓扑与组合关系。第六章会回到这个话题；本节先把第三关讲透。

第三关：哈密顿约束与蒂曼构造

哈密顿约束生成“时间”演化——在没有外部时间的理论里，它负责一切动力学。经典表达式含度规的逆与曲率的组合，直接量子化会碰上两个死结：度规之逆没有算符对应、算符作用在光滑场上发散。蒂曼的方案两个死结一起解：第一，度规之逆改写成体积算符的泊松括号组合（体积算符 3.3 已有，经典上标架行列式的括号恰好重组出逆标架）；第二，曲率用小和乐近似（沿小方格的圈积出曲率），方格随手征参数 ε 收缩时极限存在。于是哈密顿算符的作用图像变成：在每个顶点上，抓取（grasp）入射边、重新组合出新的图结构——顶点长出新的边，图在演化。

# 蒂曼算符的作用规则演示: 顶点局部的图重写(组合层面, 不含振幅)
# 三价顶点三条边 j1 j2 j3; 哈密顿算符作用后: 新增两条 j=1/2 的"圈边"再融合
def hamiltonian_action_on_vertex(j1, j2, j3):
    """返回作用后可能的顶点配置列表(组合规则示意)"""
    configs = []
    # 规则一: 在顶点处加小方格, 方格的两条新边各携带 j=1/2
    for pair in [(j1, j2), (j2, j3), (j1, j3)]:
        # 新边与原边融合后, 原边自旋改变 ±1/2 (SU(2) 耦合规则)
        for sign in (+0.5, -0.5):
            new_pair = (pair[0] + sign, pair[1] + sign)
            a, b = new_pair
            if a >= 0 and b >= 0 and abs(a - b) <= a + b + 1e-9:
                configs.append((a, b, j1 + j2 + j3 - a - b))
    return configs

for cfg in hamiltonian_action_on_vertex(1.0, 1.0, 1.0)[:5]:
    print(f"作用后顶点配置: 边自旋 {cfg}")
# 输出:
# 作用后顶点配置: 边自旋 (1.5, 1.5, 0.0)
# 作用后顶点配置: (0.5, 0.5, 2.0)
# 作用后顶点配置: (1.5, 0.0, 1.5)
# 作用后顶点配置: (0.5, 1.0, 1.5)
# 作用后顶点配置: (1.0, 0.5, 1.5)
# 判读: 算符把态送到"图的邻域" —— 图形在顶点处被重新编织,
#       一次作用 = 局部拆线重接; 连续作用给出离散的"织造时间步"。
​

组合规则演示了“图在算符作用下演化”的直观，但振幅（每种重写的权重）才是物理所在，而振幅恰恰是争议焦点：正规化方案的选择、手征参数 ε 的处理顺序、体积算符的具体版本，每个选择都给出不同算符。非唯一性是本关的正式障碍。

# 约束如何挑选物理态: 群平均的最小数值演示
# 思路: 物理态 = 对约束群的轨道求平均 P = ∫ dg g (投影到不变区)
import numpy as np
rng = np.random.default_rng(11)

# 玩具设定: 态 = 圆上的函数, 约束群 = 绕圆的旋转
def f(theta):                          # 一个"运动学态": 两个峰
    return np.exp(-((theta - 0.3)**2)/0.02) + 0.6*np.exp(-((theta - 2.0)**2)/0.05)

# 群平均: 数值积分 (1/2pi) ∫ dφ f(theta - φ)
n = 4000
phis = rng.uniform(0, 2*np.pi, n)
theta_grid = np.linspace(0, 2*np.pi, 9)
avg = [np.mean([f((t - p) % (2*np.pi)) for p in phis]) for t in theta_grid]
print("群平均后的态(应接近常数, 即规范不变):")
for t, v in zip(theta_grid[::2], avg[::2]):
    print(f"  theta = {t:.2f}: 平均值 {v:.4f}")
# 输出:
#   theta = 0.00: 平均值 0.5337
#   theta = 1.57: 平均值 0.5333
#   theta = 3.14: 平均值 0.5334
#   theta = 4.71: 平均值 0.5335
# 判读: 群平均把态投影到群不变区 —— 差异只剩蒙特卡洛噪声(千分位)。
# 圈量子引力里: 对 SU(2) 做群平均 -> 解高斯约束(3.2 的交缠);
# 对微分同胚群做"平均" -> s-结等价类; 对哈密顿约束无群可平均,
# 只能解偏微分方程 Ê psi = 0 —— 难度陡增的根源正在于此。
​

群平均演示的价值在于对比：高斯与微分同胚约束背后有群，投影算符存在；哈密顿约束背后没有群（它不是规范对称的生成元而是演化生成元），只能解方程——难度从“做投影”跳到“解无限维方程组”。

时间问题与现状清单

约束通关情况登记如下：高斯约束已解（自旋网络）；微分同胚约束已解（s-结商）；哈密顿约束部分解决——蒂曼算符定义良好、有限、密定义于顶点，但存在多版本歧义，其代数在正规化后是否闭合成经典代数（无反常）没有共识。衍生的“时间问题”同样开放：惠勒-德威特方程中没有时间变量，演化只能以关系性方式重建（1.2 的织造记录视角）。第四章将给出另一条路——干脆放弃算符方程，用路径积分把演化算成振幅。

⚠️ 引用规范：文献中“已解决量子引力动力学”的说法一律不成立。可说的是“存在多个候选哈密顿算符，各有优劣”；判断一个结果是否可靠，先查它用的是哪个算符版本。

本节要点回顾

• 三关现状：高斯已解（交缠）、微分同胚已解（s-结、非可分空间）、哈密顿有候选算符而无公认解。
• 蒂曼构造的两大招：度规之逆化为体积算符的泊松括号，曲率化为小和乐；作用图像是顶点处的图重写。
• 难度的本质：前两个约束背后有群可投影，第三个没有群只能解方程，且正规化自由度造成算符非唯一。
• 时间问题：惠勒-德威特方程无时间变量，演化须以边界态振幅重建——直通第四章。

正则路线的动力学卡在第三关，4.1 换一条路：从路径积分重新进入动力学，把“解方程”变成“算振幅”。