[论文解读] Time evolution algorithms for Matrix Product States and DMRG

Juan José García‐Ripoll|arXiv (Cornell University)|Oct 8, 2006

Cold Atom Physics and Bose-Einstein Condensates参考文献 30被引用 60

一句话总结

本文提出并基准测试了一维量子系统中矩阵乘积态（MPS）的先进时间演化算法，比较了泰勒、帕德和阿诺尔迪方法在近似时间演化算符方面的表现。结果表明，阿诺尔迪方法在中等计算资源下实现了更优的精度，使得在非绝热扫掠条件下，对一维强关联系子玻色子中费什巴赫共振动力学的模拟成为可能，原子关联可被高保真度地映射到分子态。

ABSTRACT

In this work we develop several new simulation algorithms for 1D many-body quantum mechanical systems combining the Matrix Product State variational ansatz with Taylor, Padé and Arnoldi approximations to the evolution operator. By comparing with previous techniques based on MPS and DMRG we demonstrate that the Arnoldi method is the best one, reaching extremely good accuracy with moderate resources. Finally we apply this algorithm to studying how correlations are transferred from the atomic to the molecular cloud when crossing a Feschbach resonance with two-species hard-core bosons in a 1D optical lattice.

研究动机与目标

• 开发并比较一维量子系统中矩阵乘积态（MPS）的时间演化算法的效率。
• 评估不同近似方法——泰勒、帕德和阿诺尔迪——在时间演化算符上的精度与性能。
• 将性能最佳的算法应用于研究一维强关联系子玻色子在费什巴赫共振穿越过程中的关联转移。
• 证明即使在非绝热条件下，分子态关联也能忠实反映初始原子关联。

提出的方法

• 使用矩阵乘积态（MPS）变分试探法，以有界矩阵维度D表示多体量子态。
• 应用最优投影算符将态截断到MPS流形上，同时最小化误差。
• 对时间演化算符e^(-iHt)采用泰勒和帕德展开，实现短时间步长演化。
• 提出一种基于阿诺尔迪的新方法，通过构建MPS的Krylov子空间以提高精度与稳定性。
• 结合阿诺尔迪方法与自适应时间步长及截断误差控制，管理数值误差的累积。
• 使用自旋-1/2模型验证算法，并将最优方法应用于具有费什巴赫共振耦合的一维双组分强关联系子玻色子模型。

实验结果

研究问题

• RQ1 在固定计算成本下，泰勒、帕德或阿诺尔迪方法中，哪一种时间演化方法在MPS模拟中提供最高精度？
• RQ2 截断误差在时间演化过程中如何累积？能否在MPS模拟中有效控制？
• RQ3 在费什巴赫共振扫掠后，分子态关联是否能忠实反映初始原子二体关联？
• RQ4 在一维原子到分子转变中，实现高保真度关联映射所需的最短扫掠时间是多少？
• RQ5 有限跃迁项如何影响共振扫掠过程中相干性与关联转移效率？

主要发现

• 阿诺尔迪方法在精度与稳定性方面优于泰勒和帕德近似，尤其在大时间步长下表现更优。
• 截断误差与舍入误差是所有方法中的主要限制因素，阿诺尔迪方法展现出最稳健的误差控制能力。
• 当扫掠时间T ≈ 1.5/Ω时，分子态保有约90%的初始原子二体关联强度。
• 即使在J = 0.4Ω时，分子关联仍与初始原子关联高度一致，表明其在绝热区之外也具有鲁棒性。
• 由于分子态混杂，系统在共振附近显著抑制了原子跃迁，有效“钉住”了原子，从而实现关联测量。
• 即使采用非绝热扫掠，关联映射协议仍保持有效，表明其在当前冷原子实验中具有实际可行性。