近共振光传播的量子相空间方法基准测试

Benchmarking quantum phase-space methods for near-resonant light propagation

arXiv: 2602.17660v1

论文信息

标题: Benchmarking quantum phase-space methods for near-resonant light propagation

作者: Mojdeh S. Najafabadi, Joel F. Corney, Luis Sanchez Soto, et al.

发布日期: 2026-02-19

arXiv ID: 2602.17660v1

PDF 链接: 下载 PDF

3 分钟速览

  • 研究问题：在光与原子介质近共振相互作用的背景下，如何系统性地对比 truncated Wigner approximation（TWA）和 positive P representation（PPR）两种相空间量子模拟方法的准确性、稳定性与适用范围。
  • 核心方法：利用 Jordan–Schwinger 映射将原子自由度转换为玻色子模式，在此基础上分别推导出基于 PPR 和 TWA 的随机微分方程组，并通过模拟自感应透明（SIT）光脉冲的压缩特性进行基准测试。
  • 关键结果：在无光学耗散时，TWA 与 PPR 的结果在中等相互作用强度下定量一致，但在更高原子密度或存在光库噪声时 TWA 出现明显系统性偏差（见图 1 和图 2）。
  • 主要局限：TWA 仅通过初始条件的 Wigner 分布引入量子涨落，忽略动力学过程中高阶导数项的贡献，无法先验地确定该方法在特定参数下的失效边界，必须依赖 PPR 作为对照。
  • 适合读者：从事量子光学、超冷原子物理或开放量子系统模拟的研究人员，尤其是需要选择或验证相空间数值方法的研究生与理论物理工作者。

论文背景和研究动机

量子力学中的相空间表述提供了一种绕过希尔伯特空间复杂性、以准概率分布描述量子态的框架。其中，Wigner 函数、Glauber–Sudarshan P 函数等虽然直观，却在处理非经典态或多体开放系统时面临严重障碍：有的分布高度奇异，有的产生负值，有的演化方程不满足正定扩散条件。

positive P representation（PPR）和 truncated Wigner approximation（TWA）是应对这些困难的两种主流方法。PPR 通过扩展相空间维数保证扩散矩阵正定性，但伴随数值不稳定性和更高的计算开销；TWA 直接截断高阶导数项，把量子效应完全归入初始噪声，计算效率高，但在强相互作用或长时演化下可能出现定量误差。二者在强驱动与强耗散情形下各具优势，但此前缺少在开放光–原子系统中系统性的对比基准。

本研究正是为填补这一空白而展开。论文通过 Jordan–Schwinger 映射将二能级原子系综转换为玻色子模式，使其同时适配 PPR 与 TWA 的数学结构，避免了传统集体算符描述带来的高阶项困扰。在此基础上，作者推导了两套随机微分方程，并以自感应透明（SIT）光脉冲在共振介质中的传播为具体场景，比较不同相互作用强度与不同耗散条件下三种模拟方案（PPR、PPR_DR、TWA）对压缩比预测的一致性。

核心方法和技术细节

Jordan–Schwinger 映射与模型构建

论文将每个二能级原子的泡利算符映射为两个独立的玻色子模式 b^1,b^2\hat{b}_1, \hat{b}_2：基态对应 b^1\hat{b}_1，激发态对应 b^2\hat{b}_2。激发过程相当于湮灭一个基态玻色子并产生一个激发态玻色子，具体映射关系为

σ^+=b^2†b^1,σ^−=b^1†b^2,σ^z=12(b^1†b^1−b^2†b^2)。\hat{\sigma}^+ = \hat{b}_2^\dagger \hat{b}_1,\quad \hat{\sigma}^- = \hat{b}_1^\dagger \hat{b}_2,\quad \hat{\sigma}^z = \frac{1}{2}(\hat{b}_1^\dagger\hat{b}_1 - \hat{b}_2^\dagger\hat{b}_2)。

这种映射消除了原本集体算符非正则对易关系带来的高阶相位空间项，使标准的相干态 PPR 可直接应用，也让 TWA 能够自然地推广到原子系统中。

原子与单模导引光场的相干相互作用由旋转波近似下的哈密顿量描述，包含原子自身能量 HAH_A、光场传播色散 HFH_F 和原子–场耦合 HAFH_{AF} 三部分。光学损耗则通过 Lindblad 形式的光库超算符引入，包含热光子数 nˉ\bar{n} 和损耗率 γ\gamma。

两种相空间动力学推导

利用论文表 1 中的算符对应规则，将主方程映射为 Fokker–Planck 方程，再转换为伊藤随机微分方程。

  • PPR 框架：在扩展相空间中，(αj,αj+)(\alpha_j, \alpha_j^+)、(β1n,β1n+)(\beta_{1n}, \beta_{1n}^+)、(β2n,β2n+)(\beta_{2n}, \beta_{2n}^+) 各对被视为独立变量而非复共轭。从式 (15)–(16) 可见，原子–场非线性耦合不仅贡献漂移项，还产生乘法噪声项 FαjF^{\alpha_j}、Fβ1nF^{\beta_{1n}} 等，其强度与耦合系数 gg 及动力学变量自身相关（如 −igβ2n/2\sqrt{-ig\beta_{2n}/2} 形式的系数）。光库噪声的扩散矩阵为 D(1)\mathbf{D}^{(1)}，可与相互作用噪声矩阵 D(0)\mathbf{D}^{(0)} 独立分解后拼接为总噪声矩阵，各组元噪声由独立的实 Wiener 增量驱动。

  • TWA 框架：在 Wigner 表示下，相互作用噪声的二阶项恰好抵消，只留下三阶及以上的导数项。TWA 将这些高阶项全部截断，方程 (18) 中原子部分的动力学完全失去乘法噪声，仅光学模式保留来自光库的热噪声（式 (19)）。量子涨落全部通过初态中按 Wigner 分布采样的高斯噪声项 δβ1n,δβ2n,δϕ(τ)\delta\beta_{1n}, \delta\beta_{2n}, \delta\phi(\tau) 引入，见式 (29)–(30)。

