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里耦合静电模块和波动光学:
- 静电场计算石墨烯的载流子浓度变化
- 把结果映射到波动光学模块的化学势参数
- 扫电压参数看反射率变化
实测中发现边界条件设置有个坑:石墨烯层要设为表面电流密度边界,而不是薄层。薄层近似在强耦合时会翻车,特别是当结构尺寸接近石墨烯费米波长时。
最后给个性能参考:当化学势从0.2eV调到0.6eV,在3-5THz频段反射率能从15%飙到89%。不过实际调谐速度受限于离子凝胶层的响应,用固态电解质比液态的靠谱,虽然电导率差点但稳定性好。别问我怎么知道的——都是烧了三个样品的教训。
搞超表面的兄弟记住:仿真结果和实测差个20%都算正常。微纳加工时的侧壁粗糙度、石墨烯转移时的褶皱,分分钟教你做人。所以参数优化时记得留点余量,别把结构尺寸卡死在仿真最优值上。


被折叠的 条评论
为什么被折叠?



