CN1269328C - 优化输入端的低窜扰平顶型波分复用器件 - Google Patents

优化输入端的低窜扰平顶型波分复用器件 Download PDF

Info

Publication number
CN1269328C
CN1269328C CNB021120544A CN02112054A CN1269328C CN 1269328 C CN1269328 C CN 1269328C CN B021120544 A CNB021120544 A CN B021120544A CN 02112054 A CN02112054 A CN 02112054A CN 1269328 C CN1269328 C CN 1269328C
Authority
CN
China
Prior art keywords
width
broadening
input
waveguide
quality factor
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
CNB021120544A
Other languages
English (en)
Other versions
CN1383284A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CNB021120544A priority Critical patent/CN1269328C/zh
Publication of CN1383284A publication Critical patent/CN1383284A/zh
Application granted granted Critical
Publication of CN1269328C publication Critical patent/CN1269328C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Optical Integrated Circuits (AREA)

Abstract

本发明公开了一种优化输入端的低窜扰平顶型波分复用器件及构成方法。它包括输入波导自由扩散区,刻蚀凹面光栅,输出波导阵列,多模干涉部件。在输入波导的末端接有带有预展宽输入波导的曲线型多模干涉部件。本发明的优点是:1)预展宽输入的锥形波导,能把输入波导的模场进行预展宽,使平坦化频谱获得非常陡的下降沿,从而有效的降低了窜扰;2)多模干涉区是根据自映像原理设计,通过预展宽输入波导输入产生模场激励,形成多模干涉,最终在输出得到分开一定距离的两个像,从而实现平坦化;3)通过抛物线线形的优化调节,能有效的降低带通纹波,增大1dB带宽。本发明具有较多的一般性,能够适用于AWG和EDG器件中的任何一种。

Description