两种方法的根本差别由此清晰显现：PPR 在整个演化中持续注入和转化量子噪声，而 TWA 仅在初始时刻布撒一颗 “噪声种子”，后续演化由经典非线性方程决定。

创新点和贡献

本文的主要贡献并非提出全新方法，而是在于系统化基准和公式化框架的统一。

首先，Jordan–Schwinger 映射的引入是关键创新。此前该课题组的工作中，PPR 必须依赖集体原子描述以压制高阶项，这不仅限制了模型普适性，还无法兼容 TWA。新的玻色子表述一举消除了高阶相位空间项，使同一物理系统可以同时用 PPR 和 TWA 处理，为公平对比奠定了数学基础。

其次，首次在开放光–原子系统中完成 TWA 与 PPR 的直接对比。论文不仅考察了幺正演化（无光库）的情况，还专门分析了存在光学热库（nˉ=26\bar{n}=26，κ=1×10−6 m−1\kappa=1\times 10^{-6}\,\text{m}^{-1}）时的差异，清晰地揭示了耗散对 TWA 准确性的加速侵蚀：在较低原子密度（2.65×1021 m−32.65\times 10^{21}\,\text{m}^{-3}）时，图 2(a) 中 TWA 的偏离已经相当显著，而相同密度在无耗散时的图 1(a) 中三者仍高度吻合。

第三，对 “扩散规范” 的实践性讨论。论文指出 PPR 内部（如本文 PPR 与 Drummond 和 Raymer 的 PPR_DR）因扩散矩阵因子分解方式不同而产生不同的噪声分配方案。从原理上，这种规范自由度不影响物理结果，但在数值稳定性上可能产生巨大差异。作者明确指出这种自由度可用于调优，为后续算法优化提供了明确路径。

实验结果分析

幺正演化下的对比（图 1）

在无光学耗散时，三种方法在总原子密度 ρ=2.65×1021 m−3\rho=2.65\times 10^{21}\,\text{m}^{-3} 时给出了近乎重合的压缩比曲线。随着密度提升至 3.7×1022 m−33.7\times 10^{22}\,\text{m}^{-3}，TWA 的曲线逐渐偏离两个 PPR 结果的共同包络，且 PPR 自身的统计误差条显著增宽（论文提到 PPR 需约 2×1052\times 10^5 条轨道，而 TWA 仅需 5×1035\times 10^3 条）。

这一趋势符合理论预期：TWA 的截断误差随相互作用强度累积而放大。在弱耦合短距离下，初始噪声的非线性混合足以逼近真实量子动力学；一旦非线性效应深度重塑噪声统计，缺失的高阶项便无法忽视。

含光库耗散的对比（图 2）

耗散的加入使对比格局更加分明。即便在最低密度 ρ=2.65×1021 m−3\rho=2.65\times 10^{21}\,\text{m}^{-3} 时，TWA 的压缩比已明显低于两个 PPR 结果，且偏差随传播距离单调扩大。论文解释：当热光子数 nˉ\bar{n} 不可忽略时，光库噪声与相互作用噪声形成耦合，TWA 因缺少后者在传播中的动态自洽更新而系统性低估压缩。

值得注意的是，PPR_DR 与本文 PPR 在整个参数扫描范围内保持统计一致性，说明 Jordan–Schwinger 映射的 PPR 实现与传统集体算符 PPR 在不违背后者假设（每个时空–频率格点包含足够多原子）时是等价的。但本文 PPR 的采样误差在部分区域更大，暗示噪声结构仍有优化空间。

局限与待解决问题

尽管本文完成了富有价值的基准工作，但其局限同样不可忽视。

首先，TWA 的失效边界无法先验判断。作者在结论部分坦承：“当 PPR 与 TWA 的结果重叠时，我们可以信赖预测；一旦没有独立基准，单独使用 TWA 时难以确知是否仍在有效区间”。这意味着 TWA 在实践中必须伴随 PPR 作为参照，丧失了其本应具备的 “省钱替代品” 优势。

其次，PPR 本身的可扩展性受限。相空间维度加倍导致噪声诱导的数值散逸问题在大规模模拟中尤为棘手，论文未提出针对该问题的具体解决方案。尽管 Jordan–Schwinger 映射消除了高阶项，原子–场耦合噪声仍可能驱动轨迹进入不稳定区域。

第三，原子库耗散（碰撞阻尼）缺失。论文指出这是 “未来工作的重要方向”。实际原子介质除辐射阻尼外还存在退相干和碰撞效应，这些过程会引入额外的二阶及以上噪声，对 TWA 的截断假设构成更严峻挑战。在当前框架下，尚无法评估 TWA 在包含碰撞阻尼时的表现。

最后，模拟限定在一维单模近似内，忽略横向空间效应和三维传播中的衍射耦合。这一简化虽有助于方法对比，却限制了结果向真实空心光纤等系统的直接推广。如何在保持计算可行性的同时纳入更高维度动力学，是相空间方法走向实用化必须跨越的障碍。