Comsol石墨烯二维材料 包含太赫兹德鲁得和近红外Kubo两种模型咱们今天来聊聊COMSOL里整石墨烯二维材料的活儿特别是太赫兹波段的德鲁得模型和近红外的Kubo模型。这俩模型在模拟石墨烯的光电特性时就像火锅配香油碟——绝配但用起来得看准场景下菜。先看德鲁得模型这老伙计太赫兹波段模拟它最拿手。在材料属性里填参数的时候记得把石墨烯当成二维面电流来处理。这里有个典型的代码片段sigma_drude (e**2 * tau * Ef)/(pi * hbar**2 * (1 - 1j*omega*tau))这串式子里的tau弛豫时间和Ef费米能级是灵魂参数。之前有个哥们把单位搞混了用eV当费米能级结果算出来的电导率直接飞升到外太空。特别注意当频率omega接近1/tau时这个模型就开始闹脾气不准了。转场到近红外波段就得请出Kubo公式了。这里有个隐藏关卡——温度参数T的设置。见过最离谱的错误是有人用300开尔文模拟低温石墨烯结果光电导曲线抖得跟心电图似的。正确的打开方式应该是kubo_cond (omega) (e^2/(4*hbar}))*(... % 这里展开久保公式的积分项 fermi_dirac(Ef/(2*kB*T)) ... % 费米分布处理 (hbar*omega)./(2*kB*T).*sech((hbar*omega)./(4*kB*T)).^2);这段代码里的温度T和频率omega的耦合关系是重点。有个实战技巧当模拟波长在1550nm附近时记得把网格尺寸压到纳米级不然场分布会糊成一团马赛克。Comsol石墨烯二维材料 包含太赫兹德鲁得和近红外Kubo两种模型模型切换有个经典坑位在1-3THz这个过渡区间两种模型的结果可能会打架。这时候不如把两个模型的结果叠在一起做个加权平均就像下面这样处理double hybrid_sigma alpha*drude_sigma (1-alpha)*kubo_sigma; // alpha随频率渐变 0到1过渡这么干的时候要注意过渡区的平滑性检验别搞出二阶导数不连续的幺蛾子。最后说个骚操作在频域研究里同时激活两个材料模型用if语句根据频率自动切换。不过记得在结果分析时把交界处的数据点加密不然拟合曲线会露出破绽。COMSOL的模型树里设置条件判断时用(freq 1e12)[Hz]这种写法比直接写数值更保险单位换算的坑能避开一大半。说到底模型选对事半功倍参数调准锦上添花。下次遇到石墨烯模拟不准先检查费米能级单位是不是eV再看温度设置是不是跟实验条件对得上——这两个雷区踩中的人能凑个加强连了。