-
强流脉冲电子束源是高功率微波装置的核心部件之一,主要用于产生强流相对论电子束[1-3]。随着高功率微波源在军事领域的广泛应用,对微波源的需求逐渐从提高输出功率向提高输出能量发展。在高功率微波源设备运行过程中,阴极表面的发射微点在预击穿电流的作用下,温度升高,释放其表面吸附气体以及内部滞留的气体,并在其附近加剧等离子体的形成,进而产生脉冲缩短现象,导致微波脉冲提前终止[4-6]。因此,选取低放气率阴极材料是确保高功率微波源长期稳定运行的关键[7]。SiC纳米材料是一种宽带隙半导体材料,以其具有纳米级发射尖端、大长径比、高熔点、高导热率、高击穿电压、低脱气率及优良的热、化学稳定性等优点,成为优选的强流脉冲电子束及高功率微波源的阴极材料[8-10]。本研究在高真空环境中采用长脉冲激光模拟热冲击对SiC的局部受热释气特性进行研究,评估不同激光参数对SiC表面的损伤。研究方法依托于近年来本课题组发展的激光诱导解吸附−四极杆质谱(LID-QMS)[11-13]和激光诱导烧蚀−四极杆质谱(LIA-QMS)技术[14-17]。这两种技术已经被成功应用于离线检测托卡马克装置面向等离子体壁材料表面氘、氦滞留的诊断研究。激光诱导解吸附和激光诱导烧蚀四极杆质谱实验装置的工作原理就是通过脉冲激光局部加热材料使吸附在材料中的气体成分释放,并通过QMS来监测被释放的物种成分。课题组成员Lyu Yan等[12]通过使用LID-QMS技术对W基底上的聚变燃料D滞留进行了定量究。Oelmann J等[14]使用LIA-QMS技术,对多层膜太阳能电池中的氢进行了检测分析。在本论文中,通过改变激光参数(激光能量、激光光斑尺寸及激光脉冲数)对强场阴极材料SiC样品进行辐照,使用精密真空规监测激光脉冲所引起的真空腔室内气压的变化(可换算被释放的总粒子数)通过四极杆质谱仪(QMS),监测了激光脉冲热冲击的过程中真空腔室内气压和气体成分的变化。因为强场电极材料的真空释气主要成份为N2 (比O2和CO2高约3个量级),所以在本论文中使用QMS重点研究了N2离子流强度随激光参数改变的变化趋势。为了进一步分析激光烧蚀对SiC释气行为的具体影响,通过激光共聚焦显微镜测量了烧蚀区域的具体尺寸和深度。
-
图1为实验装置示意图,脉冲解吸附激光是Nd:YAG激光器(Beamtech Nimma 650),激光波长为1064 nm,激光脉宽为750 μs,激光光束的空间轮廓近似于高斯分布,时间轮廓近似于矩阵分布。激光通过位于真空室赤道面的石英窗口耦合到真空室,石英窗口距离样品表面约为20 cm,样品被置于如图1所示的样品架上。为了避免激光能量的波动,实验是在最大平均输出激光能量600±10 mJ的基础上,通过调整由半波片和偏振立方体组成的能量衰减器来改变激光能量。激光能量是由激光能量计(Newport842-PE)测量得到。使用焦距为60 cm的石英平凸透镜对激光进行聚焦,并通过移动此平凸透镜改变激光辐照在样品表面的光斑大小。通过调节样品架与安装在二维精密移动镜架上的镜面,实现激光束在样品表面的扫描。由于激光束与样品表面法线的夹角小于±1.5°,扫描过程中激光光斑大小的变化可以忽略不计。在脉冲激光辐照后,使用精密真空规(PKR251,Pfeiffer)记录真空腔室内气压的变化,使用精密四极杆质谱仪(HAL3F-RC, HIDEN)检测激光脉冲诱导释放的气体的成分。样品表面形貌的变化通过使用激光共聚焦显微镜(MICROMESURE 2; STIL. SA, France)进行测量表征。本研究采用20×20×1 mm的SiC薄片,在进行正式实验之前,将SiC薄片放入无水乙醇中使用超声波清洗50 min以去除表面的污渍。为了避免酒精所携带的碳氢化合物污染真空腔室,使用去离子水超声清洗10 min,并将其烘干,放入腔室,准备实验。
-
如图2所展示,为了更准确地测量由激光脉冲诱导的SiC释气量变化,本研究的所有实验均在静态真空条件下进行。具体操作为,在进行激光脉冲辐照前,先关闭连接到真空泵的闸板阀门。此时,气压信号会呈线性上升,这一现象主要是由于真空室内壁的脱气所导致的。图2中的蓝色曲线是将SiC样品放入测量室后,在未进行激光辐照的情况下,关闭抽气阀门之后气压的线性上升情况。通过对这一过程进行线性拟合,得到斜率约为1.12×10−7 Pa/s。考虑到实验室真空计的精度,本研究真空计的最小气压检测信号估计为2.0×10−7 Pa。这意味着,如果激光脉冲前后气压变化小于2.0×10−7 Pa,则认为没有气体被释放。而激光辐照后观察到气压骤然增加,如图3红色曲线,随后由于真空室的脱气而再次线性上升。由激光辐照而产生的实际信号通过激光辐照后的气压信号上升值减去背景信号的上升值来获得,在进行一次激光辐照后打开真空腔室的闸板阀门,使背景压力恢复至与上次测量相同(pBG=8.5×10−7 Pa)后,再次关闭阀门,重复上次操作。为使数据更加准确,每组实验重复3次。
假设在整个实验过程中,腔室内的气体处于理想状态,通过使用式(1)可算得由于单次脉冲而引起气压变化导致腔室内物质的量的增长,再通过式(2)可得到SiC表面脱附的总粒子数,由式(3)可算得由一次激光脉冲导致SiC表面单位面积释放的粒子数。
其中,Δp是脉冲激光引起的腔室内净气压增加量,单位为Pa;Δn是脉冲激光引起的SiC表面脱附的总气体的物质的量,单位为mol;V是腔室的体积,本论文中V = 0.019441 m3;R是理想气体常数,其值为 8.314 J/(mol·K);T表示气体的温度,本论文T = 300 K;Na为阿伏伽德罗常数,其值为Na = 6.02214076×1023 mol−1;ΔN是单次激光释放的总粒子数;ΔNb为单位面积释放的粒子数,单位为cm−2;A为激光辐照在SiC表面的光斑面积,单位为cm2。
-
图3展示了在基压pBG = 8.5×10−7 Pa,恒定光斑直径d = 1 mm 的条件下,通过使用不同激光能量辐照SiC表面,进而改变辐照在表面的激光能量密度。实验测量的单位面积释放的粒子数随激光能量密度的变化如图3(a)所示,从图中可以看到,当激光能量密度在低于约6.46 J/cm2时,真空规未检测到任何气压信号,随着激光能量密度的增大,SiC表面开始有粒子被释放,但其增长速度缓慢,此时所检测到的粒子数信号是由激光解吸SiC表面所吸附的气体信号。当激光能量密度约为22.49 J/cm2时,其释放的粒子数明显增多,此时单位面积释放的粒子数约为1.43×1017 cm−2。通过激光共聚焦显微镜对样品表面进行扫描,发现SiC表面在激光能量密度大于等于22.49 J/cm2时下开始出现轻微损伤(定义该值为SiC材料的激光烧蚀/损伤阈值)。实验表明,SiC材料的激光烧蚀斑的形状是不均匀的,如图3(b)所示,且损伤的程度随激光能量密度的增加而加剧,相应的单位面积释放粒子数也随之增强。当出现SiC表面损伤之后,所测得的释气量数值不仅反映了SiC在体相中滞留的气体和表面吸附的气体的释放,还包括由于表面烧蚀而产生的团簇粒子的贡献。图3(c)给出了采用四级杆质谱仪记录的SiC表面单位面积产生的m/q = 28 amu/e离子(N2+)流与激光能量密度的关系。与图3(a)SiC表面单位面积释放的粒子数相比,这两者的变化趋势大致相似。当激光能量高于29.36 J/cm2时,单位面积产生m/q = 28 amu/e离子(N2+)流的变化趋势开始减缓。这表明在较高的激光能量密度下,可能有更多种类的气体分子被释放,与单纯的m/q = 28 amu/e离子(N2+)释放相比,这些其他气体的释放对总体气压升高的影响更为显著。以上结果说明,适当的激光能量密度可以有效诱导SiC表面的气体释放,但过高的能量密度可能导致材料损伤。
-
在固定的激光能量条件下(能量约为139 mJ),通过移动聚焦透镜的位置来改变焦点到样品表面的距离(Δz),变化的范围是从−12 cm到0 cm,研究其对SiC单位面积释放粒子数及质谱信号强度的影响。Δz的改变直接影响激光辐照在样品表面的光斑尺寸,进而影响作用于SiC表面的激光能量密度。图4 (a) 为在给定激光能量条件下,SiC表面的光斑尺寸和激光能量密度随Δz变化的规律,Δz在−12 cm ~ 0 cm之间变化时,光斑直径d的变化范围为1.495 mm ~ 0.637 mm,激光能量密度从7.92 J/cm2增加到43.61 J/cm2。如图4 (b) 给出了SiC表面单位面积释放的气体粒子数随激光光斑尺寸的变化趋势,从图中可以明显观察到,随着光斑尺寸的增大,SiC表面单位面积释放的粒子数先增大后减小,在光斑直径为0.682 mm时达到最大值为6.16×1018 cm−2。而后随着光斑尺寸的进一步增大,单位面积释放的粒子数开始下降。这表明由激光辐照释放的单位面积粒子数并不是在光斑尺寸最小处(对应激光能量密度最大值),而是存在一个特定的光斑尺寸d = 0.682 mm,此时激光能量密度为38.04 J/cm2。值得注意的是,单位面积产生m/q = 28 amu/e离子(N2+)流的变化趋势与单位面积释放的粒子数有所不同,如图4 (c) 所示,在光斑直径d = 0.645 mm时,激光能量密度为42.63 J/cm2时释放的粒子数达到最大值。
-
在此研究中,作者对SiC在激光能量密度为25.76 J/cm2的条件下施加不同数量的激光脉冲,探讨了连续热冲击(即激光脉冲数)对其粒子数释放量和质谱信号强度的影响,考虑到在研究中每一次激光脉冲都会导致烧蚀坑大小的变化,故仅计算了释放的总粒子数。实验中脉冲间隔设定为1分钟,以确保能够准确捕捉每一次激光脉冲后的信号变化。结果如图5所示,随着脉冲数的增加,烧蚀深度、释放的总粒子数和m/q = 28 amu/e离子(N2+)流都在增加,且增加速率有减缓趋势。具体如表1所示,该表展示了在不同脉冲区间里平均单次脉冲所引起的烧蚀坑深度和释放的粒子数的变化情况, 可以看到第一个脉冲所引起的烧蚀深度以及粒子数的增加是最强的,随着激光脉冲数的增加,每次脉冲所引起的深度变化和释放粒子数都在减弱。这种减弱可能是因为在非连续激光脉冲中,烧蚀深度和粒子释放量的增加并非仅由激光能量的直接累积所导致,而是与材料表面形态的演变紧密相关。随着激光的重复辐照,辐照区域的材料会反复融化和冷却,进而形成了新的表面形貌,降低表面的粗糙度,影响了激光能量的吸收效率。另外,初期激光脉冲可能已经移除了大量吸附能(Ea)较低的气体物种分子,使得在后续脉冲中释放的粒子数减少。
通过改变激光的参数,大致估算了在不同条件下单位面积释放的粒子数和释放的总粒子数。在该实验中,为了尽可能准确地测量SiC的释放粒子数,将总释放粒子数分为了两部分,即SiC自身放气及腔室的本底释气和漏气引起的粒子数增多与样品受到激光脉冲加热后的释放粒子数[18]。而对于SiC来说,其释放粒子数与材料微结构及制备处理工艺紧密相关。还有一些其他外界因素也会影响释放粒子数的多少,比如SiC样品安装时,暴露大气中,也将引起粒子数的增高;不同的表面粗糙度也会影响材料对激光能量的吸收,进而影响所测材料的释放粒子数。对于同一块样品,表面吸附的物质(通常是水汽)与周围空间的压力成一定的比例[19],因此在封闭的真空室中,当泵组向真空腔室外抽气时,SiC吸附的水汽也会减少。而实验前的处理过程(比如烘烤和清洗)的差别,也会影响碳SiC所释放的粒子数特征。
-
本文在超高真空(8.5×10−7 Pa)条件下,通过采用长脉冲激光辐照模拟热冲击的方式,研究了SiC在不同激光能量密度、不同激光光斑尺寸和不同脉冲数条件下的释气特性,研究表明,在不同激光能量密度条件下,在达到其释放粒子数的能量阈值约6.46 J/cm2后,随着激光能量密度的增加,SiC的表面单位面积脱附的粒子数也增加,当能量密度到达到22.49 J/cm2时,样品表面发生了轻微的损伤,SiC单位面积脱附的粒子数迅速增加,约1.43×1017 cm−2,因此,SiC烧蚀损伤阈值为22.49 J/cm2。在不同激光光斑尺寸条件下,存在一个最优的光斑直径d = 0.68 mm,此时脉冲激光能量密度为38.04 J/cm2,单位面积释放的粒子数约为6.16×1018 cm−2。当激光光斑直径为最小值d = 0.64 mm时,此时辐照在SiC表面的激光能量密度达到最大值43.60 J/cm2,单位面积释放的粒子数约为4.36×1018 cm−2。在激光能量密度为25.76 J/cm2下,对样品表面进行不同脉冲数的激光辐照,结果表明,随着脉冲数的增加,释放的粒子数逐渐增加,但单次脉冲引起的粒子释放和烧蚀深度逐渐减小,这可能是由于前面脉冲已经移除了大量吸附能(Ea)较低的气体物种分子,后续脉冲释放的粒子为具有高吸附能的物种(深度滞留的气体),其数目较少,另外在多次激光辐照下,样品表面产生融化,粗糙度变小,影响了样品的激光的吸收效率,使得后续脉冲粒子释放效率下降。
长脉冲激光辐照强场阴极材料碳化硅释气特性研究
The Outgassing Characteristics of Silicon Carbide as an Intense-Field Cathode Material Under Long-Pulse Laser Irradiation
-
摘要: 文章在超高真空(8.5×10−7 Pa)条件下采用波长1064 nm,脉宽750 μs的激光模拟脉冲热冲击,对碳化硅(SiC)材料的释气特性进行了实验研究。主要探讨了激光能量、激光光斑大小和激光脉冲数对SiC释气特性的影响。通过全量程真空规和四极杆质谱仪(QMS)监测了激光脉冲热冲击过程中真空腔室内气压和质谱检测的气体成分(N2+)的变化。研究表明,增加激光能量密度为6.46 J/cm2时,SiC表面开始有气体分子(粒子)释放;随着激光能量密度的进一步增加,释气量也进一步增加。当激光能量密度达到22.49 J/cm2时,SiC表面释放的气体分子(粒子)数明显增多。研究还表明,适当的激光能量密度可以有效诱导SiC材料表面的气体释放,但过高的能量密度可能导致材料损伤,实验测定的SiC损伤阈值为22.49 J/cm2。在激光能量为139 mJ的条件下,通过调节聚焦透镜焦点到SiC表面的距离改变光斑直径,单位面积释放的粒子数在光斑直径为0.68 mm处达到最大值,此时的激光能量密度为38.04 J/cm2,并不是最大的激光能量密度。这说明SiC单位面积释放的粒子数不仅与激光能量密度相关,还与激光辐照区域面积大小相关。在激光能量密度为25.76 J/cm2的条件下,随着脉冲数的增加,SiC的释气量逐渐增加,但每次脉冲引起的粒子释放量和烧蚀坑深度随脉冲数逐渐减小。其原因可能是由于激光重复辐照导致材料表面粗糙度降低,从而降低材料对激光能量的吸收效率。Abstract: In this work, a Nd:YAG laser operating at a wavelength of 1064 nm with a pulse width of 750 μs was utilized to simulate pulsed thermal shock on silicon carbide (SiC) material, aiming to investigate its outgassing characteristics under ultra-high vacuum conditions (pressure of 8.5×10−7 Pa). The study primarily delved into the influence of laser energy, laser spot size, and the number of laser pulses on SiC gas release characteristics. Throughout the laser pulse thermal shock process, chamber pressure variations and outing gas composition (N2+) changes were monitored by using a full-range vacuum gauge and a quadrupole mass spectrometer (QMS). The results showed that gas molecules (particles) began to release when the laser energy density was increased to 6.46 J/cm2, with a noticeable increase when the laser energy density further increased. Additionally, it was found that an appropriate laser energy density could effectively induce gas release from the SiC material surface, while excessively high energy density could lead to material damage, with the experimentally determined damage threshold of SiC being 22.49 J/cm2. Moreover, the laser spot size on the SiC surface was adjusted by moving the distance from the to the SiC surface with a fixed laser energy of 139 mJ to investigate the gas releasing characteristics. The maximum number of particles released per unit area is at a beam diameter of 0.68 mm, corresponding to a laser energy density of 38.04 J/cm2, which was not the maximum laser energy density. This indicates that the number of particles released per unit area of SiC is not only related to the laser energy density but also to the size of the laser irradiation area. Furthermore, under a laser energy density of 25.76 J/cm2, an increase in the number of pulses was found to incrementally elevate the gas release quantity of SiC. However, the quantity of particles released per pulse and the depth of ablation pits induced by each pulse demonstrated a gradual decline with increasing pulse count. This may be due to the decrease in material surface roughness caused by repeated laser irradiation, thereby reducing the material's absorption efficiency of laser energy.
-
Key words:
- Laser thermal shock /
- Outgassing /
- SiC /
- High-power microwave source .
-
-
图 3 SiC材料在固定光斑尺寸为1 mm时。(a)单位面积释放粒子数随激光能量密度的变化趋势,(b)在激光能量密度为22.49 J/cm2时SiC表面烧蚀2D分布,(c)单位面积产生m/q = 28 amu/e离子(N2+)流随激光能量密度的变化趋势
Figure 3. When the fixed spot size of laser in SiC material is 1 mm. (a) Areal density of outgassing gas as a function of laser fluence, (b) 2D distribution of SiC surface ablation at a laser fluence of 22.49 J/cm2, (c) areal density of ion current with m/q = 28 amu/e (N2+) as a function of laser fluence
图 4 SiC材料在激光能量为139 mJ时。(a) 光斑尺寸和激光能量密度随焦点到SiC表面距离的变化趋势, (b) SiC表面单位面积释放粒子数随光斑尺寸的变化趋势,(c) 单位面积产生m/q=28 amu/e离子(N2+)流随光斑尺寸的变化趋势
Figure 4. When the fixed energy of laser in SiC material is 139 mJ. (a) Trends in spot size and laser fluence with changes in the distance from the focus to the SiC surface, (b) areal density of outgassing gas as a function of spot size, (c) areal density of ion current with m/q = 28 amu/e (N2+) as a function of spot size
图 5 SiC材料在激光能量密度为25.76 J/cm2时。(a) SiC释放总粒子数与烧蚀深度随脉冲个数的变化趋势,(b) m/q = 28 amu/e离子(N2+)流与烧蚀深度随脉冲个数的变化趋势
Figure 5. When the fixed laser fluence in SiC material is 25.76 J/cm2. (a) Total number of particles as a function of pules number, (b) ion current with m/q = 28 amu/e (N2+) as a function of pules number
表 1 不同脉冲区间下平均单次激光脉冲对SiC深度变化及粒子释放数
Table 1. Average depth change and number of particles released from SiC per laser pulse across different pulse intervals
脉冲区间 平均单次脉冲深度变化/μm 平均单次脉冲释放粒子数 1 2.14 3.65×1015 2−3 0.781 2.14×1015 4−5 0.776 1.30×1015 6−9 0.396 1.21×1015 21−30 0.230 4.70×1014 -
[1] Schamiloglu E. High power microwave sources and applications; proceedings of the 2004 IEEE MTT-S International Microwave Symposium Digest (IEEE Cat No 04CH37535), F, 2004 [C]. IEEE [2] Xun T, Zhao Y, Yang H, et al. Developments of pulsed electron beam sources for high-power microwave applications[J]. IEEE Access,2020,8:101351−101358 doi: 10.1109/ACCESS.2020.2998088 [3] 荀涛, 杨汉武, 张军, 等. 高性能强流脉冲电子束源关键技术研究[J]. 强激光与粒子束,2020,32(2):7 (in Chinese) Xun T, Yang H W, Zhang J, et al. Development of high performance, high-current pulsed electron beam sources[J]. High power laser and particle beam,2020,32(2):7 [4] Agee F. Evolution of pulse shortening research in narrow band, high power microwave sources[J]. IEEE transactions on plasma science,1998,26(3):235−245 doi: 10.1109/27.700749 [5] 李姝敏, 李永东, 刘震. 相对论返波管中击穿现象粒子模拟[J]. 强激光与粒子束,2017,29(6):16−19 (in Chinese) Li S M, Li Y D, Liu Z. Particle-in-cell simulation of field breakdown in a relativistic backward wave oscillator[J]. High power laser and particle beam,2017,29(6):16−19 [6] 张晓微, 肖仁珍, 陈昌华, 等. 一种提高相对论返波管功率容量的方法[J]. 强激光与粒子束,2011,23(11):3069−3072. (in chinese) doi: 10.3788/HPLPB20112311.3069 Zhang X, Xiao R, Chen C, et al. A method for improving power capability of relativistic backward wave oscillator[J]. High power laser and particle beam,2011,23(11):3069−3072 doi: 10.3788/HPLPB20112311.3069 [7] Krasik Y, Gleizer Z, Yarmolich D, et al. Characterization of the plasma on dielectric fiber (velvet) cathodes[J]. Journal of Applied Physics,2005,98(9):093308 doi: 10.1063/1.2126788 [8] Jang C, Kim T, Lee S, et al. Low-resistance ohmic contacts to SiC nanowires and their applications to field-effect transistors[J]. Nanotechnology,2008,19(34):345203 doi: 10.1088/0957-4484/19/34/345203 [9] Xun T, Zhao X, Li G, et al. High-current, pulsed electron beam sources with SiC nanowire cathodes[J]. Vacuum,2016,125:81−84 doi: 10.1016/j.vacuum.2015.11.023 [10] Wang L, Li C, Yang Y, et al. Large-scale growth of well-aligned SiC tower-like nanowire arrays and their field emission properties[J]. ACS Appl Mater Interfaces,2015,7(1):526−533 doi: 10.1021/am506678x [11] Zlobinski M, Philipps V, Schweer B, et al. Laser induced desorption as tritium retention diagnostic method in ITER[J]. Fusion Engineering and Design,2011,86(6-8):1332−1335 doi: 10.1016/j.fusengdes.2011.02.030 [12] Lyu Y, Li C, Wu D, et al. Characterization on deuterium retention in tungsten target using spatially resolved laser induced desorption-quadrupole mass spectroscopy[J]. Physica Scripta,2021,96(12):124040 doi: 10.1088/1402-4896/ac2c24 [13] Yehia-Alexe S, Groza A, Serbanescu M, et al. Considerations on hydrogen isotopes release from thin films by laser induced ablation and laser induced desorption techniques[J]. Spectrochimica Acta Part B-Atomic Spectroscopy,2023,208:106774 doi: 10.1016/j.sab.2023.106774 [14] Oelmann J, Gierse N, Li C, et al. Depth-resolved sample composition analysis using laser-induced ablation-quadrupole mass spectrometry and laser-induced breakdown spectroscopy[J]. Spectrochimica Acta Part B: Atomic Spectroscopy,2018,144:38−45 doi: 10.1016/j.sab.2018.03.009 [15] Oelmann J, Li C, Brezinsek S, et al. Depth resolved analysis of hydrogen in W7-X graphite components using laser-induced ablation-quadrupole mass spectrometry (LIA-QMS)[J]. Nuclear Materials and Energy,2019,18:153−158 doi: 10.1016/j.nme.2018.12.019 [16] Zhao D, Yi R, Oelmann J, et al. Ex situ analysis of W7-X divertor plasma-facing components by picosecond laser diagnostics[J]. Physica Scripta,2020(T171):014018 [17] Lyu Y, He Z, Wang X, et al. Characterization of helium retention in the inhomogeneous co-deposited layers using a long pulse laser induced ablation-quadrupole mass spectroscopy[J]. Nuclear Materials and Energy,2022,33:101268 doi: 10.1016/j.nme.2022.101268 [18] Li A, Fan Y, Qian B, et al. Outgassing rate analysis of a velvet cathode and a carbon fiber cathode[J]. Journal of Applied Physics,2017,122(18):185901 doi: 10.1063/1.4996649 [19] Shen Y, Zhang H, Xia L, et al. Vacuum outgassing behavior of carbon nanotube cathode with high-intensity pulsed electron emission[J]. Plasma Science and Technology,2015,17(2):129−133 doi: 10.1088/1009-0630/17/2/06 -