[论文解读] Quantum Friction: Cooling Quantum Systems with Unitary Time Evolution

Aurel Bulgac, Michael McNeil Forbes+2 more|arXiv (Cornell University)|May 29, 2013

Advanced Thermodynamics and Statistical Mechanics参考文献 4被引用 5

一句话总结

本文提出了一种局域量子摩擦——一种幺正、无记忆的方法，通过在哈密顿量中添加时间反演对称性破缺势能，实现对量子多体系统的冷却。通过实现 $O(N^2ackepsilon N)$ 的标度，该方法显著加速了时间依赖密度泛函理论（TDDFT）模拟中的基态制备，相较于传统 $O(N^3)$ 的对角化方法，在具有 $N \sim 10^5$ 到 $10^6$ 个单体态的大型费米子系统中表现更优。

ABSTRACT

We introduce a type of quantum dissipation -- local quantum friction -- by adding to the Hamiltonian a local potential that breaks time-reversal invariance so as to cool the system. Unlike the Kossakowski-Lindblad master equation, local quantum friction directly effects unitary evolution of the wavefunctions rather than the density matrix: it may thus be used to cool fermionic many-body systems with thousands of wavefunctions that must remain orthogonal. In addition to providing an efficient way to simulate quantum dissipation and non-equilibrium dynamics, local quantum friction coupled with adiabatic state preparation significantly speeds up many-body simulations, making the solution of the time-dependent Schrödinger equation significantly simpler than the solution of its stationary counterpart.

研究动机与目标

• 解决在强关联费米子系统的大规模时间依赖密度泛函理论（TDDFT）模拟中准备初始态的计算瓶颈问题。
• 克服在具有 $N \sim 10^5$ 到 $10^6$ 个单体态的系统中，通过迭代对角化求解自洽基态时 $O(N^3)$ 的高昂计算成本。
• 开发一种无需密度矩阵演化、无记忆的幺正替代方法，以替代Kossakowski-Lindblad主方程，用于建模量子耗散。
• 通过简化时间依赖薛定谔方程的求解，相比其定态形式，实现对核物理、超冷原子和电子结构中非平衡动力学的高效模拟。
• 展示通过绝热开关与局域量子摩擦相结合，可在三维大尺度系统（如 $16 \times 16 \times 64$ 晶格）中实现基态制备的可行性。

提出的方法

• 引入一个局域势能 $U_t$，通过显式破坏时间反演不变性，实现波函数的幺正演化以诱导冷却。
• 将在时间依赖薛定谔方程中应用局域量子摩擦势能，保持多体态的正交性，且无需演化密度矩阵。
• 将局域量子摩擦与绝热态制备相结合，实现从简单初始态到目标基态的渐进演化。
• 使用修改后的薛定谔方程进行实时演化，以准备TDDFT的初始构型，避免对大型 $N \times N$ 哈密顿矩阵进行自洽对角化。
• 在大规模计算平台（如NCCS的UW Hyak和Titan）上实现该方法，采用GPU加速求解耦合PDE系统。
• 在具有 $32,768 \times 32,768$ 个单体哈密顿矩阵和 $15,322$ 个耦合PDE的三维简并费米气体系统上验证该方法，仅用24块GPU在约1小时内实现收敛。

实验结果

研究问题

• RQ1 能否设计一种幺正、无记忆的方法，在不演化密度矩阵的前提下，模拟强关联费米子系统中的量子耗散？
• RQ2 在引入时间反演对称性破缺的局域势能后，是否能实现多体波函数的高效冷却，同时保持正交性？
• RQ3 与迭代对角化相比，局域量子摩擦是否能显著降低TDDFT模拟中初始态制备的计算成本？
• RQ4 在大规模系统中，使用局域量子摩擦的实时演化在求解时间依赖薛定谔方程方面，与定态方法相比能提升多少性能？
• RQ5 该方法是否可扩展至模拟具有 $N \sim 10^5$ 到 $10^6$ 个单体态的系统（如重核或超冷费米气体）中的非平衡动力学？

主要发现

• 局域量子摩擦通过时间反演对称性破缺势能，实现了费米子多体系统的幺正冷却，避免了密度矩阵演化和内存密集型主方程。
• 该方法通过用实时演化替代迭代对角化，将基态制备的计算成本从 $O(N^3)$ 降低至 $O(N^2\backepsilon N)$。
• 在使用24块GPU的UW Hyak集群上，仅用约1小时即完成了对具有 $16 \times 16 \times 64$ 晶格点和 $32,768$ 个单体态的三维简并费米气体系统的基态制备。
• 该方法在不同晶格尺寸（$24^3$、$32^3$、$40^3$）下均表现出红外收敛，证实了其数值稳定性和可扩展性。
• 绝热开关与局域量子摩擦的结合，相比纯绝热演化，显著缩短了切换时间，加速了收敛过程。
• 该方法可扩展至未来百亿亿次级系统（如NCCS的Titan，拥有18,680块GPU），支持在 $100^3$ 晶格上模拟数百万个耦合PDE。