CN107870446A - 一种将倏逝波转化为行波的方法 - Google Patents

一种将倏逝波转化为行波的方法 Download PDF

Info

Publication number
CN107870446A
CN107870446A CN201711016184.0A CN201711016184A CN107870446A CN 107870446 A CN107870446 A CN 107870446A CN 201711016184 A CN201711016184 A CN 201711016184A CN 107870446 A CN107870446 A CN 107870446A
Authority
CN
China
Prior art keywords
wave
nano brick
converted
nano
traveling wave
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
CN201711016184.0A
Other languages
English (en)
Other versions
CN107870446B (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.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201711016184.0A priority Critical patent/CN107870446B/zh
Publication of CN107870446A publication Critical patent/CN107870446A/zh
Application granted granted Critical
Publication of CN107870446B publication Critical patent/CN107870446B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/58Optics for apodization or superresolution; Optical synthetic aperture systems
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012Optical design, e.g. procedures, algorithms, optimisation routines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供一种将倏逝波转化为行波的方法,利用偏振独立型纳米砖阵列构成的超表面光栅实现,通过对纳米砖长度和宽度尺寸参数的调整,设计相位梯度并构造超表面光栅,可以对分别沿纳米砖长轴和短轴方向偏振入射的光波实现不同的横向波矢改变量,进而将倏逝波转化为行波传递至远场,实现超分辨成像。这种纳米砖阵列构成的超表面具有高度集成、透射率高、加工工艺相对简单等突出优势,可应用于生物医学显微成像、光学光刻、超高密度光存储等领域。

Description

一种将倏逝波转化为行波的方法
技术领域
本发明属于光学成像领域,具体涉及一种利用偏振独立型纳米砖阵列构成的超表面材料将倏逝波转化为行波的方法。
背景技术
依据傅里叶光学成像理论,来自物体的光波可以被分解为具有不同空间频率的单色平面波,其中横向空间频率代表了物体的细节信息,横向空间频率越大,其相应的横向波矢也就越大。这样根据横向波矢的大小把光波分为可以传播到成像面的行波以及传播距离为波长量级的倏逝波。传统的光学成像系统因为在成像面只收集了行波而受到衍射极限的限制,无法获得物体更小的细节,而在现代科技许多领域,如生物医学显微成像、光学光刻、超高密度光存储等,希望获得超出衍射极限的物体成像,这就要求在成像面不仅要收集行波,同时也要收集倏逝波。目前在近场实现超分辨成像已有相关报道,也有通过逐点扫描的方式实现倏逝波转化为行波的超分辨成像技术(光子隧道显微镜),但是如何实现远场并行超分辨成像仍然是一个难题,其主要难点在于如何在远场收集到倏逝波所携带的物体高频细节信息。因此,有必要设计一种能简单讲倏逝波转化为行波的方法。
发明内容
针对现有技术存在的问题,本发明为解决现有技术中存在的问题采用的技术方案如下:
一种将倏逝波转化为行波的方法,其特征在于:包括如下步骤:
步骤1、根据入射光波的口径确定纳米砖阵列的大小x*x,根据所需观察的图像高频细节设定不同偏振方向入射光波的波矢改变量Δk,由此写出不同偏振方向的相位分布矩阵;
步骤2、由不同偏振方向入射光波的波矢改变量Δk,计算确定纳米砖阵列的单元结构周期;
步骤3、在电磁仿真软件COMSOL中进行模型的构建,分别对纳米砖单元结构进行一维和二维尺寸优化扫描,设置的尺寸扫描范围不超过步骤2中所得到的纳米砖单元结构周期;
步骤3.1、进行一维尺寸扫描时,将纳米砖的长度和宽度均用同一参数表示,根据仿真得到的相位分布与透过率曲线,确定纳米砖的最优高度;
步骤3.2、进行二维尺寸扫描时,将纳米砖的长度和宽度用两个不同的参数表示,扫描范围与一维扫描保持相同,将得到的相位分布以及透过率数据导入到matlab程序中,选取出与相位分布矩阵中的相位调制值最接近的纳米砖,以此构成实现倏逝波转化为行波的超表面材料。
所述步骤2中因一个周期内纳米砖阵列相位变化为2π,故将波矢改变量Δk与纳米砖阵列周期Λ的关系表示为:
Δk×Λ=2π (1)
根据纳米砖阵列的排布方式,进而可以得到单元结构周期cellsize为:
cellstze=Λ/x (2)
所述步骤3.1中进行一维扫描时,以得到的透过率曲线在扫描范围内保持较高水平且相对均匀,同时相位分布曲线在略大于0到2π范围内变化为标准确定纳米砖最优高度。
所述纳米砖阵列基底材料为透明光学玻璃材质,纳米砖材料为高折射率材料。
所述纳米砖阵列基底材料为二氧化硅,纳米砖材料为硅。
本发明提供了一种将倏逝波转化为行波的方法,利用偏振独立型纳米砖阵列构成的超表面光栅实现,通过对纳米砖长度和宽度尺寸参数的调整,设计相位梯度并构造超表面光栅,可以对分别沿纳米砖长轴和短轴方向偏振入射的光波实现不同的横向波矢改变量,进而将倏逝波转化为行波传递至远场,实现超分辨成像,这种纳米砖阵列构成的超表面具有高度集成、透射率高、加工工艺相对简单等突出优势,可应用于生物医学显微成像、光学光刻、超高密度光存储等领域。
与现有技术中基于金属银构造的亚波长光栅实现倏逝波转化为行波的技术相比,本发明将倏逝波转化为行波的方法具有如下优点:
1、本发明涉及的纳米砖采用低损的电介质作为材料,避免了金属对光波的极大损耗,提高了倏逝波转化为行波的效率;
2、相比于亚波长光栅,偏振独立型纳米砖阵列加工工艺较为简单。
附图说明
图1是本发明中实现倏逝波转化为行波的偏振独立型纳米砖阵列几何外形示意图;
图2是本发明实施例中对纳米砖单元结构进行一维扫描后的透过率图:图2(a)以及相位分布图:图2(b);
图3是本发明实施例中以纳米砖单元结构二维扫描结果为基础筛选的16个纳米砖结果图:四列纳米砖分别由正方形、三角形、加号、圆形表示;
图4是本发明中实施例一种4*4的纳米砖阵列仿真几何结构图,图4(a)、(b)分别为XOY方向和YOZ方形仿真结构视图,;
图5是本发明实施例中仿真不同波矢分量的x方向偏振倏逝波入射时的电场分布图;
图6是本发明实施例中仿真y方向偏振倏逝波入射时的电场分布图;
其中:1-完美匹配层A,2-空气层A,3-入射面,4-SiO2基底,5-Si纳米砖阵列,6-空气层B,7-完美匹配层B。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步具体的说明:
如图1所示,为本发明中实现倏逝波转化为行波功能的偏振独立型纳米砖阵列几何外形示意图,通过选用不同尺寸纳米砖以特定方式排布成阵列为入射光波提供不同波矢改变量,从而将空间波矢分量在一定范围内的倏逝波转化为行波出射。
以一种4*4偏振独立型纳米砖阵列构成的超表面材料为例,设定对x方向偏振入射的光波提供2.5k0横向波矢改变量,让y方向偏振入射的光波不受调制直接透过。一个周期内纳米砖阵列相位变化为2π,则波矢改变量Δk与纳米砖阵列周期Λ的关系可以表达为:
Δk×Λ=2π (3)
由于纳米砖阵列按4*4的方式排布,因此单元结构周期为:
cellstze=Λ/4=63.5nm (4)
在得到纳米砖单元结构周期后,对其分别进行一维和二维尺寸优化扫描。一维扫描即在建立几何模型时将纳米砖长、宽用同一个参数表示,其目的是为了确定纳米砖的高度,每经过一次扫描对纳米砖高度做出调整,通过对比相位分布以及透过率曲线来决定单元结构的最优高度。扫描范围为8nm~60nm,在对比不同高度的仿真结果后,确定纳米砖单元结构高度为340nm,其相应的透过率与相位分布如图2所示,其中图2(a)代表透过率分布图像,图2(b)代表相位分布图像。可以看出,在高度为340nm时,透过率在所有扫描尺寸上均保持在较高水平且相对均匀,同时透过率在0~2π之间变化并具有一定的余量,这样可以在一定程度上消除二维扫描中的计算误差,固选定纳米砖单元结构最优高度为340nm。
在得到纳米砖高度的基础上对单元结构进行二维尺寸优化扫描,即分别用两个不同的参数表示纳米砖的长度和宽度,扫描范围依然为4-60nm。将扫描得到的透过率以及相位分布数据导入到matlab程序中进行纳米砖的选取,运行程序后得到与理想相位调制值最为接近的16个纳米砖,具体选取过程为将记载相位与透过率分布的txt文件导入matlab后,先对数据进行简单的处理,使相位值全部在0~2π范围内。根据纳米砖单元结构的周期计算可知,φx(x,y)与φy(x,y)的相位变化梯度均为π/2,选取的纳米砖的相位调制作用Φ应尽可能接近相位分布矩阵中的相位值Φ0(0,90°,180°,270°),从这一点出发,在程序中构造一个偏移量deviation=abs(Φ-Φ0),选取出具有最小偏移量且透过率符合要求的的那一组数据,在这里我们要求选取纳米砖对两种偏振方向光波的透过率均不低于70%,这样一组数据对应的就是我们需要的16个纳米砖。运行结果如图3所示,从图3中可以直观的看出选择的纳米砖相位调制作用与理想相位调制作用的符合程度较好,用来构建纳米砖阵列是可行的。16个纳米砖相关参数如表1所示,其中Lx与Ly分别为纳米砖的长宽尺寸,φx与φy分别为纳米砖对x与y偏振方向入射光波的相位调制大小,Tx与Ty则分别为纳米砖在x与y偏振方向上的透过率。
表1
表1中16个纳米砖在表格中从左到右分别对应图3中每列从下到上的各纳米砖参数。
本实施例中设计的纳米砖阵列预期实现的功能是对x方向偏振入射的光波提供2.5k0波矢改变量,而让y方向偏振入射的光波直接透过,结合单元结构周期为63.5nm,即意味着φx(x,y)的相位变化梯度为π/2,而φy(x,y)没有梯度,相位调制大小全部为0。据此可以写出φx(x,y)和φy(x,y)的相位分布矩阵分别为:
根据相位分布矩阵从16个纳米砖中进行选取构造纳米砖阵列便得到如图4所示的仿真几何结构,通过对波长进行扫描确定透过率最高的波长作为工作波长,得到其中入射光波长为635nm,从空气端入射,入射端口距离纳米砖距离为80nm。
入射端口光波相应的电模式场为:
其中ky为真空中波矢k0在y方向的分量。另一方面,由于dφx(x,y)/dy=π/2,纳米砖阵列提供的波矢改变量Δk会作用在ky上,因此入射光波在空间各方向波矢分量为:
倏逝波在经过纳米砖阵列后,y方向的空间波矢分量变为k′y=ky+Δk,得到出射光波在YOZ面内偏转角θt的正弦为sin(θt)=k′y/k0。在入射端口输入不同波矢分量的倏逝波进行仿真,得到的电场分布图如图5所示,其中图5(a)ky=-3.5k0;图5(b)ky=-3.1k0;图5(c)ky=-2.5k0;图5(d)ky=-1.9k0。计算出四种情况的理论出射角分别为-90°,-37°,0,37°,比较发现仿真结果中的实际出射角度与理论计算值符合的很好。接下来在入射端口输入y方向偏振倏逝波,仿真后的电场分布如图6所示。从图6中可以看出,在出射端并没有观察到明显的光波,说明入射的倏逝波在经过纳米砖阵列时没有受到相位调制作用,在传播过程中就衰减了。综上,本发明成功实现了一种利用偏振独立型纳米砖阵列将倏逝波转化为行波的方法。
本发明的保护范围并不限于上述的实施例,显然,本领域的技术人员可以对本发明进行各种改动和变形而不脱离本发明的范围和精神。倘若这些改动和变形属于本发明权利要求及其等同技术的范围内,则本发明的意图也包含这些改动和变形在内。

Claims (5)

1.一种将倏逝波转化为行波的方法,其特征在于:包括如下步骤:
步骤1、根据入射光波的口径确定纳米砖阵列的大小x*x,根据所需观察的图像高频细节设定不同偏振方向入射光波的波矢改变量Δk,由此写出不同偏振方向的相位分布矩阵;
步骤2、由不同偏振方向入射光波的波矢改变量Δk,计算确定纳米砖阵列的单元结构周期;
步骤3、在电磁仿真软件COMSOL中进行模型的构建,分别对纳米砖单元结构进行一维和二维尺寸优化扫描,设置的尺寸扫描范围不超过步骤2中所得到的纳米砖单元结构周期;
步骤3.1、进行一维尺寸扫描时,将纳米砖的长度和宽度均用同一参数表示,根据仿真得到的相位分布与透过率曲线,确定纳米砖的最优高度;
步骤3.2、进行二维尺寸扫描时,将纳米砖的长度和宽度用两个不同的参数表示,扫描范围与一维扫描保持相同,将得到的相位分布以及透过率数据导入到matlab程序中,选取出与相位分布矩阵中的相位调制值最接近的纳米砖,以此构成实现倏逝波转化为行波的超表面材料。
2.如权利要求1所述的一种将倏逝波转化为行波的方法,其特征在于:所述步骤2中因一个周期内纳米砖阵列相位变化为2π,故将波矢改变量Δk与纳米砖阵列周期Λ的关系表示为:
Δk×Λ=2π (1)
根据纳米砖阵列的排布方式,进而可以得到单元结构周期cellsize为:
cellsize=Λ/x (2)。
3.如权利要求1所述的一种将倏逝波转化为行波的方法,其特征在于:所述步骤3.1中进行一维扫描时,以得到的透过率曲线在扫描范围内保持较高水平且相对均匀,同时相位分布曲线在略大于0到2π范围内变化为标准确定纳米砖最优高度。
4.如权利要求1所述的一种将倏逝波转化为行波的方法,其特征在于:所述纳米砖阵列基底材料为透明光学玻璃材质,纳米砖材料为高折射率材料。
5.如权利要求4所述的一种将倏逝波转化为行波的方法,其特征在于:所述纳米砖阵列基底材料为二氧化硅,纳米砖材料为硅。
CN201711016184.0A 2017-10-25 2017-10-25 一种将倏逝波转化为行波的方法 Expired - Fee Related CN107870446B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711016184.0A CN107870446B (zh) 2017-10-25 2017-10-25 一种将倏逝波转化为行波的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711016184.0A CN107870446B (zh) 2017-10-25 2017-10-25 一种将倏逝波转化为行波的方法

Publications (2)

Publication Number Publication Date
CN107870446A true CN107870446A (zh) 2018-04-03
CN107870446B CN107870446B (zh) 2019-06-11

Family

ID=61752618

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711016184.0A Expired - Fee Related CN107870446B (zh) 2017-10-25 2017-10-25 一种将倏逝波转化为行波的方法

Country Status (1)

Country Link
CN (1) CN107870446B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031660A (zh) * 2018-09-04 2018-12-18 武汉大学 消单色像差大视场双级联变焦透镜及其设计方法
CN109164536A (zh) * 2018-09-06 2019-01-08 武汉大学 基于超表面材料的智能光功率分配器件
CN109239851A (zh) * 2018-08-17 2019-01-18 武汉邮电科学研究院有限公司 一种基于线偏振的光纤耦合器及其制造与使用方法
CN109634086A (zh) * 2019-02-01 2019-04-16 清华大学 基于惠更斯超表面的复振幅调制方法及装置
CN110109198A (zh) * 2019-05-22 2019-08-09 清华大学 一种渐变折射率的超构材料及其纠缠光子对产生系统
CN110609345A (zh) * 2019-09-24 2019-12-24 武汉大学 基于微纳半波片的振幅型光栅的制作方法
CN113433606A (zh) * 2021-06-22 2021-09-24 武汉大学 实现片上波前整形的类金属线结构和非对称传输的应用
CN114325886A (zh) * 2022-02-23 2022-04-12 深圳迈塔兰斯科技有限公司 超表面及其设计方法、装置及电子设备
CN114384612A (zh) * 2022-01-12 2022-04-22 深圳迈塔兰斯科技有限公司 超表面单元、具有其的相位可调超表面和光学系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204422813U (zh) * 2015-02-11 2015-06-24 武汉大学 一种透射式硅纳米阵列光分束器
CN104777545A (zh) * 2015-05-05 2015-07-15 武汉大学 一种硅纳米砖阵列偏振分光器
CN105807598A (zh) * 2016-05-27 2016-07-27 武汉大学 应用于傅里叶彩色全息的硅纳米砖阵列结构及其设计方法
CN105843025A (zh) * 2016-05-27 2016-08-10 武汉大学 应用于傅里叶三维全息的硅纳米砖阵列结构及其设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204422813U (zh) * 2015-02-11 2015-06-24 武汉大学 一种透射式硅纳米阵列光分束器
CN104777545A (zh) * 2015-05-05 2015-07-15 武汉大学 一种硅纳米砖阵列偏振分光器
CN105807598A (zh) * 2016-05-27 2016-07-27 武汉大学 应用于傅里叶彩色全息的硅纳米砖阵列结构及其设计方法
CN105843025A (zh) * 2016-05-27 2016-08-10 武汉大学 应用于傅里叶三维全息的硅纳米砖阵列结构及其设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GUOXING ZHENG ETAL: "Ultracompact high-efficiency polarising beam splitter based on silicon nanobrick arrays", 《OPTICS EXPRESS》 *
ZHONGYI GUO ETAL: "Dielectric metasurface based high-efficiency polarization splitters", 《RSC ADVANCES》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109239851B (zh) * 2018-08-17 2020-07-28 武汉邮电科学研究院有限公司 一种基于线偏振的光纤耦合器及其制造与使用方法
CN109239851A (zh) * 2018-08-17 2019-01-18 武汉邮电科学研究院有限公司 一种基于线偏振的光纤耦合器及其制造与使用方法
CN109031660A (zh) * 2018-09-04 2018-12-18 武汉大学 消单色像差大视场双级联变焦透镜及其设计方法
CN109164536A (zh) * 2018-09-06 2019-01-08 武汉大学 基于超表面材料的智能光功率分配器件
CN109164536B (zh) * 2018-09-06 2020-08-07 武汉大学 基于超表面材料的智能光功率分配器件
CN109634086A (zh) * 2019-02-01 2019-04-16 清华大学 基于惠更斯超表面的复振幅调制方法及装置
CN109634086B (zh) * 2019-02-01 2020-06-02 清华大学 基于惠更斯超表面的复振幅调制方法及装置
CN110109198A (zh) * 2019-05-22 2019-08-09 清华大学 一种渐变折射率的超构材料及其纠缠光子对产生系统
CN110609345A (zh) * 2019-09-24 2019-12-24 武汉大学 基于微纳半波片的振幅型光栅的制作方法
CN110609345B (zh) * 2019-09-24 2020-10-30 武汉大学 基于微纳半波片的振幅型光栅的制作方法
CN113433606A (zh) * 2021-06-22 2021-09-24 武汉大学 实现片上波前整形的类金属线结构和非对称传输的应用
CN113433606B (zh) * 2021-06-22 2022-03-04 武汉大学 实现片上波前整形的类金属线结构和非对称传输的应用
CN114384612A (zh) * 2022-01-12 2022-04-22 深圳迈塔兰斯科技有限公司 超表面单元、具有其的相位可调超表面和光学系统
CN114384612B (zh) * 2022-01-12 2024-02-02 深圳迈塔兰斯科技有限公司 超表面单元、具有其的相位可调超表面和光学系统
CN114325886A (zh) * 2022-02-23 2022-04-12 深圳迈塔兰斯科技有限公司 超表面及其设计方法、装置及电子设备
CN114325886B (zh) * 2022-02-23 2023-08-25 深圳迈塔兰斯科技有限公司 超表面及其设计方法、装置及电子设备

Also Published As

Publication number Publication date
CN107870446B (zh) 2019-06-11

Similar Documents

Publication Publication Date Title
CN107870446B (zh) 一种将倏逝波转化为行波的方法
CN204422813U (zh) 一种透射式硅纳米阵列光分束器
CN109814195A (zh) 基于偏振的多功能超表面结构、超表面元件及加密方法
CN110426758A (zh) 基于超表面的光学薄膜材料、构建方法及其应用
CN108490509A (zh) 低深宽比的电介质几何相位超表面材料及其结构优化方法
Dietrich et al. Elevating optical activity: Efficient on-edge lithography of three-dimensional starfish metamaterial
CN110780366A (zh) 一种抗温度干扰的光学超表面聚焦成像元件及其方法
CN110531458B (zh) 一种可实现非互易性功能的超表面
CN210155355U (zh) 一种对光场振幅与相位均可调的超表面结构
CN105629462A (zh) 一种采用超构表面实现中红外频段隐身的方法
CN112882133B (zh) 一种产生艾里光束的几何相位型超透镜
CN112578490A (zh) 一种3d打印的低折射率大角度偏折稀疏光栅
CN111856746B (zh) 基于叠层超表面的双档振幅光栅的设计方法
CN110568540A (zh) 一种具有双图显示功能的微纳波片阵列及其构建方法
CN112684602A (zh) 用于实现近场自旋角动量复用的超表面材料的设计方法
CN110794662B (zh) 消除零级光的振幅型超表面计算全息片的设计方法
Yang et al. Multiparameter controllable chiral optical patterns
CN112130231A (zh) 产生偏振阶数可调的柱矢量光束的超表面系统及构造方法
CN113655547B (zh) 一种分辨率可调的超透镜阵列及实现方法
CN105467600B (zh) 非线性介质微螺旋器件
Yan et al. Creation of isotropic super-resolved magnetization with steerable orientation
CN110568527B (zh) 一种共振型soi超表面及其在纳米印刷术中的应用
Benedetti et al. Numerical tailoring of linear response from plasmonic nano-resonators grown on a layer of polystyrene spheres
CN106896615A (zh) 非线性螺旋相位器件
CN111722420B (zh) 一种基于超表面的光自旋角动量空间模式转换器

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190611

Termination date: 20211025