使用光子晶体工程光流

2019年7月25日

1980年,来自贝尔通讯研究的Eli Yablonovitch思考了如何在特定频率范围内减少半导体激光器中的损耗。他在透明的介质中切成周期性的圆孔,并观察到它不允许频率范围导致损失通过。Yablonovitch发现,这些结构与带有传导和价带的半导体相似,并将其命名光子晶体(与普林斯顿大学的Sajeev John一起)。让我们讨论光子晶体控制光的三个例子。

编者注:此博客文章于2020年11月18日更新,以反映更新的建模功能。

引入光子晶体

光子晶体是光学结构,定期重复以实现不同的应用,例如弯曲,滤波和镜子。光子晶体的最简单配置之一是建立两种不同的介电材料,这些介电材料会在单个方向上定期重复(+/--y该模型的方向),也称为多层膜。如参考文献第4章所述,我们基于两层(蓝色和绿色)的相对介电常数对1D光子晶体进行了三种不同的分析。1。

  1. 情况1:\ epsilon_ {1}= 13和\ epsilon_ {2}= 13
  2. 案例2:\ epsilon_ {1}= 13和\ epsilon_ {2}= 12
  3. 案例3:\ epsilon_ {1}= 13和\ epsilon_ {2}= 1

为了简单起见,与其对单个方向重复的多层膜重复的实际3D几何形状进行建模,不如对2D横截面的几何形状进行建模。这Floquet周期性条件应用于+/-y方向,而连续性周期性条件应用于+/-X光子晶体的方向。这连续性周期性条件允许将晶胞扩展到+/-中的无穷大。X通过将相同的解决方案应用于源和目的地边界来方向。而Floquet周期性条件允许在+/-中无限重复单位电池y通过在源边界和目的地边界之间应用适当的相位变化来方向。

特征频率分析是用从中扫除的波矢量(KY)进行的\ frac { - \ pi} {a}\ frac {\ pi} {a},其中a是层的厚度(该模型中的200 nm)。特征频率分析在100 thz左右进行,并搜索更高数量的本征频(在此模型中8)。可以使用代表特征频率的“ FREQ”表达式为每种情况绘制分散关系。

从分散关系,可以推断出第一种情况,而在材料属性被认为是均匀的情况下,模式沿着光线为(\ Omega(k)= \ frac {ck} {\ sqrt {\ epsilon}})。在第二种情况下,相对介电常数的轻微变化导致在105 THz和107 THz之间形成光子带隙。对于第三种情况,由于相对介电常数的显着差异,光子带隙进一步扩散,拒绝了114 THz和192 THz之间的频率,最终是作为滤波器的工作。

1D光子晶体模型作为3D和2D几何形状的可视化。
1D光子晶体的实际3D几何形状和Comsol多物理学中的2D几何形状的示意图。

3个并排图显示了三种一维光子晶体模型的分散分析。
1D光子晶体的三种情况的分散分析。

带有光子晶体的弯曲光

如下图所示,当GAAS支柱以周期性的方式排列时,该设备具有角度弯曲的能力(在这种情况下为90°),并且还充当频带的滤波器(也称为A光子带隙)。

光子晶体的示意图。
光子晶体的示意图。

在该中建模该光子晶体comsolMultiphysics®软件和附加波光学模块,横向电(TE)波(在z方向)用于使用左边界传播散射边界条件,幅度为1 v/m。其余的边界分配散射边界条件没有任何事件场。当我们用不同的波长为入射光源运行模型时,我们将获得透射和反射的图表,如下所示。

光子晶体模型的归一化透射率图。
光子晶体的归一化透射率。

左:以1000 nm波长为单位。右:在700 nm波长处的停止条带。

建模光子晶体纤维

