太赫兹频域超表面-石墨烯混合结构表面阻抗预测

摘要

精准控制石墨烯基混合超表面表面阻抗对下一代太赫兹（THz）光电子学至关重要。本工作推导方形贴片超表面闭式解析阻抗表达式，结合石墨烯本征片阻抗建立混合阻抗整体模型。全波仿真验证解析预测在宽频段准确性。研究贴片尺寸和石墨烯化学势对复合阻抗影响，为利用工程表面阻抗的吸收体等THz器件提供设计准则。

引言

人工阻抗超表面因其独特操控电磁波能力近年受广泛关注，应用于全息、高分辨率成像、地毯隐身、宽带吸收等。然而常规超表面固有色散特性限制其窄带工作。石墨烯凭借电压可调电导率提供动态平台，克服此限制，实现THz和光谱段等离子体和超表面特性实时控制。因此，混合超表面-石墨烯结构成为可调、紧凑THz器件有力候选。

文献中已开发超表面和石墨烯片等效表面阻抗解析模型，但能准确捕捉方形贴片超表面与石墨烯片相互作用的统一模型仍缺失。本研究推导综合解析框架并经全波仿真验证，填补这一空白。

方法

方形贴片与石墨烯片阻抗

典型超表面-石墨烯吸收体由薄图案层置于金属接地面上方，中间介质间隔厚度h隔开。单元为边长D、窄槽宽w的方形贴片（图1a）。在D≪λ条件下可构建传输线等效电路（图1b）。单元输入阻抗表达为：

1/Z_in = 1/Z₁ + 1/Z_mg = 1/(jZ_h tan(k_zh h)) + 1/Z_mg (1)

其中Z_h和k_zh为基底阻抗和传播常数。正入射总吸收率由反射系数S₁₁得出：

A(ω)=1-|S₁₁|²=1-|(Z_in-120π)/(Z_in+120π)|² (2)

式(1)使从模拟S参数提取混合表面阻抗Z_mg成为可能，揭示贴片几何与石墨烯电导率相互作用。

方形贴片阻抗

电致密完全导电贴片阵列网格阻抗为：

Z_m = (D/w)Z_s - j(η_eff)/(2α) (3)

有效波阻抗η_eff = √(μ₀/ε₀ε_eff)，几何因子

α = (k_eff D/π) ln[1/sin(πw/2D)] (4)

其中ε_eff≈(ε_r+1)/2，k_eff=k₀√ε_eff。这些表达式适用于贴片周期远小于工作波长时。

验证解析模型，在Ansoft HFSS中全波仿真：电阻片Z_s=35 Ω/sq，间隔层20 µm，相对介电常数ε_r=3.2(1-j0.045)。从模拟输入阻抗提取贴片阻抗与解析预测紧密匹配，确认模型准确性。

石墨烯片阻抗

石墨烯视为无限薄导电表面。其面电导率由Kubo公式描述：

σ_g = (j e² k_B T)/(πħ²(ω+j/τ)) [μ_c/(k_B T) + 2 ln(e^{-μ_c/k_B T}+1)] + (j e²)/(4πħ) ln[(2|μ_c|-ħ(ω+j/τ))/(2|μ_c|+ħ(ω+j/τ))] (5)

假设室温（T=300 K）、弛豫时间τ=0.1 ps，片阻抗Z_g=1/σ_g=R_g+jX_g。图4显示电阻和电抗随化学势μ_c单调递减，固定μ_c下在0.2-6 THz带近乎平坦。

结果与讨论

不同于传统假设混合阻抗仅为贴片与石墨烯阻抗并联，仿真揭示更复杂相互作用。图5对比不同化学势（w=19 µm）下Z_mg实部虚部解析与模拟值。解析趋势捕捉整体频率依赖，但绝对值差异显著，尤其低μ_c下石墨烯阻抗主导时。

为进一步阐明尺寸依赖性，固定μ_c=0.4 eV下贴片宽度从17变至19.5 µm提取Z_mg（图7）。实部随w增加至0.31 THz先减后增，反映容性耦合与石墨烯片电阻平衡。虚部遵循相似趋势，确认几何与材料响应相互作用。

3 THz电流密度分布（图8）显示较大贴片支持较低表面电流，符合式(11)给出贴片长度与电流幅值反比关系。HFSS积分电流值佐证此趋势，强化阻抗变化物理解释。

结论

我们开发并验证了预测THz频域超表面-石墨烯混合结构表面阻抗的综合解析框架。方形贴片阻抗推导公式结合Kubo基石墨烯片阻抗，使吸收体、天线等器件快速设计无需详尽全波仿真。研究阐明贴片尺寸和石墨烯化学势如何联合影响有效阻抗，为可调THz元件提供清晰设计杠杆。

缩略语

HFSS

高频结构仿真

TEM

横向电磁波

THz

太赫兹