《角分辨光电子能谱数据的快速拼接算法设计》由会员分享，可在线阅读，更多相关《角分辨光电子能谱数据的快速拼接算法设计（7页珍藏版）》请在金锄头文库上搜索。
1、角分辨光电子能谱数据的快速拼接算法设计 贾小文 陆军军事交通学院基础部 摘 要： 根据角分辨光电子能谱的实验特点, 提出快速拼接算法:首先生成均匀分布的动量点, 然后逆用角度-动量转换公式获取对应的角度值, 最后在角度空间进行插值.由于避开了普通插值的三角剖分和几何排序等复杂耗时过程, 快速拼接算法能将运算速度提高 2 个数量级, 插值的效果和传统插值方法相当.结合并行运算, 还可进一步提高数据处理速度.关键词： 角分辨光电子能谱; 数据拼接; 网格化数据; 作者简介：贾小文 (1983-) , 男, 陕西神木人, 陆军军事交通学院基础部讲师, 博士, 研究方向为凝聚态超导物理、实验数据分析.
2、收稿日期：2017-06-27Fast stitching algorithm for angle-resolved photoemission spectroscopyJIA Xiao-wen General Course Department, Army Military Transportation University; Abstract： Based on the characteristics of angle-resolved photoemission spectroscopy, a fast stitching algorithm was proposed.Firstly,
3、uniformly distributed points in momentum space were generated, then angles for each momentum point were calculated using the angle-momentum conversion formula inversely, finally the values for the points in momentum space by interpolation in angle space were obtained.Because the fast stitching algor
4、ithm did not involve complicated time-consuming processes such as triangulation and geometric ordering, the computation speed increased 2 orders of magnitude compared to that of the traditional algorithm, and the interpolation quality did not degrade.Combined with parallel computing, computation spe
5、ed could be further improved.Keyword： angle-resolved photoemission spectroscopy; data stitching; gridded matrix data; Received： 2017-06-271887 年, 赫兹发现当入射光照射到金属表面时, 会有电子从金属表面逸出.1905 年, 爱因斯坦提出了“光电子”的概念, 从量子力学的角度完美地解释了光电效应.由于光电子携带了材料内部电子结构的信息, 因此可以通过探测光电子来研究材料的电子结构, 角分辨光电子能谱就是该原理的具体应用1-3.实验上, 通过电子能量分析器
6、探测出射的光电子, 得到数据为某一方向立体角内电子的统计数.从实验得到的光电子统计数是按照角度分布的.为了方便研究, 需要将角度转换为动量 (实际上是动量在 kx-ky面上的投影) , 经过转换后的数据为光电子统计数随动量 (k x, ky) 的分布.在角度空间分布于一条线上的光电子 (目前主流探测器都是这种模式) , 转换到动量空间的分布是一条弯曲的曲线.因此要获得某个动量平面内所有光电子的分布, 就要求对各个不同角度下 (通常这个角度就是测量样品的取向) 的测量数据进行拼接4.在拼接数据时, 传统的方法先由角度换算为动量, 然后在动量空间进行均匀插值, 这种方法存在耗时较大的缺点.由于角分
7、辨光电子能谱实验要求极度清洁的样品表面, 如果耗时过长, 样品清洁度将会下降, 从而导致实验失败.因此无论是从数据处理效率还是从开展实验的角度来看, 都需要对传统的方法进行改进.本文充分考虑了角分辨光电子能谱实验数据的相邻倾角数据基本类似的特点, 提出了快速数据拼接算法, 该算法首先生成均匀分布的动量点, 然后逆用角度动量转换公式获取对应的角度值, 最后在角度空间进行插值.由于避开了普通插值的三角剖分和几何排序等复杂过程, 可将运算速度提高 2 个数量级.1 角分辨光电子能谱实验原理及数据的拼接角分辨光电子能谱实验的基本原理是光电效应, 如图 1 所示.探测器收集电子的窗口是 1 条狭缝, 1
8、 次收集 1 个狭长立体角空间的所有电子.原始数据为网格化的矩阵, 横坐标是动能, 纵坐标是角度, 数据表示光电子的谱重, 如图 2 所示.处理数据时, 需要将角度转换成动量, 转换的公式取决于样品与电子探测器相对方位的布局和定义.以文献4为例, 分析器相对于样品的角度定义取图 3 所示定义:将与狭缝平行的方向定义为偏转角度 , 与狭缝垂直的方向定义为倾角 , 样品转过的角度定义为旋转角 4-5.动量 (k x和 ky面的投影) 和实验测量角度的转换公式为式 (1) 中, k 0表示出射电子的动量, 由式 (2) 决定:从式 (1) 中可以看到, 探测器探测到的动量实际上是水平方向的动量, 因
9、此角分辨光电子能谱经常用于研究电子结构二维性较好的材料.图 1 角分辨光电子能谱实验基本原理示意图 下载原图图 2 典型角分辨光电子能谱原始数据图谱 下载原图图 3 实验中样品、探测器狭缝取向与动量换算几何关系 下载原图为了得到不同动量区域的光电子信息, 需要不断地改变狭缝相对于样品的倾角, 并将这些数据拼接在一起.在拼接时, 图 2 中的原始数据在角度空间是均匀分布的, 但是利用式 (1) 转换到动量空间 (其实是 kz=C 的某个平面) 后, (k x, ky) 就不再是均匀的网格分布了, 有些动量点有数据, 有些动量点则没有数据, 如图 4 所示.图 4 (a) 所示为铁基超导体 FeS
10、e0.3Te0.7在 20K 下费米能级处谱重随动量的分布 (光子能量 22.218eV, 光电子能谱仪为 Scienta 公司的 R4000) , 可以看到图 4 (a) 的数据分布并不均匀, 有些点有数据, 有些点没有数据 (相当于谱重等于 0) , 是非均匀分布的数据图.出现这种情况的原因有:1) 角度到动量的转换算法式 (1) 是非线性的, 连续均匀的角度分布换算到动量空间就不再是连续均匀分布了;2) 探测器在倾角方向分辨率太大 (倾角人为设定, 一般是 12, 远大于探测器狭缝方向的角度分辨率) , 数据呈现明显的条状分布.这种非均匀分布的数据给后续的分析带来不便, 实践中需要将图
11、4 (a) 的结果通过插值拼接成分布连续、光滑的数据, 如图 4 (b) 所示.拼接的过程中需要解决数据不均匀分布的问题, 即非均匀数据网格化.图 4 费米能级处谱重随动量的分布 下载原图2 非均匀数据网格化非均匀数据网格化是常见的数据处理方法, 有非常成熟的数值算法, 其过程是先对自变量空间进行三角剖分, 再利用三角剖分结果重新抽样.很多软件 (如Igor, Matlab 等) 都有相应的函数完成上述运算.这里介绍 Igor 下的做法6-7.首先创建 N3 维的矩阵 M, 对每个倾角下的谱应用式 (1) 计算 (k x, ky) , 分别记录在 M 的第 1 列和第 2 列, M 的第 3
12、列存放对应的测量值.然后利用ImageInterpolate 命令进行插值运算:以上做法是目前角分辨光电子能谱实验中数据拼接的传统做法, 图 4 (b) 的结果就是利用该方法计算得到的, 其缺点是运算时间过长, 通常需要 1min 左右 (取决于谱的数量) , 甚至更长.耗时来源于三角剖分过程及以三角剖分为条件的插值几何排序过程.3 快速数据拼接算法三角剖分法适用的范围广, 但某些应用场合效率较低.实际上, 角分辨光电子能谱实验数据有自己的特点:数据的自变量并非完全随机分布, 相邻狭缝之间除了倾角不同, 其他的实验条件完全相同.根据这个特点构建了快速拼接算法:逆用式 (1) 从动量值反过来推算
13、角度, 利用角度值在原始数据中插值, 将插值作为该动量处的结果.快速拼接算法的优点是:1) 在动量空间, 可以预先定义均匀分布的网格数据, 省去了重新抽样的过程;2) 在角度空间, 由于数据的分布是完全均匀的, 可直接通过线性插值来获取插值结果, 省去了三角剖分的计算过程.快速拼接算法的复杂度为 O (n) , 因此运算速度非常快.实验结果表明:快速拼接算法要比传统算法运算速度提高 2 个数量级.快速拼接算法的详细描述如下:1) 创建矩阵 M (kx, ky) .2) 取出 (k x, ky) , 逆用式 (1) 计算 (, ) .3) 对原始测量数据通过偏转角、倾角进行排序, 找出相邻倾角
14、1和 2, 使得 1 2, min max. min和 max为原始谱最小角和最大角加上偏转角 .4) 设 1对应的原始数据是 data1, 2对应的原始数据是 data2, 利用双线性插值分别获取 2 个值 z1=data1 (0) () , z 2=data2 (0) () .这里 0 表示费米能级.5) 通过线性插值得到 (, ) 处数据:6) 将 M 中 (k x, ky) 处相应的值设置为 z.下面是算法具体实现 (Igor 代码) :4 实验结果以图 4 所示数据拼接为例, 对上述算法进行测试.测试数据一共有 77 张谱, 每张谱在角度方向的数据长度为 360.图 5 是快速拼接算
15、法给出的插值结果.可以看到图 5 (a) 和图 4 (b) 的结果几乎没有区别.为了进一步确认快速拼接算法的可靠程度, 图 5 (b) 给出了沿不同方向切面ProfileLine, 方向如图 5 (a) 中红色虚线箭头所指的动量分布曲线 (Momentum distribution curves, MDC) , 为了方便比较, 曲线都进行了归一化.可以看出快速拼接算法在曲线数据特征如峰位、半高全宽等方面与传统算法是一致的.尽管如此, 如果仔细观察, 传统算法和快速拼接算法还是存在一些微小的差异.出现这种现象的原因主要是插值机制的区别.Igor 下传统算法使用的泰森多边形插值 (Voronoi
16、interpolation) 是平滑插值方法, 而本文中的快速拼接算法则完全采用线性插值.另外, 由于计算过程不同, 计算机的舍入误差也会给插值结果带来一些差别.图 5 快速拼接算法结果 下载原图表 1 给出了 Igor Pro (Ver 6.37) 下二者所用时间的对比.从表 1 可以看出快速拼接算法的运算速度要比传统算法提高了近 200 倍.利用多线程技术, 还可以将快速拼接算法并行分配到多个处理器 (多线程) , 能够进一步将运算速度提高到 300 多倍.多线程实现的方法如下 (Igor 代码) :1) 将渲染函数 render_FS 声明为线程安全函数.2) 利用 multithread 指令调用渲染函数对 M_FS 赋值.快速拼接算法 (多线程) 和传统三角剖分算法运行消耗时间比较如表 1 所示.表 1 快速拼接算法、多线程和传统三角剖分算法运行消耗时间比较 下载原表 5 讨论原始数据谱的数量
