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

仿真周报表面等离激元线光栅作者姬应科报告日期2017-9-25摘要模拟一束平面波入射到介电基底上的金属线光栅上，计算了不同入射角度情况下零阶和一阶衍射（透射和反射中）的衍射效率。通过施加额外的的背景散射场实现利用端口和散射边界完成散射，透射以及反射功率的积分求和。从二维到三维，利用周期性边界条件实现周期性光学微结构的模拟。从而大大减少计算所用的空间和时间。全局定义全局设置名称PlasmonicGratingPML.mph路径E:\姬应科--组会\研三组会\2017-9-25组会\Plasmonic_Grating_PML.mphCOMSOL版本COMSOL5.3(开发版本:223)使用的产品COMSOLMultiphysicsRFModule参数1参数名称表达式值描述na11空气折射率nb1.21.2底衬折射率d400[nm]4E−7

m光栅常数lam0441[nm]4.41E−7

m真空波长f0c_const/lam06.798E14

1/s频率alpha30[deg]0.5236

rad入射角变量Variables1选择几何实体层整个模型名称表达式单位描述alpha_pasin(sin(alpha)+lam0/(d*na))rad反射衍射角,m=1alpha_masin(sin(alpha)-lam0/(d*na))rad反射衍射角,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=-1Model1定义变量Variables2选择几何实体层整个模型名称表达式单位描述kaemw.k0rad/mPropagationconstant,airkaxka*sin(alpha)rad/mkxforincidentwavekay-ka*cos(alpha)rad/mkyforincidentwavekapxka*sin(alpha_p)rad/mkxforreflecteddiffractedwave,m=1kapyka*cos(alpha_p)rad/mkyforreflecteddiffractedwave,m=1kamxka*sin(alpha_m)rad/mkxforreflecteddiffractedwave,m=-1kamyka*cos(alpha_m)rad/mkyforreflecteddiffractedwave,m=-1kbnb*emw.k0rad/mPropagationconstant,dielectrickbxkb*sin(beta)rad/mkxforrefractedwavekby-kb*cos(beta)rad/mkyforrefractedwavekbpxkb*sin(beta_p)rad/mkxforrefracteddiffractedwave,m=1kbpy-kb*cos(beta_p)rad/mkyforrefracteddiffractedwave,m=-1kbmxkb*sin(beta_m)rad/mkxforrefracteddiffractedwave,m=-1kbmy-kb*cos(beta_m)rad/mkyforrefracteddiffractedwave,m=-1组件耦合积分1耦合类型积分算子名称intop1积分2耦合类型积分算子名称intop2选择BoundarySystem1Geometry1Geometry1单位长度单位m角度单位deg材料AirAir选择几何实体层域选择域3–4,6DielectricDielectric选择几何实体层域选择域1–2,5ElectromagneticWavesElectromagneticWavesEquationsWaveEquation,Electric1方程PerfectElectricConductor1方程WaveEquation,Electric2方程Port1方程Port2方程PeriodicCondition1方程ElectromagneticWaves2ElectromagneticWaves2EquationsWaveEquation,Electric1方程PerfectElectricConductor1方程WaveEquation,Electric2方程PeriodicCondition1方程ScatteringBoundaryCondition1方程Mesh1Mesh1Study1计算信息计算时间25sCPUIntel(R)Core(TM)i7-4790KCPU@4.00GHz,4

内核数操作系统Windows7FrequencyDomain频率(GHz)f0研究设置描述值包含几何非线性关物理场和变量选择物理场接口离散化ElectromagneticWaves(emw)physicsElectromagneticWaves2(emw2)physics网格选择几何网格Geometry1(geom1)mesh1研究扩展描述值辅助扫描开扫描类型指定组合助参数alpha参数值列表range(0,0.07853981633974483,1.4922565104551517)单位ResultsDataSetsSolution1解描述值解Solver1组件SavePointGeometry1数据集:Solution1CutLine2D1数据描述值数据集Solution1线数据描述值线定义方法两点点{{0,-6e-7},{4e-7,-6e-7}}高级描述值空间变量cln1x法向变量{cln1nx,cln1ny}数据集:CutLine2D1DerivedValuesLineIntegration1数据描述值数据集CutLine2D1表达式表达式单位描述-emw2.PoavyW/m-emw2.Poavy输出计算位置Table1TablesTable1LineIntegration1(-emw2.Poavy)Table1alpha-emw2.Poavy(W/m)0.00000.667060.0785400.674040.157080.596030.235620.535050.314160.488780.392700.450360.471240.416960.549780.386760.628320.358390.706860.331010.785400.304290.863940.278030.942480.251801.02100.225071.09960.197601.17810.169171.25660.139441.33520.108341.41370.0758121.49230.038372Table2Table2Column1Column20.00000.663870.0785400.670830.157080.593900.235620.533500.314160.487460.392700.449160.471240.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.039204绘图组1DPlotGroup5PlasmonicWireGrating