优化输入端的低窜扰平顶型波分复用器件
技术领域
本发明属于光通讯波分复用领域,特别涉及一种具有低窜扰、宽带通、低带通纹波、大的1dB带宽的波分复用器件。
背景技术
波分复用/解复用是现代光通讯技术的核心技术之一,它是指把不同波长的光变作复合光及将复合光中各组成波长的光分离出来。波分复用器件就是用来实现波分复用/解复用的器件。现在很有发展潜力的波分复用/解复用器件主要有阵列波导光栅(AWG)和刻蚀衍射光栅(EDG)两种。
传统的波分复用/解复用器件的带通是高斯型,这使得实际中心波长偏离设计中心时,造成透过率大大下降。这必然增加对光源波长严格控制的要求,同时环境温度改变、材料的老化、以及器件刻蚀过程里存在的公差等外界因素都会对器件性能造成大的损害。平坦的频谱响应能够有效的消除这些因素对器件性能造成的影响。
已有一些方法用来实现波分复用/解复用器件频谱的平坦化,例如:
美国专利NO.5412744提出使用Y分支波导实现频谱平坦化,由于Y分支尖角的模场泄漏会造成大的损耗。
K.Okamoto等发表了一篇题为“Eight-Channel flat spectral responseArrayed-waveguide multiplexer with asymmetrical Mach-Aehnder filters”,IEEEPhot.Tech.Lett.,vol.8,no.3,march 1996,pp.373-374.文章涉及一种使用M-Z干涉仪来实现平坦化,但是这种方法会造成大的器件尺寸,大的插损,同时还有可能引起模斑变形。
Y.P.Ho等发表过一篇文章,“Flat channel-passband-wavelength multiplexingand demultiplexing devices by multiple Rowland circles”,IEEEPhot.Tech.Lett.,(9),pp.342-344,1997中提到使用双罗兰圆来实现平坦化。
A.Rigney等发表了令一篇文章,“Double-phased array for a flattened spectralresponse”,Proc.23rd ECOC,Edinburgh,UK,pp.79-82,Sept.1997,中提到使用双相位阵列来实现平坦化,该方法和上一种方法类似,平坦化效果相对有限,而且复杂度增大。
M.R.Amersfoort等发表了令一篇题为“Phased-arrayed wavelengthdemultiplexer with flattened wavelength response”,Electron.Lett.,1994,30,(4),pp.300-302的文章提出使用多模波导输出来实现平坦化,但这种方法只适用于波分复用/解复用器直接与探测器相连的情况。
T.Chiba等在OECC(2000),13B2-2里提到使用interleaver来实现信号复共轭实现平坦化,但是这种方法使器件尺寸过大,不宜于集成,同时价格也相对昂贵。
M.R.Amersfoort等发表了一篇题为,”Passband broadening of integratedarrayed waveguide filter using multimode interference couplers”,Electron.Lett.,vol.32(1996),pp.1661-1662,提出使用MMI耦合器实现平坦化,由于MMI结构简单、工艺容差大,可以方便的用于实现频谱平坦化。但是单纯的使用MMI平坦化的同时会带来大的窜扰、插损和带通纹波等不利因素。
发明内容
本发明的目的是研制一种把带间窜扰降到很低值,把带通纹波减小到接近于零的优化输入端的低窜扰平顶型波分复用器件及构成方法。
本发明采用的技术方案如下:
它的结构为输入波导,其后接有自由扩散区,自由扩散区后接有刻蚀凹面光栅,输入波导侧面有输出波导阵列,在输入波导的末端与自由扩散区之间接有多模干涉部件。输入波导末端与多模干涉部件前端之间接有带有预展宽输入波导,且多模干涉部件的多模干涉区的边界为抛物线结构;预展宽输入波导的初始端宽度d1大于输入波导的宽度,末端宽度d2大于初始端宽度d1,长度 L 1 = 3 8 L πeff = n r 2 λ ( d 2 + C ) ( d 1 + C ) , 其中包层折射率nc,芯层折射率nr,中心波长λ,Lπeff为拍长,波导等效厚度 C = λ π · ( n c n r ) 2 σ ( n r 2 - n c 2 ) ( - 1 / 2 ) , 对横电偏振σ的值是0,对横磁偏振σ的值是1;
1)预展宽区域尺寸的确定:
预展宽的最终目的是要出现高阶模次,改变模场分布,因此,通过BPM模拟分析这一过程里的模场分布,找到使得展宽模场边缘下降最陡的结构参数,
在任一截面上预展宽区域的宽度可以表示为: d ( z ) = d 2 + ( d 1 - d 2 ) ( 1 - z L 1 ) 其中z是传输方向,d1、d2分别为预展宽区域初始端宽度和末端宽度,L1是预展宽区域的长度, C = λ π · ( n c n r ) 2 σ ( n r 2 - n c 2 ) ( - 1 / 2 ) , 其中对横电偏振σ的值是0,对横磁偏振σ的值是1;
则可以表示预展宽区域的设计长度为: L 1 = 3 8 L πeff = n r 2 λ ( d 2 + C ) ( d 1 + C ) 其中包层折射率nc,芯层折射率nr,中心波长λ,Lπeff为拍长;
预展宽区域的三个参数d1,d2,L1中,长度L1是初始端宽度d1和末端宽度d2的函数,因此,只需要确定预展宽区域的宽度d1和d2,先固定初始端宽度d1,找出末端宽度d2,依据计算光谱和插入损耗的方法,分别找出不同末端宽度d2和品质因数、插入损耗的关系,作出曲线,随着末端的宽度d2的增大,品质因数和插入损耗都会随之增大,能得到品质因数大,而插入损耗小的末端宽度d2
2)曲线型多模干涉区为抛物线型干涉区,其结构设计:
初始干涉区宽度d3的优化其主要衡量指标是平坦化频谱的品质因数,品质因数大小为1dB带宽与33dB带宽的比值,在给定的预展宽末端宽度d2下,改变初始干涉区宽度d3,做出d3随品质因数的变化关系曲线,能找到一个最优的初始干涉区宽度d3,使得品质因数具有最大值;
抛物线线形和干涉区长度L2的优化调节能消除带通纹波,在上述干涉区初始宽度d3已经确定的情况下,抛物线必然经过其边缘点,同时其对称轴必然为器件中心对称轴,因此,决定抛物线线形的三个参数只剩下了一个,通过改变该参数,能得到最小的带通纹波,另外能在干涉区成像点附近找到一个使带通纹波最小,1dB带宽最大的干涉区长度L2
本发明具有的有益的效果是:
1)预展宽输入的锥形波导,能把输入波导的模场进行预展宽,使平坦化频谱获得非常陡的下降沿,从而有效的降低了窜扰;
2)多模干涉区是根据自映像原理设计,通过预展宽输入波导输入产生模场激励,形成多模干涉,最终在输出得到分开一定距离的两个像,从而实现平坦化。
3)通过抛物线线形的优化调节,能有效的降低带通纹波,增大1dB带宽。本发明具有较多的一般性,能够适用于AWG和EDG器件中的任何一种。
附图说明
图1是本发明的结构原理示意图;
图2是图1中输入波导末端带有预展宽输入的抛物线型MMI结构的放大图;
图3是带有预展宽输入的MMI结构;
图4是普通最优化MMI结构;
图5是图3,图4两种结构用于平坦化的频谱比较(证明预展宽结构的低窜扰特性);
图6是使用图3结构得到的低窜扰特性频谱;
图7是使用抛物线结构的MMI对带通纹波的消除和1dB带宽的增大;
图8是频带间隔200GHz下的最优频谱响应;
图9是频带间隔100GHz下的最优频谱响应;
图10是品质因数,1dB带宽随预展宽宽度变化的关系曲线;
图11是插入损耗随预展宽宽度变化的关系曲线;
图12是品质因数,1dB带宽随多模干涉区初始宽度变化的关系曲线。
具体实施方式
本发明具有较多的一般性,适用于AWG和EDG器件中的任何一种,为了方便说明,下面仅以EDG器件为例。
本发明所有结构设计基于硅基底的二氧化硅波导,包层折射率nc=1.445,芯层折射率nr=1.454,中心波长λ0=1.55μm。
如图1所示本发明包括输入波导1,自由扩散区2,刻蚀凹面光栅3,输出波导阵列4。复合光由输入波导1入射,进入自由扩散区2发散传输,然后经过刻蚀凹面光栅3的衍射作用,不同的波长会聚到相应的输出波导端口4。在输入波导1末端接有带预展宽输入的抛物线型多模干涉部件5,用来实现平坦化。
如图2所示,给出抛物线型多模干涉部件5的放大图,预展宽输入波导6是本发明的第一个改进。
预展宽区域尺寸的确定:
预展宽的最终目的是要出现高阶模次,改变模场分布,因此,能通过BPM模拟分析这一过程里的模场分布,找到使得展宽模场边缘下降最陡的结构参数。
在任一截面上预展宽区域的宽度能表示为: d ( z ) = d 2 + ( d 1 - d 2 ) ( 1 - z L 1 ) - - - ( 1 ) 其中z是传输方向,d1、d2分别为预展宽区域初始端宽度和末端宽度,L1是预展宽区域的长度。
如果表示 C = λ π · ( n c n r ) 2 σ ( n r 2 - n c 2 ) ( - 1 / 2 ) , 对横电偏振σ的值是0,对横磁偏振σ的值是1            (2)则可以表示预展宽区域的设计长度为: L 1 = 3 8 L πeff = n r 2 λ ( d 2 + C ) ( d 1 + C )
其中Lπeff表不拍长。
从上面的分析能看出,预展宽区域的三个参数d1、d2和L1中,长度L1是初始端宽度d1和末端宽度d2的函数,因此,只需要确定预展宽区域的宽度d1、d2
例如,固定初始端宽度为d1=6μm,找出合适的末端宽度d2。过程依据下面提到的计算光谱和插入损耗的方法,分别找出不同宽度d2和品质因数、插入损耗的关系,作出曲线如图10和图11所示,可见随着d2的增大,品质因数和插入损耗都会随之增大,而目标是要得到品质因数尽量大,插入损耗相对小的设计结构,所以应根据实际带宽需要,综合两者选出最为合适的d2大小。
我们将带有预展宽输入波导的MMI(如图3)和宽度最优化以后的MMI(如图4)的频谱做了比较,并将它们的频谱响应画在图5中。
其中图3各参数取值为:
d0=5μm;d1=6μm;d2=10μm;d3=25μm;L1=400μm;L2=475μm。
图4各参数取值为:
d0=5μm;d3′=18μm;L2′=160μm。
图5中虚线为通常的高斯频谱响应,点划线为使用普通最优化MMI频谱响应,实线则为使用带有预展宽输入的结构以后的频谱响应曲线。从图中能看到,使用改进的带有预展宽输入的MMI在20dB以下能有效降低窜扰,从图6可以看到,使用该结构的EDG频谱窜扰能降低到50dB左右,这是通常的平坦化难以达到的效果。
如图2所示,抛物线结构8是本发明的第二个重要改进,抛物线干涉区的结构设计如下:
初始干涉区宽度d3的优化其主要衡量指标是平坦化频谱的品质因数,其大小为1dB带宽与33dB带宽的比值。在给定的预展宽末端宽度d2下,改变初始干涉区宽度d3,做出其随品质因数的变化关系曲线,则相应能找到一个最优的初始干涉区宽度d3,此时品质因数具有最大值,此时对应的d3大小为设计最佳值。对应的搜索过程实例如图12所示,此时的d2大小为12μm,可以看到当d3≈24μm时,具有最大的品质因数。
抛物线线形和干涉区长度L2的优化调节能有效地消除带通纹波,在上述初始干涉区宽度d3已经确定的情况下,抛物线必然经过其边缘点,同时其对称轴必然为器件中心对称轴,因此,决定抛物线线形的三个参数只剩下了一个。通过改变该参数,能得到相对最小的带通纹波,另外在干涉区成像点附近进一步优化长度L2,同样能找到更加好的干涉区长度L2,使得带通纹波最小,1dB带宽最大。
图7中,实线为使用抛物线结构后的频谱响应,点划线为同样尺寸下使用图3结构的最优频谱响应,从图中可以看到,使用改进的抛物线结构8可以有效的降通纹波;另外从图中还可以看到,使用了抛物线结构8还使得1dB带宽明显增大,从而进一步的改善了平坦化性能。
图2中多模干涉区7是根据自映像原理设计,通过预展宽输入波导6输入激励模场激励,形成多模干涉,最终在输出得到分开一定距离的两个像,从而实现平坦化。
图8和图9分别给出了适用于频带间隔为200GHz和100GHz情况下的两种最优化频谱响应。
其中图8工作于Δf=200GHz,并在输出端用15μm的tapper接收能量。其具体参数为:
d0=5μm;d1=6μm;d2=18μm;d3=40μm;d4≈64μm;L1=400μm;L2=775μm。
其得到的性能参数如下:
1dB带宽:74.5600%;品质因数F=0.8543;插入损耗约4dB;窜扰:≤-120dB;带通纹波:≤0.1dB。
由于输出端用15μm的tapper接收能量,这样做使得器件尺寸相应增大,适用于通道数较少的情况,例如200GHz,其频谱相应已经接近于理想的矩形响应,其品质因数有85%左右,带通纹波也已经接近于零。
图9工作于Δf=100GHz,并直接用6μm的输出波导来接收能量。其具体结构参数如下:
d0=5μm;d1=6μm;d2=12μm;d3=24μm;d4≈31μm;L1=400μm;L2=550μm。
其得到的性能参数如下:
1dB带宽:53.4000%;品质因数F=0.7744;插入损耗约4.2dB;窜扰:≤-70dB;带通纹波:≤0.3dB。
此时由于输出端是直接使用宽度为6μm的输出波导来接收能量,这样做使得结构紧凑,适用于多通道的密集波分复用,图示为频带间隔为100GHz时的频谱响应,在保持了相对优越的平坦化频谱的同时其插损也保持了较小的水平。
器件的频谱响应计算如下:
通过BPM计算出经过MMI区的出射场分布,把其表示为E0,对得到的进入自由扩散区的初始场E0作快速傅立叶变换,其变换结果记为A(u)=FFT(E0),根据角谱理论,一个任意分布的空间模场的传播可以看作是其角谱以平面波形式在空间各点的叠加过程。因此,在空间任意一点,我们可以计算出其场值为:
U ( x , z ) = ∫ - ∞ + ∞ A ( u ) exp [ ikz 1 - ( λu ) 2 ] exp [ i 2 π ( ux ) ] du - - - ( 1 )
再由基尔霍夫衍射公式可以得到之后的输出模场为:
E ′ ( x ′ z ′ ) = 1 2 ( n λ ) 1 2 ∫ Gratings U ( x , z ) r 0 ( cos θ i + cos θ d ) exp ( - ikr 0 ) ds - - - ( 2 )
其中r0为输出波导到光栅上任意一点之间的距离,n为折射率,λ为波长,θi、θd分别为光栅齿面上的入射角和衍射角。
假设输出波导本征模为E2(x),经过基尔霍夫衍射计算得到的输出波导端模场为E1(x),那么耦合进入波导中的频谱响应为:
T ( u ) = | ∫ - ∞ + ∞ E 1 ( x - u ) · E 2 * ( x ) dx | 2

