Non-stationary Statistics and Energetics of Brownian Motion under Stochastic Harmonic Confinement

技术摘要：随机谐振约束下布朗运动的非平稳统计与能量学

问题陈述
本研究调查了在陷阱刚度随时间随机波动的谐振陷阱中，布朗粒子的统计动力学和热力学性质。尽管先前的研究已经探索了具有波动参数（退火无序）的系统（例如由平稳 Ornstein–Uhlenbeck 过程建模的系统），但这些模型通常允许负刚度值，这对应于物理上不切实际且瞬时不稳定的势能。作者通过将陷阱刚度 $k(t)$ 建模为布朗运动的平方（$k(t) = a^2(t)$）来解决这一局限性。这种公式确保了刚度始终保持严格非负，同时引入了一个均值刚度随时间线性增加的非平稳过程。主要目标是在随机热力学框架内，表征由此产生的非平稳位置涨落以及相关的能量交换（功与热）。

方法论
该系统由耦合的 Langevin 方程描述，其中粒子的位置 $x(t)$ 受热噪声和随时间变化的刚度 $k(t)$ 驱动，而 $k(t)$ 本身则作为布朗过程 $a(t)$ 的平方进行演化。作者分析了两种机制下的系统：

1. 零温情况 ($T=0$)： 粒子与热浴脱耦，动力学仅由随机刚度驱动。
2. 有限温情况 ($T \neq 0$)： 粒子与热库耦合。

为了求解动力学方程，作者采用了 Feynman-Kac 形式体系。他们推导了一个关于积分刚度泛函的后向 Fokker-Planck 方程。通过利用高斯假设，他们将该形式体系中产生的非线性 Riccati 方程转化为一个二阶线性微分方程，并对其进行了精确求解。

对于有限温度情况，作者利用了 次级化方法 (method of subordination)。他们将具有随机刚度的 Fokker-Planck 方程映射为一个带有随机时间变量的标准扩散方程。这使得通过对刚度轨迹系综进行平均，从而计算出位置矩的精确解析表达式成为可能。

热力学量（特别是对系统所做的平均功以及与热浴交换的平均热量）是利用随机热力学第一定律推导出来的。作者通过对内能的全导数并对热噪声和刚度实现进行平均，计算了功和热的速率。

核心贡献与结果

• 精确概率密度函数 (PDF)： 在缺乏热噪声的情况下 ($T=0$)，作者推导了粒子位置 $P(x,t)$ 的精确解析表达式。结果表示为一个包含 Heaviside $\theta$ 函数和指数项的无穷级数。数值模拟证实了与该解析解的高度一致性。
• 位置矩： 研究提供了在存在热噪声时，任意时刻前四阶位置矩 ($\langle x^n \rangle$) 的精确解析表达式。
  • 奇数阶矩： 第一阶和第三阶矩在长时间极限下呈指数级衰减至零。
  • 偶数阶矩： 第二阶和第四阶矩表现出代数衰减，在长时间极限下按 $t^{-1/2}$ 的比例缩放。
  • 物理阐释： 这种衰减行为归因于刚度的非平稳性质；由于平均刚度 $\langle k(t) \rangle = t$ 随时间线性增加，约束力随之增强，从而抑制了位置涨落。
• 热力学：
  • 功： 由波动陷阱对粒子所做的平均功呈非线性增长，在长时间极限下按 $t^{1/2}$ 比例缩放。这与由平稳无序驱动的模型形成对比，在后者中，功通常呈线性增长。
  • 热： 与热浴交换的平均热量也按 $t^{1/2}$ 比例缩放。
  • 内能： 内能的变化 $\Delta U$ 在长时间极限下趋于一个定值 $-\sigma_0 x_0^2 / 2$，这代表了由于将粒子从陷阱中心位移而注入的初始能量。
  • 平衡： 结果证实，在长时间极限下，热耗散速率与功注入速率平衡，符合热力学第一定律。

意义与主张
本文声称提供了一个受乘性退火无序驱动的、在物理上严格遵守刚度非负性的 Langevin 系统实例。其主要意义在于证明，将刚度建模为布朗运动的平方会导致一个非平稳过程，其渐近行为与之前研究的平稳模型（例如使用 Ornstein–Uhlenbeck 过程的模型）截然不同。

作者强调，刚度波动的非平稳性质从根本上改变了系统的能量学，导致功和热呈现亚线性 ($t^{1/2}$) 增长，而非在平稳无序模型中观察到的线性增长。这些理论预测得到了耦合 Langevin 方程数值模拟的严格支持。这项工作有助于深化对无序系统中非平衡热力学的理解，并暗示了其在模拟具有波动代谢策略的生态群落以及通过调制刚度的光学捕获胶体颗粒实现的实验应用方面的潜力。