[论文解读] Analytic influence functionals for numerical Feynman integrals in most open quantum systems

Nikesh S. Dattani, Felix A. Pollock+1 more|arXiv (Cornell University)|Mar 20, 2012

Spectroscopy and Quantum Chemical Studies参考文献 34被引用 7

一句话总结

本文提出了开放量子系统中离散化影响泛函的完全解析公式，用有限且可解析处理的阻尼响应函数 α(t) 的积分替代传统上数值不稳定的谱密度的非正常积分。对于一大类物理上相关的谱分布——包括洛伦兹-德鲁德和幂律形式——当 α(t) 为指数函数之和时，这些积分可得到闭合表达式，从而实现更快、更精确且更严谨的数值费曼路径积分计算。

ABSTRACT

Fully analytic formulas, which do not involve any numerical integration, are derived for the discretized influence functionals of a very extensive assortment of spectral distributions. For Feynman integrals derived using the Trotter splitting and Strang splitting, we present general formulas for the discretized influence functionals in terms of proper integrals of the bath response function. When an analytic expression exists for the bath response function, these integrals can almost always be evaluated analytically. In cases where these proper integrals cannot be integrated analytically, numerically computing them is much faster and less error-prone than calculating the discretized influence functionals in the traditional way, which involves numerically calculating integrals whose bounds are both infinite. As an example, we present the analytic discretized influence functional for a bath response function of the form α(t)=\sum_{j}^{K}p_{j}e^{Ω_{j}t}, which is a natural form for many spectral distribution functions (including the very popular Lorentz-Drude/Debye function), and for other spectral distribution functions it is a form that is easily obtainable by a least-squares fit . Evaluating our analytic formulas for this example case is much faster and easier to implement than numerically calculating the discretized influence funcitonal in the traditional way. In the appendix we provide analytic expressions for p_{j} and Ω_{j} for a variety of spectral distribution forms, and as a second example we provide the analytic bath response function and analytic influence functionals for spectral distributions of the form J(ω)\proptoω^{s}e^{-(ω/ω_{c})^{q}}. The value of the analytic expression for this bath response function extends beyond its use for calculating Feynman integrals.

研究动机与目标

• 消除传统费曼路径积分计算中开放量子系统所涉及的非正常积分带来的计算负担与数值不稳定性。
• 推导以阻尼响应函数 α(t) 而非谱密度 J(ω) 表示的离散化影响泛函（DIF）的一般解析表达式。
• 为一大类物理上相关的谱分布实现 DIF 的精确、快速且鲁棒的数值评估。
• 提供开源的 MATLAB 和 Mathematica 代码，以实现解析 DIF 的即用型部署。
• 通过确保环境耦合模型的一致性与解析精确性，支持近似量子动力学方法（如 HEOM 和 Nakajima-Zwanzig 方程）的基准测试。

提出的方法

• 推导以阻尼响应函数 α(t) 的正常积分表示的离散化影响泛函的一般表达式，替代 J(ω) 的无限区间积分。
• 聚焦于 α(t) = ∑ⱼ pⱼ e^{Ωⱼ t} 的情形，该形式自然源于洛伦兹-德鲁德和德拜函数等谱密度。
• 证明当 α(t) 为解析函数时，由此产生的 DIF 积分可因有限积分区间而实现精确或高效率的数值求解。
• 给出以谱分布参数（如 J(ω) = ω^s e^{-(ω/ω_c)^q}）和逆温度 β 表示的 pⱼ 和 Ωⱼ 的显式解析表达式。
• 使用帕德逼近（如 [N-1/N]）表示费米-狄拉克与玻色-爱因斯坦分布，以实现 α(t) 分解的高精度。
• 实现并发布开源的 MATLAB 和 Mathematica 代码，可直接从谱参数计算解析 DIF。

实验结果

研究问题

• RQ1 费曼路径积分中的离散化影响泛函能否以阻尼响应函数 α(t) 的有限、可解析求解积分形式表达？
• RQ2 在哪些物理上相关的谱分布中，α(t) 取 ∑ⱼ pⱼ e^{Ωⱼ t} 的形式，从而实现 DIF 的解析求解？
• RQ3 将 J(ω) 的无限区间积分替换为 α(t) 的有限积分，如何提升数值费曼积分计算的速度、精度与稳定性？
• RQ4 解析 DIF 在多大程度上可支持近似开放量子系统动力学方法（如 HEOM 和 Nakajima-Zwanzig 方程）的基准测试？
• RQ5 解析 DIF 是否能通过标准数值软件高效计算并实际部署？

主要发现

• 本文推导出以阻尼响应函数 α(t) 的正常积分表示的离散化影响泛函的一般公式，避免了对 J(ω) 的数值不稳定的非正常积分的依赖。
• 对于洛伦兹-德鲁德/德拜及幂律形式等谱分布（J(ω) ∝ ω^s e^{-(ω/ω_c)^q}），对应的 α(t) 取 ∑ⱼ pⱼ e^{Ωⱼ t} 的形式，从而实现完全解析的 DIF。
• 解析 DIF 在计算上更优：其计算速度更快，且比依赖无限区间数值积分的传统方法更不易产生数值误差。
• 为三类主要谱分布（包括通过帕德逼近表示的费米-狄拉克与玻色-爱因斯坦形式）推导出 pⱼ 和 Ωⱼ 的显式解析表达式。
• 该方法可对具有复杂谱密度（包括幂律形式 J(ω) ∝ ω^s e^{-(ω/ω_c)^q}）的系统实现精确的 DIF 计算，推导出以不完全伽马函数和三角函数项表示的 α(t) 的解析表达式。
• 作者提供了开源的 MATLAB 和 Mathematica 代码，可直接从谱参数计算解析 DIF，确保广泛的可复现性与可用性。