Learning Forum

[202310263]

老师您好，我在仿真复刻论文的时候，使用FDTD画能带图时遇到如图所示的问题，参数如下：基地介电平板折射率为1，光子晶体为折射率8.9，半径0.2a的圆柱，a是晶格常数，这里取a=0.5，但是spectrum返回值总是为空。

[Guilin Sun]

这种情况一般是因为仿真时间太短，扫描没有提到数据造成的。你可以打开已经仿真过的扫描文件看看。

网上例子所用的仿真时间是针对该例子的，你现在用如此高的折射率，仿真时间需要适当增加。但也不是要求很长时间，因为计算能带只需要找到谐振波长即可，一般autoshutoff 数值在0.95左右就可以了，你还要根据最后的能带结构是否有丢失点，是否有额外点造成不易找到能带等适当调整仿真时间。能带仿真是少数几个靠仿真时间结束仿真的例子之一，不需要达到autshutoff min 1e-5。

[202310263]

好的，谢谢老师，按照您的指导从1000增加到10000之后仍然存在一些问题，而且能带图也有一条的一段缺失，此外，在画各条线的值时存在autoshutoff大部分大于1的情况，想请问一下老师是什么问题

[Guilin Sun]

很明显分析组中的脚本出问题了。请到分析组脚本处，点击run analysis，再看看是什么问题。估计是不是里面不小心把额外字符加里面了。

在autoshutoff大部分大于1如果很接近1 ，可能问题不大，如果远大于1 ，可能是要发散了。你需要找如何修正发散的方法，例如参考这个帖子：

Ansys Insight: 修改PML哪些参数可以有效抑制FDTD发散

Ansys Insight: FDTD 和 varFDTD 常见的发散原因以及处理

 

[202310263]

老师您好，谢谢您给的建议，后来根据老师的建议重新调整运行之后，脚本并没有任何问题，将步长稍作调整之后能带图也总是缺少某些东西，看起来很不对劲，想请问老师是什么问题，还是矩形光子晶体，折射率为8.9，考虑过是不是因为点源问题，将点源增加到20之后能带图也没有什么变化，之前也考虑过是不是带宽的问题，也做修改之后，将带宽调整为50-700之后才能得到较为完整的图。以下是（0，0）-（0，0.4725）-（0.4725，0.4725）-（0，0）的路径，模式为TE偏振光，频率范围是50-700THz。并且当论文中做归一化时：纵坐标为wa/2*pi*c，也有一点数值不对应的问题。