Mean first passage time of chiral active Brownian particles

技术摘要：手性活性布朗粒子的平均首次通过时间

问题陈述
手性活性布朗粒子（CABPs）是具有内在旋转动力学的自推进粒子，这种动力学诱导其产生圆形轨迹。虽然标准活性布朗粒子（ABPs）表现出平移持久性，但 CABPs 打破了时间反演对称性和宇称对称性，从根本上改变了其力学响应和输运行为。理解这些粒子如何探索受限环境并定位目标，对于表征从运动细菌到合成微泳体等物理和生物系统中的输运与搜索效率至关重要。尽管兴趣日益增长，但关于 CABPs 平均首次通过时间（MFPT）的定量理论仍然稀缺。现有文献主要关注被动或非手性活性粒子，往往忽略了内在手性，或者仅针对特定的狭窄逃逸场景，而未将其推广到多样的几何构型。本研究旨在填补这一空白，阐明手性如何调节各种一维和二维受限域中的逃逸动力学。

方法论
作者采用基于向后福克 - 普朗克方程（backward Fokker–Planck equation）的理论框架来支配 MFPT，$T(x, q)$，定义为从位置 $x$ 出发、取向为 $q$ 的粒子逃逸出该区域所需的预期时间。控制方程为：
$$\mathcal{L}^\dagger[T] = -1$$
其中 $\mathcal{L}^\dagger$ 是伴随算子，包含了游动速度（$U_s$）、平移扩散系数（$D_x$）、角速度（$\Omega$）和旋转扩散系数（$D_R$）。

本研究考察了一维（1D）和二维（2D）两种典型几何构型：

1. 一维区间：
   • 对称区间 $[-L, L]$，两端均为吸收边界。
   • 非对称区间 $[-L, L]$，左侧为反射边界，右侧为吸收边界。
2. 二维圆盘：
   • 具有单个吸收弧的圆形区域。
   • 具有两个对称吸收弧的圆形区域。

分析采用了双重方法：

• 渐近分析： 在高手性机制下（$\chi \to \infty$，其中 $\chi = \Omega L / U_s$），作者采用匹配渐近展开法。这涉及将解分离为体区域和出口附近的边界层，推导闭式渐近解以捕捉 MFPT 随 chirality 的标度律。
• 数值模拟： 对于有限手性和任意参数，通过 FreeFem++ 使用有限元法（FEM）数值求解控制偏微分方程。结果针对特定情况与布朗动力学（BD）模拟进行了验证。

主要贡献与结果

• 高 chirality 机制（1D）：

  • 对于具有两个吸收边界的区间，作者推导出了复合渐近解。在高 chirality 极限下，主导阶 MFPT 收敛于被动布朗粒子（PBPs）的 MFPT。
  • 对于从中点出发的粒子，对被动结果的首次修正出现在 $O(1/\chi^{3/2})$ 阶，而对于偏离中心的粒子或非对齐取向，则出现在 $O(1/\chi)$ 阶。
  • 分析表明，手性引入了依赖于取向的修正，形成了边界层，粒子的旋转将其推离或拉向出口，导致对位置的非单调依赖。

• 有限 chirality 机制（1D）：

  • 两个吸收边界： 对于初始取向垂直于出口的粒子，MFPT 随 chirality 呈现非单调行为。存在一个最优的中间 chirality 值，使逃逸时间最小化。这一最优值源于游动时间尺度（$L/U_s$）与旋转时间尺度（$1/\Omega$）相当的平衡，使得 chirality 能够在粒子失去持久性之前将其导向出口。
  • 反射 - 吸收边界： 当粒子初始取向指向反射壁时，chirality 显著增强了逃逸。在低 chirality 下，粒子被持续的推进力困在壁面；中间 chirality 将其重新定向远离壁面，使 MFPT 低于被动布朗粒子。然而，在极高 chirality 下，MFPT 再次增加，因为定向运动受到抑制，最终收敛于被动布朗粒子极限。

• 二维圆盘几何：

  • 单个吸收弧： 在高佩克莱特数（$Pe$）下，MFPT 在中间 chirality 值（$\chi \approx 1$）处呈现尖锐的最小值。这一最优 chirality 在很大程度上独立于弧的大小。在低 $Pe$ 下，由于扩散占主导地位，非单调行为不太明显。
  • 两个对称吸收弧： 与一维对称情况类似，对于高活性粒子（$Pe=10$）观察到了非单调行为。然而，出现了一个关键区别：在单弧几何中，当$\chi \to 0$ 时，CABPs 的逃逸速度慢于被动布朗粒子（被反射边界困住）；而在双弧几何中，当 $\chi \to 0$ 时，CABPs 的逃逸速度快于被动布朗粒子，因为双出口允许粒子沿边界高效滑动至出口。

意义与主张
本文确立了 chirality 作为受限系统中输运和搜索的稳健控制参数。核心发现是存在一个最优 chirality，其处于中间值，能使 MFPT 最小化。这一现象在高活性粒子（大佩克莱特数）以及具有反射边界的几何构型中尤为显著，在这些构型中，标准活性粒子可能会被困住。

作者声称，他们的工作提供了对典型受限环境中 CABPs 逃逸动力学的全面表征。他们证明，虽然中间 chirality 通过减轻困陷并将粒子导向目标来增强逃逸效率，但过高的 chirality 是有害的，因为它抑制了持久游动，导致系统退化为类被动扩散。

该研究将 chirality 不仅仅视为一种生物特征，而是将其作为工程合成活性系统的可调参数。结果表明，通过控制旋转速率（例如通过磁场或几何受限），可以优化靶向药物递送、微流体输运和环境传感的搜索策略。这项工作为理解 chirality 与边界几何之间的相互作用如何产生非平凡的输运行为提供了定量基础，将先前的狭窄逃逸理论扩展到了更广泛的参数范围和几何构型。