CN103744034B - 一种提高cpt原子磁力仪灵敏度和绝对精度的差分方法 - Google Patents

一种提高cpt原子磁力仪灵敏度和绝对精度的差分方法 Download PDF

Info

Publication number
CN103744034B
CN103744034B CN201310753433.XA CN201310753433A CN103744034B CN 103744034 B CN103744034 B CN 103744034B CN 201310753433 A CN201310753433 A CN 201310753433A CN 103744034 B CN103744034 B CN 103744034B
Authority
CN
China
Prior art keywords
laser
magnetic field
polarization
hand circular
regulate
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.)
Expired - Fee Related
Application number
CN201310753433.XA
Other languages
English (en)
Other versions
CN103744034A (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.)
China Aero Geophysical Survey & Remote Sensing Center For Land And Resources
Hangzhou Mutual Inductance Information Technology Co ltd
Original Assignee
China Aero Geophysical Survey & Remote Sensing Center For Land And Resources
Hangzhou Mutual Inductance Information 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 China Aero Geophysical Survey & Remote Sensing Center For Land And Resources, Hangzhou Mutual Inductance Information Technology Co ltd filed Critical China Aero Geophysical Survey & Remote Sensing Center For Land And Resources
Priority to CN201310753433.XA priority Critical patent/CN103744034B/zh
Publication of CN103744034A publication Critical patent/CN103744034A/zh
Application granted granted Critical
Publication of CN103744034B publication Critical patent/CN103744034B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开一种提高CPT原子磁力仪灵敏度和绝对精度的差分方法。该方法是恒温下将弱磁场探头置于待测磁场中;调节激光器控制电路,保持激光波长稳定,且保持激光波长与碱金属原子基态能级共振;将二分之一波片、偏振分束棱镜垂直置于激光,调节二分之一波片光轴与激光偏振的角度;将四分之一波片垂直于两束平行激光;调节四分之一波片光轴与激光偏振的角度,平行激光的偏振分别变为左旋圆偏振和右旋圆偏振;两个光电转换器件对左旋圆偏振激光、右旋圆偏振激光探测,将探测到的光强值经减法器相减,经数据采集处理设备采集、信号处理,得到待测磁场的磁感应强度。本发明的优点是操作简单、灵敏度高、精度高。

Description

