惰性平面壁附近自扩散泳扁椭球Janus颗粒的极限环现象学：长宽比与活性/被动半球泳动迁移率比驱动的复杂Hopf分岔

《Soft Matter》：Limit cycle phenomenology of self-diffusiophoretic oblate Janus spheroids near an inert planar wallOpen Access

　　

研究人员研究惰性平面壁附近自扩散泳（self-diffusiophoresis）扁椭球Janus颗粒的通量预设模型（prescribed-flux model），考察长宽比（aspect ratio）以及两个半球表面上泳动迁移率（phoretic mobility）之比的影响。该工作将既有研究扩展至远离球形的形状（far-from-spherical shapes），并揭示颗粒形状对壁面效应（wall effects）的重要性。研究人员从动力系统（dynamical systems）视角展示这类系统的研究价值。在由长宽比与活性/被动半表面（active/passive half-surfaces）上泳动迁移率之比构成的参数平面中，相图（phase diagram）的显著特征包括超临界与亚临界Hopf分岔（supercritical and subcritical Hopf bifurcation），二者在低分辨率（low resolution）下看似平滑合并。对连接二者的新型分岔结构（novel bifurcation structure）的阐明表明，这类模型具有意外复杂的极限集行为（limit set behavior）。

研究背景与问题：微游泳者在受限环境中的运动具有基础与应用意义，实验已显示壁面等限制效应可显著影响其运动。自扩散泳Janus胶体颗粒具有活性表面，可产生或消耗溶质，溶质与颗粒表面相互作用驱动滑移速度并导致颗粒运动。Squirmer模型可预先规定表面滑移速度并考虑流体颗粒-壁相互作用，但不能恰当处理壁诱导的溶质浓度扰动及由此引起的滑移速度变化；更精细模型需处理溶质扩散，而在未采用完整化学动力学处理时，通量预设模型（prescribed-flux model，预先指定表面溶质通量和泳动迁移率及颗粒形状）是可用的模型类别。已有球形Janus颗粒研究考察活性区域大小及活性区与惰性区泳动迁移率比，但远离球形形状下的壁效应和极限集行为仍需系统映射，尤其需要从动力系统角度关注分岔。为此，研究人员将研究扩展至扁椭球（oblate spheroid）Janus颗粒，考察形状与迁移率比对极限环（limit cycle，相空间中孤立周期轨道）和分岔的影响。

研究内容、结论与意义：研究人员建立通量预设模型，将问题约化为二维自治动力系统，在长宽比a/b（0.3–0.8）与活性/被动半球泳动迁移率比μ_a/μ_p（1.0–10.0）参数平面绘制相图，分析固定点、稳定或不稳定极限环以及固定点不稳定性。研究发现相图中存在超临界和亚临界Hopf分岔（Hopf bifurcation，焦点稳定性改变并产生或吸收极限环），二者在低分辨率下似乎平滑合并；进一步揭示连接二者的新型分岔结构，表明该类模型具有意外复杂的极限集行为。该工作将自泳颗粒近壁研究推广至高度非球形颗粒，显示颗粒形状对壁效应具有定性影响，并可能首次详细映射微游泳者模型中的此类Hopf分岔；其意义在于从动力系统角度揭示丰富分岔结构，并为传统取向下的问题提供分析工具。

关键技术方法（不超过250字）：本研究无样本队列来源；方法为理论建模与数值分析。采用通量预设模型描述活性/被动半球上的溶质通量J_a与泳动迁移率μ_a、μ_p；在准稳态（quasi-steady-state）下将控制方程约化为Laplace方程与Stokes方程；利用Lorentz互易定理由滑移速度求颗粒平动与转动速度；用边界元法（boundary element method，BEM）求解Laplace与Stokes边值问题，粒子表面离散为512个六节点三角形单元、1026个网格点，壁面边界通过壁约束Green函数处理，线性系统用LAPACK DGESV直接求解；在参数网格上构建相平面速度场，固定点用SciPy optimize.root求解，轨迹用SciPy RK45积分，并用修正Poincaré映射定位极限环。

