MPS–HEOM 混合方法给出无序分子极化激元达到热力学极限所需的最小分子数 NT，并揭示声子时间尺度调控暗态激活
核心概要
作者提出矩阵乘积态与层级运动方程（MPS–HEOM）混合方法，对静态与动态无序下的分子极化激元动力学进行数值精确模拟，并引入收敛尺度 NT（光子动力学达到热力学极限所需的分子数），发现动态无序比静态无序要求更大的 NT，且 NT 随声子浴趋于马尔可夫而呈现先增后减的翻转行为，其微观根源是声子时间尺度调控亮态到暗态的能量转移与集体行为的抑制。
FIG. 1. (a) Schematic of a molecular ensemble in an optical
· 第 2 页深度剖析
作者构建了 MPS–HEOM 混合框架，将 HEOM 的辅助密度算符表示为矩阵乘积态，从而在统一张量网络内同时处理非马尔可夫振动弛豫与马尔可夫腔损耗和外部驱动，可模拟多达约 100 个二能级体系的中间耦合区，计算代价随体系大小线性增长。 此前非微扰 HEOM 因指数标度只能处理少数发射体，而能处理热力学极限的方法往往忽略耗散与泵浦；该工作把两者结合，使无序 TC 与 HTC 模型从少数发射体到宏观极限的动力学可在同一框架内数值精确求解。 文中给出 MPS 张量网络结构含 2+N(K+1) 个格点，采用 TDVP 时间演化并 GPU 加速，报告约 60 倍加速、N=10 时 TC 每步约 0.3 秒、HTC 每步约 2 秒，并给出 HEOM 截断参数（动态无序 Nb=4、L=5；静态无序 Nb=8）的收敛检验。
作者引入收敛尺度 NT，即光子动力学达到热力学极限所需的分子数，用相邻分子数之间动力学的时间归一化 RMSE（阈值 10^-4）判定，并据此定量回答集体极化激元体系达到热力学极限所需的最小体系大小。 此前实验用宏观样品（一般 >10^5 分子）而多数理论方法适用于 <20 分子的体系，二者之间的鸿沟缺乏定量标尺；NT 提供了连接少数发射体与宏观区间的可操作判据。 以实验可测的平均光子数 ⟨a†a⟩ 为观测量，参数对齐 BODIPY-Br 类有机极化激元（ωc=ω0=2.0，κ=Γ↓=0.02，Q=100），在无无序时得到 NT=3，并系统扫描无序强度 σ 与浴特征频率 γ。
作者发现光–物质耦合无序比频率无序要求更大的 NT，且动态无序总体上比静态无序更苛刻；NT 随浴趋于马尔可夫（γ 增大）先增后减，呈类似 Kramers 翻转的非单调行为，静态（非均匀）极限在该参数区间给出 NT 的下界。 此前对无序如何影响热力学极限收敛缺乏定量刻画；该工作把 NT 与无序类型、无序强度和声子时间尺度直接联系起来，并指出声子时间尺度同时控制集体行为的破坏与 NT 的增长。 基于图 2 中 NT 随 σ 与 γ 的扫描，以及图 3(d) 中 NT 随观测时间的演化：在 γ=0.08 时 NT 先升、在 t≈40 附近达峰后下降并饱和，而在静态无序与 γ=0.02 时 NT 几乎不随时间变化。
作者把上述趋势归因于无序对集体光–物质动力学的抑制：频率无序通过物质哈密顿量耦合亮态与暗态，而耦合无序直接破坏集体对称性使腔模与暗态耦合、暗态变为光学活性的灰态；暗态布居的翻转行为可由二阶微扰 Fermi 黄金规则解释，其转移速率主要由噪声谱在拉比劈裂处的谱重决定。 该工作给出无序激活非集体自由度的微观机制，并将其确立为支配无序光–物质体系热力学收敛的关键机制，同时指出声子时间尺度匹配对亮–暗转移的重要性。 补充材料给出亮–暗耦合矩阵元与 FGR 速率推导，得到 k_B→D=(N−1)/N·J(Ω_R)coth(βΩ_R/2)，在热力学极限趋于 J(Ω_R)coth(βΩ_R/2)；微扰理论复现了数值精确结果中的翻转行为，但在小 γ 区低估转移速率。
启示与展望
该结果面向研究分子极化激元集体行为的理论与从头模拟工作者：在 TC/HTC 模型、Debye-Drude 振动谱密度、单激发初始条件与所设腔损耗/泵浦速率下，NT 可作为选择模拟体系大小的定量依据，并提示可通过声子工程选择性增强把激发导入暗态流形的耗散通道。对实验工作者而言，它给出无序与声子时间尺度如何影响集体行为何时失效的判断框架。
本次载入的文本为不完整读取，正文图 2、图 3 的具体数值曲线与补充材料图 S1–S4 只能通过文字描述获知，无法核对图中细节；NT 的定量取值依赖收敛阈值 10^-4、观测窗口 tmax 与 HEOM 截断参数的选择，不同设定下的数值可能不同；微扰 FGR 在小 γ 区低估转移速率，说明该区间的定量解释仍以数值结果为准；此外，结论建立在特定谱密度形式与参数对齐单一实验体系之上，向其他分子体系与多模腔情形的适用性仍是开放问题。