一种提高CPT原子磁力仪灵敏度和绝对精度的差分方法
技术领域
本发明属于弱磁测量技术领域,涉及一种提高CPT原子磁力仪灵敏度和绝对精度的差分方法,用于提高CPT(CoherentPopulationTrapping)原子磁力仪的灵敏度和绝对精度。
背景技术
在现代科技中,微弱磁场的探测是一项非常重要的技术。目前,原子磁力仪作为主要的微弱磁场探测仪器之一,种类主要有质子磁力仪、光泵磁力仪、CPT磁力仪等等,其中有的灵敏度较高,但绝对精度较差,有的具有高绝对精度,但灵敏度较低,无法同时保证磁力仪的灵敏度和绝对精度性能。现有的应用于原子磁力仪的技术中,没有可以同时提高灵敏度和绝对精度的技术方法,传统的差分方法由于使用含有磁场信息和不含有磁场信息的两个信号相减,只能提高原子磁力仪的灵敏度,无法提高其绝对精度。
发明内容
为了克服现有原子磁力仪无法同时保证高灵敏度和高绝对精度的问题,本发明提供了一种用于同时提高CPT原子磁力仪灵敏度和绝对精度的差分方法。这种差分方法通过使用两个均含有磁场信息的信号相减,可以同时优化CPT原子磁力仪的灵敏度和绝对精度性能。
本发明的工作机理是:在待测磁场下,碱金属原子能级将发生分裂,分裂情况与待测磁场大小相关。一束圆偏振激光通过碱金属原子气体时,若激光频率与分裂的碱金属原子基态能级完全共振,则碱金属原子对激光的吸收率最低,测量透射光强大小与激光频率的关系可以得到待测磁场大小。受非磁场因素的影响,在同一待测磁场下,由偏振正交的两束圆偏振光测量到的透射光强与激光频率的关系存在差异,若将这两个信号进行差分,不仅可以消除共模噪声,提高仪器灵敏度,还可以去除非磁场因素对测量磁场的影响,提高仪器的绝对精度。
CPT原子磁力仪包括激光光源、弱磁场探头、信号探测器;激光光源、弱磁场探头、信号探测器通过激光光路连接;
所述的激光光源由半导体激光器、激光器控制电路、二分之一波片、偏振分束棱镜、四分之一波片构成;
激光器控制电路控制半导体激光器发射出激光,激光依次透过二分之一波片、偏振分束棱镜、四分之一波片、弱磁场探头;
所述的激光控制电路由电流源、温度控制器、微波源、信号发生器构成,其中电流源与温度控制器直接控制半导体激光器,信号发生器的输出信号用于调制微波源,经过调制后的微波信号由微波源输出给半导体激光器。
所述的弱磁场探头由包含碱金属饱和蒸汽的玻璃泡构成;
所述的信号探测器由两个光电转换器件、减法器、数据采集处理设备构成;
两个光电转换器件采集透过弱磁场探头的光束的电信号,然后传递给减法器、数据采集处理设备处理。
本发明对CPT原子磁力仪的调节方法具体是:
步骤(1).保持恒温状态下,将弱磁场探头置于待测磁场中;
步骤(2).调节激光光源:
调节激光器控制电路中电流源、温度控制器、微波源、信号发生器,保持半导体激光器发射出的激光波长稳定,且使得半导体激光器发射出的激光波长与弱磁场探头中碱金属原子基态能级共振;将二分之一波片、偏振分束棱镜均垂直置于激光光束方向;激光波长均在所选的二分之一波片与偏振分束棱镜的波长范围内;调节二分之一波片光轴与激光偏振的角度,使得由半导体激光器发射出的一束激光经偏振分束棱镜分开后形成的两束平行激光光强相等;将四分之一波片垂直于两束平行激光的传播方向;激光波长均在所选的四分之一波片的波长范围内;调节四分之一波片光轴与激光偏振的角度,使得两束平行激光的偏振分别变为左旋圆偏振和右旋圆偏振;调节弱磁场探头的位置,使左旋圆偏振激光、右旋圆偏振激光透过弱磁场探头;
步骤(3).调节信号探测器:
两个光电转换器件分别对透过弱磁场探头的左旋圆偏振激光、右旋圆偏振激光进行探测,得到光电转换器件探测到的光强值与信号发生器输出信号频率的关系,关系如公式(1)、公式(2):
Y 1 = k π [ υ ( X - f 0 + Δ f ) 2 + υ 2 ] + N c ( X ) - - - ( 1 ) ;
Y 2 = k π [ υ ( X - f 0 - Δ f ) 2 + υ 2 ] + N c ( X ) - - - ( 2 ) ;
其中,Y1与Y2为两个光电转换器件探测到的光强值,X为信号发生器输出信号频率,π是圆周率,k为比例系数,是定值,υ为信号线宽,对于恒定温度下的包含碱金属饱和蒸汽的玻璃泡9,υ也为定值,Nc(X)为共模噪声,Δf为非磁场因素引起的碱金属原子的能级分裂,f0为由待测磁场引起的碱金属原子的能级分裂,其关系如公式(3):
f0=γ·B(3);
其中,γ为定值。
将光电转换器件探测到的光强值经过减法器相减,得到两个光电转换器件探测到的光强差值ΔY,如公式(4):
ΔY=Y1-Y2(4);
然后经数据采集处理设备采集、信号处理,根据公式(1)、(2)、(4),求解ΔY=0时X对应的值f0,最后根据公式(3)得到待测磁场的磁感应强度B。
传统的CPT磁力仪的方法是利用一束激光得到公式(1)或者公式(2)中的关系,通过求Y1或Y2为极值时X对应的值f0,根据公式(3)得到待测磁场的磁感应强度B。根据公式(1)或公式(2)可知,Y1或Y2为极值时X对应的值f0不仅受到共模噪声Nc(X)的影响,灵敏度较差,还受到非磁场因素引起的碱金属原子的能级分裂Δf的影响,绝对精度也较差。
利用本发明方法后,通过求解ΔY=0得到X对应的值f0,根据公式(3)得到待测磁场的磁感应强度B。根据公式(4)可知,ΔY中不包含共模噪声Nc(X)的影响,同时,求解ΔY=0得到X对应值f0的方法也避免了非磁场因素引起的碱金属原子的能级分裂Δf的影响。因此,本发明同时提高了灵敏度和绝对精度。
对于传统的CPT磁力仪,可以不采用信号发生器调制微波源的方式,但对于本发明方法,需要使用信号发生器调制微波源的方式。
本发明的优点是:一、操作简单,仅需操作CPT原子磁力仪中的激光光源、弱磁场探头和信号探测器三个部分;二、提高灵敏度,这是由于差分方法去除了共模噪声;三、提高绝对精度,这是由于差分方法去除了非磁场因素对测量磁场的影响。
附图说明
图1为本发明的流程示意图;
图2为本发明的详细流程示意图;
图3为本发明激光控制电路的流程示意图。
具体实施方式
下面结合附图对本发明做进一步的分析。
在待测磁场下,碱金属原子能级将发生分裂,分裂情况与待测磁场大小相关。一束圆偏振激光通过碱金属原子气体时,若激光频率与分裂的碱金属原子基态能级完全共振,则碱金属原子对激光的吸收率最低,测量透射光强大小与激光频率的关系可以得到待测磁场大小。受非磁场因素的影响,在同一待测磁场下,由偏振正交的两束圆偏振光测量到的透射光强与激光频率的关系存在差异,若将这两个信号进行差分,不仅可以去除共模噪声,提高仪器灵敏度,还可以去除非磁场因素对测量磁场的影响,提高仪器的绝对精度。
如图1所示,CPT原子磁力仪包括激光光源1、弱磁场探头2、信号探测器3;激光光源1、弱磁场探头2、信号探测器3通过激光光路连接;
如图2所示,所述的激光光源1由半导体激光器4、激光器控制电路5、二分之一波片6、偏振分束棱镜7、四分之一波片8构成;
如图3所示,激光器控制电路5控制半导体激光器4发射出激光,激光依次透过二分之一波片6、偏振分束棱镜7、四分之一波片8、弱磁场探头2;
所述的激光控制电路5由电流源13、温度控制器14、微波源15、信号发生器16构成,其中电流源13与温度控制器14直接控制半导体激光器4,信号发生器16的输出信号用于调制微波源15,经过调制后的微波信号由微波源输出给半导体激光器4。
所述的弱磁场探头2由包含碱金属饱和蒸汽的玻璃泡9构成;
所述的信号探测器3由两个光电转换器件10、减法器11、数据采集处理设备12构成;
两个光电转换器件10采集透过弱磁场探头2的光束的电信号,然后传递给减法器11、数据采集处理设备12处理。
具体调节CPT原子磁力仪的差分方法是:
实施例中弱磁场探头2碱金属原子采用铷-87原子,铷原子饱和蒸汽的玻璃泡9尺寸为Φ15×20mm,将其周围均匀加热且恒温至45℃,并置于待测磁场中。在使用过程中,先开启激光器控制电路5,其中电流源13使用美国Agilent公司生产的型号为B2912A的电流源,温度控制器14使用美国Thorlab公司生产的型号为TED200C的温度控制器,微波源15使用美国Agilent公司生产的型号为E8267D的微波源,信号发生器16使用中国普源公司生产的型号为DG4162的信号发生器,调节激光二极管的电流为1.3mA,温度为60℃,使半导体激光器4的波长稳定至794.984nm;将适用波长为794.984nm的二分之一波片6和偏振分束棱镜7垂直置于激光光束方向,调节二分之一波片6的相对角度至43°,使得由偏振分束棱镜7分开的两束平行激光光强相等;将适用波长为794.984nm的四分之一波片8垂直置于两束平行激光方向,调节四分之一波片8的相对角度至84°,使得两束偏振正交的线偏振平行激光变为两束偏振正交的圆偏振激光;调节铷原子饱和蒸汽的玻璃泡9的位置,使两束平行圆偏振激光纵向通过铷原子饱和蒸汽的玻璃泡9,并使用两个高灵敏度硅光电二极管10分别对两束平行圆偏振激光进行探测,将两个高灵敏度硅光电二极管10探测采集到的信号输入给减法器11相减,相减后的信号输入给数据采集处理设备12中锁相放大器进行相敏检测,并将锁相放大器的输出信号输入计算机进行采集、处理和最终输出。
上述提到的高灵敏度硅光电二极管10为光电转换器件10。
上述实施例并非是对于本发明的限制,本发明并非仅限于上述实施例,只要符合本发明要求,均属于本发明的保护范围。