踏板纤维可以通过芯的高折射率引导光,而光子晶体纤维(PCF)由微结构光纤组成,这些光纤通过索引引导或带隙限制引导光。在这篇博客文章中,我们专注于索引引导的PCF,其中PCF的核心由覆层材料制成,并被充气孔包围。我们假设气孔的半径为0.3*螺距,其中螺距是相邻孔之间的距离。为了量化色散关系,假设波传播在平面方向上发生,则使用该光子晶体纤维的2D横截面进行模式分析。

索引引导的光子晶体的示意图。
索引引导光子晶体的示意图。

为了绘制色散图(有效折射率与归一化波长),我们执行模式分析,并在孔半径上从0.23 UM到4.69 UM进行参数扫描。为了检测基本和高阶模式,搜索模式的数量增加到50。然后,挑战从总共50个确定的模式中正确识别基本和高阶模式。识别这些模式的一种方法是在核心区域进行集成,以针对不同的有效模式指数(或有效的屈光指标)进行集成。

比较集成在核心区域中的功率和光子晶体纤维的有效模式指数编号的图。
集成在核心区域的功率与4.65-UM孔半径和15.5-UM螺距的有效模式索引编号。

有两种方法可以滤除不必要的模式并仅捕获有意义的基本和高阶模式:

  1. 将过滤器应用于电源;例如,然后是所需的有效折射率ewfd.neff*(intop1(ewfd.poavz)> p_threshold), 在哪里p_threshold是消除不必要模式的力量
  2. 观察基本和高阶模式的有效模式索引编号,以及是否重复

在这种情况下,我们观察到有效模式索引编号40和45之间重复了基本模式,并且在有效模式索引编号20和25之间重复了高阶模式。随附的两个过滤器绘制了分散图,以删除两个滤波器图。不必要的模式。色散图与图4,第9章中的第9章匹配参考。1

在ComsolMultiphysics®中模拟的光子晶体纤维的色散图。
分散图:有效折射率与归一化波长(LDA0/PITCE)。

另一个例子,分析了微结构光纤的泄漏模式。进行模式分析以量化不同模式的复杂有效折射率,例如HE,TE01和TM01样模式。模式分析使用模式搜索方法作为地区并明确提到了折射率的真实和虚构部分的最小和最大值。完美匹配的层(PML)以及散射边界条件(SBC)用于截断模拟域。

以Comsol多物理学模型的TE01样模式下的微结构光纤的图。
微结构光纤的TE01状模式。

分析光子带隙

建模带隙的另一种方法是制定特征值问题,如图所示光子晶体模型的带隙分析。在这种情况下,对GAAS支柱的定期排列进行了建模,其中将支柱彼此等距。与我们在第一个示例中所做的那样,我们没有对GAAS支柱进行建模,而是仅建模一个单位单元格并应用一个Floquet周期性边界条件,如下所示。

光子带隙的阵列到单位近似的示意图。
与单位近似的数组,并在+/- x和+/- y处应用的周期性边界条件。

在波矢量上进行辅助扫描k从0到0.5,以评估(1,1)方向的分散关系。另外,GAAS材料的折射率被认为是频率依赖性的。从分散关系可以看出,在波段3和4之间的(1,1)方向上,没有EM波传播,也称为光子晶体的带隙。

显示光子带隙模拟的分散关系的图。
当波矢量在(1,1)方向上从0到0.5变化时,分散关系。

关于光子晶体的最终想法

无论您执行不同频率的参数扫描,分析模式还是求解特征值的参数扫描,都有不同的方法来建模光子晶体设备。光子晶体可以用作过滤器和工具来设计光线,这在一边特别有用设计光子集成电路。此外,在执行模式分析时,集成核心区域的功率并过滤特定模式索引编号可以帮助从其他不必要的模式中划出基本和高阶模式。最后,一个单位单元模型以及Floquet可以对周期性边界条件进行建模以执行带隙分析。

下一步

了解有关在comsol®软件中建模光子晶体的专门功能的更多信息:

参考

  1. J.D. Joannopoulus,R.D。Meade和J.N.温恩,光子晶体(造型光流),普林斯顿大学出版社,2008年。

评论(15)

