CN104111235A - 一种测量二维薄膜材料复折射率谱的方法 - Google Patents
一种测量二维薄膜材料复折射率谱的方法 Download PDFInfo
- Publication number
- CN104111235A CN104111235A CN201410331238.2A CN201410331238A CN104111235A CN 104111235 A CN104111235 A CN 104111235A CN 201410331238 A CN201410331238 A CN 201410331238A CN 104111235 A CN104111235 A CN 104111235A
- Authority
- CN
- China
- Prior art keywords
- film material
- dimensional film
- phi
- spectrum
- contrast
- 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
Links
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种测量二维薄膜材料复折射率谱的方法,包括:测量所述样品的反射光谱,包括覆盖有二维材料区域的反射光谱和没有覆盖二维材料区域的反射光谱,利用二位薄膜材料在衬底上的衬度公式,得到二维薄膜材料的衬度谱;改变透明介质厚度,得到不同透明介质厚度下的衬度谱;在固定照明波长下通过多层薄膜反射光强的理论公式拟合衬度与透明介质厚度的关系获得二维薄膜材料的复折射率参数;改变照明波长,最终获得二维薄膜材料的复折射率谱。该方法具有空间分辨率高、测量范围宽、重复性好、成本低等优点。避免了传统椭偏仪和filmetrics测量方法不适用于新型小尺度二维材料的缺点,在以二维薄膜材料为功能单元的新型器件设计与制备方面有广泛的应用前景。
Description
技术领域
本发明提供了一种直接测量二维薄膜材料折射率谱的方法,该方法可以应用于小尺度材料。涉及新型二维薄膜材料物理及相关器件研究领域。
背景技术
二维薄膜材料通常指单层或者少层(通常少于10层原子层)原子构成的薄膜材料,例如:单层硫化钼(MoS2)等。由于其新奇的物理性质,二维薄膜材料已经成为物理、化学、生物以及材料科学领域的一个研究热点。折射率是描述电磁波在介质中如何传播的一个重要的物理参量。与MoS2体材料不同,单层MoS2是直接禁带半导体。这种晶体结构和能带结构的变化使得单层硫化钼的折射率与其体材料有很大不同。迄今为止,采用微机械剥离法或化学气相沉积(CVD)法制备的单层MoS2尺寸小(通常在几-几十微米量级),不能用传统的方法(例如:椭偏仪,或Filmetrics等)对其折射率进行测量。
发明内容
本发明的目的是提供一种测量二维薄膜材料复折射率谱的方法,该方法可以应用于小尺度材料。
其中所用的二维薄膜材料尺度只要大于2μm×2μm即可。
本发明提供的技术方案如下:
一种测量二维薄膜材料复折射率谱的方法,其特征是,包括如下步骤:
1)在衬底上生长二维薄膜材料或者转移二维薄膜材料到衬底上,保证样品干净;
2)将带有二维薄膜材料的衬底置于光学显微镜(如图1)之下,采用白光光源(如卤素灯)反射照明方式,结合光纤耦合输出对目标区域进行光谱测量,分别获得样品上有二维薄膜材料和没有二维薄膜材料处的反射光谱;
3)利用二位薄膜材料在衬底上的衬度公式对获得的反射光谱进行处理,得到二维薄膜材料衬度谱,衬度公式如下:
式中,I1(λ)和I0(λ)分别表示有二维薄膜材料处和没有二维薄膜材料处的反射光强度;
4)用具有不同透明介质厚度的衬底重复步骤1)至3)多次(重复的次数越多则测量结果精度越高),得到不同入射波长下衬度与透明介质厚度的关系曲线;
5)在固定入射波长下,采用多层薄膜干涉的理论模型对所述的样品结构(例如:空气/单层MoS2/SiO2/Si)进行计算,得到以二维薄膜材料折射率实部、虚部和透明介质层厚度为自变量的衬度的表达式;
6)用该表达式对测得的衬度与介质厚度的关系曲线进行拟合,得到该波长下二维薄膜材料的复折射率;改变入射波长,得到整个测量波段下该二维薄膜材料的复折射率谱。
优选的:
所述的测量二维薄膜材料复折射率谱的方法,其特征是,所用衬底为具有透明介质、非透明介质双层结构的固态衬底。
所述的测量二维薄膜材料复折射率谱的方法,其特征是,步骤2)中,通过光纤直径来控制所测区域的大小可以。
所述的测量二维薄膜材料复折射率谱的方法,其特征是,步骤4)中,重复执行步骤1)至3)至少20次。
所述的测量二维薄膜材料复折射率谱的方法,其特征是,改变衬底透明介质的厚度来获得多种不同厚度衬底所对应的衬度谱数据,利用多层薄膜反射光强的理论公式对该数据进行拟合,得到最终的复折射率谱,公式如下:
其中
是相应材料的复折射率,实数部分ni为折射率;虚数部分的ki决定于光波在吸收性介质中传播时的衰减,叫做吸收系数;di是相应材料的厚度(见图2中所示),λ是入射光波长。
所述的测量二维薄膜材料复折射率谱的方法,其特征是,所用的二维薄膜材料通过微机械剥离或CVD方法制备在SiO2/Si衬底上。
与现有技术相比,本发明的有益效果是:
本发明所述的方法具有空间分辨率高、测量范围宽、重复性好、成本低等优点。该方法相比传统的椭偏仪测量方法具有一定优势。首先椭偏仪测量一般需要材料具有较大的面积,尺度至少约100μm。而很多二维薄膜材料目前较难获得这么大面积的单晶,例如单层MoS2,采用微机械剥离法或CVD方法制备的单层单晶MoS2尺寸通常在几-几十微米量级。除此之外,从椭偏仪的直接测量结果要得到折射率等数据,一般需要选用合适的色散关系模型,而对于新型二维材料,相关模型可靠性有待验证(现在的商用椭偏仪折射率测量一般需要样品厚度大于5nm才能获得比较准确的结果)。本发明所述方法空间分辨率可达2μm(空间分辨率由物镜放大倍数和在像平面上采集光谱的光纤直径共同决定。物镜的放大倍数越大,分辨率越高,采集光纤的直径越细,分辨率越高),完全可以用来测量新型二维单晶薄膜材料。同时不需要借助色散关系等模型,结果可靠性高。并且本发明所述方法对于仪器要求低,不需要高精度的对准或移动器件,使得该方法比椭偏仪更加简单,重复性更高。对于Filmetrics方法,其测量的样品只能是具有宏观尺度的体材料,而无法对纳米材料进行微区测量。并且由于其无法利用入射光的偏振信息,很难准确测量出非常薄的材料(如二维材料)的折射率信息(现在的商用Filmetrics仪器通常需要材料厚度大于20nm时才能获得准确的折射率信息)。在以二维薄膜材料为功能单元的新型器件设计与制备方面有广泛的应用前景。
附图说明
图1.本发明(一种测量二维薄膜材料复折射率的新方法)的测量光路图。其中,1-光谱仪;2-采集光纤;3-像平面;4-反射镜;5-半透半反镜;6-物镜;7-样品;8-白光光源。
图2.本发明所测量的样品结构示意图。图中di为对应材料的厚度,为对应材料的复折射率。
图3.本发明采用的单层MoS2的光学显微镜照片。
图4.测量得到的280nm SiO2/Si衬底上的单层MoS2的反射谱。
图5.测量得到的280nmSiO2/Si衬底上的没有单层MoS2处的反射谱。
图6.从图4中相应数据计算得到的280nm SiO2/Si衬底上的单层MoS2的衬度。
图7.衬度与SiO2厚度及照明波长的关系图。
图8.本发明得到的单层MoS2的复折射率谱。
具体实施方式
本发明方法,通过分别测量二维薄膜材料和衬底反射谱,得到衬度谱并利用公式加以拟合,最终得到二维薄膜材料在整个照明波段各个波长所对应的折射率,包括实部和虚部。本方法适用于任何二维薄膜材料,尤其是不能被传统方法测量的小尺度材料。
本发明所用样品为二维薄膜材料长在固态衬底上。所用衬底为具有透明介质、非透明介质双层结构的任意固态衬底。测量所述样品的反射光谱(覆盖有二维材料区域的反射光谱和没有覆盖二维材料区域的反射光谱),利用公式(1),得到二维薄膜材料的衬度谱。改变透明介质厚度,得到不同透明介质厚度下的衬度谱。通过在固定照明波长下通过公式(2)拟合衬度与透明介质厚度的关系获得二维薄膜材料的复折射率参数。改变照明波长,最终获得二维薄膜材料的复折射率谱。这种脱离椭偏仪的折射率测量方法,可以保持高空间分辨率的同时获得宽谱范围的折射率参数。在以二维薄膜材料为功能单元的新型器件设计与制备方面有广泛的应用前景。
以单层MoS2/SiO2/Si样品为例,主要实验过程如图1所示,样品置于显微镜物镜下方,用卤钨灯对样品进行照明。调节样品与物镜距离,使目镜中可以看到清晰的图像。将收集光纤置于像平面位置,分别收集有MoS2区域的反射光谱和没有MoS2区域的反射光,测量结果导入光谱仪中得到反射光谱。步骤如下:
通过湿法腐蚀SiO2,获得不同厚度的SiO2层。
1)用CVD方法在不同厚度的SiO2/Si衬底上生长单层MoS2。
2)将单层MoS2/SiO2/Si样品置于光学显微测量***之下,采用白光光源(卤素灯)反射照明方式,结合光纤耦合输出光谱测量,分别获得样品上有单层MoS2和没有单层MoS2处的反射光谱。
3)利用公式对获得的反射光谱进行处理,得到二维薄膜材料衬度谱。式中,和分别表示有无MoS2处的反射光强度。
4)用具有不同SiO2厚度的衬底重复步骤1),2)和3)。例如,采用25组不同SiO2厚度的衬底。在此基础上,得到不同入射波长下衬度与介质厚度的关系。
5)在固定入射波长下,采用多层薄膜干涉的理论模型对实验中的样品结构(例如:空气/单层MoS2/SiO2/Si)进行计算,得到以二维薄膜材料折射率实部、虚部、和透明介质层厚度为自变量的衬度的表达式。
6)用该表达式对测得的衬度与介质厚度的关系曲线进行拟合,得到该波长下二维薄膜材料的复折射率。改变入射波长,可以得到整个测量波段下该二维薄膜材料的复折射率谱。
实施例1
下面以单层MoS2为例对本发明所述的方法进行介绍。
在SiO2/Si衬底上生长或转移二维薄膜材料。衬底上的SiO2的厚度从100nm变化到400nm,在此范围中选取约25个不同厚度,从而保证最后的测量精度。
在光学显微镜反射照明光路中,保持照明光功率稳定,分别测量衬底上有MoS2和没有MoS2处的反射光谱,如图2所示。从这两个光谱在同一波长下的值可以计算得到相应的衬度C(λ)值。所用的公式如下:
式中,和分别表示有无MoS2处的反射光强度。
采用上述方法测出不同厚度SiO2衬底上单层MoS2衬度谱。在此基础上,固定入射波长,得到单层MoS2衬度和SiO2厚度关系曲线。利用多层薄膜反射光强理论公式(2)结合公式(3)可以得到单层MoS2衬度与SiO2厚度关系的公式。
将公式(3)中的和用公式(2)相应光强参量表示,则得到不同波长下单层MoS2衬度值与SiO2厚度的理论关系。用该理论关系在固定入射波长下,对测得的单层MoS2衬度和SiO2厚度关系曲线进行拟合即可得到该波长下单层MoS2折射率的实部和虚部。对入射光波段所有光波长进行拟合即可得到单层MoS2折射率谱。
图3是本发明采用的单层MoS2的光学显微镜照片。图4是测量得到的280nm SiO2/Si衬底上的单层MoS2的反射谱。图5是测量得到的280nmSiO2/Si衬底上的没有单层MoS2处的反射谱。图6为从图4中相应数据计算得到的280nm SiO2/Si衬底上的单层MoS2的衬度。图7是衬度与SiO2厚度及照明波长的关系图。图8是本发明得到的单层MoS2的复折射率谱。
Claims (6)
1.一种测量二维薄膜材料复折射率谱的方法,其特征是,包括如下步骤:
1)在衬底上生长二维薄膜材料或者转移二维薄膜材料到衬底上,保证样品干净;
2)将带有二维薄膜材料的衬底置于光学显微镜之下,采用白光光源反射照明方式,结合光纤耦合输出对目标区域进行光谱测量,分别获得样品上有二维薄膜材料和没有二维薄膜材料处的反射光谱;
3)利用二位薄膜材料在衬底上的衬度公式对获得的反射光谱进行处理,得到二维薄膜材料衬度谱,衬度公式如下:
式中,I1(λ)和I0(λ)分别表示有二维薄膜材料处和没有二维薄膜材料处的反射光强度;
4)用具有不同透明介质厚度的衬底重复步骤1)至3)多次,得到不同入射波长下衬度与透明介质厚度的关系曲线;
5)在固定入射波长下,采用多层薄膜干涉的理论模型对所述的样品结构进行计算,得到以二维薄膜材料折射率实部、虚部和透明介质层厚度为自变量的衬度的表达式;
6)用该表达式对测得的衬度与介质厚度的关系曲线进行拟合,得到该波长下二维薄膜材料的复折射率;改变入射波长,得到整个测量波段下该二维薄膜材料的复折射率谱。
2.如权利要求1所述的测量二维薄膜材料复折射率谱的方法,其特征是,所用衬底为具有透明介质、非透明介质双层结构的固态衬底。
3.如权利要求1所述的测量二维薄膜材料复折射率谱的方法,其特征是,步骤2)中,通过光纤直径来控制所测区域的大小可以。
4.如权利要求1所述的测量二维薄膜材料复折射率谱的方法,其特征是,步骤4)中,重复执行步骤1)至3)至少20次。
5.如权利要求1所述的测量二维薄膜材料复折射率谱的方法,其特征是,改变衬底透明介质的厚度来获得多种不同厚度衬底所对应的衬度谱数据,利用多层薄膜反射光强的理论公式对该数据进行拟合,得到最终的复折射率谱,公式如下:
其中
是相应材料的复折射率,实数部分ni为折射率;虚数部分的ki决定于光波在吸收性介质中传播时的衰减,叫做吸收系数;di是相应材料的厚度,λ是入射光波长。
6.如权利要求1所述的测量二维薄膜材料复折射率谱的方法,其特征是,所用的二维薄膜材料通过微机械剥离或CVD方法制备在SiO2/Si衬底上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410331238.2A CN104111235B (zh) | 2014-07-11 | 2014-07-11 | 一种测量二维薄膜材料复折射率谱的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410331238.2A CN104111235B (zh) | 2014-07-11 | 2014-07-11 | 一种测量二维薄膜材料复折射率谱的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104111235A true CN104111235A (zh) | 2014-10-22 |
CN104111235B CN104111235B (zh) | 2016-10-05 |
Family
ID=51708111
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410331238.2A Expired - Fee Related CN104111235B (zh) | 2014-07-11 | 2014-07-11 | 一种测量二维薄膜材料复折射率谱的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104111235B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106017338A (zh) * | 2016-05-26 | 2016-10-12 | 国家纳米科学中心 | 一种确定薄膜连续性的临界厚度的方法 |
CN106770039A (zh) * | 2017-03-10 | 2017-05-31 | 厦门大学嘉庚学院 | 一种复折射率测量装置及其测量方法 |
CN109211837A (zh) * | 2018-08-21 | 2019-01-15 | 厦门大学嘉庚学院 | 一种液体吸收介质的复折射率测量方法 |
CN109883347A (zh) * | 2019-03-14 | 2019-06-14 | 北京大学 | 一种基于光学三倍频的测量二维材料中应变张量的装置及方法 |
CN110057401A (zh) * | 2019-04-19 | 2019-07-26 | 华中科技大学 | 一种透明超薄膜折射率及厚度测量方法 |
CN110187434A (zh) * | 2019-05-16 | 2019-08-30 | 南京星隐科技发展有限公司 | 光波导及电磁波传输装置 |
CN114858751A (zh) * | 2022-02-27 | 2022-08-05 | 西北工业大学 | 一种测量二维材料复折射率的方法 |
CN115980902A (zh) * | 2022-12-15 | 2023-04-18 | 华中科技大学 | 一种基于二维硒化亚锡的可调波长偏振器及其设计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1405550A (zh) * | 2002-10-31 | 2003-03-26 | 中国科学院上海技术物理研究所 | 测量光学薄膜等效折射率及物理厚度的设备和方法 |
CN101865827A (zh) * | 2010-06-11 | 2010-10-20 | 山东大学 | 一种磁光椭偏测量装置及测量方法 |
CN102362146A (zh) * | 2009-03-23 | 2012-02-22 | Snu精度株式会社 | 反射度分布曲线的建模方法及应用该方法的厚度检测方法以及厚度检测反射仪 |
US20130038883A1 (en) * | 2011-08-12 | 2013-02-14 | Otsuka Electronics Co., Ltd. | Optical characteristic measuring apparatus and optical characteristic measuring method |
CN102967581A (zh) * | 2012-10-18 | 2013-03-13 | 中国人民解放军电子工程学院 | 薄膜夹及具有该薄膜夹的膜系材料红外复折射率测定装置 |
-
2014
- 2014-07-11 CN CN201410331238.2A patent/CN104111235B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1405550A (zh) * | 2002-10-31 | 2003-03-26 | 中国科学院上海技术物理研究所 | 测量光学薄膜等效折射率及物理厚度的设备和方法 |
CN102362146A (zh) * | 2009-03-23 | 2012-02-22 | Snu精度株式会社 | 反射度分布曲线的建模方法及应用该方法的厚度检测方法以及厚度检测反射仪 |
CN101865827A (zh) * | 2010-06-11 | 2010-10-20 | 山东大学 | 一种磁光椭偏测量装置及测量方法 |
US20130038883A1 (en) * | 2011-08-12 | 2013-02-14 | Otsuka Electronics Co., Ltd. | Optical characteristic measuring apparatus and optical characteristic measuring method |
CN102967581A (zh) * | 2012-10-18 | 2013-03-13 | 中国人民解放军电子工程学院 | 薄膜夹及具有该薄膜夹的膜系材料红外复折射率测定装置 |
Non-Patent Citations (3)
Title |
---|
R.ULRICH ET AL.: "Measurement of Thin Film Parameters with a Prism Coupler", 《APPLIED OPTICS》 * |
傅竹西等: "ZnO薄膜光学常数测量", 《发光学报》 * |
魏斌等: "偶氮金属螯合物薄膜的光学常数和吸收光谱", 《光学学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106017338B (zh) * | 2016-05-26 | 2018-08-24 | 国家纳米科学中心 | 一种确定薄膜连续性的临界厚度的方法 |
CN106017338A (zh) * | 2016-05-26 | 2016-10-12 | 国家纳米科学中心 | 一种确定薄膜连续性的临界厚度的方法 |
CN106770039B (zh) * | 2017-03-10 | 2023-04-21 | 厦门大学嘉庚学院 | 一种复折射率测量装置及其测量方法 |
CN106770039A (zh) * | 2017-03-10 | 2017-05-31 | 厦门大学嘉庚学院 | 一种复折射率测量装置及其测量方法 |
CN109211837A (zh) * | 2018-08-21 | 2019-01-15 | 厦门大学嘉庚学院 | 一种液体吸收介质的复折射率测量方法 |
CN109211837B (zh) * | 2018-08-21 | 2020-12-25 | 厦门大学嘉庚学院 | 一种液体吸收介质的复折射率测量方法 |
CN109883347A (zh) * | 2019-03-14 | 2019-06-14 | 北京大学 | 一种基于光学三倍频的测量二维材料中应变张量的装置及方法 |
CN110057401A (zh) * | 2019-04-19 | 2019-07-26 | 华中科技大学 | 一种透明超薄膜折射率及厚度测量方法 |
CN110187434A (zh) * | 2019-05-16 | 2019-08-30 | 南京星隐科技发展有限公司 | 光波导及电磁波传输装置 |
CN114858751A (zh) * | 2022-02-27 | 2022-08-05 | 西北工业大学 | 一种测量二维材料复折射率的方法 |
CN114858751B (zh) * | 2022-02-27 | 2024-05-10 | 西北工业大学 | 一种测量二维材料复折射率的方法 |
CN115980902A (zh) * | 2022-12-15 | 2023-04-18 | 华中科技大学 | 一种基于二维硒化亚锡的可调波长偏振器及其设计方法 |
CN115980902B (zh) * | 2022-12-15 | 2024-01-16 | 华中科技大学 | 一种基于二维硒化亚锡的可调波长偏振器及其设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104111235B (zh) | 2016-10-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104111235B (zh) | 一种测量二维薄膜材料复折射率谱的方法 | |
CN110398479B (zh) | 一种基于光学芯片基底的显微光谱测量装置及方法 | |
US10145785B2 (en) | Optical element rotation type Mueller-matrix ellipsometer and method for measuring Mueller-matrix of sample using the same | |
KR101656436B1 (ko) | 박막 웨이퍼의 막두께 분포 측정 방법 | |
CN107270822B (zh) | 测定多孔薄膜厚度和孔隙率的方法 | |
US11022752B2 (en) | Optical fibers having metallic micro/nano-structure on end-facet, and fabrication method, and application method thereof | |
KR20130018553A (ko) | 막 두께 측정 장치 | |
Goyal et al. | Porous photonic crystal structure for sensing applications | |
CN107389653B (zh) | 确定含有单层石墨烯区的石墨烯样品堆垛次序的方法 | |
CN106054300B (zh) | 一种co2气体检测用双通道红外滤光片及其制备方法 | |
CN107121080A (zh) | 一种测量有序多孔纳米薄膜厚度的方法 | |
Karacali et al. | Novel design of porous silicon based sensor for reliable and feasible chemical gas vapor detection | |
Xie et al. | Dual-degree-of-freedom multiplexed metasensor based on quasi-BICs for boosting broadband trace isomer detection by THz molecular fingerprint | |
Valenta et al. | Radiometric calibration of optical microscopy and microspectroscopy apparata over a broad spectral range using a special thin-film luminescence standard | |
CN108332674B (zh) | 单光谱求取多孔膜厚度和孔隙率的方法 | |
Baranov et al. | Broadband antireflective coatings based on two-dimensional arrays of subwavelength nanopores | |
CN101881672B (zh) | 一种反射式光纤温度传感器 | |
Khanikar et al. | A fiber optic refractive index sensor with a high index ZnO overlay | |
CN101697020B (zh) | 一种双通道光谱能量调谐滤光片 | |
CN112198140B (zh) | 基于多层介质膜干涉滤光片的液体折射率测量方法及其装置 | |
Mann et al. | Spectral imaging, reflectivity measurements, and modeling of iridescent butterfly scale structures | |
Guo et al. | Single point Linnik white-light spectral microscopic interferometer for surface measurement | |
CN105628673B (zh) | 测试复合硅衬底上过渡金属硫族化合物样品层数的方法 | |
Damin et al. | Application of scanning angle Raman spectroscopy for determining the location of buried polymer interfaces with tens of nanometer precision | |
Geng et al. | A scanning surface plasmon resonance sensor based on the phase shift algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
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: 20161005 Termination date: 20200711 |