[论文解读] A Wave Function approach to dissipative processes

Yvan Castin, Jean Dalibard+1 more|ArXiv.org|May 26, 2008

Quantum Mechanics and Applications参考文献 2被引用 4

一句话总结

本文提出一种蒙特卡洛波函数（MCWF）方法，用于模拟与马尔可夫 reservoir 耦合的耗散量子系统，通过随机量子跃迁和非厄米哈密顿量模拟开放量子动力学。该方法在形式上等价于主方程，但对大系统具有计算优势，并能提供对单个量子轨迹的直接物理洞察，将模拟结果与单个系统的实际实验测量联系起来。

ABSTRACT

Due to growing interest in quantum measurement, control and feedback, we reproduce a manuscript from 1992, presenting a simple physical and mathematical derivation of stochastic differential equations for wave functions of probed quantum systems. V. P. Belavkins seminal quantum filtering theory with similar equations, developed in the 1980es, was not known to the authors at the time of writing of the present manuscript.

研究动机与目标

• 开发一种适用于马尔可夫环境下开放量子系统的随机波函数方法，避免完整密度矩阵演化带来的计算开销。
• 建立 MCWF 方法与标准主方程方法在弛豫过程中的等价性。
• 提供一个物理直观的框架，其中量子跃迁用于建模耗散与涨落，将模拟与实际的单系统测量序列联系起来。
• 实现对具有大希尔伯特空间（N ≫ 1）的复杂量子光学与原子物理系统的高效模拟。
• 将模拟的量子轨迹与实验可观测量（如单系统测量中的量子噪声）联系起来。

提出的方法

• 系统在非厄米哈密顿量 H = H_S - iħ/2 ∑ₘ Cₘ†Cₘ 作用下演化，该哈密顿量在小时间步长 δt 内引起非幺正且规范长度减小的演化。
• 在每个时间步，以概率 δpₘ = δt ⟨ϕ(t)|Cₘ†Cₘ|ϕ(t)⟩ 发生量子跃迁，对应于算符 Cₘ 作用于波函数。
• 若发生跃迁，则波函数按 |ϕ(t+δt)⟩ = Cₘ|ϕ(t)⟩ / ||Cₘ|ϕ(t)⟩|| 更新，并随后归一化。
• 在连续时间极限情况下，该方法被重述为包含零均值、方差为 √Δt 的高斯噪声项 Δζₘ 的随机微分方程。
• 通过考虑大量无穷小跃迁分量 Dₘ,ε，推导出连续极限，得到与伊藤形式等价的随机方程。
• 通过证明波函数轨迹的系综平均精确重现主方程动力学，验证了该方法的有效性。

实验结果

研究问题

• RQ1 随机波函数方法能否重现由主方程控制的耗散量子系统动力学？
• RQ2 如何利用量子跃迁对开放量子系统中的耗散与涨落进行建模？
• RQ3 MCWF 方法与单个量子系统实际实验测量序列之间有何关联？
• RQ4 对于大希尔伯特空间，MCWF 方法在计算上与标准主方程方法相比有何优势？
• RQ5 在时间步长极小且跃迁通道数量众多的极限下，离散跃迁过程能否映射为连续随机微分方程？

主要发现

• MCWF 方法在形式上等价于主方程方法，波函数轨迹的系综平均可得到相同的约化密度矩阵演化。
• 该方法将计算成本从密度矩阵的 O(N²) 降低至波函数的 O(N)，使大希尔伯特空间维数 N ≫ 1 的系统模拟更加高效。
• 量子跃迁与物理测量过程直接关联：每条轨迹对应于在连续观测下单个量子系统可能的历史。
• MCWF 模拟中可观测量 A 所观测到的量子噪声，与真实实验中单系统测量所预期的量子噪声完全一致。
• 在时间步长趋于无穷小且跃迁通道数量众多的极限下，离散跃迁过程收敛为连续的伊藤型随机微分方程。
• 通过引入复跃迁算符 Dₘ,ε（其中 ε = ±1, ±i），实现连续极限，所得方程与文献中已知的随机波函数方程完全一致。