24小时热门版块排行榜    

查看: 9601  |  回复: 3

huaer0

新虫 (小有名气)

[交流] [Nanophotonics]手把手教你怎样用Comsol和Matlab提取共振所产生的相位突变已有2人参与

Hello大家好,模数哥最新Comsol仿真交流帖,详情请点击超链接:Nanophotonics-几篇Nature子刊里的法诺共振

[本文原载于微信公众号“Comsol科研助手”,提供Comsol仿真咨询与帮助。敬请查询关注。]

共振是自然界里一种常见的现象,往往表现为某些物理量的增强。如往期例文Nanophotonics- Comsol仿真石墨烯Graphene超材料metamaterial中的完美吸收器功能,石墨烯条带上表面电场的谐振增强发挥了巨大作用。共振也存在于人的思维和意识,比如模数哥关于仿真的作用的观点模数哥谈Comsol仿真系列1——仿真有什么用,就引起了很多朋友的共鸣。

今天模数哥要跟大家分享的,是一种特殊的共振,法诺共振(Fano Resonance),指具有明显的反对称性线性的一种共振形式,其命名是为了纪念法诺物理学家 Ugo Fano 在 1961 年对原子自电离态的反对称共振给出了合理的解释。在光学领域,我们所讨论的法诺共振,一般对应反对称的谱线。这种具有反对称特点的谱线,同样非常sharp,类似于之前讨论的电磁感应透明xxx,所以同样在研究慢光、传感、开关等方面有较大前景。感兴趣的朋友可以参考2010年Nature Materials上的一篇综述文章The Fano resonance in plasmonic nano structures and meta materials(作者为Plasmonic界大牛,美国Rice University的Halas和Nordlander夫妇,他们今年四月下旬在南京参加了第五届等离激元光子学前沿国际会议,小伙伴们有木有去学习呢?)

在电磁感应透明现象中,我们需要一个亮模和一个暗模进行叠加,形成谱线非常尖锐的EIT现象;在法诺共振系统中,我们同样需要一个窄谱(一般为电偶极子谐振)和一个宽谱(如电多极子)的相互作用。比较经典的结构如下图所示巨石阵结构(Nano Lett. 9, 1663–1667 (2009), Fig.2):
[img][Nanophotonics]手把手教你怎样用Comsol和Matlab提取共振所产生的相位突变[/img]
   对于x方向偏振的入射电磁波来说,三个金属纳米棒均产生同向的电偶极子谐振,整个结构形成了一个更大的电偶极子,所以对应一个比较对称而平坦的消光谱(红线);而当电场的偏振变为y方向时,右侧纳米棒虽然能激发出电偶极子,但是左侧两个纳米棒的电偶极子恰好相反,形成了所谓的dark mode,所以他们共同作用形成了消光谱上的谷(蓝线)。a图是实测的结果,b图是FDTD的仿真结果。c图从上到下,展示了入射光偏振方向由y偏振变为x偏振所产生的光谱变化。e, f则分别为780nm入射光在两种偏振下所激发出的电荷分布。