Claims (3)

1.一种提高CPT原子磁力仪灵敏度和绝对精度的差分方法,其特征在于该方法包括以下步骤:
步骤(1).保持恒温状态下,将弱磁场探头置于待测磁场中;
步骤(2).调节CPT原子磁力仪中的激光光源:
激光光源由半导体激光器、激光器控制电路、二分之一波片、偏振分束棱镜、四分之一波片构成;
首先调节激光器控制电路中电流源、温度控制器、微波源、信号发生器,保持半导体激光器发射出的激光波长稳定,且保持半导体激光器发射出的激光波长与弱磁场探头中碱金属原子基态能级共振;将二分之一波片、偏振分束棱镜均垂直置于激光光束方向;调节二分之一波片光轴与激光偏振的角度,使得由半导体激光器发射出的一束激光经偏振分束棱镜分开后形成的两束平行激光光强相等;将四分之一波片垂直于两束平行激光的传播方向;调节四分之一波片光轴与激光偏振的角度,使得两束平行激光的偏振分别变为左旋圆偏振和右旋圆偏振;调节弱磁场探头的位置,使左旋圆偏振激光、右旋圆偏振激光透过弱磁场探头;
步骤(3).调节CPT原子磁力仪中的信号探测器:
信号探测器由两个光电转换器件、减法器、数据采集处理设备构成;
两个光电转换器件分别对透过弱磁场探头的左旋圆偏振激光、右旋圆偏振激光进行探测,得到光电转换器件探测到的光强值与信号发生器输出信号频率的关系,关系如公式(1)、公式(2):
Y 1 = k π [ υ ( X - f 0 + Δ f ) 2 + υ 2 ] + N c ( X ) - - - ( 1 ) ;
Y 2 = k π [ υ ( X - f 0 - Δ f ) 2 + υ 2 ] + N c ( X ) - - - ( 2 ) ;
其中,Y1与Y2为两个光电转换器件探测到的光强值,X为信号发生器输出信号频率,π是圆周率,k为比例系数,是定值,υ为信号线宽,对于恒定温度下的包含碱金属饱和蒸汽的玻璃泡,υ也为定值,Nc(X)为共模噪声,Δf为非磁场因素引起的碱金属原子的能级分裂,f0为由待测磁场引起的碱金属原子的能级分裂,其关系如公式(3):
f0=γ·B(3);其中,γ为定值;
将光电转换器件探测到的光强值经过减法器相减,得到两个光电转换器件探测到的光强差值ΔY,如公式(4):
ΔY=Y1-Y2(4);
然后经数据采集处理设备采集、信号处理,根据公式(1)、(2)、(4),求解ΔY=0时X对应的值f0,最后根据公式(3)得到待测磁场的磁感应强度B。
2.如权利要求1所述的一种提高CPT原子磁力仪灵敏度和绝对精度的差分方法,其特征在于弱磁场探头由包含碱金属饱和蒸汽的玻璃泡构成。
3.如权利要求1所述的一种提高CPT原子磁力仪灵敏度和绝对精度的差分方法,其特征在于激光波长在二分之一波片、偏振分束棱镜、四分之一波片的波长范围内。
CN201310753433.XA 2013-12-30 2013-12-30 一种提高cpt原子磁力仪灵敏度和绝对精度的差分方法 Expired - Fee Related CN103744034B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310753433.XA CN103744034B (zh) 2013-12-30 2013-12-30 一种提高cpt原子磁力仪灵敏度和绝对精度的差分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310753433.XA CN103744034B (zh) 2013-12-30 2013-12-30 一种提高cpt原子磁力仪灵敏度和绝对精度的差分方法

