石墨烯表面电导率（局域近似）

COMSOL案例石墨烯可重构超反射镜 可提供详细讲解 可重构超透镜石墨烯这玩意儿在超表面领域简直就是外挂级材料。今天咱们来聊怎么用COMSOL整活——做个电压调谐的超反射镜。别被超表面吓到说白了就是通过纳米结构设计让光场跳舞。先上硬菜石墨烯的电导率模型得整明白。掏出Python算个基本参数import numpy as np def graphene_sigma(mu, gamma, T, omega): kB 8.617e-5 # eV/K e 1.602e-19 hbar 6.582e-16 # eV·s sigma_intra (1j * e**2 * kB * T) / (np.pi * hbar**2 * (omega 1j*2*gamma)) sigma_inter (e**2 / (4*hbar)) * (0.5 (1/np.pi)*np.arctan((hbar*omega - 2*mu)/(2*kB*T)) - 1j/(2*np.pi)*np.log((hbar*omega 2*mu)**2/((hbar*omega - 2*mu)**2 (2*kB*T)**2)))) return sigma_intra sigma_inter # 示例参数 mu 0.4 # 化学势0.4eV print(graphene_sigma(mu, 0.1e12, 300, 2*np.pi*1e12)) # 太赫兹波段这段代码的关键在于同时考虑带内intra和带间inter跃迁。带内项主导低频响应带间项在可见光更明显。调mu就像拧水龙头——化学势一变电导率立马翻脸。COMSOL建模时建议用边界偏微分方程PDE来刻画石墨烯层。在材料定义里塞进这个电导率模型记得把化学势mu设置为全局参数方便后面扫参。结构方面推荐C形开口环比方形结构更容易产生磁共振% COMSOL LiveLink脚本片段 model.component(comp1).geom.create(c1, Circle); model.component(comp1).geom(c1).set(r, 150e-9); model.component(comp1).geom.create(r1, Rectangle); model.component(comp1).geom(r1).set(size, [200e-9, 50e-9]); model.component(comp1).geom(c1).set(pos, [-100e-9, 0]);这种C形结构在2μm周期下对10μm中红外波的相位调控效果拔群。重点在于让结构尺寸小于波长但又不能太小导致加工翻车——玩的就是纳米尺度上的微操。COMSOL案例石墨烯可重构超反射镜 可提供详细讲解 可重构超透镜电压调谐才是灵魂所在。在COMSOL里耦合静电模块和波动光学静电场计算石墨烯的载流子浓度变化把结果映射到波动光学模块的化学势参数扫电压参数看反射率变化实测中发现边界条件设置有个坑石墨烯层要设为表面电流密度边界而不是薄层。薄层近似在强耦合时会翻车特别是当结构尺寸接近石墨烯费米波长时。最后给个性能参考当化学势从0.2eV调到0.6eV在3-5THz频段反射率能从15%飙到89%。不过实际调谐速度受限于离子凝胶层的响应用固态电解质比液态的靠谱虽然电导率差点但稳定性好。别问我怎么知道的——都是烧了三个样品的教训。搞超表面的兄弟记住仿真结果和实测差个20%都算正常。微纳加工时的侧壁粗糙度、石墨烯转移时的褶皱分分钟教你做人。所以参数优化时记得留点余量别把结构尺寸卡死在仿真最优值上。