后来有人把纳米颗粒数做成两个的,见2011年Nature Materials的Fano- resonant asymmetric metamaterials for ultrasensitive spectroscopy and identification of molecular monolayers (下面左图)。也有做成4个的,见2013年Nature Nanotechnology,A subwavelength plasmonic metamolecule exhibiting magnetic- based optical Fano resonance (下面右图)。有趣的是,上述的几篇关于法诺共振的工作,作者单位都是在美国德州。
[img][Nanophotonics]手把手教你怎样用Comsol和Matlab提取共振所产生的相位突变-1[/img]
今天模数哥要跟大家具体介绍的,是2013年Scientific Reports上面的一篇理论文章,A dynamically reconfigurable Fano metamaterial through graphene tuning for switching and sensing applications,DOI: 10.1038/ srep02105。这篇工作设计的结构也很简单,如下图所示,即两个相同周期的金属光栅和石墨烯光栅嵌套在一起,是一种hybrid grating with metal and graphene。
[img][Nanophotonics]手把手教你怎样用Comsol和Matlab提取共振所产生的相位突变-2[/img]
前面介绍的几篇法诺共振的文章,都是通过设计特殊的导体结构分别激发电偶极子(窄谱)、电多极子(宽谱),进而叠加产生反对称型的光谱。而在这篇工作中,主要是利用金属和石墨烯两种材料对电磁波的不同响应来产生反对称型(法诺)共振谱线的。下图分别是单纯的金属光栅和石墨烯光栅所产生的透射谱线,与金属光栅相比,石墨烯由于具有非常高的载流子迁移速度,所以其透射谷非常尖锐,带宽很窄,品质因子很高(原文Fig.2 a & b)。
[img][Nanophotonics]手把手教你怎样用Comsol和Matlab提取共振所产生的相位突变-3[/img]
对于既有金属又有石墨烯的hybrid grating,通过优化金属光栅的尺寸以及调节石墨烯的费米能级,我们可以将两个谱线的位置调到一起,宽谱与窄谱的相互作用,便产生了反对称型的法诺谱线,如下面左图。为了分别使金属光栅和石墨烯光栅的谱线谷位置重合,我们再次见证了石墨烯强大的tunable特性(原文Fig.2 c):
[img][Nanophotonics]手把手教你怎样用Comsol和Matlab提取共振所产生的相位突变-4[/img]
对于透射谱线的提取,之前的帖子已经介绍了很多,即通过利用port的S参数。今天主要介绍如何提取上面右图的Transmission Phase,模数哥认为这是此文的一大亮点,即研究了超材料的界面对入射电磁波相位的影响(不同频率的电磁波通过这层光栅结构所产生的相位变化)。
[img][Nanophotonics]手把手教你怎样用Comsol和Matlab提取共振所产生的相位突变-5[/img]
如上图所示,设入射Port1,出射Port2,因为S参数是透射、反射系数,此处S21 = Eout/Ein,则我们可以通过arg(ewfd.S21),得出出射端对入射端的相位变化。我们需要清楚的是,这里所得的相位,包含了三个部分,即入射端材料里的光程,出射端材料里的光程,还有界面所造成的相位变化。所以界面相位Phase_surf = arg(ewfd.S21) + k0*n_in*h_in + k0*n_out*h_out (注意此处为何用+号)。

到这里就结束了吗?抱歉还没有,此时你只能画出上面的右图,是不是离文献原图还差得挺远?此时我们需要做的,是对相位的修正。因为相位信息满足2*pi的周期性,所以上图看似的不连续,其实在相位上算是连续的。我们可以将数据导入到Matlab中,使用unwrap函数,来修正抹平,便可以得到原图了。

几点点评:

1,Matlab中求复数角度的函数是angle,而在Comsol中函数是arg。虽然Comsol的底层语言是Matlab,但是有些函数的语法还是有些不一样。这一点,在我们使用Comsol With Matlab的时候要稍加注意。

2,Matlab里面的unwrap函数,不知道大家用的多不多?尤其是在处理相位信号等方面,还是很强大的,感兴趣的朋友可以打开Matla去help一下。

3,这篇工作的Fig. 4,Fig.5,都是对结构参数、石墨烯性质等参数的扫描。均是由Fig. 2c变化而来,细节不表。比较有趣的是Fig. 3,即画出光栅界面上的电荷分布。这个电荷该怎么求呢?模数哥抛砖引玉,不妨让让大家现自己好好想想,翻翻《电动力学》。有需要的朋友可以后台联系。


声明:

1,本公众号提供Comsol的付费咨询,有需要的朋友请后台咨询细节。

2,本文旨在介绍Comsol强大的仿真能力,引用截图和文献若涉及侵权,请后台联系,我们将进行修改或者删除。
回复此楼

