CN114815323A - 一种光子晶体调频方法和系统 - Google Patents

一种光子晶体调频方法和系统 Download PDF

Info

Publication number
CN114815323A
CN114815323A CN202210619697.5A CN202210619697A CN114815323A CN 114815323 A CN114815323 A CN 114815323A CN 202210619697 A CN202210619697 A CN 202210619697A CN 114815323 A CN114815323 A CN 114815323A
Authority
CN
China
Prior art keywords
photonic crystal
frequency modulation
crystal frequency
structure device
preprocessing
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
CN202210619697.5A
Other languages
English (en)
Other versions
CN114815323B (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.)
Suzhou Dianguangbo Industrial Intelligent Technology Co ltd
Original Assignee
Suzhou Dianguangbo Industrial Intelligent Technology Co ltd
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 Suzhou Dianguangbo Industrial Intelligent Technology Co ltd filed Critical Suzhou Dianguangbo Industrial Intelligent Technology Co ltd
Priority to CN202210619697.5A priority Critical patent/CN114815323B/zh
Publication of CN114815323A publication Critical patent/CN114815323A/zh
Application granted granted Critical
Publication of CN114815323B publication Critical patent/CN114815323B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
    • G02F1/01—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour 
    • G02F1/0102—Constructional details, not otherwise provided for in this subclass
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B1/00—Optical elements characterised by the material of which they are made; Optical coatings for optical elements
    • G02B1/002—Optical elements characterised by the material of which they are made; Optical coatings for optical elements made of materials engineered to provide properties not available in nature, e.g. metamaterials
    • G02B1/005—Optical elements characterised by the material of which they are made; Optical coatings for optical elements made of materials engineered to provide properties not available in nature, e.g. metamaterials made of photonic crystals or photonic band gap materials
    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics

