Variational Multi-Gaussian Quantum Trajectories

技术摘要：变分多高斯量子轨迹

问题陈述
对开放量子多体系统（特别是相互作用玻色子系统）进行理论研究面临着巨大的挑战，其原因在于局部希尔伯特空间是无界的。随着占据数和系统规模的增加，精确模拟会迅速变得难以处理。虽然诸如截断维格纳近似（TWA）和高斯近似等半经典方法具有较高的效率，但在相互作用和非高斯相关性变得显著时，这些方法会失效。此外，尽管量子轨迹为描述在连续监测条件下进行的随机演化提供了一个强大的框架，但其精确模拟同样受限于复杂的多体问题。因此，需要一种既能高效模拟随机量子轨迹，又能捕捉超越半经典极限的非高斯物理特性的方法。

方法论：随机时间相关变分原理 (sTDVP)
作者提出了一种随机变分多高斯 (sVMG) 框架。该方法的核心是将时间相关变分原理 (TDVP) 进行推广，以适应随机量子态扩散过程。

1. 随机动力学： 系统由林德布拉德（Lindblad）主方程控制，该方程通过异测向解构（量子态扩散）被“解构”为量子轨迹。这导致了一个 Itô 随机薛定谔方程 (SSE)。为了应用变分微积分，作者将其转换为等效的 Stratonovich SSE，从而保留了将动力学投影到变分流形上所需的普通链式法则。
2. 变分拟设 (Ansatz)： 波函数 $|\psi\rangle$ 被近似为 $N_c$ 个高斯波包的相干叠加：
   $$\langle x | \psi[\theta] \rangle = \sum_{k=1}^{N_c} G(x; \theta_k)$$
   每个高斯分量都由一组实参数 $\theta$ 参数化（包括平均位置/动量、协方差矩阵和权重）。对于大型系统，协方差矩阵可以通过“低秩”分解进行压缩，以减少参数数量。
3. 运动方程： 通过将波函数的精确随机增量投影到变分流形的切空间上，作者推导出了 sTDVP 方程。这产生了一组关于参数 $\theta$ 的随机微分方程：
   $$d\theta = \text{Re}(S^{-1}) \text{Re} \left( V dt + \sum_m Y^{(m)} \circ dZ_m^* \right)$$
   其中 $S$ 是量子几何张量，$V$ 和 $Y^{(m)}$ 是漂移项和噪声力向量。至关重要的是，由于拟设由高斯函数组成，用于组装 $S, V$ 和 $Y^{(m)}$ 的积分（涉及多项式加权高斯函数）可以完全解析地完成，利用维克定理 (Wick's theorem) 即可实现。

核心贡献

• TDVP 的推广： 该工作将时间相关变分原理从单纯的幺正或林德布拉德演化扩展到了由连续测量驱动的随机动力学（量子轨迹）。
• 解析组装： 多高斯拟设允许解析地推导运动方程，从而实现高效的数值积分，无需对波函数本身进行蒙特卡洛采样。
• 受控非高斯性： 该方法通过改变高斯分量的数量 ($N_c$) 来系统地控制波函数的非高斯性。随着 $N_c$ 的增加，近似过程会趋向于精确动力学。

结果与基准测试
该方法在驱动-耗散 Bose-Hubbard 模型上进行了基准测试和应用：

1. Bose-Hubbard 二聚体 ($N=2$)：

   • 使用较少分量（$N_c=4$）的 sVMG 方法能够忠实地重现精确的单个轨迹和纠缠动力学，即使在半经典近似（TWA、高斯轨迹近似）失效的机制下也是如此。
   • 观测量的相对误差随 $N_c$ 的增加而系统性下降。
   • 该方法成功捕捉到了对称性破缺相变和在热力学极限（大 $N$ 标度）下的纠缠相变，并在非高斯性减弱处与平均场预测相匹配。

2. 扩展 Bose-Hubbard 点阵（1D 和 2D）：

   • 该框架被应用于 1D 链和 2D 点阵，以研究对称性破缺转变（棋盘格图案）。
   • 维度依赖性： 在 2D 点阵中，sVMG 结果显示出明显的对称性破缺相，并伴有棋盘格子间明显的区分。相比之下，1D 链表现出由畴壁分隔的局部类棋盘格区域，导致平均序参量被冲淡。这表明对称性破缺相在 2D 中存在，但在 1D 中由于增强的涨落而被抑制。
   • 该方法利用 $N_c=4$ 且秩 $K=2$ 的低秩拟设捕捉到了这些不同的行为，展示了超越暴力精确对角化能力的扩展性。

意义
论文声称 sVMG 框架为模拟相互作用玻色子系统的随机量子动力学提供了一个通用且高效的工具。其主要意义在于能够弥合计算成本较低的半经典方法与难以处理的精确模拟之间的差距。通过捕捉单个轨迹上的非高斯相关性，该方法使得研究在标准高斯或截断维格纳近似不足的机制下的现象（如纠缠动力学和对称性破缺相变）成为可能。其能够解析地组装运动方程，使其成为探索开放多体系统中非平衡相态的一种实用方法。