Publications (2)

Publication Number Publication Date
CN103744034A CN103744034A (zh) 2014-04-23
CN103744034B true CN103744034B (zh) 2016-01-20

Family

ID=50501069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310753433.XA Expired - Fee Related CN103744034B (zh) 2013-12-30 2013-12-30 一种提高cpt原子磁力仪灵敏度和绝对精度的差分方法

Country Status (1)

Country Link
CN (1) CN103744034B (zh)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104698410B (zh) * 2015-03-02 2019-03-01 北京大学 用于磁力仪的原子磁传感器及消除磁力仪探测盲区的方法
CN104698404B (zh) * 2015-03-02 2018-07-17 北京大学 一种用于全光光泵磁力仪的原子磁传感器
CN104677508B (zh) * 2015-03-16 2017-09-12 北京航空航天大学 一种基于圆偏振探测光的原子自旋进动检测方法及装置
CN105068025B (zh) * 2015-07-16 2017-08-11 山西大学 基于电磁感应透明效应测量微弱磁场场强的方法及装置
CN105699919B (zh) * 2016-03-01 2018-07-17 中国科学院武汉物理与数学研究所 一种差分探测相干布居囚禁磁强计的实现方法
CN106019179B (zh) * 2016-05-19 2018-09-18 北京航天控制仪器研究所 一种cpt磁力仪快速锁频方法
CN106073751B (zh) * 2016-05-30 2019-06-04 中国科学院电工研究所 一种便携式脑磁检测装置
CN105866716B (zh) * 2016-06-23 2018-08-10 梁尚清 一种新型全光型激光光泵磁力仪及其实现方法
CN106595626A (zh) * 2016-12-19 2017-04-26 北京航空航天大学 一种核磁共振陀螺仪椭圆偏振光检测系统
CN107656219A (zh) * 2017-08-02 2018-02-02 兰州空间技术物理研究所 一种铷原子磁力仪
CN110045430B (zh) * 2018-01-17 2021-03-30 杭州电子科技大学 一种地磁日变监测方法
CN110045309B (zh) * 2018-01-17 2021-03-30 杭州电子科技大学 一种新型的光泵磁力仪精度自检测方法
CN108614224B (zh) * 2018-04-03 2020-09-18 北京航天控制仪器研究所 一种用于cpt磁力仪的气室工作温度自动标定系统及方法
CN108693490A (zh) * 2018-04-25 2018-10-23 中国科学技术大学 磁场强度测量装置及方法
CN109342980A (zh) * 2018-10-31 2019-02-15 浙江工业大学 基于椭圆光的单光Mx原子磁力仪
CN109839606B (zh) * 2019-03-13 2020-11-06 中国科学院武汉物理与数学研究所 一种新型原子磁强计装置及检测方法
CN110649923B (zh) * 2019-11-13 2021-08-17 浙江大学 一种双频探测相干布居囚禁原子钟及其工作方法
CN111030712B (zh) * 2019-11-27 2021-08-06 山东航天电子技术研究所 一种透地磁场通信装置
CN111007443B (zh) * 2019-11-28 2022-07-29 北京航天控制仪器研究所 一种采用cpt磁力仪进行弱磁场测量的方法及cpt磁力仪
CN111610470B (zh) * 2020-05-09 2022-08-16 杭州电子科技大学 一种新型射频原子磁力仪及其实现方法
CN112304442B (zh) * 2020-10-28 2022-09-02 中国科学院国家授时中心 双调制cpt差分探测方法及系统
CN113655413B (zh) * 2021-07-21 2023-07-28 杭州电子科技大学 一种光纤式原子磁力仪探头及其调节方法
CN113791370A (zh) * 2021-08-12 2021-12-14 北京量子信息科学研究院 磁强计和磁场强度确定方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6472869B1 (en) * 2001-06-18 2002-10-29 United States Of America As Represented By The Secretary Of The Air Force Diode laser-pumped magnetometer
CN101441253A (zh) * 2008-12-02 2009-05-27 浙江大学 高灵敏度原子磁力仪

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7038450B2 (en) * 2002-10-16 2006-05-02 Trustees Of Princeton University High sensitivity atomic magnetometer and methods for using same
US8605282B2 (en) * 2009-05-08 2013-12-10 Ben-Gurion University Of The Negev Research And Development Authority Method and apparatus for high precision spectroscopy
JP5446731B2 (ja) * 2009-10-29 2014-03-19 セイコーエプソン株式会社 磁場測定装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6472869B1 (en) * 2001-06-18 2002-10-29 United States Of America As Represented By The Secretary Of The Air Force Diode laser-pumped magnetometer
CN101441253A (zh) * 2008-12-02 2009-05-27 浙江大学 高灵敏度原子磁力仪

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
全光学高灵敏度铷原子磁力仪的研究;李曙光 等;《物理学报》;20100228;第59卷(第2期);877-882 *
原子磁力仪的研究;李曙光;《中国博士论文全文数据库》;20090401;1-104 *