Landscapes

  • Physics & Mathematics (AREA)
  • Nonlinear Science (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Optical Integrated Circuits (AREA)

Abstract

本申请公开了一种光子晶体调频方法和系统，通过构建光子晶体调频结构装置；基于所述光子晶体调频结构装置进行预处理，获取预处理结果；基于所述预处理结果，通过傅里叶级数展开法进行计算，获取计算结果；基于所述计算结果获得衰减常数随频率的变化曲线,完成光子晶体调频。本申请通过调整金银合金光子晶体中金的含量，即可连续调控金属光子晶体禁带的范围，通过改变金银合金中金比例以达到调频目的的方法比现有方法要简单的多，且成本更低。并减小了光在传播阶段的衰减损耗。

Description

一种光子晶体调频方法和系统
技术领域
本申请属于光子晶体调频领域，具体涉及一种光子晶体调频方法和系统。
背景技术
光子晶体作为光电集成线路发展的关键，一直以来受到人们的广泛关注。光子晶体的一个重要特征是光子禁带，落在禁带频率范围内的电磁波将无法通过光子晶体传播。光子晶体最大的特点是能够控制光波或电磁波的流动，由此可以带来许许多新的应用。如何有规律的控制光子禁带范围是近年来科学家们的研究热点之一。
根据光子晶体的基本理论，改变光子禁带可以通过两种途径实现，一是调整组成光子晶体介质的介电常数ε，二是调整结构的晶格常数，所谓晶格常数就是周期性结构一个基本单元的标度，是一个长度量。当一个器件制造好以后，改变晶格常数即改变光子晶体结构的大小是很困难的方案，因此我们选择的是改变介质的介电常数。
但现有技术缺点是虽然目前计算理论已经成熟，但是液体并不能满足任意介质常数的要求，无法连续调控金属光子晶体禁带的范围；同时也限制了介质的材料，有些场合下需要固体材质且特定介质的液体配制成本过高。而且在金属光子晶体中填充液体介质的方法，即使在更高成本的情况下达到调控禁带的目的，也会造成光在传播阶段出现更大的衰减损耗。
发明内容
本申请提出了一种光子晶体调频方法和系统，通过构建光子晶体调频结构装置并通过傅里叶级数展开法进行计算，获取计算结果；基于所述计算结果获得衰减常数随频率的变化曲线。
为实现上述目的，本申请提供了如下方案：
一种光子晶体调频方法，包括以下步骤：
构建光子晶体调频结构装置；
基于所述光子晶体调频结构装置进行预处理，获取预处理结果；
基于所述预处理结果，通过傅里叶级数展开法进行计算，获取计算结果；
基于所述计算结果获得衰减常数随频率的变化曲线，完成光子晶体调频。
优选的，所述光子晶体调频结构装置包括：圆柱型金银合金介质、空气和PML层。
优选的，所述预处理包括：光从左侧射入所述光子晶体调频结构装置，在PML层之间传播。
优选的，所述傅里叶级数展开法进行计算的方法包括：
依据金属的波长和介电常数的参数及空气的介电常数，通过傅里叶级数展开法计算出各个频率点对应的衰减常数和相位常数，计算参数包括：空气介质常数、属介质的介电常数、计算空间、周期和PML层的厚度。
为了实现上述技术目的，本申请还提出了一种光子晶体调频系统，包括：
光子晶体调频结构装置、光子晶体调频预处理模块、光子晶体调频数据处理模块和光子晶体调频数据处理模块；
所述光子晶体调频结构装置用于构建光子晶体调频结构装置；
所述光子晶体调频预处理模块用于基于所述光子晶体调频结构装置进行预处理，获取预处理结果；
所述光子晶体调频数据处理模块用于基于所述预处理结果，通过傅里叶级数展开法进行计算，获取计算结果；
所述光子晶体调频图像绘制模块用于基于所述计算结果获得衰减常数随频率的变化曲线，完成光子晶体调频。
优选的，所述光子晶体调频结构装置中构建所述光子晶体调频结构装置包括：圆柱型金银合金介质、空气和PML层。
优选的，所述光子晶体调频预处理模块中所述预处理包括：光从左侧射入所述光子晶体调频结构装置，在PML层之间传播。
优选的，所述光子晶体调频数据处理模块中所述傅里叶级数展开法进行计算的方法包括：
依据金属的波长和介电常数的参数及空气的介电常数，通过傅里叶级数展开法计算出各个频率点对应的衰减常数和相位常数，计算参数包括：空气介质常数、属介质的介电常数、计算空间、周期和PML层的厚度。
本申请的有益效果为：本申请公开了一种光子晶体调频方法和系统，本申请通过调整金银合金光子晶体中金的含量，即可连续调控金属光子晶体禁带的范围，通过改变金银合金中金比例以达到调频目的的方法比现有方法要简单的多，且成本更低。并减小了光在传播阶段的衰减损耗。
附图说明
为了更清楚地说明本申请的技术方案，下面对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例方法流程示意图；
图2为本申请实施例二维金属光子晶体的结构图示意图；
图3为本申请实施例金属光子晶体的结构图示意图；
图4为本申请实施例晶格常数示意图；
图5为本申请实施例矩形棒离散圆棒示意图；
图6为本申请实施例不同比例金属纳米棒阵列的相位常数随频率的变化曲线示意图；
图7为本申请实施例不同比例金属纳米棒阵列的衰减常数随频率的变化曲线示意图；
图8为本申请实施例波长为359nm银链式结构的近场传播分布示意图；
图9为本申请实施例波长为359nm银链式结构的沿y轴的频率场的分布示意图；
图10为本申请实施例波长为383nm Au:Ag(1:4)链式结构的近场传播分布示意图；
图11为本申请实施例波长为383nm Au:Ag(1:4)链式结构的沿y轴的频率场的分布示意图；
图12为本申请实施例波长为530nm Au:Ag(1:2)链式结构的近场传播分布示意图；
图13为本申请实施例波长为530nm Au:Ag(1:2)链式结构的沿y轴的频率场的分布示意图；
图14为本申请实施例波长为288nm纯金链式结构的近场传播分布示意图；
图15为本申请实施例波长为288nm纯金链式结构的沿y轴的频率场的分布示意图；
图16为本申请实施例系统组成示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
为使本申请的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本申请作进一步详细的说明。
实施例1
如图1所示，一种光子晶体调频方法，具体包括以下内容：
构建光子晶体调频结构装置；
基于光子晶体调频结构装置进行预处理，获取预处理结果；
基于预处理结果，通过傅里叶级数展开法进行计算，获取计算结果；
基于计算结果获得衰减常数随频率的变化曲线。
其中，构建光子晶体调频结构装置包括：圆柱型金银合金介质、空气和PML层。
其中，预处理包括：光从左侧射入，在PML层之间传播。
其中，傅里叶级数展开法进行计算的方法包括：
依据金属的波长和介电常数的参数及空气的介电常数，计算出各个频率点对应的衰减常数和相位常数，计算参数包括：基于圆柱的底面半径r,空气介质常数，金属介质的介电常数ε，计算空间Λ，周期h，PML层的厚度。
然后，通过使用傅里叶级数展开法，计算出各个频率点对应的衰减常数和相位常数，然后绘制成衰减常数和相位常数随频率的变化曲线。
最后，得到最终的输出结果，实现最终效果，减小光在传播阶段的衰减损耗，使光的传输效率更高。
在本实施例中，具体的操作步骤如下所示：
构建光子晶体调频结构装置具体包括：圆柱型金银合金介质、空气和PML层。
如图2所示，1表示圆柱型金银合金介质，其底面半径r＝55nm，圆柱个数为无数个，圆柱在垂直方向为无限长；
2表示周期h＝55nm；
3表示介质为空气；
4表示为PML层；
5表示PML层的厚度为55nm；
6表示虚拟周期A为3.3ume；
7表示光从左侧垂直入射。
如图3所示，光从左侧射入，在两层PML层之间传播；
然后依据金属的波长和介电常数的参数及空气的介电常数，使用傅里叶级数展开法，计算出各个频率点对应的衰减常数和相位常数，然后绘制成衰减常数和相位常数随频率的变化曲线。
如图4所示，图中的点，横坐标是频率，纵坐标分别是是衰减常数和相位常数，
在z方向上，散射单元无限长，相互平行分布。假设底面圆半径为r的圆柱是纯电介质，其介电常数为ε。具有这种结构的数组沿y方向重复相同的配置，假设周期为Λ，然后用该数组位于0≤y≤Λ的单元格来近似初始结构。假设我们引导了E(Ez，Hx，Hy)波和H(Hz，Ex，Ey)波的传播。表面等离子体是发生在金属-介质界面上，仅由H波(垂直于纳米线截面的磁场分量)激发的大量自由电子集体振动。首先，我们进行描述提法的核心内容。对于H波，麦克斯韦方程为:
Figure BDA0003675823740000081
Figure BDA0003675823740000082
Figure BDA0003675823740000083
Dy＝ε0ε(y)Ey (4)
其中H′z＝(μ0/ε0)1/2Hz，G(y)＝[1+iσ(y)]-1是拉伸坐标变量，代表假设的完美匹配层(PML)，其中σ(y)＝σmax(1-y/ω)d是电导率函数。在传播场的虚拟周期中，把电场和磁场近似为看作截断的傅立叶级数：
Figure BDA0003675823740000084
Figure BDA0003675823740000085
Figure BDA0003675823740000086
Figure BDA0003675823740000087
式中Dy是电位移向量，
Figure BDA0003675823740000088
将式(2)代入式(1)，根据傅里叶基底的正交性质进行计算。有些H模有场分量(Hz，Ex，Ey)的特点，我们将Hz定为主导场。垂直于芯包层边界的电场存在不连续的问题，为了解决这个问题，这里的求解我们参考李的分解规则，我们将Ey场进行傅立叶级数展开。推导出确定傅里叶系数{hz，m(x)}和{ey，m(x)}的一组线性方程为：
Figure BDA0003675823740000091
Figure BDA0003675823740000092
hz(x)＝[hz，-M...hz，0...hz，M]T (11)
ey(x)＝[ey，-M...ey，0...ey，M]T (12)
Figure BDA0003675823740000093
Figure BDA0003675823740000094
Figure BDA0003675823740000095
Figure BDA0003675823740000096
Figure BDA0003675823740000097
式中δmm′为单位矩阵；上标T是向量的转置符号；k0是自由空间中的波数；ε(y)为0≤y≤Λ周期内沿y轴的介电常数。矩阵G的特征值kn＝ξn 2(n＝1，2，3，...，2M+1)确定了传播常数ξn和假设波导中引导模式。特征向量Pn，确定了辐射模式的场分布。式(8)(9)的解为
Figure BDA0003675823740000101
Figure BDA0003675823740000102
Figure BDA0003675823740000103
P＝[p1 p2...p2M+1] (21)
u(±)(x)＝[exp(±ik0ξnx)δnn′] (22)
B＝[ξn δmm′] (23)
a(x)＝[a(+)(x)a(-)(x)] (24)
Figure BDA0003675823740000104
其中an (±)(x)表示前向波和后向波的n阶模的振幅。u(x)，P，和a(±)(x)的元素按引导和辐射模式的顺序重新排列，Re[ξn]递减。
紧接着，如图5所示，将每个圆柱细分为一定数量的平行矩形柱，然后把x方向周期链的散射单元替换为分层平行平面波导的级联连接。计算每段波导，方程(1)-(8)的解均有式(18)-(25)进行求解。我们将截面两侧的傅里叶系数相等，使Hz和Ey符合在每个台阶不连续处的边界条件。这将使散射矩阵Sj在界面x＝xj处有如下关系：
Figure BDA0003675823740000111
Figure BDA0003675823740000112
散射矩阵Sj+1在x＝xj+1处的定义方法与式(25)相同，另外我们还应当考虑在距离xj+1-xj处的模态传播，如下：
Figure BDA0003675823740000113
和
Figure BDA0003675823740000114
Figure BDA0003675823740000115
Figure BDA0003675823740000116
Figure BDA0003675823740000117
Figure BDA0003675823740000118
Figure BDA0003675823740000119
其中工为单位矩阵；u(+)(x)由式(22)定义。方程(29)-(31)表示波导在每一个不连续处的散射矩阵的递归关系。当光在波导中传播过程中不断变化时，过渡段近似为大量的阶跃不连续。假如在x＝xj(j＝1，2，...，N)处沿波导存在，(N-1)倍递归过程会使整个系统的广义散射矩阵SN[33]为:
Figure BDA0003675823740000121
我们在每一个周期h内，按照以上关系计算x方向上的广义散射矩阵，其中N代表着单元格内细分小矩形的数量。根据计算出的结果，得到符合以下关系[33]的转移矩阵K:
Figure BDA0003675823740000122
所以，第k模的传播常数γk＝βk+iαk由下式确定：
γk＝-ilogχk/h (37)
其中χk表示的是转移矩阵KN。根据Floquet定理，第k阶导模的传播常数可以由以下关系确定:
χk＝exp(iβkh) (38)
如图6、图7所示，我们改变金银合金介质中金银比例的大小，发现随着金含量的增加，衰减和截止频率范围逐渐减小，即禁带范围逐渐减小，拥有连续变化的特征。随着增加合金中金的比例，相位常数和衰减常数随频率变化的曲线图变得越来越平滑，且含金量越多的混合型链式结构的衰减系数越小。这意味着在传播频率范围内，能量的损耗越小，传输功率越高。
其中，h/λ：频率、βh/2π:相位常数、αh/2π:衰减常数、Ag纯银、Au纯金、Au:Ag＝1:4金银比例为1比4、Au:Ag＝1:2金银比例为1比2、Au:Ag2:1金银比例为2比1。
A点-D点频率分别为0.1532，0.1437，0.1037和0.1909，分别是图8、图10、图12和图14的频率。
如表1所示，表1为不同比例金属链式结构的传播特性表，之后我们测试了由以上结构中40个圆柱组成的有限链式结构，得到不同比例金属的链式结构传播特性表并且举出4个频率点的近场传播分布和沿y轴的电场的分布情况。
表1
Figure BDA0003675823740000131
如图8-图15所示，通过设置圆柱的底面半径r,空气介质常数，金属介质的介电常数ε，计算空间Λ，周期h，PML层的厚度参数，通过傅里叶级数展开公式，即可计算出每个h/λ对应的βh/2π和αh/2π，然后我们通过计算各个频率点近场传播分布，通过三维图结果判断波的传播状态，划分出传播状态频率范围，衰减和截止状态频率范围。
实施例2
如图16所示，为了更好的实现上述技术内容，本实施例还提出了一种光子晶体调频系统，包括：光子晶体调频结构装置、光子晶体调频预处理模块、光子晶体调频数据处理模块和光子晶体调频数据处理模块；
光子晶体调频结构装置用于构建光子晶体调频结构装置；
光子晶体调频预处理模块用于基于光子晶体调频结构装置进行预处理，获取预处理结果；
光子晶体调频数据处理模块用于基于预处理结果，通过傅里叶级数展开法进行计算，获取计算结果；
光子晶体调频图像绘制模块用于基于计算结果获得衰减常数随频率的变化曲线。
在本实施例中，光子晶体调频结构装置中构建光子晶体调频结构装置包括：圆柱型金银合金介质、空气和PML层。
在本实施例中，光子晶体调频预处理模块中预处理包括：光从左侧射入光子晶体调频结构装置，在PML层之间传播。
在本实施例中，光子晶体调频数据处理模块中傅里叶级数展开法进行计算的方法包括：
依据金属的波长和介电常数的参数及空气的介电常数，计算出各个频率点对应的衰减常数和相位常数，计算参数包括：空气介质常数、属介质的介电常数、计算空间、周期和PML层的厚度。
以上所述的实施例仅是对本申请优选方式进行的描述，并非对本申请的范围进行限定，在不脱离本申请设计精神的前提下，本领域普通技术人员对本申请的技术方案做出的各种变形和改进，均应落入本申请权利要求书确定的保护范围内。

Claims (8)

1.一种光子晶体调频方法，其特征在于，
构建光子晶体调频结构装置；
对所述光子晶体调频结构装置进行预处理，获取预处理结果；
基于所述预处理结果，通过傅里叶级数展开法进行计算，获取计算结果；
基于所述计算结果获得衰减常数随频率的变化曲线，完成光子晶体调频。
2.如权利要求1所述一种光子晶体调频方法，其特征在于，
所述光子晶体调频结构装置包括：圆柱型金银合金介质、空气和PML层。
3.如权利要求1所述一种光子晶体调频方法，其特征在于，
所述预处理包括：光从左侧射入所述光子晶体调频结构装置，在PML层之间传播。
4.如权利要求1所述一种光子晶体调频方法，其特征在于，
所述傅里叶级数展开法进行计算的方法包括：
依据金属的波长和介电常数的参数及空气的介电常数，通过傅里叶级数展开法计算出各个频率点对应的衰减常数和相位常数，计算参数包括：空气介质常数、属介质的介电常数、计算空间、周期和PML层的厚度。
5.一种光子晶体调频系统，其特征在于，
包括：光子晶体调频结构装置、光子晶体调频预处理模块、光子晶体调频数据处理模块和光子晶体调频数据处理模块；
所述光子晶体调频结构装置用于构建光子晶体调频结构装置；
所述光子晶体调频预处理模块用于基于所述光子晶体调频结构装置进行预处理，获取预处理结果；
所述光子晶体调频数据处理模块用于基于所述预处理结果，通过傅里叶级数展开法进行计算，获取计算结果；
所述光子晶体调频图像绘制模块用于基于所述计算结果获得衰减常数随频率的变化曲线，完成光子晶体调频。
6.如权利要求5所述一种光子晶体调频系统，其特征在于，
所述光子晶体调频结构装置中构建所述光子晶体调频结构装置包括：圆柱型金银合金介质、空气和PML层。
7.如权利要求5所述一种光子晶体调频系统，其特征在于，
所述光子晶体调频预处理模块中所述预处理包括：光从左侧射入所述光子晶体调频结构装置，在PML层之间传播。
8.如权利要求5所述一种光子晶体调频系统，其特征在于，
所述光子晶体调频数据处理模块中所述傅里叶级数展开法进行计算的方法包括：
依据金属的波长和介电常数的参数及空气的介电常数，通过傅里叶级数展开法计算出各个频率点对应的衰减常数和相位常数，计算参数包括：空气介质常数、属介质的介电常数、计算空间、周期和PML层的厚度。
CN202210619697.5A 2022-06-02 2022-06-02 一种光子晶体调频方法和系统 Active CN114815323B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210619697.5A CN114815323B (zh) 2022-06-02 2022-06-02 一种光子晶体调频方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210619697.5A CN114815323B (zh) 2022-06-02 2022-06-02 一种光子晶体调频方法和系统

Publications (2)

Publication Number Publication Date
CN114815323A true CN114815323A (zh) 2022-07-29
CN114815323B CN114815323B (zh) 2023-04-14

Family

ID=82519986

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210619697.5A Active CN114815323B (zh) 2022-06-02 2022-06-02 一种光子晶体调频方法和系统

Country Status (1)

Country Link
CN (1) CN114815323B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1996100A (zh) * 2006-12-30 2007-07-11 中国科学院光电技术研究所 温度调控一维光子晶体产生动态光子晶体的方法
CN101408623A (zh) * 2008-01-23 2009-04-15 北京航空航天大学 宽带综合孔径上变频成像系统
CN102393571A (zh) * 2011-11-09 2012-03-28 南开大学 高速光子晶体波导太赫兹调制器
CN103235361A (zh) * 2013-04-27 2013-08-07 江苏大学 一种二维等离子体光子晶体带隙控制方法
WO2018126841A1 (zh) * 2017-01-08 2018-07-12 复旦大学 一种焦距大幅可调的光子晶体透镜及其设计方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1996100A (zh) * 2006-12-30 2007-07-11 中国科学院光电技术研究所 温度调控一维光子晶体产生动态光子晶体的方法
CN101408623A (zh) * 2008-01-23 2009-04-15 北京航空航天大学 宽带综合孔径上变频成像系统
CN102393571A (zh) * 2011-11-09 2012-03-28 南开大学 高速光子晶体波导太赫兹调制器
CN103235361A (zh) * 2013-04-27 2013-08-07 江苏大学 一种二维等离子体光子晶体带隙控制方法
WO2018126841A1 (zh) * 2017-01-08 2018-07-12 复旦大学 一种焦距大幅可调的光子晶体透镜及其设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DAN ZHANG等: "Numerical Analysis of Leaky Modes in Two-Dimensional Photonic Crystal Waveguides Using Fourier Series Expansion Method with Perfectly Matched Layer", 《IEICE TRANS. ELECTRON.》 *
ISRAEL DE LEON等: "Fourier analysis of reflection and refraction in two-dimensional photonic crystals", 《PHYSICAL REVIEW B》 *

Also Published As

Publication number Publication date
CN114815323B (zh) 2023-04-14

Similar Documents

Publication Publication Date Title
Ziolkowski Propagation in and scattering from a matched metamaterial having a zero index of refraction
Yasumoto et al. Modeling of photonic crystals by multilayered periodic arrays of circular cylinders
Sigalas et al. Waveguides in three-dimensional metallic photonic band-gap materials
Matekovits et al. Controlling the bandlimits of TE-surface wave propagation along a modulated microstrip-line-based high impedance surface
Ishizaka et al. Three-dimensional finite-element solutions for crossing slot-waveguides with finite core-height
Borzdov Frequency domain wave-splitting techniques for plane stratified bianisotropic media
Zhu et al. Wideband absorbing plasmonic structures via profile optimization based on genetic algorithm
Silveirinha et al. Nonresonant structured material with extreme effective parameters
CN114815323B (zh) 一种光子晶体调频方法和系统
Ahmed et al. A time-domain adjoint variable method for materials with dispersive constitutive parameters
Zong et al. Mixed propagation modes in three Bragg propagation periods of variable chain structures
Baibak et al. Surface electromagnetic waves in finite semiconductor-dielectric periodic structure in an external magnetic field
Zhang et al. Numerical analysis of leaky modes in two-dimensional photonic crystal waveguides using Fourier series expansion method with perfectly matched layer
Tayeb et al. Compact directive antennas using metamaterials
Rodríguez‐Esquerre et al. Frequency‐dependent envelope finite‐element time‐domain analysis of dispersion materials
Zhang et al. PSO-algorithm-assisted design of compact SSPP transmission line
Silveirinha et al. Efficient calculation of the band structure of artificial materials with cylindrical metallic inclusions
Choubani et al. Analysis and design of electromagnetic band gap structures with stratified and inhomogeneous media
Váry et al. Propagation of surface plasmon polaritons through gradient index and periodic structures
Thevenot et al. FDTD to analyze complex PBG structures in the reciprocal space
Yasumoto et al. Analysis of post-wall waveguides using a model of two-dimensional photonic crystal waveguides
Jamali et al. Implementation of Periodic Boundary Conditions in FDTD Algorithm
CN117559141B (zh) 一种基于拓扑优化的透明宽带超材料吸波体及其制备方法
Fallahi et al. Diffraction analysis of periodic structures using a transmission-line formulation: principles and applications
Archemashvili et al. Numerical Analysis of Dielectric Post-Wall Waveguides

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