石墨烯光子晶体仿真

这个仿真旷日良久，现在还没有一个定论。
目前，石墨烯在COMSOL中通过一个导电平面表示，三维结构，四壁采用的是完美电导体，上下是散射边界条件，尺寸是外面正方形5um和中间的圆3um。我们计算本征频率，得到的电场分布图很有规律，相关的仿真文件及更进一步的解释在坚果云上面。

当四周采用完美电导体的时候，本征频率普标偏高，起步在3.077THz；同样的结构用散射边界条件，本征频率就成为了2.458THz，电场分布几乎没什么区别。
采用了周期性边界条件后，奇怪的就是没有规律了。

关于石墨烯的表示方法，COMSOL方便的是使用一个电流平面，在ewfd里增加一个表面电流密度，然后定义ewfd.Ex*sigma，这里的sigma是-e_const^2*Ef*i/((pi*hbar_const^2)*(ewfd.omega-i/tao))。

说明文档：https://www.jianguoyun.com/p/DWjWE-QQtNTQBhitqz4
仿真例程https://www.jianguoyun.com/p/DZcb5lMQtNTQBhipqz4