非周期光子纳米结构中的近场驱动辐射热动力学

论文信息：

M. Prado, A. Manjavacas, F.A. Pinheiro, W.J.M. Kort-Kamp, Near-field-driven Radiative Thermal Dynamics in Aperiodic Photonic Nanostructures,arXiv:2512.18974(2025).

论文链接：

https://doi.org/10.48550/arXiv.2512.18974

研究背景

在纳米光子学中，调控结构关联性是操纵光-物质相互作用的核心。确定性非周期系统，如Vogel螺旋，通过单一参数（发散角）即可连续调节其结构有序度，从而支持丰富的电磁模式，如临界局域态与多重分形场分布，为新颖的光学调控提供了平台。然而，现有研究多集中于支持自由空间光传播的体系，对于由表面束缚模式（如等离激元与声子极化激元）主导的体系——特别是近场辐射热传递过程——结构关联性的影响尚未被揭示。该过程在纳米间距下由倏逝波主导，热流可远超黑体极限，是研究关联效应与能量输运动力学的理想载体。本研究旨在填补这一空白，探究确定性非周期序如何重塑近场热传递的模态谱与瞬态弛豫动力学。

研究内容

本研究以确定性非周期 Vogel螺旋阵列作为模型体系，深入探究了结构关联性对纳米颗粒系统中近场辐射热传递动力学的调控机制。研究核心在于阐明，如何通过单一几何参数——发散角α，实现对体系结构有序度的连续调节，并最终揭示这种可控的非周期序如何直接影响以表面极化激元为主导的热能弛豫路径与时间尺度。为此，研究构建了基于涨落电动力学框架和严格多体散射理论的理论模型，对由碳化硅纳米颗粒组成的多种阵列构型进行了系统的数值计算与对比分析。

图1:两个沃格尔螺旋中热化动力学的示意图，其中中心粒子最初的温度高于阵列中的所有其他粒子。

研究首先建立了适用于描述纳米颗粒近场辐射热交换的微观理论框架。在偶极近似下，每个半径远小于热波长（R< <λ< pan>T）的球形SiC颗粒被建模为一个具有频率依赖极化率α(ω)的局域电偶极子，其响应由材料的光学声子极化激元共振主导。颗粒体系与真空环境的热涨落作为驱动源，通过自洽的线性散射方程耦合，该方程完全计入了多体相互作用与多次散射效应，并以紧凑的矩阵形式表达了“穿饰”后的偶极矩与局域电场。应用涨落-耗散定理关联涨落源的统计性质后，可推导出颗粒i吸收的净辐射功率Pi的表达式，该表达式最终表述为与所有颗粒温度偏差相关的线性求和，其系数由耦合核函数fij(ω)在频谱上的积分决定，这一耦合核函数本质上编码了体系所有电磁散射信息。

在热动力学描述上，研究采用了本征模形式体系。假设电磁响应的时间尺度远快于温度演化（绝热近似），并在环境温度T0附近对玻色-爱因斯坦函数进行线性化，可将体系的温度弛豫过程转化为一个线性微分方程组。该系统的动力学完全由矩阵H=Γ-1F掌控，其中Γ是对角热容矩阵，F是对称的辐射耦合矩阵。矩阵H的本征谱提供了对热化过程的完整分解：每一个本征值λμ对应一个特征衰减率，其相应的本征向量ΔTμ则描述了与之对应的空间温度偏差模式。最慢的模式（μ=1，最小λ）对应整个系综作为一个整体与环境热浴的最终平衡，而更快的模式（μ＞1）则反映了由近场相互作用介导的、具有复杂空间结构的局域弛豫过程。通过求解该本征问题，研究可以将任意初始温度分布分解为这些本征模式的叠加，从而清晰解析热流在阵列中扩散的瞬态过程。

图2:半径为R=25nm、粒子间距为d1的N=1024的SiC纳米粒子方形阵列在不同初始条件下的热化动力学。彩色曲线显示了所选纳米粒子的温度随时间的演化，该粒子初始温度为ΔT=50K，而所有其他粒子最初与温度为T0=300K的环境处于热平衡状态。插图显示了初始热粒子的位置:蓝色(红色)表示离阵列几何中心最近(最远)的点，绿色点位于中间距离处。灰色曲线对应于所有纳米粒子最初都处于ΔT=50/1024K的情况，箭头表示蓝色初始条件下(d1=5R)的粒子间热化时间tTh。