Also Published As

Publication number Publication date
CN103744034A (zh) 2014-04-23

Similar Documents

Publication Publication Date Title
CN103744034B (zh) 一种提高cpt原子磁力仪灵敏度和绝对精度的差分方法
US10288701B2 (en) Optically pumped atomic magnetometer and magnetic sensing method
US9244137B2 (en) Optical pumping magnetometer and magnetic sensing method
CN104062608B (zh) 一种serf原子自旋磁强计光位移消除方法
CN105866716B (zh) 一种新型全光型激光光泵磁力仪及其实现方法
CN110401492A (zh) 一种基于量子效应的无线电调幅信号接收方法及调幅量子接收机
Luhmann Jr et al. Chapter 3: microwave diagnostics
JP2017026405A (ja) 光ポンピング磁力計および磁気センシング方法
CN104714110A (zh) 基于电磁诱导透明效应测量高频微波场强的装置和方法
CN104698410A (zh) 用于磁力仪的原子磁传感器及消除磁力仪探测盲区的方法
CN107085130A (zh) 采用主动温度补偿的偏振不敏感电流和磁场传感器
CN102928110B (zh) 一种非接触式原子气室温度测量装置及测量方法
CN108535668B (zh) 一种补偿激光原子磁力计磁屏蔽罩内部剩余磁场的方法
CN106226713A (zh) 一种serf原子磁强计的光频移抑制方法
CN111492252B (zh) 用于测量在电能传输导体中的电场和/或磁场的装置
CN103969604A (zh) 射频原子磁力仪及其测量核磁共振信号的方法
CN105929458A (zh) 航空磁场矢量检测装置及监测方法
CN110045309B (zh) 一种新型的光泵磁力仪精度自检测方法
CN107389315A (zh) 光器件频响测量方法及测量装置
CN109631959A (zh) 一种基于光纤Sagnac干涉的原子自旋进动检测装置信号抗干扰能力优化方法
CN102129050B (zh) 基于光栅调制相位比较等离子体磁场测量方法与装置
CN110207822A (zh) 高灵敏度光学时延估计系统、方法及介质
CN203859921U (zh) 一种测量等离子体电流分布的极化ece诊断系统
US8947663B2 (en) Dual-modulation faraday rotation spectroscopy
Akiyama et al. Recent Progress on dispersion interferometers for nuclear fusion and low-temperature plasmas

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

Granted publication date: 20160120

Termination date: 20171230

CF01 Termination of patent right due to non-payment of annual fee