First passage time for an underdamped harmonic oscillator and application to the power of an information engine

技术摘要：欠阻尼谐振子的首次通过时间及其在信息机中的应用

问题陈述
首次通过时间（First Passage Time, FPT）分布的估计是统计物理学中的一个基本问题，其应用范围涵盖从化学反应速率到信息机等多个领域。虽然对于由朗之万方程描述的过阻尼系统，FPT 分布已得到充分研究，但欠阻尼机制下的情况仍然具有挑战性。在欠阻尼系统中，惯性引入了位移的二阶导数，使得随机运动方程呈现非马尔可夫特性。因此，计算 FPT 概率在历史上一直是一项“艰巨的任务”，闭式解通常仅在特定极限下（例如长时极限下的自由粒子）才可用。本文旨在解决布朗谐振子在有限品质因子 $Q$ 下跨越位置阈值 $x_B$ 的 FPT 分布 $P(t_{fp})$ 的计算问题。

方法论
作者提出了一个统一的理论框架，结合了三种不同的方法，以覆盖全范围的时间尺度和品质因子：

1. 概率密度函数（PDF）的分解： 总概率密度函数 $P(t_{fp})$ 被分解为三个贡献部分：

   • $P_I \delta(t_{fp})$：代表瞬时触发，即初始位置满足 $x(0) \geq x_B$。
   • $P_{II}(t_{fp})$：描述第一个振荡周期内的动力学过程（$0 < t_{fp} \leq 2\pi$）。
   • $P_{III}(t_{fp})$：捕捉长时行为（$t_{fp} > 2\pi$）。

2. 短时机制（$t_{fp} \leq 2\pi$）： 对于低阻尼情况（$Q \gg 1$），作者采用了哈密顿近似（Hamiltonian approximation）。他们假设总能量 $E$ 在一个周期内近似为常数。

   • 瞬时触发（$P_I$）： 使用初始位置的玻尔兹曼分布进行计算，得出 $P_I = \frac{1}{2}\text{erfc}(\sqrt{B})$，其中 $B = \frac{1}{2}x_B^2$ 是能垒。
   • 第一个周期（$P_{II}$）： 动力学被视为确定性的正弦运动 $x(t) = x_E \cos(\theta_0 + t)$。该概率通过对初始能量 $E$ 和相位 $\theta_0$ 的联合分布进行积分得出。这产生了一个平台期并随后快速衰减，其衰减由在时间 $t$ 处可达到的最大能量 $E^\dagger(t)$ 决定。

3. 长时机制（$t_{fp} > 2\pi$）： 对于超过一个振荡周期的时刻，系统的行为受热噪声引起的能量扩散支配。

   • 高 $Q$ 极限： 该问题被映射为跨越能量势垒的克拉默斯（Kramers）逃逸问题。通过基于随时间变化的能量阈值 $E^\dagger(t)$ 计算逃逸率 $\Gamma(t)$，最终收敛至常数速率 $\Gamma_B$。
   • 有限 $Q$： 作者通过计算与随机方程相关的克拉默斯微分算子的最慢特征值 $-\lambda_1(x_B)$ 来确定逃逸率 $\Gamma$，该算子带有位于 $x_B$ 处的吸收边界条件。
   • 由此产生的尾部 $P_{III}$ 遵循指数衰减 $e^{-\Gamma_B t_{fp}}$。

主要结果
论文给出了一个仅由能量势垒 $B$ 和品质因子 $Q$ 参数化的 FPT 分布闭式解析表达式（公式 2）。

• 实验验证： 理论预测通过使用室温下欠阻尼微悬臂梁（$Q \approx 7$）进行的实验得到了测试。该实验以 2 MHz 的采样率进行皮米级分辨率的位置测量。
• 一致性： 实验数据在三个不同机制下与理论模型表现出“极佳的一致性”：
  1. 瞬时触发的概率（$P_I$）。
  2. 分布在 $0 < t < \pi$ 期间的平台值（$P_{II}$）。
  3. 长时尾部的指数衰减速率（$P_{III}$）。
• 信息机应用： 推导出的 FPT 分布被用于估算信息机的功率。在这种设置中，当振子跨越阈值 $x_B$ 并导致势阱发生位移时，便提取出功。平均功率 $\langle P \rangle$ 被计算为平均提取功与平均 FPT 及停留时间之和的比值。理论功率曲线与实验数据相匹配，识别出了一个最优工作点，其中阈值 $x_B \approx 2$ 且位移 $L \approx 1$。

意义与主张
作者声称这项工作通过提供一个覆盖全范围品质因子的解，为欠阻尼系统的 FPT 问题提供了“新见解”，弥合了短时哈密顿动力学与长时扩散逃逸之间的鸿沟。

• 理论进展： 论文指出，所推导的闭式解（公式 2）代表了欠阻尼系统 FPT 统计学领域的一项进展。文中提到，此前已知的唯一关于欠阻尼系统的闭式解仅限于长时极限下的自由粒子（随机加速过程）。
• 与克拉默斯理论的区别： 作者强调，他们的描述与标准的克拉默斯速率理论不同，后者通常适用于跨越大型势垒的扩散过程。在他们的分布中，指数尾部仅在显著大于一个振荡周期的时刻才会出现，而短时行为则由惯性振荡主导。
• 实用价值： 该方法能够精确估算具有惯性的信息机的功率，而这在以往由于缺乏完整的分布知识而难以量化。作者表示，这为将 FPT 统计学应用于其他涉及欠阻尼动力学的领域（如随机热力学）铺平了道路。

论文最后指出，关于品质因子影响的完整计算细节和数值模拟已在随附论文 [22] 中提供。