发表评论
乐动体育app
加载...
丹红周
丹红周
2019年7月29日

Shweta Mittal
Shweta Mittal
2019年12月17日

我们如何在COMSOL中绘制这种分散和有效的索引曲线?

乌塔姆·帕尔
乌塔姆·帕尔
2019年12月19日 comsol员工

嗨,Shweta,我建议您通过模型//www.dvdachetez.com/model/photonic-crystal-fiber-blog-model-75371,我们在其中创建了具有有效折射率的分散曲线。

Lazyalways
Lazyalways
2020年2月25日

如何计算透射图?

乌塔姆·帕尔
乌塔姆·帕尔
2020年2月27日 comsol员工

嗨,Lazyalways,您可以使用散射变量S21或Prebuild变量EWBE.TPORT_2评估透射图,如您在光学环谐振器模型中所见://www.dvdachetez.com/model/22221

雷切尔·达西(Rachel Darthy R)
雷切尔·达西(Rachel Darthy R)
2021年3月31日

如何计算等频轮廓

乌塔姆·帕尔
乌塔姆·帕尔
2021年4月6日 comsol员工

嗨,雷切尔,
是的,您可以绘制ISO频率轮廓图。
如果您看到Sonic Crystal模型>>研究1,则可以运行KX和KY的“所有组合”扫描。
然后,您执行全局评估>> ABS(ACPR.FREQ),它将为您提供一个表格。
您可以执行此表的“表面”图,并将X轴作为KX,Y轴作为KY和Z轴作为ABS(ACRP.FREQ)。
相同的方法适用于光学应用。对于任何进一步的后续问题,请在support@comsol.com

雷切尔·达西(Rachel Darthy R)
雷切尔·达西(Rachel Darthy R)
2021年5月12日

感谢您的答复。
全局评估后使用表面表面绘图时,错误消息弹出。

乌塔姆·帕尔
乌塔姆·帕尔
2021年5月13日 comsol员工

亲爱的雷切尔,请在support@comsol.com,我们将研究您的模型。

Subrahmanyam Gantasala
Subrahmanyam Gantasala
2021年6月17日

嗨,乌塔姆·帕特尔(Uttam Patel),
我已经尝试了同样的绘制等音轮廓。但是,在选择表表面图的同时弹出的错误是“表中的数据没有定义表面”。你能帮我么?

乌塔姆·帕尔
乌塔姆·帕尔
2021年8月6日 comsol员工

嗨,Subrahmanyam,最好的方法是接触support@comsol.com与您的模型。我们会研究一下。

Sanjay Kumawat
Sanjay Kumawat
2021年6月27日

你好,
我可以计算comsol中的1D语音晶体的带盖吗?如果是,请告诉我如何?我一直在使用周期性的连续性BC和周期性的BC进行尝试,但无法区分P波和S波的带隙,因此获得随机的分散曲线。

乌塔姆·帕尔
乌塔姆·帕尔
2021年8月6日 comsol员工

嗨,Sanjay,您有兴趣通过三种情况进行一维光子晶体模型(//www.dvdachetez.com/model/89741)。

查尔斯·兰博(Charles Rambo)
查尔斯·兰博(Charles Rambo)
2021年10月18日

我遵循为光子晶体波导提供的指南。它不包括生成动画的说明,这是如何完成的?另外,您如何修改该模型以考虑入射光的极化?具体而言,引入了圆极化。

乌塔姆·帕尔
乌塔姆·帕尔
2021年10月29日 comsol员工

嗨查尔斯,
1.当人们使用频域模拟时,然后在导出>>动画中,序列类型应为“动态数据扩展”。因此,动画将取决于相位。您也可以浏览此视频(但这是与时间有关的动画):https://www.comsol.co.in/video/create-export-animations-comsol-multiphysics
2.极化定义为端口条件下的输入,例如(1,0,0)x极化,而(1,-j,0)是圆极化。

探索comsol乐动体育赛事播报博客