Claims (2)

1.优化输入端的低窜扰平顶型波分复用器件,它的结构为输入波导(1),其后接有自由扩散区(2),自由扩散区(2)后接有刻蚀凹面光栅(3),输入波导(1)侧面有输出波导阵列(4),在输入波导(1)的末端与自由扩散区(2)之间接有多模干涉部件(5),其特征在于:输入波导(1)末端与多模干涉部件(5)前端之间接有带有预展宽输入波导(6),且多模干涉部件(5)的多模干涉区(7)的边界为抛物线结构(8);预展宽输入波导(6)的初始端宽度d1大于输入波导(1)的宽度,末端宽度d2大于初始端宽度d1,长度 L 1 = 3 8 L πeff = n r 2 λ ( d 2 + C ) ( d 1 + C ) 其中包层折射率nc,芯层折射率nr,中心波长λ,Lπeff为拍长,波导等效厚度 C = λ π · ( n c n r ) 2 σ ( n r 2 - n c 2 ) ( - 1 / 2 ) , 对横电偏振σ的值是0,对横磁偏振σ的值是1。
2.用于权利要求1所述的优化输入端的低窜扰平顶型波分复用器件的构成方法,其特征在于:
1)预展宽区域尺寸的确定:
预展宽的最终目的是要出现高阶模次,改变模场分布,因此,通过BPM模拟分析这一过程里的模场分布,找到使得展宽模场边缘下降最陡的结构参数,
在任一截面上预展宽区域的宽度可以表示为: d ( z ) = d 2 + ( d 1 - d 2 ) ( 1 - z L 1 ) 其中z是传输方向,d1、d2分别为预展宽区域初始端宽度和末端宽度,L1是预展宽区域的长度,
C = λ π · ( n c n r ) 2 σ ( n r 2 - n c 2 ) ( - 1 / 2 ) , 其中对横电偏振σ的值是0,对横磁偏振σ的值是1;
则可以表示预展宽区域的设计长度为: L 1 = 3 8 L πeff = n r 2 λ ( d 2 + C ) ( d 1 + C ) 其中包层折射率nc,芯层折射率nr,中心波长λ,Lπeff为拍长;
预展宽区域的三个参数d1,d2,L1中,长度L1是初始端宽度d1和末端宽度d2的函数,因此,只需要确定预展宽区域的宽度d1和d2,先固定初始端宽度d1,找出末端宽度d2,依据计算光谱和插入损耗的方法,分别找出不同末端宽度d2和品质因数、插入损耗的关系,作出曲线,随着末端的宽度d2的增大,品质因数和插入损耗都会随之增大,能得到品质因数大,而插入损耗小的末端宽度d2
2)曲线型多模干涉区为抛物线型干涉区,其结构设计:
初始干涉区宽度d3的优化其主要衡量指标是平坦化频谱的品质因数,品质因数大小为1dB带宽与33dB带宽的比值,在给定的预展宽末端宽度d2下,改变初始干涉区宽度d3,做出d3随品质因数的变化关系曲线,能找到一个最优的初始干涉区宽度d3,使得品质因数具有最大值;
抛物线线形和干涉区长度L2的优化调节能消除带通纹波,在上述干涉区初始宽度d3已经确定的情况下,抛物线必然经过其边缘点,同时其对称轴必然为器件中心对称轴,因此,决定抛物线线形的三个参数只剩下了一个,通过改变该参数,能得到最小的带通纹波,另外能在干涉区成像点附近找到一个使带通纹波最小,1dB带宽最大的干涉区长度L2
CNB021120544A 2002-06-11 2002-06-11 优化输入端的低窜扰平顶型波分复用器件 Expired - Fee Related CN1269328C (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB021120544A CN1269328C (zh) 2002-06-11 2002-06-11 优化输入端的低窜扰平顶型波分复用器件

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB021120544A CN1269328C (zh) 2002-06-11 2002-06-11 优化输入端的低窜扰平顶型波分复用器件

Publications (2)

Publication Number Publication Date
CN1383284A CN1383284A (zh) 2002-12-04
CN1269328C true CN1269328C (zh) 2006-08-09

Family

ID=4741868

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB021120544A Expired - Fee Related CN1269328C (zh) 2002-06-11 2002-06-11 优化输入端的低窜扰平顶型波分复用器件

Country Status (1)

Country Link
CN (1) CN1269328C (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103278926A (zh) * 2013-04-22 2013-09-04 天津工业大学 一种阵列波导光栅输出光谱平坦化的方法
CN114637072B (zh) * 2022-01-18 2023-10-13 浙江大学 一种浅刻蚀多模干涉耦合的多通道平顶型波分复用接收器

Also Published As

Publication number Publication date
CN1383284A (zh) 2002-12-04

Similar Documents

Publication Publication Date Title
CN1130047C (zh) 具有平坦的光谱响应的低损耗的阵列式波导多路分用器
KR100333900B1 (ko) 모드모양 변환기, 그 제작 방법 및 이를 구비한 집적광학 소자
JP3349950B2 (ja) 波長分波回路
CN1202634C (zh) 阵列波导光栅类型的光信号多路复用器/多路分离器
WO2017075977A1 (zh) 通道损耗均匀且通带响应平坦的蚀刻衍射光栅器件
JPH116928A (ja) アレイ導波路格子型波長合分波器
JPH09189817A (ja) 光分岐器用光導波路
CN108469651B (zh) 损耗均匀的阵列波导光栅路由器
US6836591B2 (en) Arrayed waveguide grating type optical multiplexer/demultiplexer and optical waveguide circuit
CN1269328C (zh) 优化输入端的低窜扰平顶型波分复用器件
CN2552017Y (zh) 优化输入端的高平坦低窜扰波分复用器件
CN1129021C (zh) 声光器件
US6728435B2 (en) Arrayed waveguide grating type optical multiplexer/demultiplexer and optical waveguide circuit
CN1246715C (zh) 基于多个子光栅的平顶型蚀刻衍射光栅波分复用器件
CN1142650C (zh) 采用空气槽渐变输出波导的通带平坦化波分复用器件
CN1492617A (zh) 阵列式波导多路复用/解复用器
CN2585502Y (zh) 采用空气槽渐变输出波导的宽通带低串扰波分复用器件
CN1190677C (zh) 加宽阵列波导光栅通带响应的方法和设备
CN1308432A (zh) 多路复用/多路分路光线路
CN1732397A (zh) 阵列波导格子型光合分波电路
CN1236337C (zh) 三点聚焦实现的平顶型通带波分复用器件
CN2579094Y (zh) 平顶型通带波分复用器件
CN1453952A (zh) 平坦型波分复用器件
CN2588675Y (zh) 平坦化波分复用器件
CN1466300A (zh) 平顶型通带波分复用器件

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20060809