表面等离激元线光栅comsol仿真.docx

1、表面等离激元线光栅comsol仿真仿真周报表面等离激元线光栅作者姬应科报告日期2017-9-25 摘要模拟一束平面波入射到介电基底上的金属线光栅上，计算了不同入射角度情况下零阶和一阶衍射（透射和反射中）的衍射效率。通过施加额外的的背景散射场实现利用端口和散射边界完成散射，透射以及反射功率的积分求和。从二维到三维，利用周期性边界条件实现周期性光学微结构的模拟。从而大大减少计算所用的空间和时间。1全局定义全局设置名称Plasmonic Grating PML.mph路径E:姬应科-组会研三组会2017-9-25 组会Plasmonic_Grating_PML.mphCOMSOL 版本COMSOL

2、5.3 (开发版本: 223)使用的产品COMSOL MultiphysicsRF Module1.1参数 1参数名称表达式值描述na11空气折射率nb1.21.2底衬折射率d400nm4E7m光栅常数lam0441nm4.41E7m真空波长 f0c_const/lam06.798E141/s频率alpha30deg0.5236rad入射角1.2变量1.2.1Variables 1选择几何实体层整个模型名称表达式单位描述alpha_pasin(sin(alpha) + lam0/(d*na)rad反射衍射角, m=1alpha_masin(sin(alpha) - lam0/(d*na)rad

3、反射衍射角, m=-1betaasin(na*sin(alpha)/nb)rad折射角beta_pasin(na/nb*sin(alpha) + lam0/(d*nb)rad折射衍射角 m=1beta_masin(na/nb*sin(alpha) - lam0/(d*nb)rad折射衍射角, m=-12Model 12.1定义2.1.1变量Variables 2选择几何实体层整个模型名称表达式单位描述kaemw.k0rad/mPropagation constant, airkaxka*sin(alpha)rad/mkx for incident wavekay-ka*cos(alpha)ra

4、d/mky for incident wavekapxka*sin(alpha_p)rad/mkx for reflected diffracted wave, m=1kapyka*cos(alpha_p)rad/mky for reflected diffracted wave, m=1kamxka*sin(alpha_m)rad/mkx for reflected diffracted wave, m=-1kamyka*cos(alpha_m)rad/mky for reflected diffracted wave, m=-1kbnb*emw.k0rad/mPropagation con

5、stant, dielectrickbxkb*sin(beta)rad/mkx for refracted wavekby-kb*cos(beta)rad/mky for refracted wavekbpxkb*sin(beta_p)rad/mkx for refracted diffracted wave, m=1kbpy-kb*cos(beta_p)rad/mky for refracted diffracted wave, m=-1kbmxkb*sin(beta_m)rad/mkx for refracted diffracted wave, m=-1kbmy-kb*cos(beta_

6、m)rad/mky for refracted diffracted wave, m=-12.1.2组件耦合积分 1耦合类型积分算子名称intop1积分 2耦合类型积分算子名称intop22.1.3选择Boundary System 12.2Geometry 1Geometry 1单位长度单位m角度单位deg2.3材料2.3.1AirAir选择几何实体层域选择域 34, 62.3.2DielectricDielectric选择几何实体层域选择域 12, 52.4Electromagnetic WavesElectromagnetic WavesEquations2.4.1Wave Equati

7、on, Electric 1方程2.4.2Perfect Electric Conductor 1方程2.4.3Wave Equation, Electric 2方程2.4.4Port 1方程2.4.5Port 2方程2.4.6Periodic Condition 1方程2.5Electromagnetic Waves 2Electromagnetic Waves 2Equations2.5.1Wave Equation, Electric 1方程2.5.2Perfect Electric Conductor 1方程2.5.3Wave Equation, Electric 2方程2.5.4Pe

8、riodic Condition 1方程2.5.5Scattering Boundary Condition 1方程2.6Mesh 1Mesh 13Study 1计算信息计算时间25 sCPUIntel(R) Core(TM) i7-4790K CPU 4.00GHz, 4内核数操作系统Windows 73.1Frequency Domain频率 (GHz)f0研究设置描述值包含几何非线性关物理场和变量选择物理场接口离散化Electromagnetic Waves (emw)physicsElectromagnetic Waves 2 (emw2)physics网格选择几何网格Geometry

9、 1 (geom1)mesh1研究扩展描述值辅助扫描开扫描类型指定组合助参数alpha参数值列表range(0, 0.07853981633974483, 1.4922565104551517)单位4Results4.1Data Sets4.1.1Solution 1解描述值解Solver 1组件Save Point Geometry 1数据集: Solution 14.1.2Cut Line 2D 1数据描述值数据集Solution 1线数据描述值线定义方法两点点0, -6e-7, 4e-7, -6e-7高级描述值空间变量cln1x法向变量cln1nx, cln1ny数据集: Cut Lin

10、e 2D 14.2Derived Values4.2.1Line Integration 1数据描述值数据集Cut Line 2D 1表达式表达式单位描述-emw2.PoavyW/m-emw2.Poavy输出计算位置Table 14.3Tables4.3.1Table 1Line Integration 1 (-emw2.Poavy)Table 1alpha-emw2.Poavy (W/m)0.00000.667060.0785400.674040.157080.596030.235620.535050.314160.488780.392700.450360.471240.416960.549

11、780.386760.628320.358390.706860.331010.785400.304290.863940.278030.942480.251801.02100.225071.09960.197601.17810.169171.25660.139441.33520.108341.41370.0758121.49230.0383724.3.2Table 2Table 2Column 1Column 20.00000.663870.0785400.670830.157080.593900.235620.533500.314160.487460.392700.449160.471240.

12、415700.549780.385440.628320.357240.706860.330320.785400.304060.863940.277990.942480.251791.02100.225081.09960.197631.17810.169331.25660.139671.33520.108581.41370.0755311.49230.0392044.4绘图组4.4.11D Plot Group 5 Plasmonic Wire Grating4.4.2Electric Field (emw)20: alpha=1.4923, freq=6.798E14 表面: Electric field norm (V/m)4.4.3Electric Field (emw2)1: alpha=0, freq=6.798E14 表面: Electric field norm (V/m)4.4.4一维绘图组 8线图: 积分 1 (W/m)