为了量化这一瞬态过程，研究定义了“颗粒间热化时间”tTh，它表征从初始非平衡态出发，直至阵列内部温度达到均匀（此后仅剩整体衰减模式主导）所需的时间。这一时间尺度是评估不同几何结构调控热传递效率的关键指标。作为基准，研究首先分析了周期性正方阵列的热化行为。结果表明，当热量初始沉积于单个颗粒时，其热化动力学强烈依赖于该颗粒在阵列中的位置，中心位置由于被最近邻包围而热化最快。这一差异源于不同初始条件对同一组固定热本征模的投影不同。无论初始热点位于何处，最终所有颗粒都会达到一个均匀的中间温度，随后整个阵列以速率λ1整体弛豫。若初始热量均匀分布，则系统将直接进入这一最终衰减阶段，这验证了λ1模式的普适性。

随后，研究聚焦于四种具有不同发散角α的Vogel螺旋阵列：黄金角螺旋、τ螺旋、μ螺旋和π螺旋。这些螺旋由无理数生成，其结构介于短程关联强烈的准周期图案与长程无序的随机分布之间。通过保持平均最近邻间距 1>=5R不变，研究重点考察了其结构无序度的变化，该无序度通过粒子间距的标准差σd1量化。分析显示，从GA到π螺旋，σd1单调增加，意味着结构有序度逐步降低，局部环境均匀性变差。

图3:(a) 所考虑的不同Vogel螺旋的示意图。(b-e) 分别为GA(b)、τ(c)、μ(d)和π(e)螺旋的粒子间距离分布。橙色和蓝色点分别对应具有N=100和 N=1024纳米粒子的系统。(f)(b)-(e)中Vogel螺旋的粒子间距离标准差。

计算这些螺旋在中心颗粒被加热的初始条件下的温度演化，揭示了结构关联性对tTh的显著调控作用。尽管所有构型的最终整体衰减率λ1相近，但达到内部均匀化的瞬态路径截然不同。高度有序的GA螺旋展现出最平滑、最快速的内部热化（tTh≈0.074s），其动力学与周期阵列中心加热情形类似。随着发散角变化至π螺旋，结构无序度增加，热化过程变得缓慢且曲线呈现出更多起伏特征（tTh≈0.63s），这反映了中心颗粒所处局部环境的不规则性增强，可供高效近场耦合的最近邻数量相对减少，导致热流向外扩散的路径受阻、效率下降。

图4:与图2相同条件下，GA(a)、τ(b)、μ(c)和π(d) 螺旋在 1> =5R时的热化动力学。插图显示了相应螺旋几何形状。

为确立确定性非周期序的独特优势，研究进一步与完全随机的阵列进行了对比。在约束相同平均间距和最小间距的条件下，生成大量随机构型并统计其热化时间。研究发现，尽管某些随机实现的热化时间可能与μ螺旋相当，但随机系综的平均表现以及其最佳实现，均无法超越GA或τ螺旋所能达到的最快热化速度。这明确证实，具有精心设计的非周期关联性的结构，在优化近场辐射热传递动力学方面，能够超越无关联的完全随机体系。

最后，通过系统改变Vogel螺旋阵列的尺寸（从N=16到N=1024）并考察不同初始加热位置（中心或边缘），研究建立了结构无序度参数σd1与热化时间tTh之间普适的、单调递增的定量关系。该关系在不同的系统尺度下均保持稳健，清晰地表明：更高的结构无序度直接导致更慢的颗粒间热化。Vogel螺旋的核心价值在于，它通过一个简单的几何参数α，为实现从“有序”到“无序”的连续、可预测插值提供了理想平台，从而为近场热输运工程提供了一个前所未有的控制旋钮。

图5:对于具有(a)N=16、(b)N=64、(c)N=100、(d)N=256和(e)N=1024 纳米颗粒的阵列，热化时间tTh作为粒子间距离σd1标准差函数的关系。蓝色(红色)点对应于热量在阵列中心(边缘)的初始条件。对于随机阵列，tTh是在 214/N 个独立配置上进行平均，误差条对应于标准差。

结论与展望

综上所述，本研究证实了通过调控Vogel螺旋的发散角可连续改变其结构关联度，从而主动操控近场辐射热传递的弛豫动力学。较强的非周期序（如黄金角螺旋）表现出更快的颗粒间热化过程，其性能优于周期阵列与完全无序体系。这揭示了结构关联性是设计纳米尺度热输运行为的关键自由度。展望未来，需进一步考察实际制造中的尺寸与位置无序对模式寿命的影响，并将研究拓展至金属纳米颗粒等不同材料体系。该框架为开发基于确定性非周期序的动态热光子器件（如可编程热管理界面与纳米热开关）奠定了理论基础。