湖北大学吴小虎MTP：表面电导剪切重构近场热光子隧穿

论文息：Lu Yang, Zihao Liu, Leyong Jiang, and Xiaohu Wu, “Surface-conductivity shear rewires near-field thermal photon tunneling,” Materials Today Physics, accepted for publication (2026).

研究背景

近场辐射热传递的本质并不仅仅是“热辐射变强”，而是大量具有不同频率和面内动量的倏逝电磁通道在纳米间隙中发生隧穿。过去利用各向异性材料调控这些通道时，最常见的做法是直接旋转发射端和接收端，使两侧极化激元模式发生匹配或失配。但这种“整体扭转”实际上改变的是两个端体之间的相对取向。一个更基础的问题是：如果发射端和接收端始终保持整体对齐，仅仅改变每个端体内部二维导电层与各向异性基底之间的相对方向，能否重新组织热光子通道？这一问题把热扭转调控从端体之间的模式匹配，进一步推进到单个复合界面的张量结构设计。

研究内容

本文构建了由两个完全相同复合端体组成的近场辐射体系，每个端体均由一层各向异性二维导电层和半无限各向异性介质基底组成，中间由真空间隙隔开。二维导电层采用具有黑磷代表性参数的各向异性面电导模型，基底选用双轴α-MoO₃。与传统扭转结构不同，两侧二维导电层始终相互平行，两侧基底也保持相同取向；调控参数θ只描述每个端体内部二维导电层与基底之间的局域错轴。因此，改变θ不会额外引入发射端与接收端之间的整体失配，而是从内部改变每一个复合界面的电磁响应。

图1. 表面电导剪切调控近场辐射热传递的结构示意图。两个完全相同的复合端体由各向异性二维导电层和各向异性基底组成，中间由宽度为d的真空间隙隔开。二维导电层的主轴定义实验室坐标系。二维导电层保持固定，而两个各向异性基底均相对于相邻二维导电层的主轴旋转相同角度θ。上、下端体温度分别为T1和T2，同时两端体在真空间隙两侧始终保持整体对齐。

图2. （a）基于黑磷参数的二维导电层模型中，主轴方向表面电导虚部Im(σxx)和Im(σyy)随角频率的变化，其中实线表示Im(σxx)，虚线表示Im(σyy)，单位为mS。载流子浓度分别为n = 5、25和50 × 10¹² cm⁻²。蓝色、青绿色和紫色分别表示逐渐增大的载流子浓度，上方横坐标给出了对应的光子能量。（b）α-MoO₃主轴方向介电常数实部Re(εxx)、Re(εyy)和Re(εzz)。蓝色、青绿色和紫色分别表示x、y和z方向的响应，阴影区域表示相应主轴方向选择性的Reststrahlen带。

这一局域错轴会产生一个关键的张量效应。当二维导电层本身满足σx ≠ σy，而基底满足εx ≠ εy时，中间扭转角使面电导张量与介电张量不再共轴，两者无法通过同一个坐标旋转同时对角化。此时在基底坐标系中出现横向面电导分量σxy = (σy − σx)sinθcosθ，从而形成互易的“表面电导剪切”。这种剪切并不是简单把原有极化激元色散整体旋转，而是将二维材料的电子各向异性与α-MoO₃的声子各向异性进行重新混合，改变复合界面的杂化等离激元—声子极化激元结构。

图3. 在T0 = 300 K时，线性化倏逝热传递系数H随真空间隙宽度d的变化。对于复合端体，两侧均处于共轴构型，即θL = θU = 0°，并在整体上保持对齐。实线表示两个完全相同的BP/α-MoO₃复合端体，载流子浓度分别为n = 5、25和50 × 10¹² cm⁻²。粉色虚线和灰色点线分别表示n = 50 × 10¹² cm⁻²时仅含BP二维导电层的对照结构和仅含α-MoO₃基底的对照结构。水平虚线表示黑体热传递系数hBB = 6.12 W m⁻² K⁻¹。

图4. 两个完全相同的BP/α-MoO₃复合端体在d = 10 nm、T0 = 300 K条件下的光谱倏逝热传递系数hω。（a）—（d）分别对应局域二维导电层—基底扭转角θ = 0°、30°、60°和90°。蓝色、青绿色和红色曲线分别表示载流子浓度n = 5、25和50 × 10¹² cm⁻²。

