FDTD中反射率仿真
本案例以WO3/W薄膜为例,介绍FDTD中反射率测量的主要过程。软件版本为Lumerical的FDTD Solutions 2020a。

下面介绍主要步骤:
- 导入材料属性
由于WO3材料在FDTD的材料库中没有内置,需要自己查找并导入,对于一般材料可以从如下网站中查找相应折射率:
http://www.ioffe.ru/SVA/NSM/nk/
点击Materials,

在弹出的窗口中,点击Add,在弹出的选项框中选择Sampled 3D data。

点击Import,在Select File中选择折射率文件,导入即可。
[图片上传失败...(image-2f59ec-1587692644459)]
- 添加结构
结构添加通过Structures添加,选择Rectangle,添加WO3薄膜,

设置结构参数,可根据图中参数进行修改,

材料选择刚才导入的WO3。

其他两项设置默认即可,不用修改。
相同的思路,添加W层,WO3是在W上方,因此设置的时候Z方向的数值应契合好。相关设置如下:

W材料是FDTD内置的,可以直接从中选取。

- 设置FDTD Region
点击Simulation,添加FDTD。由于是薄膜结构,X和Y方向是无限延伸的,因此X,Y方向设置成周期性边界条件,周期性边界条件的情况下PML可以设置成steep
angle,以提高计算速度。

结构尺寸上按需求设置,X,Y方向没有特别要求,Z方向最小值应设置在W层中,最大值在WO3上方,同时远离结构至少半个波长。其他不需要特别修改,默认即可。

- 添加Mesh
点击Simulation右方的下拉菜单,添加Mesh。建模设计过程没有明显的先后次序要求,Mesh什么时候设置都可以。网格有两种添加方案,一是设置Mesh的尺寸,二是基于结构,这里主要考虑WO3/W薄膜反射率,可以把mesh加到WO3上。


- 光源设置
点击Sources,加入光源。光源选择平面波,入射方向设置z负方向。光源尺寸大于FDTD尺寸(如果没有大于会在仿真时自动拓展),光源放置在结构上方,具体可以对比结构尺寸和光源尺寸查看。

随后设置入射光波长范围,这里设置成0.4~0.9。

- 监视器设置
主要设置反射率监视器。
监视器类型选择Frequency-domain field and power。

监视器选择Z-normal,XY平面尺寸大于FDTD尺寸,大了无所谓,超出的部分不会记录数据,主要是设置Z方向的位置。反射率监视器在光源和FDTD上边界中间。

在General选项中,点击Set global monitor settings,

Frequency points数量设置为50,这个数值越大,谱线越光滑。

- 检查运行内存
点击Check,选择Check simulation and memory requirements,

保证mesh和running
simulation不超过电脑内存大小,否则无法计算,需要降低精度才能计算。

- 运行
上一步检查没有问题的话,点击Check旁边的Run,运行即可,等待运行结束。

一般会提前结束,提前结束的条件是Autoshutoff的数值满足FDTD中设置的停止条件。

- 结果分析
右键点击R,可以查看该监视器平面的电场磁场等信息,选择T,查看反射率信息。

反射率曲线如下图所示

另外,FDTD仿真得到的计算结果,也可以导入Matlab中进行进一步的数据处理。
欲将反射光谱导入Matlab,可以在script prompt处输入以下语句:
R0=getresult('R','T');
R=R0.T;
lambda=R0.lambda;
matlabsave('WO3.mat',R,lambda);
原文链接:FDTD测反射率
网友评论