» 收录本帖的淘贴专辑推荐

淘贴集

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

欢迎关注微信公众号:Comsol科研助手,提供仿真咨询。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qli1028

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
非常感谢。

大致上看,unwrap 之后相位应该还是有个负号的差别

是因为comsol用的exp(jwt)?
2楼2018-05-27 22:44:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huaer0

新虫 (小有名气)

引用回帖:
2楼: Originally posted by qli1028 at 2018-05-27 18:44:27
非常感谢。

大致上看,unwrap 之后相位应该还是有个负号的差别

是因为comsol用的exp(jwt)?

嗯的。哈哈,欢迎你关注我的公众号。
欢迎关注微信公众号:Comsol科研助手,提供仿真咨询。
3楼2018-05-28 21:24:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

songguo215

捐助贵宾 (初入文坛)

学习了,谢谢~
4楼2018-10-31 11:35:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huaer0 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金评审时,评审专家还回去看申请人代表作之外的文章吗? +8 孤单12站 2024-06-26 12/600 2024-06-26 21:02 by 孤单12站
[基金申请] 自然基金变化?? +15 今晚推荐22 2024-06-21 25/1250 2024-06-26 19:40 by 朋友别哭010
[基金申请] 基你太美 +13 lyfbangong 2024-06-24 15/750 2024-06-26 17:12 by lyfbangong
[硕博家园] 联培博士文章第一单位署名问题交流~ +9 橙成成c 2024-06-23 21/1050 2024-06-26 16:45 by 橙成成c
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[硕博家园] 博士该不该读? +8 L1009225316 2024-06-25 8/400 2024-06-26 16:26 by huixiong0627
[硕博家园] 怎么带研究生? +24 豆豆小小2 2024-06-20 42/2100 2024-06-26 16:13 by wangyu1986w
[基金申请] 博后面上今天下午会公布吗?大家有无消息? +9 地理学1995 2024-06-24 13/650 2024-06-26 16:02 by kyukitu
[基金申请] 青年基金E02口青基去年几个函评专家? +6 他山攻玉之石 2024-06-25 9/450 2024-06-26 15:09 by 他山攻玉之石
[基金申请] 省基金对国自然申请是有利还是有弊 +6 zycumt04 2024-06-26 7/350 2024-06-26 15:02 by 潇湘之迷
[基金申请] 要持续整治滥发“帽子”、“牌子”之风 +6 babu2015 2024-06-25 6/300 2024-06-26 14:52 by felicity6056
[基金申请] 今天能不能出来名单 +8 地理学1995 2024-06-25 10/500 2024-06-26 09:46 by msjy
[硕博家园] 申博雅思求助 +6 刘嘉研究究 2024-06-24 7/350 2024-06-26 07:08 by 小李醒yy
[基金申请] 基金申请书名称有变化 +5 xuel2011 2024-06-25 7/350 2024-06-26 00:07 by 老虎当猫养
[考博] 申博好难 +6 自强不息a?a 2024-06-24 6/300 2024-06-25 23:02 by 考研吒儿
[基金申请] 博后面上到底什么时候出结果??? +7 爱学lsy 2024-06-24 7/350 2024-06-25 19:24 by 暴走的蒜泥
[第一性原理] Vasp 版权问题 10+4 竹叶青9 2024-06-22 5/250 2024-06-25 14:58 by 无所谓109
[考博] 申请25博士,可以提前进组做科研助理 +6 逐梦途中w 2024-06-22 8/400 2024-06-24 21:49 by keyaner23
[博后之家] 在国内某高校做全职博士后2年,现在找到新的单位,出站或退站对新工作有什么影响? +10 nxplfcc 2024-06-20 10/500 2024-06-22 07:52 by 徐长安
[论文投稿] ACS 编辑的意见 10+3 哈哈妞1993 2024-06-20 3/150 2024-06-21 17:06 by 投个论文
信息提示
请填处理意见