数值结果首先揭示了载流子浓度与真空间隙之间非常明显的协同效应。在300 K、d = 10 nm时，载流子浓度n = 5×10¹² cm⁻²对应的热传递系数达到5.70×10⁴ W m⁻² K⁻¹，约为黑体热传递系数的9.32×10³倍；随着载流子浓度增加，10 nm间隙下的总热传递反而下降。但当间隙增加到100 nm时，这一排序完全反转，高载流子浓度对应更高的热传递系数。其物理原因不是简单的“载流子越多，等离激元越强”，而是真空间隙本身具有动量滤波作用：随着距离增大，大波矢通道被指数压制，不同载流子浓度所对应的频率—动量相空间因此被重新筛选。

图5. 低载流子浓度下的动量空间能量传输系数（ETC）重构。图中给出了n = 5 × 10¹² cm⁻²、d = 10 nm条件下，局域扭转角θ分别为0°、30°、60°和90°时的p偏振投影倏逝能量传输系数ξ(ω,kx,ky)，各角度由上至下排列。不同列对应光谱热传递系数hω中的四个代表性特征。频率单位为10¹⁴ rad s⁻¹，波矢单位为100k0。多色实线和青色虚线分别表示简化复合界面极化激元和孤立BP二维导电层的等频轮廓（IFCs）。

图6. 高载流子浓度下的动量空间能量传输系数（ETC）重构。图中给出了n = 50 × 10¹² cm⁻²、d = 10 nm条件下，局域扭转角θ分别为0°、30°、60°和90°时的p偏振投影倏逝能量传输系数ξ(ω,kx,ky)，各角度由上至下排列。不同列对应光谱热传递系数hω中与高载流子浓度相关的四个代表性特征。频率单位为10¹⁴ rad s⁻¹，波矢单位为100k0。多色实线和青色虚线分别表示简化复合界面极化激元和孤立BP二维导电层的等频轮廓（IFCs）。

更重要的是，局域扭转对积分热传递系数的影响虽然只有几个百分点，却能显著重构谱空间和动量空间中的热光子隧穿通道。随着θ变化，不同谱峰并不是整体平移，而是出现选择性的上移和下移；在动量空间中，能量传输系数分布则发生明显的位移、倾斜、偏转和对称性重构。θ = 0°和90°时，σxy消失，kx和ky方向的镜面对称性恢复；在30°和60°等中间角度，有限的表面电导剪切混合两个面内方向，使高透射通道偏离晶体主轴并发生倾斜。由此形成一个清晰的层级调控图景：载流子浓度决定可参与热输运的频率和动量区域，真空间隙筛选其中能够隧穿的部分，而局域张量错轴进一步决定这些剩余通道在动量空间中的方向、对称性和衰减结构。

图7. 随频率变化的简化对称轴偏转诊断量。（a）n = 5 × 10¹² cm⁻²；（b）n = 50 × 10¹² cm⁻²。不同曲线对应θ = 10°至90°、以10°为步长变化的局域二维导电层—基底扭转角。α相对于二维导电层x轴定义，并折叠到[-45°，45°)区间。青色虚线表示θ = 90°，此时σxy = 0，折叠后的诊断量α = 0°。

图8. 载流子浓度对简化对称轴偏转诊断量的调控。（a）θ = 30°；（b）θ = 60°。颜色表示基于BP参数的二维导电层载流子浓度，其范围为5 × 10¹²至50 × 10¹² cm⁻²。α相对于二维导电层x轴定义，并折叠到-45°，45°)区间。

结论与展望

这项工作的核心并不是获得一个特别大的“扭转开关比”，而是重新定义了近场热辐射中的扭转自由度。传统热扭转主要依赖发射端和接收端整体旋转，从而改变跨间隙模式匹配；本文则将调控位置推进到单个复合界面内部，通过二维面电导张量与基底介电张量的非共轴排列产生互易表面电导剪切。其结果是，即使两端在宏观上完全对齐，热光子的频谱位置、动量方向、通道拓扑和镜面对称性仍可被重新编程。尤其值得注意的是，积分HTC变化较小并不意味着电磁响应弱，而是不同动量通道在积分后发生补偿。文章因此提示，评价近场热调控不能只看一个标量热传递系数，还应关注其背后的频率—动量通道结构。更进一步，载流子浓度、真空间隙和局域张量错轴形成了三个层级清晰、功能不同的调控维度，为可编程各向异性热光子学和热扭转学提供了一种新的局域张量设计思路。