6，comsol超表面仿真：带积分号的石墨烯模型构建与分析

在上一篇文章中，那个石墨烯电导率公式不带积分号，有的文章中石墨烯电导率是带积分号的。比如今天这篇文章《Tunable terahertz band-stop filter based on self-gated graphene monolayers with antidot arrays 》。

Tunable terahertz band-stop filter based on self-gated graphene monolayers with antidot arrays.pdf

文中给出石墨烯电导率公式如下

看到积分号，首先回会想到matlab。但其实不借助matlab，完全在comsol中也能处理这种积分号问题。观察上面两个积分式，发现积分上限为无穷，这说明被积分的在无穷远处可能是收敛于0的。

下图是（2）式中 被积分函数 与epsilon的关系。（uc=0eV,gamma=0.1meV）

 下图是（3）式中 被积分函数 与epsilon的关系。

看到上面两幅图，就不用怕那个积分号了。

下图是我复现结果与原文对比

文章模型如下

在相对介电常数为epsilon1和epsilon2的介质基板，在epsilon1的介质板上下有石墨烯，石墨烯不是完整的，而是中间有个圆孔。太赫兹波垂直照下来，发现有一个很窄的吸收峰如上图透射光谱dB所示。上图是周期性结构的一个单元。