研究结果：
A 固定点：通过线性化相平面速度并分析矩阵特征值，将固定点分类为节点、鞍点或焦点；固定点处相平面速度散度等于矩阵迹，可定量表征不稳定性。相图显示在所研究参数范围内固定点唯一且总为焦点，颜色和等值线给出稳定性；接近白色区域时固定点移向壁而特征值实部不趋零，构成有壁条件下特有的分岔。
B 极限环必须环绕固定点：通过指标理论论证，极限环必须环绕至少一个固定点，固定点与其外围极限环的稳定性向外交替。
C Hopf分岔：通过固定点稳定性映射和零散度等值线定位Hopf分岔；当a/b?0.455时为超临界分岔并产生稳定极限环，0.455?a/b?0.68时为亚临界分岔并产生不稳定极限环，a/b?0.68时无法清晰分辨；简单Hopf分岔产生的极限环增长后撞壁消失。
D 连接Hopf分岔：通过检查a/b≈0.45附近发现同时存在两个极限环的微小区域，并用示意与数值结果说明：超临界Hopf产生的稳定极限环可沿S–W撞壁消失，或沿W–F与收缩的不稳定极限环相互湮灭；不稳定极限环从壁在W–U分离，沿F–Z1被固定点吸收而构成亚临界Hopf；W–F为折叠（鞍结）分岔；WIF区域内有两个极限环，W–F上单一周期轨道由稳定与不稳定极限环合并，既非线稳也非线不稳，吸引内部点并排斥外部点，其上方消失。
E 核、帽和杯：通过相平面全局结构分析，定义核（core）为固定点及其所有外围极限环内部或上的强不变区域，并讨论其外部。1 核：核的向前和向后极限集均包含于核内，核按其最外层元素稳定性表现为吸引或排斥。2 帽：通过相图数值示例得出，在足够大μ_a/μ_p下出现帽（cap），帽内点双向流向壁，动态隔离核与无穷。3 杯：通过相图示例提出在足够小μ_a/μ_p下出现杯（cup），杯内点双向流向无穷，动态隔离核与壁；杯不能被数值积分直接观测，识别不完全确定但可能性很高。4 无帽无杯：通过相图示例得出中间μ_a/μ_p范围存在既无帽也无杯的结构，核不与壁或无穷断开，固定点任意小邻域内有些点连无穷、有些点连壁，在初值敏感意义下呈混沌；该区域存在性似乎可靠。5 全局分岔：通过相图识别帽或杯的出现/消失为全局分岔，其非局部且比Hopf分岔更难分析。
F 涌现时间尺度与取向扩散：通过固定点散度估计颗粒向吸引固定点沉降的新时间尺度t_fp，在参数平面显著区域至少为(10^?2–10^?1)t_slip；无壁颗粒速度约为量级一，即v_p～b/t_slip；由Stokes转动阻力和Einstein关系估计水中300 K球体取向扩散系数D≈0.2 s^?1(R/μm)^?3，以b替换R后可小于10^?1，取向扩散不必成为严重问题。
G 通过修正Poincaré映射定位极限环：通过定义Θ=θ?θ₀及半线L（h=h₀, Θ≥0），构造修正Poincaré（第一返回）映射ΔΘ(Θ)，以零点定位极限环，以零点处斜率定量表征不稳定性（负斜率表示稳定），并可清晰显示两个极限环；a/b=0.45面板显示两个极限环，a/b=0.44顶部曲线接近折叠分岔W–F。1 极限环如何消失：通过数值积分和相图得出，亚临界Hopf产生的极限环撞壁消失；超临界情形极限环增长很高，a/b=0.40时网格延至h_max=20b仍继续增长，但在稍大μ_a/μ_p=1.75时观察到从不稳定焦点连向壁的轨道，且两极限环在WIF内有界，支持稳定极限环也以撞壁方式消失而非增长至无穷。

讨论与结论部分总结：A 与球体研究比较：Uspal et al.研究了非常相似的球体模型；参数对应关系为本文μ对应其?b，μ_a/μ_p对应β?1，取向角范围90°至?90°对应0–π且方向相反；其半活性情形为χ=0。其相图在小μ_a/μ_p区域与本文有相似之处但无固定点，减小β后的相平面速度与本文较大μ_a/μ_p有某些相似，但其固定点与本文固定点并不相同，因为本文固定点在a/b≈0.73时消失入壁。本文固定点对应文献中的sliding（颗粒以恒定高度平行壁稳定运动），±90°固定点对应hovering（对称性使平行壁速度为零）；二者在Uspal et al.参数区域均出现，但hovering在本文未出现。B 未来方向：尚未定位帽和杯出现/消失的分岔；尤其对杯，需要与闭式远场近似无缝匹配的方法，以便不显式积分确定流动性质或把轨道跟踪到所需h，并用于追踪增长至大h的极限环。C 分岔是核心：主要发现是简单自泳颗粒近壁表现出意外复杂的分岔结构，即使从纯动力系统角度看也很有意义；实验研究此类系统最有趣的方式可能是瞄准分岔，因为名义相同颗粒群体不可避免存在微小差异，这些差异会被非线性放大，导致定性不同且易于区分的长期行为。

研究结论部分翻译：研究人员的主要发现是，简单自泳颗粒在壁附近表现出意外复杂的分岔结构，即使从纯动力系统角度看也很有意义。实验研究此类系统最有趣的方式可能是瞄准分岔。名义相同的颗粒群体不可避免存在微小差异；可以说，这些微小差异会被非线性放大，导致定性不同且易于区分的长期行为。