CN114296156A - 基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统 - Google Patents

基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统 Download PDF

Info

Publication number
CN114296156A
CN114296156A CN202111648512.5A CN202111648512A CN114296156A CN 114296156 A CN114296156 A CN 114296156A CN 202111648512 A CN202111648512 A CN 202111648512A CN 114296156 A CN114296156 A CN 114296156A
Authority
CN
China
Prior art keywords
magneto
optical
periodic structure
interface
composite layered
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111648512.5A
Other languages
English (en)
Other versions
CN114296156B (zh
Inventor
曾然
兰阳
黄宇航
徐静俏
李浩珍
杨淑娜
李齐良
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN202111648512.5A priority Critical patent/CN114296156B/zh
Publication of CN114296156A publication Critical patent/CN114296156A/zh
Application granted granted Critical
Publication of CN114296156B publication Critical patent/CN114296156B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

本发明涉及基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统,其中,光学拓扑转换方法,包括以下步骤:(1)建立磁光材料与石墨烯复合层状周期结构模型;(2)分析并计算复合层状周期结构的介电张量;(3)计算并确定复合层状周期结构中磁场和电场表达式以及对应边界条件;(4)计算复合层状周期结构的传输矩阵;(5)由Bloch定理导出无限周期结构的色散关系方程;(6)基于色散关系方程,调节外加磁场得到等频面的变化,实现光学拓扑转换。本发明不局限于现有的金属‑电介质传统双曲材料,拓宽了材料选择的范围,选择了色散双曲特性更加明显的材料结构组成方式,更容易能够实现光学拓扑转换,分析新的拓扑结构特性。

Description

基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方 法及系统
技术领域
本发明属于光电信息领域,具体涉及基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统。
背景技术
随着纳米制造技术和表征技术的发展,光学超材料越来越受到人们的关注。其中,一类高度各向异性的超材料呈现双曲型色散,导致光的传播行为发生了巨大的变化。电磁波在材料中的传播是由介电张量和磁导率张量共同定义,控制介电常数和磁导率便能对任意偏振和非偏振光的色散进行完全的调节。此外,通过对介质介电常数和磁导率的控制,可以针对性设计材料的色散和阻抗,特别是实现双曲材料和自由空间之间的阻抗匹配。阻抗匹配可防止界面处的任何光反射,并允许有效的光耦合。
但近年来,主要聚焦于常见的金属-电介质构成的双曲材料,对于使用其他人造超材料等不同的层状结构以及考虑耗散存在等情况下等频面拓扑结构的转换研究还有待开发。
发明内容
基于现有技术中存在的上述不足,本发明的目的是提供基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统,具体使用一种新型的材料结构即使用电各向异性的石墨烯薄层和磁光材料来构造层状周期材料,通过调节外加磁场强度、各向异性材料层在晶胞中的厚度比等,实现光学拓扑转换,分析其双曲色散特性。
为了达到上述发明目的,本发明采用以下技术方案:
基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法,包括以下步骤:
(1)建立磁光材料与石墨烯复合层状周期结构模型;
(2)分析并计算复合层状周期结构的介电张量;
(3)计算并确定复合层状周期结构中磁场和电场表达式以及对应边界条件;
(4)计算复合层状周期结构的传输矩阵;
(5)由Bloch定理导出无限周期结构的色散关系方程;
(6)基于色散关系方程,调节外加磁场得到等频面的变化,实现光学拓扑转换。
作为优选方案,所述步骤(1)中,建立的磁光材料与石墨烯复合层状周期结构模型为:磁光材料层与石墨烯层交替排列的层状周期结构;
其中,磁光材料层的厚度为d1,磁导率为μ1,介电常数切向分量用εMO表征,εz1为介电常数在z方向分量;石墨烯层的厚度为d2,磁导率为μ2,介电常数切向分量εGr,εz2为介电常数z方向分量,层状周期结构的单位晶胞厚度为h=d1+d2;真空中介电常数为ε0。ωp为等离子体频率,ωp=7.85×1011Hz,z轴为层叠加的方向,磁光材料光轴与z轴重合,即外加磁场方向为z轴。
作为优选方案,所述步骤(2)中,磁光材料的介电张量为:
Figure BDA0003445851580000021
石墨烯材料的介电张量为:
Figure BDA0003445851580000022
其中,磁光材料的介电常数切向分量用εMO表征,εz1为磁光材料的介电常数在z方向分量;石墨烯材料的介电常数切向分量εGr,εz2为石墨烯材料的介电常数z方向分量;
Figure BDA0003445851580000031
ωp为等离子体频率,ωp=7.85×1011Hz;真空中介电常数为ε0
作为优选方案,所述步骤(3)包括:
由于磁光材料和石墨烯材料都为电各向异性材料,只在TM极化下存在双曲特性,TM极化根据边界条件采用磁场
Figure BDA0003445851580000032
以及电场
Figure BDA0003445851580000033
的连续性;区域1为磁光材料层,区域2为石墨烯层;
k1z为区域一中波矢在z方向分量,k2z为区域二中波矢在z方向分量,kx为波矢在x方向分量。
Figure BDA0003445851580000034
Figure BDA0003445851580000035
分别为磁场在区域一左侧x负方向和正方向分量,
Figure BDA0003445851580000036
Figure BDA0003445851580000037
分别为磁场在区域一右侧x负方向和正方向分量。H1y,H2y,H3y分别为区域一、二、三种磁场y方向分量;E1x、E2x为电场在区域一、二中x方向分量;
在区域1、2中磁场表达式为:
Figure BDA0003445851580000038
Figure BDA0003445851580000039
由安培定律得到切向电场表达式为:
Figure BDA00034458515800000310
Figure BDA00034458515800000311
其中,σ为石墨烯的电导率,其表达式为:
Figure BDA0003445851580000041
其中,e为单位电荷,ξ为电子能量,
Figure BDA0003445851580000048
为约化普尔克常量fd(ξ)={exp[(ξ-μc)/(KBT)]+1}-1为费米狄克拉分布,ω为角频率,KB为玻尔兹曼常数,T为温度,μc为化学势,τ为弛豫时间;
其中,1-和1+分别表示两种材料层交界面的左侧和右侧;界面1为区域一磁光材料和区域二石墨烯层交界面;界面2为区域二石墨烯层和区域三磁光材料层交界面;界面3为区域三磁光材料和材料下一个周期的石墨烯层分界面;
由电磁场边界条件可知,在界面1处,电场和磁场的切向分量守恒:
Figure BDA0003445851580000042
Figure BDA0003445851580000043
作为优选方案,所述步骤(4),包括:
Figure BDA0003445851580000044
Figure BDA0003445851580000045
的解为:
Figure BDA0003445851580000046
Figure BDA0003445851580000047
根据区域1、2中磁场关系,将界面2左侧和界面1左侧之间的关系用矩阵M表示;若从界面1左侧到界面2左侧的关系可用传输矩阵中元素M11表示,得到从界面2左侧到界面3左侧的关系用传输矩阵中元素M22表示;
则电磁波传播经过一个由两种材料层组成的周期(从界面1左侧传播到界面3左侧)的传输矩阵为:
Figure BDA0003445851580000051
其中,
Figure BDA0003445851580000052
Figure BDA0003445851580000053
Figure BDA0003445851580000054
Figure BDA0003445851580000055
同理在整个材料结构中有:
Figure BDA0003445851580000056
Figure BDA0003445851580000057
Figure BDA0003445851580000058
分别为磁场在区域三左侧x负方向和正方向分量。
作为优选方案,所述步骤(5),包括:
在双层周期结构中,根据Bloch定理,
Figure BDA0003445851580000059
具有周期性,则M矩阵本征值为
Figure BDA00034458515800000510
其中,kz为有效波矢量在z方向分量,则有:
Figure BDA00034458515800000511
将界面3振幅换算为界面1振幅:
Figure BDA00034458515800000512
得到无限周期结构的色散关系:
Figure BDA00034458515800000513
将矩阵元素带入后得到:
Figure BDA0003445851580000061
其中,在TM模式下,波矢量在材料层中的表示为:
Figure BDA0003445851580000062
Figure BDA0003445851580000063
其中,c为光速。
作为优选方案,所述步骤(6),包括:
以ωp作为归一化频率,令kxc/ωp、kc/ωp为平面坐标,ω/ωp为纵坐标,绘制拓扑转换图,在亚波长尺度内研究复合层状周期结构;其中,ωp=2πc/200h,εMO=14,d1=fh,d2=(1-f)h,
Figure BDA0003445851580000064
ε0为真空中介电常数,τ-1=1×10-3s-1
首先,令ωcp=2不变,比较填充因子f=0.4、0.6时,等频面拓扑转换情况;
当取f=0.4时,在ω/ωp等于0.299和0.3125时出现新的等频面。
当取f=0.6时,在ω/ωp等于0.4002时出现新的等频面,实现光学拓扑转换得到新的拓扑结构。
本发明基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换系统,应用如上任一项所述的光学拓扑转换方法,所述光学拓扑转换系统包括:
模型构建模块,用于建立磁光材料与石墨烯复合层状周期结构模型;
分析计算模块,用于分析并计算复合层状周期结构的介电张量,还用于计算并确定复合层状周期结构中磁场和电场表达式以及对应边界条件,还用于计算复合层状周期结构的传输矩阵;
色散关系方程导出模块,用于根据Bloch定理导出无限周期结构的色散关系方程;
光学拓扑转换模块,用于基于色散关系方程,调节外加磁场得到等频面的变化,实现光学拓扑转换。
本发明与现有技术相比,有益效果是:
(1)本发明不局限于现有的金属-电介质传统双曲材料,拓宽了材料选择的范围,选择了色散双曲特性更加明显的材料结构组成方式,更容易能够实现光学拓扑转换,分析新的拓扑结构特性。
(2)本发明主要基于传输矩阵和Bloch定理进行设计,更加容易计算推导得到准确的色散关系,并且计算过程误差更小。
(3)本发明的材料结构能够通过人为调节外部磁场电场,以及材料占空比等实现不同情况下的等频面变换,使最终拓扑转换结果受材料自身影响的局限性减弱,参数调整方案更加灵活。
附图说明
图1是基于磁光材料与石墨烯构造的新型层状周期结构的光学拓扑转换方法的流程图;
图2是本发明实施例的复合层状周期结构的结构示意图;
图3是本发明实施例的填充因子为0.4的等频面拓扑转换情况示意图;
图4是本发明实施例的填充因子为0.6的等频面拓扑转换情况示意图;
图5、6、7是本发明实施例得到的新的拓扑结构的示意图。
具体实施方式
以下通过具体实施例对本发明的技术方案作进一步解释说明。
本发明实施例提供一种新型的磁光材料与石墨烯复合层状周期材料,该材料结构可以通过调节外加磁场强度、各向异性材料层在晶胞中的厚度比等,分析其双曲色散特性,实现光学拓扑转换。
具体地,如图1所示,基于磁光材料与石墨烯构造的新型层状周期结构的光学拓扑转换方法,包括以下步骤:
第一步:建立复合层状周期结构模型。
建立磁光材料层与石墨烯结构层交替排列的复合层状周期结构,如图2所示。
设磁光材料介电常数切向分量为εMO,层厚度为d1,磁导率为μ1,εz1为介电常数在z方向分量。石墨烯层的材料参数中磁导率为μ2,介电常数切向分量为εGr,层厚度为d2,εz2为介电常数在z方向分量,单位晶胞厚度(即一个周期层厚度)为h=d1+d2。ωp为等离子体频率,ωp=7.85×1011Hz,z为材料层叠加的方向,磁光材料光轴与z轴重合,即外加磁场方向为z轴。由于材料的对称性,只在xz平面内讨论电磁波在其中的传播特性,yz平面的情况与其大致相同。其中用″1-″和″1+″分别表示两种材料层交界面的左侧和右侧。
第二步:分析并计算复合结构的介电张量。
磁光材料的介电张量为:
Figure BDA0003445851580000081
石墨烯材料的介电张量为:
Figure BDA0003445851580000091
第三步:计算并确定复合层状周期结构中磁场和电场表达式以及对应边界条件。
由于磁光材料和石墨烯材料都为电各向异性材料,只在TM极化下存在双曲特性,此处只考虑TM传播的情况。与TE极化不同,TM极化根据边界条件需要采用磁场
Figure BDA0003445851580000092
以及电场
Figure BDA0003445851580000093
的连续性。设k1z为区域一中波矢在z方向分量,k2z为区域二中波矢在z方向分量,kx为波矢在x方向分量。
在区域1,2中磁场表达式:
Figure BDA0003445851580000094
Figure BDA0003445851580000095
由安培定律得到切向电场表达式为:
Figure BDA0003445851580000096
Figure BDA0003445851580000097
其中,σ(ω)为石墨烯的电导率,其表达式为:
Figure BDA0003445851580000098
式中,e为单位电荷,ξ为电子能量,
Figure BDA0003445851580000099
为约化普尔克常量fd(ξ)={exp[(ξ-μc)/(KBT)]+1}-1为费米狄克拉分布,ω为角频率,KB为玻尔兹曼常数,T为温度,μc为化学势,τ为弛豫时间。
其中,1-和1+分别表示两种材料层交界面的左侧和右侧;界面1为区域一磁光材料和区域二石墨烯层交界面;界面2为区域二石墨烯层和区域三磁光材料层交界面;界面3为区域三磁光材料和材料下一个周期的石墨烯层分界面;
由电磁场边界条件可知,在界面1处,电场和磁场的切向分量守恒:
Figure BDA0003445851580000101
Figure BDA0003445851580000102
第四步:计算复合层状周期结构的传输矩阵。
Figure BDA0003445851580000103
Figure BDA0003445851580000104
的解为:
Figure BDA0003445851580000105
Figure BDA0003445851580000106
由上可得,根据区域1,2中磁场关系,可以将界面2左侧和界面1左侧之间的关系用矩阵M表示,若从界面1处左侧到界面2左侧的关系可用传输矩阵中元素M11表示,同样,可以得到从界面2左侧到界面3左侧的关系用传输矩阵中元素M22表示。
则电磁波传播经过一个由两种材料层组成的周期(从界面1左侧传播到界面3左侧)的传输矩阵为:
Figure BDA0003445851580000107
其中:
Figure BDA0003445851580000108
Figure BDA0003445851580000111
Figure BDA0003445851580000112
Figure BDA0003445851580000113
同理,在整个材料结构中有:
Figure BDA0003445851580000114
第五步:由Bloch定理导出无限周期结构的色散关系方程。
在双层周期结构中,根据Bloch定理
Figure BDA0003445851580000115
也具有周期性,则M矩阵本征值为
Figure BDA0003445851580000116
其中kz为有效波矢量在z方向分量,则有:
Figure BDA0003445851580000117
将界面3振幅换算为界面1振幅:
Figure BDA0003445851580000118
综上,得到无限周期结构的色散关系:
Figure BDA0003445851580000119
将矩阵元素带入后得到:
Figure BDA00034458515800001110
其中,在TM模式下,波矢量在材料层中的表示为:
Figure BDA00034458515800001111
Figure BDA0003445851580000121
第六步:基于色散关系方程,调节外加磁场以及材料填充因子得到等频面的变化,实现光学拓扑转换,找到新的拓扑转换。
具体地,以ωp作为归一化频率,令kxc/ωp,kc/ωp为平面坐标,ω/ωp为纵坐标,绘制拓扑转换图,在亚波长尺度内研究复合结构,层厚度取值远远小于入射波长,覆盖的石墨烯很薄,厚度忽略不计。
其中,ωp=2πc/200h,εMO=14,d1=fh,d2=(1-f)h,
Figure BDA0003445851580000122
ε0为真空中介电常数,τ-1=1×10-3s-1,磁光材料作为电各向异性材料,只在TM极化下具有双曲特性,不考虑外加磁场强度变化的影响。过程如下:
(1)首先令ω/ωp=2不变,比较填充因子f=0.4,0.6时,等频面拓扑转换情况如图3、4所示。
(2)当取f=0.6时,在ω/ωp等于0.4002时出现新的等频面,即新的拓扑结构,如图5所示。
(3)当取f=0.4时,在ω/ωp等于0.299和0.3125时出现新的等频面,如图6和7所示。
本发明实施例还提供基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换系统,包括模型构建模块、分析计算模块、色散关系方程导出模块和光学拓扑转换模块;
其中,模型构建模块用于建立磁光材料与石墨烯复合层状周期结构模型;具体地,建立磁光材料层与石墨烯结构层交替排列的复合层状周期结构,如图2所示。
分析计算模块,用于分析并计算复合层状周期结构的介电张量,还用于计算并确定复合层状周期结构中磁场和电场表达式以及对应边界条件,还用于计算复合层状周期结构的传输矩阵。
具体地,设磁光材料介电常数切向分量为εMO,层厚度为d1,磁导率为μ1,εz1为介电常数在z方向分量。石墨烯层的材料参数中磁导率为μ2,介电常数切向分量为εGr,层厚度为d2,εz2为介电常数在z方向分量,单位晶胞厚度(即一个周期层厚度)为h=d1+d2,ωp为等离子体频率,ωp=7.85×1011Hz,z为材料层叠加的方向,磁光材料光轴与z轴重合,即外加磁场方向为z轴。由于材料的对称性,只在xz平面内讨论电磁波在其中的传播特性,yz平面的情况与其大致相同。其中用″1-″和″1+″分别表示两种材料层交界面的左侧和右侧。
(a)分析并计算复合结构的介电张量。
磁光材料的介电张量为:
Figure BDA0003445851580000131
石墨烯材料的介电张量为:
Figure BDA0003445851580000132
(b)计算并确定复合层状周期结构中磁场和电场表达式以及对应边界条件。
由于磁光材料和石墨烯材料都为电各向异性材料,只在TM极化下存在双曲特性,此处只考虑TM传播的情况。与TE极化不同,TM极化根据边界条件需要采用磁场
Figure BDA0003445851580000133
以及电场
Figure BDA0003445851580000134
的连续性。设k1z为区域一中波矢在z方向分量,k2z为区域二中波矢在z方向分量,kx为波矢在x方向分量。
在区域1,2中磁场表达式:
Figure BDA0003445851580000141
Figure BDA0003445851580000142
由安培定律得到切向电场表达式为:
Figure BDA0003445851580000143
Figure BDA0003445851580000144
其中,σ(ω)为石墨烯的电导率,其表达式为:
Figure BDA0003445851580000145
式中,e为单位电荷,ξ为电子能量,
Figure BDA0003445851580000146
为约化普尔克常量fd(ξ)={exp[(ξ-μc)/(KBT)]+1}-1为费米狄克拉分布,ω为角频率,KB为玻尔兹曼常数,T为温度,μc为化学势,τ为弛豫时间。
其中,1-和1+分别表示两种材料层交界面的左侧和右侧;界面1为区域一磁光材料和区域二石墨烯层交界面;界面2为区域二石墨烯层和区域三磁光材料层交界面;界面3为区域三磁光材料和材料下一个周期的石墨烯层分界面。
由电磁场边界条件可知,在界面1处,电场和磁场的切向分量守恒:
Figure BDA0003445851580000147
Figure BDA0003445851580000148
第四步:计算复合层状周期结构的传输矩阵。
Figure BDA0003445851580000149
Figure BDA00034458515800001410
的解为:
Figure BDA00034458515800001411
Figure BDA0003445851580000151
Figure BDA0003445851580000152
由上可得,根据区域1,2中磁场关系,可以将界面2左侧和界面1左侧之间的关系用矩阵M表示,若从界面1处左侧到界面2左侧的关系可用传输矩阵中元素M11表示,同样,可以得到从界面2左侧到界面3左侧的关系用传输矩阵中元素M22表示。
则电磁波传播经过一个由两种材料层组成的周期(从界面1左侧传播到界面3左侧)的传输矩阵为:
Figure BDA0003445851580000153
其中:
Figure BDA0003445851580000154
Figure BDA0003445851580000155
Figure BDA0003445851580000156
Figure BDA0003445851580000157
同理,在整个材料结构中有:
Figure BDA0003445851580000158
色散关系方程导出模块用于根据Bloch定理导出无限周期结构的色散关系方程。具体地,在双层周期结构中,根据Bloch定理
Figure BDA00034458515800001510
也具有周期性,则M矩阵本征值为
Figure BDA0003445851580000159
其中kz为有效波矢量在z方向分量,则有:
Figure BDA0003445851580000161
将界面3振幅换算为界面1振幅:
Figure BDA0003445851580000162
综上,得到无限周期结构的色散关系:
Figure BDA0003445851580000163
将矩阵元素带入后得到:
Figure BDA0003445851580000164
其中,在TM模式下,波矢量在材料层中的表示为:
Figure BDA0003445851580000165
Figure BDA0003445851580000166
其中,c为光速。
光学拓扑转换模块用于基于色散关系方程,调节外加磁场得到等频面的变化,实现光学拓扑转换,找到新的拓扑转换。
具体地,以ωp作为归一化频率,令kxc/ωp,kc/ωp为平面坐标,ω/ωp为纵坐标,绘制拓扑转换图,在亚波长尺度内研究复合结构,层厚度取值远远小于入射波长,覆盖的石墨烯很薄,厚度忽略不计。
其中,ωp=2πc/200h,εMO=14,d1=fh,d2=(1-f)h,
Figure BDA0003445851580000167
ε0为真空中介电常数,τ-1=1×10-3S-1,磁光材料作为电各向异性材料,只在TM极化下具有双曲特性,不考虑外加磁场强度变化的影响。过程如下:
首先令ω/ωp=2不变,比较填充因子f=0.4,0.6时,等频面拓扑转换情况如图3、4所示。
当取f=0.6时,在ω/ωp等于0.4002时出现新的等频面,即新的拓扑结构,如图5所示。
当取f=0.4时,在ω/ωp等于0.299和0.3125时出现新的等频面,如图6和7所示。
以上所述仅是对本发明的两种方式实施例及原理进行了详细说明,对本领域的普通技术人员而言,依据本发明提供的思想,在具体实施方式上会有改变之处,而这些改变也应视为本发明的保护范围。

Claims (8)

1.基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法,其特征在于,包括以下步骤:
(1)建立磁光材料与石墨烯复合层状周期结构模型;
(2)分析并计算复合层状周期结构的介电张量;
(3)计算并确定复合层状周期结构中磁场和电场表达式以及对应边界条件;
(4)计算复合层状周期结构的传输矩阵;
(5)由Bloch定理导出无限周期结构的色散关系方程;
(6)基于色散关系方程,调节外加磁场得到等频面的变化,实现光学拓扑转换。
2.根据权利要求1所述的基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法,其特征在于,所述步骤(1)中,建立的磁光材料与石墨烯复合层状周期结构模型为:磁光材料层与石墨烯层交替排列的层状周期结构;
其中,磁光材料层的厚度为d1,磁导率为μ1;石墨烯层的厚度为d2,磁导率为2,层状周期结构的单位晶胞厚度为h=d1+d2;z轴为层叠加的方向,磁光材料光轴与z轴重合,即外加磁场方向为z轴。
3.根据权利要求2所述的基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法,其特征在于,所述步骤(2)中,磁光材料的介电张量为:
Figure FDA0003445851570000011
石墨烯材料的介电张量为:
Figure FDA0003445851570000021
其中,磁光材料的介电常数切向分量用εMO表征,εZ1为磁光材料的介电常数在z方向分量;石墨烯材料的介电常数切向分量εGr,εz2为石墨烯材料的介电常数z方向分量;
Figure FDA0003445851570000022
ωp为等离子体频率,ωp=7.85×1011Hz;真空中介电常数为ε0
4.根据权利要求3所述的基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法,其特征在于,所述步骤(3)包括:
由于磁光材料和石墨烯材料都为电各向异性材料,只在TM极化下存在双曲特性,TM极化根据边界条件采用磁场
Figure FDA0003445851570000023
以及电场
Figure FDA0003445851570000024
的连续性;区域1为磁光材料层,区域2为石墨烯层;k1z为区域一中波矢在z方向分量,k2z为区域二中波矢在z方向分量,kx为波矢在x方向分量;
Figure FDA0003445851570000025
Figure FDA0003445851570000026
分别为磁场在区域一左侧x负方向和正方向分量,
Figure FDA0003445851570000027
Figure FDA0003445851570000028
分别为磁场在区域一右侧x负方向和正方向分量,在其他区域同理;H1y,H2y,H3y分别为区域一、二、三种磁场y方向分量;E1x、E2x为电场在区域一、二中x方向分量;
在区域1、2中磁场表达式为:
Figure FDA0003445851570000029
Figure FDA00034458515700000210
由安培定律得到切向电场表达式为:
Figure FDA00034458515700000211
Figure FDA00034458515700000212
其中,σ(ω)为石墨烯的电导率,其表达式为:
Figure FDA0003445851570000031
其中,e为单位电荷,ξ为电子能量,
Figure FDA0003445851570000037
为约化普尔克常量fd(ξ)={exp[(ξ-μc)/(KBT)]+1}-1为费米狄克拉分布,ω为角频率,KB为玻尔兹曼常数,T为温度,μc为化学势,τ为弛豫时间;
其中,1-和1+分别表示两种材料层交界面的左侧和右侧;界面1为区域一磁光材料和区域二石墨烯层交界面;界面2为区域二石墨烯层和区域三磁光材料层交界面;界面3为区域三磁光材料和材料下一个周期的石墨烯层分界面;
由电磁场边界条件可知,在界面1处,电场和磁场的切向分量守恒:
Figure FDA0003445851570000032
5.根据权利要求4所述的基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法,其特征在于,所述步骤(4),包括:
Figure FDA0003445851570000033
Figure FDA0003445851570000034
的解为:
Figure FDA0003445851570000035
Figure FDA0003445851570000036
根据区域1、2中磁场关系,将界面2左侧和界面1左侧之间的关系用矩阵M表示;若从界面1左侧到界面2左侧的关系可用传输矩阵中元素M11表示,得到从界面2左侧到界面3左侧的关系用传输矩阵中元素M22表示;
则电磁波传播经过一个由两种材料层组成的周期(从界面1左侧传播到界面3左侧)的传输矩阵为:
Figure FDA0003445851570000041
其中,
Figure FDA0003445851570000042
Figure FDA0003445851570000043
Figure FDA0003445851570000044
Figure FDA0003445851570000045
同理在整个材料结构中有:
Figure FDA0003445851570000046
Figure FDA0003445851570000047
Figure FDA0003445851570000048
分别为磁场在区域三左侧x负方向和正方向分量。
6.根据权利要求5所述的基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法,其特征在于,所述步骤(5),包括:
在双层周期结构中,根据Bloch定理,
Figure FDA0003445851570000049
具有周期性,则M矩阵本征值为
Figure FDA00034458515700000410
其中,kz为有效波矢量在z方向分量,则有:
Figure FDA00034458515700000411
将界面3振幅换算为界面1振幅:
Figure FDA00034458515700000412
得到无限周期结构的色散关系:
Figure FDA0003445851570000051
将矩阵元素带入后得到:
Figure FDA0003445851570000052
其中,在TM模式下,波矢量在材料层中的表示为:
Figure FDA0003445851570000053
Figure FDA0003445851570000054
其中,c为光速。
7.根据权利要求6所述的基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法,其特征在于,所述步骤(6),包括:
以ωp作为归一化频率,令kxc/ωp、kc/ωp为平面坐标,ω/ωp为纵坐标,绘制拓扑转换图,在亚波长尺度内研究复合层状周期结构;其中,ωp=2πc/200h,εMO=14,d1=fh,d2=(1-f)h,
Figure FDA0003445851570000055
ε0为真空中介电常数,τ-1=1×10-3s-1
首先,令ωcp=2不变,比较填充因子f=0.4、0.6时,等频面拓扑转换情况;
当取f=0.4时,在ω/ωp等于0.299和0.3125时出现新的等频面,实现光学拓扑转换得到新的拓扑结构;
当取f=0.6时,在ω/ωp等于0.4002时出现新的等频面,实现光学拓扑转换得到新的拓扑结构。
8.基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换系统,应用如权利要求1-7任一项所述的光学拓扑转换方法,其特征在于,包括:
模型构建模块,用于建立磁光材料与石墨烯复合层状周期结构模型;
分析计算模块,用于分析并计算复合层状周期结构的介电张量,还用于计算并确定复合层状周期结构中磁场和电场表达式以及对应边界条件,还用于计算复合层状周期结构的传输矩阵;
色散关系方程导出模块,用于根据Bloch定理导出无限周期结构的色散关系方程;
光学拓扑转换模块,用于基于色散关系方程,调节外加磁场得到等频面的变化,实现光学拓扑转换。
CN202111648512.5A 2021-12-30 2021-12-30 基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统 Active CN114296156B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111648512.5A CN114296156B (zh) 2021-12-30 2021-12-30 基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111648512.5A CN114296156B (zh) 2021-12-30 2021-12-30 基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统

Publications (2)

Publication Number Publication Date
CN114296156A true CN114296156A (zh) 2022-04-08
CN114296156B CN114296156B (zh) 2023-05-05

Family

ID=80974376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111648512.5A Active CN114296156B (zh) 2021-12-30 2021-12-30 基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统

Country Status (1)

Country Link
CN (1) CN114296156B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102552932A (zh) * 2012-02-09 2012-07-11 哈尔滨工业大学 一种氧化石墨烯双靶向药物载体材料的制备方法和负载的药物
WO2018041179A1 (zh) * 2016-08-31 2018-03-08 深圳大学 无泄漏磁光材料空隙波导磁表面快波光二极管
JP2020057648A (ja) * 2018-09-28 2020-04-09 株式会社カネカ 異方性グラファイト複合体の製造方法
CN111125900A (zh) * 2019-12-19 2020-05-08 杭州电子科技大学 一种计算各向异性铁氧体三层结构Casimir作用力的方法及系统
US20210344168A1 (en) * 2020-04-30 2021-11-04 Zhengzhou University Of Aeronautics Nanolaser based on depth-subwavelength graphene-dielectric hyperbolic dispersive cavity

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102552932A (zh) * 2012-02-09 2012-07-11 哈尔滨工业大学 一种氧化石墨烯双靶向药物载体材料的制备方法和负载的药物
WO2018041179A1 (zh) * 2016-08-31 2018-03-08 深圳大学 无泄漏磁光材料空隙波导磁表面快波光二极管
JP2020057648A (ja) * 2018-09-28 2020-04-09 株式会社カネカ 異方性グラファイト複合体の製造方法
CN111125900A (zh) * 2019-12-19 2020-05-08 杭州电子科技大学 一种计算各向异性铁氧体三层结构Casimir作用力的方法及系统
US20210344168A1 (en) * 2020-04-30 2021-11-04 Zhengzhou University Of Aeronautics Nanolaser based on depth-subwavelength graphene-dielectric hyperbolic dispersive cavity

Also Published As

Publication number Publication date
CN114296156B (zh) 2023-05-05

Similar Documents

Publication Publication Date Title
Hong et al. Electromagnetic interference shielding anisotropy enhanced by CFRP laminated structures
He et al. Finite-element-based generalized impedance boundary condition for modeling plasmonic nanostructures
CN106650093B (zh) 电磁器件分析中电磁-热-应力三场去耦合计算方法
CN110232250B (zh) 多层拓扑绝缘体结构的反射克尔极化偏转的计算方法
Maslovski et al. Equivalent circuit model of radiative heat transfer
Liu et al. Spectral numerical mode matching method for metasurfaces
Al Achkar et al. Effective electromagnetic properties of woven fiber composites for shielding applications
Quan et al. Synergistically enhanced flexibility, mechanical strength and microwave absorption performances of TPE-based hybrid films via thermally assisted homogeneous separation technology
Yi et al. Phase separation-regulated fabrication of MXene/PVA cryogel sensor with effective electromagnetic interference shielding
CN114296156A (zh) 基于磁光材料与石墨烯复合层状周期结构的光学拓扑转换方法及系统
Ma et al. Damage detection of carbon fiber reinforced polymer composite materials based on one-dimensional multi-scale residual convolution neural network
Fan et al. Convolutional neural network method for damage detection of CFRP in electrical impedance tomography
Jia et al. Thin dielectric sheet approximation-based characteristic mode formulations for lossless dielectric coated conducting bodies
CN109948266B (zh) 基于陈绝缘体-手征介质界面的Kerr极化偏转分析方法
Tsang et al. Energy conservation for reflectivity and transmissivity at a very rough surface
CN111950142A (zh) 一种各向异性材料与电介质构造双层周期结构的光学拓扑转换的分析方法
Ren et al. Efficient electromagnetic modeling of multidomain planar layered medium by surface integral equation
CN114297859A (zh) 基于磁光材料与覆盖石墨烯层的电介质复合层状周期结构的光学拓扑转换方法及系统
Miao et al. Superscattering of underwater sound via deep learning approach
Feng et al. Direct-splitting-based CN-FDTD for modeling 2D material nanostructure problems
CN111950143B (zh) 一种双曲超构材料空腔结构的Casimir-Polder效应分析方法及系统
Kishk Electromagnetic waves: propagation in complex matter
Zhang et al. Electromagnetic scattering from three-dimensional bianisotropic objects using hybrid finite element-boundary integral method
Fang et al. Study of electromagnetic interference shielding effectiveness of multilayer graphene films by Monte Carlo method
Dey et al. Method of moment analysis of carbon nanotubes embedded in a lossy dielectric slab using a multilayer dyadic green’s function

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant