Preserving Mass Shell Condition in the Stochastic Optimal Control Derivation of the Dirac Equation

技术摘要：在狄拉克方程的随机最优控制推导中保持质壳条件

问题陈述
狄拉克方程是相对论量子力学中自旋-1/2粒子的基石，它确保了每个旋量分量都满足克莱因-戈尔登方程，从而满足相对论质壳条件（$E^2 = p^2c^2 + m^2c^4$）。虽然随机力学已成功恢复了薛定谔方程和克莱因-戈尔登方程，但通过随机最优控制（SOC）推导狄拉克方程在历史上一直难以实现。此前的尝试（包括作者的前期工作 [10]）依赖于对相对论动能项（$-mc\sqrt{u^\mu u_\mu}$）和哈密顿-雅可比-贝尔曼（HJB）方程进行线性化。虽然这种线性化促进了类狄拉克结构的出现，但未能保持质壳条件；具体而言，由此产生的因子化方程遗漏了共轭算符中必要的 $+mc$ 项，导致无法为单个分量恢复克莱因-戈尔登方程。

方法论
本研究通过保留标准的非线性相对论平方根动能项，而非对其进行线性化，重新构建了狄拉克方程的 SOC 推导。其方法如下：

1. 协变拉格朗日量表述： 作者定义了一个处于外部电磁场中的单粒子相对论带电粒子（自旋-1/2）的洛伦兹协变拉格朗日量。该拉格朗日量包括：

   • 标准非线性动能项：$-mc\sqrt{w^\mu w_\mu}$。
   • 最小电磁耦合：$-\epsilon e A_\mu w^\mu$。
   • 协变自旋-场耦合（泡利项）：$-\frac{e\hbar}{2m}\sigma_{\mu\nu}F^{\mu\nu}$。

2. 标量化用于标量 HJB： 由于自旋-场项 $\sigma_{\mu\nu}F^{\mu\nu}$ 是矩阵值的，因此不能直接插入标量 HJB 方程。为了解决这个问题，作者引入了一个“标量化旋量” $\chi$，它是一个与时空无关的常数、自旋-场矩阵（$\sigma_{\mu\nu}F^{\mu\nu}$）的特征旋量（满足 $\sigma_{\mu\nu}F^{\mu\nu}\chi = f\chi$）。作者通过双线性投影 $P_\chi[M] = \frac{\chi^\dagger S M \chi}{\chi^\dagger S \chi}$ 定义了标量化拉格朗日量 $L_\chi$。这使得推导被限制在允许存在此类特征旋量的电磁背景（例如均匀磁场）中，并将系统固定在特定的自旋扇区。

3. 复数 SOC 与 HJB 推导： 该问题被置于复数 SOC 理论框架下，其中粒子在四维复时空中运动。值函数 $J$ 满足 HJB 方程。通过应用“弱条件” $w^\mu w_\mu = c^2$（仅在微分后应用，而非在线性化动能项时应用），作者推导出了依赖于值函数梯度的最优控制策略 $w_\mu$。

4. 向狄拉克方程的转换：

   • 利用 Hopf-Cole 映射（$J = i\epsilon\hbar \ln \phi$）对 HJB 方程进行变换，将非线性标量 HJB 方程转换为关于标量场 $\phi$ 的线性二阶方程。
   • 该线性方程被证明等价于一个“平方狄拉克”方程：$(i\hbar\gamma^\nu\partial_\nu - e\gamma^\nu A_\nu - mc)(i\hbar\gamma^\mu\partial_\mu - e\gamma^\mu A_\mu + mc)\phi = 0$。
   • 随后构造狄拉克旋量 $\psi$，即 $\psi = (i\hbar\gamma^\mu\partial_\mu - e\gamma^\mu A_\mu + mc)\phi \chi$。
   • 在平稳极限（$\partial_\tau \to 0$）下，这产生标准的狄拉克方程：$(i\hbar\gamma^\mu\partial_\mu - e\gamma^\mu A_\mu - mc)\psi = 0$。

5. 质壳分析： 作者分析了经典（$\hbar \to 0$）和量子极限下的质壳条件。他们证明了量子修正后的质壳关系包含了 Itô 二阶变差（量子势）和泡利项的贡献。至关重要的是，他们展示了当应用 Itô 协变修正时，物理双线性观测量 $\langle \pi_\mu \pi^\mu \rangle_{\text{phys}}$ 能精确恢复质壳条件，从而抵消量子势项并留下泡利不变量。

关键结果

• 质壳条件的恢复： 与之前的线性化推导不同，这种方法保持了相对论质壳条件。推导出的方程能够自然地分解为狄拉克算符及其共轭，确保解满足克莱因-戈尔登方程。
• 自旋扇区限制： 该推导适用于存在时空无关的 $\sigma_{\mu\nu}F^{\mu\nu}$ 特征旋量的背景。这允许将矩阵值自旋相互作用在固定的自旋扇区内视为标量。
• 数值验证（狄拉克-朗道问题）： 该理论使用均匀磁场中的电子（狄拉克-朗道问题）进行了数值测试。
  • 作用量极小化： 随机模拟确认，平均随机作用量在解析推导的最优漂移 $w^\star$ 处达到局部极小值。
  • 定量一致性： 在最优漂移下计算的随机作用量各分量（动能、电磁和自旋）与解析狄拉克-朗道值在 95% 置信区间内相匹配。
  • 局部最优性： 扰动分析表明，偏离最优漂移会导致期望作用量呈二次方增加，从而证实了局部极小性。
• 几何解释： 作者可视化了复动量空间，展示了量子质壳是如何通过微小的量子修正对经典弱质壳（$\pi_\mu \pi^\mu = m^2c^2$）进行的形变，且这些修正保持有界且呈振荡状态。

意义与主张
本文声称提供了一个“具有相对论一致性的 SOC 狄拉克方程推导”，解决了早期线性化方法在质壳条件方面的不一致问题。通过保留非线性动能项并利用标量化处理自旋-场耦合，作者在不牺牲基本能量-动量关系的前提下，架起了随机最优控制与相对论量子力学之间的桥梁。

这项工作强调，狄拉克方程可以从随机力学（特别是复数 SOC）的第一性原理中推导出来，其中 $\hbar \to 0$ 极限恢复了正确的原时相对论哈密顿-雅可比方程。作者将该理论定位为一个鲁棒的框架，其中谱和波函数的恢复被构造成一个高维优化问题，这可能适用于未来的变分量子算法（VQE, QAOA）。本文并非提出新的实验装置，而是为相对论自旋-1/2粒子的 SOC 框架提供了理论和数值验证。