CST仿真eit电磁诱导透明(包括石墨烯的建模) EIT石墨烯电磁诱导透明案例搞EIT仿真的都知道传统金属结构虽然经典但石墨烯的可调性才是现在的香饽饽——靠栅压就能调费米能级相当于给器件装了个电控遥控器在传感器、慢光器件里简直是开挂般的存在。最近用CST撸了个石墨烯EIT的实操案例踩了几个坑分享下实打实的细节少讲理论多讲干活的步骤。首先得解决最核心的问题CST里怎么建石墨烯别找内置材料库了根本没有现成的石墨烯选项必须用表面阻抗模型。石墨烯的表面阻抗和费米能级、频率直接挂钩公式是 \( Zs \frac{j\pi\hbar^2 f}{e^2 EF} \)CST里得把这个公式怼进材料设置里用VBA脚本参数化会比手动输入靠谱10倍避免手敲公式出错。先上一段自动创建石墨烯材料和结构的脚本 一键生成石墨烯表面结构材料 Sub CreateGrapheneEITStruct() 1. 创建自定义石墨烯表面阻抗材料 Dim grapMat As Material Set grapMat Materials.Add(Graphene_Tunable) grapMat.SurfaceImpedance.Active True 开启表面阻抗模型 核心写入石墨烯表面阻抗公式参数化费米能级E_F Dim E_F As Double E_F 0.6 初始费米能级设为0.6eV后续可直接改这个值调谐 Dim hbar As Double, e As Double hbar 1.0545718e-34 约化普朗克常数 e 1.602176634e-19 电子电荷量 把公式转成CST支持的表达式语法f是CST内置的频率变量 grapMat.SurfaceImpedance.Zs j * pi * CStr(hbar^2) * f / ( CStr(e^2) * CStr(E_F) ) grapMat.SurfaceImpedance.ZsUnit Ohm/square 阻抗单位是每平方欧姆 2. 创建EIT核心结构石墨烯条带谐振环 石墨烯条带亮模谐振器 Dim strip As Object Set strip Objects.AddRectangle(Graphene_Strip, , , 0, 0, 0, 20e-6, 2e-6, 0) strip.Material grapMat.Name strip.Thickness 0 表面模型必须设厚度为0不然会触发体材料求解bug 石墨烯谐振环暗模谐振器和条带间距1μm Dim ring As Object Set ring Objects.AddRectangle(Graphene_Ring_Outer, , , 0, 3e-6, 0, 18e-6, 18e-6, 0) Dim ringInner As Object Set ringInner Objects.AddRectangle(Graphene_Ring_Inner, , , 0, 4e-6, 0, 16e-6, 16e-6, 0) Objects.Subtract Graphene_Ring, Graphene_Ring_Outer, Graphene_Ring_Inner Objects(Graphene_Ring).Material grapMat.Name Objects(Graphene_Ring).Thickness 0 End Sub这段脚本的核心坑点得拎出来说为啥Thickness要设为0因为表面阻抗模型是二维面模型设成实际的0.34nm厚度会触发CST的体材料求解逻辑超薄材料的网格划分会直接崩掉仿真结果全是噪声。表面阻抗公式里的f是CST内置的频率变量求解器会自动在每个频率点计算对应的阻抗完美匹配石墨烯的色散特性比手动输入固定阻抗靠谱太多。接下来是仿真设置EIT是窄带效应用频域求解器效率最高别用时域求解器浪费时间。再上一段一键配置仿真的脚本 一键配置EIT仿真参数 Sub SetupSimulation() 1. 设置频率扫描范围EIT窗口通常在THz波段这里设0.8-1.2THz Solver.FrequencyRange.SetRange 0.8e12, 1.2e12 Solver.FrequencyRange.SamplingPoints 2001 密采点才能看清窄带透明窗口 2. 边界条件自由空间仿真用Open (add space)避免边界反射干扰 Boundaries.SetAll Open (add space) Boundaries.SetSpaceDistance 10e-6 边界到结构留10μm空间足够消反射 3. 激励源平面波垂直入射偏振方向平行于石墨烯条带 Dim planeWave As Object Set planeWave Excitations.AddPlaneWave planeWave.Direction Array(0, 1, 0) 沿y轴入射 planeWave.Polarization Array(1, 0, 0) x方向偏振耦合条带的亮模 4. 求解器收敛精度EIT是弱耦合效应精度设到-40dB足够捕捉窄带峰 Solver.ConvergenceCriterion -40 End Sub这里要注意平面波的偏振方向必须和亮模谐振器条带的方向一致不然激不起来亮模根本出不来EIT透明窗口——我一开始就因为偏振设反了扫了半小时全是平的透射谱差点砸键盘。CST仿真eit电磁诱导透明(包括石墨烯的建模) EIT石墨烯电磁诱导透明案例跑完仿真看结果透射谱里会在~1THz附近出现一个窄窄的透明窗口这就是EIT效应亮模条带和暗模谐振环的谐振频率接近发生相消干涉大部分能量直接透过去只有极窄的频段被吸收。最能体现石墨烯优势的是可调性改费米能级就能挪EIT窗口的位置。用这段脚本自动扫不同费米能级 参数化扫描费米能级自动生成多组透射曲线 Sub ScanFermiLevel() Dim E_F_list As Variant E_F_list Array(0.4, 0.5, 0.6, 0.7) 扫0.4到0.7eV的费米能级 Dim grapMat As Material Set grapMat Materials(Graphene_Tunable) For Each E_F In E_F_list 更新石墨烯表面阻抗的费米能级参数 grapMat.SurfaceImpedance.Zs j * pi * 1.0545718e-34^2 * f / (1.602176634e-19^2 * CStr(E_F) ) 跑仿真存结果 Solver.Run Results.ExportToFile Transmission_EF CStr(E_F) eV.txt, 1D Results\Transmission Next E_F End Sub扫出来的曲线会明显看到费米能级从0.4eV升到0.7eVEIT透明窗口直接从0.9THz蓝移到1.1THz——这就是石墨烯电控调谐的核心金属结构得拆了重加工石墨烯加个电压就搞定这才是未来器件该有的样子。最后补两个踩过的坑求解器别选时域EIT是窄带效应时域求解器要扫很长时间才能收敛频域求解器直接针对窄带频段计算效率至少高3倍。网格划分要加加密区石墨烯结构是亚微米级的必须在结构周围设局部加密网格不然谐振峰直接被磨平透明窗口连影子都看不到。总的来说CST搞石墨烯EIT的核心就是把表面阻抗模型玩明白结构上保证亮模和暗模的耦合剩下的就是用脚本把重复操作自动化省下来的时间多扫几组参数不香吗要是感兴趣还能把结构改成哑铃型、多环阵列看看EIT窗口的分裂情况可玩性拉满。