一、概述
4W1A-X射线成像实验站主要进行晶体形貌学和X射线高分辨成像实验研究,为用户提供实验所需要的同步辐射“白光”和单色X射线、旋转样品台、以及相应的探测仪器。实验站所用光源是4W1单周期扭摆磁铁,实验站与光源的距离是43米,通过4W1A光束线将同步辐射X射线从光源引入实验站。因为4W1A光束线是北京同步辐射装置最长的一条光束线,所以束线末端的X射线成像实验站可以提供准直性良好的X射线光束,为高分辨成像和相位衬度成像研究准备了较好的实验条件。除白光形貌实验和微米分辨成像实验外,实验站还建成了X射线纳米分辨三维成像显微镜,利用标准样品,纳米成像设备可获得的空间分辨率最高可达30纳米(实际样品可获得的空间分辨率随样品特性而异),可为生命科学、能源科学、材料科学、微电子产业、微加工技术和纳米科技等众多学科领域提供先进的研究手段。
二、光源
4W1A光束线的光源是单周期扭摆磁铁(wiggler),光源磁铁磁场强度是1.8T,因为该光源磁铁是电子储存环第4弧区的第一块wiggler扭摆磁铁,所以简称为4W1扭摆器。电子束团经过磁场时,受到洛仑兹力作用产生圆周加速运动,从电子轨道的切线方向发射同步辐射X射线。在同步辐射专用运行模式下(2.5GeV@250mA)4W1的光源参数(按耦合系数10%计算)为:
水平方向束团尺寸:sx = 1.23 mm
垂直方向束团尺寸:sy = 0.42 mm
水平方向角散度:sx¢ = 0.190 mrad
垂直方向角散度:sy¢ = 0.131 mrad
4W1光源理论计算的光谱曲线如图1所示(43m@25mm×15mm pinhole)。
图1为理论计算的4W1光谱曲线(42m@45mm×15mm pinhole)。
图1. 4W1A光源的理论光谱曲线和经铍窗、铝窗吸收后的光谱曲线
三、光束线
4W1扭摆器同时为X射线成像实验站和荧光实验站提供同步辐射X射线,因而从4W1扭摆器引出两条光束线,4W1A光束线为成像站输送同步辐射X射线,4W1B光束线为荧光站输送同步辐射X射线。因为电子储存环中具有极高的真空度,所以光束线由真空管道连接而成,通过铍窗和离子泵在管道中形成逐级差分真空,在保持电子储存环高真空的条件下,将同步辐射X射线从电子储存环引入相应的实验站。4W1A光束线设有4个200微米厚度的铍窗和出口处一个10微米厚度的铝窗,经过铍窗、铝窗的吸收衰减后,到达实验站的光子通量约为:7.5×1012 Phs/sec/0.1%BW(E=8keV,2.5GeV@250mA,25mm×15mm pinhole)。
4W1A光束线目前主要有三种实验模式:
1.纳米成像实验模式
纳米成像模式时光束线光学布置示意图如图2所示,光束线接收的同步辐射X射线依次经过准直镜(平面压弯镜)、双晶单色器和聚焦镜(柱面压弯镜)后汇聚到实验站精密四刀狭缝处,为后续的波带片全场成像设备提供照明X射线。
图2. 4W1A光束线纳米成像光路示意图
2.白光实验模式
通过向下移动准直镜和单色器第二晶体、向上移动聚焦镜,可以将准直镜、单色器第二晶体和聚焦镜移出光路,从而将光束线接收的宽谱X射线(白光)直接引入到实验站。在白光模式下,可以开展晶体白光形貌术、辐照测试实验及其它相关实验研究。
3.微米分辨成像实验模式
在白光实验模式下,通过实验站的双晶单色器,将束线出口的白光X射线单色化,选取特定能量的单色X射线照射在样品上,经过样品吸收后,在探测器上获得样品内部微观结构的二维投影像。由于4W1光源的高亮度、高准直性、发光点尺寸小以及超长的物源距离等优点,因此能获得高分辨率的图像,比如采用10X镜头的探测器,可获得约1mm的空间分辨率。
在微米分辨成像实验模式下,实验站还可开展相位衬度成像方法,包括衍射增强成像和同轴相位衬度成像。
四、 实验站
用户可以在实验站开展白光形貌术实验,微米分辨成像实验和纳米分辨全场成像实验以及相应的CT实验。实验站拥有以下设备:
1.白光形貌术实验设备
白光形貌实验装置;暗室;X射线专用胶片;显影定影液等。
2.X射线微米分辨成像及CT实验装置
微米分辨成像及CT实验装置由双晶单色器、五维样品平移旋转台和X射线高分辨探测器构成,如3图所示。通过双晶单色器,将束线出口的白光X射线单色化,选取特定能量的单色X射线照射在样品上,经过样品吸收后,在探测器上获得样品内部微观结构的二维投影像。利用三维图像采集软件,X射线探测器配合样品旋转台自动拍摄样品在不同角度下的二维投影像,通过重建软件,可以实现微米分辨三维成像。
实验站常用探测器为一套可切换放大倍数的镜头耦合成像探测器,通过选择不同放大倍数的物镜,可以获得相应的成像视场和分辨率进行成像,微米分辨成像的具体成像参数见表1所示。
图4. X射线纳米分辨三维成像显微镜
表1. X射线成像实验站样品处相关参数
成像模式 |
能量范围 |
光通量(photons/s) |
空间分辨率 |
光斑尺寸(H斑尺) |
微米分辨成像 |
8~26 keV |
~1012@ 8 keV |
10 μm |
13mm×13mm |
4 μm |
10mm×5mm | |||
2 μm |
5mm×2.5mm | |||
1 μm |
2mm×1mm | |||
纳米分辨成像 |
8keV |
~108@ 8 keV |
50 nm |
15μm×15μm |
100 nm |
65μm×65μm |
五、研究对象和范围
① 晶体材料缺陷形成机制研究,
② 原位观察晶体相变,
③ 原位观察材料在加载或加温等外界环境条件下的动态行为,
④ 生物医学材料、复合材料等各种材料的内部结构,
⑤ 相位衬度三维成像原理和方法研究,
⑥ 纳米分辨三维成像研究。
六、创新研究成果
形貌成像站自2001年起开展X射线相位衬度成像方法实验研究,利用相位传播成像方法获得分辨率和衬度明显优于传统吸收像的成像结果[1]。2003年10月在衍射增强成像实验方法研究上取得进展,利用两块与同步辐射偏振面垂直的晶体拍摄了清晰的昆虫图像[2,3]。2004年将衍射增强成像方法与计算机断层成像方法相结合,成功地重建了苍蝇和豚鼠耳蜗的三维像[4,5]。2005年利用针孔成像原理分析了衍射增强成像过程,提出了普遍的衍射增强成像方程,并基于衍射增强成像方法,提出了两个折射率导数的重建算法公式[6,7,8],并且结合美国人Dilmania在2000年提出的一个折射率导数重建算法公式,提出折射率梯度重建算法公式[9]。2006年清华大学工程物理系相位衬度成像组与我们合作,提出了折射率的重建算法公式[10]。2007年提出衍射增强CT成像数据采集新方法,只须将分析晶体固定于摇摆曲线腰位,就可获得重建折射率导数的全部投影数据[11],为了减小生物样品的曝光剂量,率先将迭代重建算法引入相位衬度CT[12,13]。图5至图12为基于北京同步辐射装置形貌成像站获得的部分成像结果。
(a) (b) (c)
图5. 同步辐射多种成像方法比较,(a) 豚鼠耳蜗吸收衬度成像,(b) 豚鼠耳蜗同轴相位衬度成像,(c) 豚鼠耳蜗衍射增强相位衬度成像。豚鼠耳蜗对声音非常灵敏,豚鼠耳蜗与人耳蜗相似,是耳科专家研究听力的好材料。豚鼠耳蜗样品由首都医科大学提供。
(a) (b)
图6. 鱼成像比较,(a) 吸收衬度成像,(b) 衍射增强相位衬度成像。从两幅图像的比较,可以显示相位衬度成像可以获得比吸收衬度成像高得多的衬度。
图7. 大鼠肝血管的衍射增强成像。这幅图像表明,相位衬度成像对生物组织中的管道系统比较敏感。大鼠肝样品由中日友好医院提供。
(a) (b)
图8. 利用相位衬度CT方法重建生物样品,显示了相位衬度成像方法的发展潜力。(a)豚鼠耳蜗的三维重建像,(b)苍蝇的三维重建像。
(a) (b) (c) (d)
图9. 北京同步辐射成像组提出折射率梯度重建方法。以上各图分别为环氧树脂中气泡的 (a) 投影像,(b) 折射率梯度Z轴一维分量模的重建,(c) 折射率梯度X-Y平面折射率梯度二维矢量模的重建,(d) 折射率梯度三维矢量模的重建。
图10. 甲虫折射率三维梯度模重建。甲虫样品由中科院动物所提供。
(a) (b) (c) (d) (e)
图11. 北京同步辐射成像组提出:从一套360度的投影数据重建样品多种物理量的方法。以上各图分别为塑料齿轮的 (a)吸收系数的重建、(b)折射率的重建、(c)沿轴的折射率导数的重建、(d)沿轴的折射率导数的重建、(e)折射率梯度模的重建。
图12. 迭代重建算法应用于相位衬度CT,仅用18幅投影像实现环氧树脂中气泡的折射率重建。(a)环氧树脂中气泡的左腰投影像和右腰投影像,(b) 折射率三维重建显示,(c) 折射率断层解析重建像,(d) 折射率断层迭代重建像。比较(c)和(d),可以看出迭代重建像要远好于解析重建像。
参考文献:
[1] Li Gang 2001 Presentation at 2001 User Conference of BSRF 【黎刚,BSRF 2001 用户会大会报告,大连】,科技日报,国内首次实现硬X 射线相位衬度成像,2001 年7 月12 日,第一版.
[2] 朱佩平,王寯越,袁清习,田玉莲,黄万霞,黎刚,胡天斗,姜晓明,吴自玉.“两块晶体衍射增强成像方法研究”,物理学报,54 (1): 58 (2005).
[3]Wanxia Huang, Qingxi Yuan, Peiping Zhu, Junyue Wang, Hang Shu, Bo Chen, Tiandou Hu, Ziyu Wu. “Quantitative comparison between two geometrical arrangements for Diffraction Enhanced Imaging”, Spectrochimica Acta Part B, 62: 654 (2007).
[4]袁清习,王寯越,朱佩平,黄万霞,舒航,胡天斗,吴自玉,同步辐射硬X射线衍射增强峰位成像CT研究,高能物理与核物理, 29(10): 1023 (2005).
[5] 舒航,朱佩平,王寯越,高欣,伊红霞,刘波,袁清习,黄万霞,罗述谦,高秀来,吴自玉,方守贤. “衍射增强成像方法在计算机断层成像中的应用”,物理学报,55(3): 1099 (2006).
[6] Peiping Zhu, Junyue Wang, Qingxi Yuan, Wanxia Huang, Hang Shu and Ziyu Wu,“Recent Theoretical and Experimental Progress in Hard X-ray Phase Contrast Imaging and CT at BSRF”, Proc. 8th Int. Conf. X-ray Microscopy (invited talk), IPAP Conf. Series 7: 407 (2006).
[7] P.P. Zhu, J.Y. Wang, Q.X. Yuan, W.X. Huang, H. Shu, B. Gao, T.D. Hu and Z.Y. Wu, “Computed tomography algorithm based on DEI set-up”,Appl. Phys. Lett. 87, 264101, (2005).
[8] P. Zhu, Q. Yuan, W. Huang, J. Wang, H. Shu, B. Chen, Y. Liu, E. Li, Z. Wu, Diffraction Enhanced Imaging: A Simple Model, J. Phys. D Appl. Phys., 39: 4142 (2006).
[9] Junyue Wang, Peiping Zhu, Qingxi Yuan, Wanxia Huang, Hang Shu, Bo Chen, Tiandou Hu, Ziyu Wu, “Reconstruction of the refractive index gradient by X-ray diffraction enhanced Computed Tomography”. Physics in Medicine and Biology. 51 (2006) 3391.
[10] Z.F. Huang, K.J. Kang, Z. Li, P.P. Zhu, Q.X. Yuan, W.X. Huang, J.Y. Wang, D. Zhang, A.M. Yu, “Direct computed tomographic reconstruction for directional-derivative projections of computed tomography of diffraction enhanced imaging”, Appl. Phys. Lett. 89: 041124 (2006).
[11] MinWang, Pei-Ping Zhu, Kai Zhang, Xiao-Fang Hu, Wan-Xia Huang, Yu-Wan Cen, Qing-Xi Yuan, Xiao-Liu Yu and Jun-Yue Wang, “A new method to extract refraction angle in diffraction enhanced imaging computed tomography”. J. Phys. D: Appl. Phys. 40:1 (2007).
[12] Y J Liu, P P Zhu1, B Chen, J Y Wang, Q X Yuan,W X Huang, H Shu, E R Li, X S Liu, K Zhang, H Mingand Z Y Wu. “A new iterative algorithm to reconstruct the refractive index”. Physics in Medicine and Biology. 52 (2007) L5.
[13] 张凯,朱佩平,黄万霞,袁清习,刘力,袁斌,王寯越,舒航,陈博,刘宜晋,李恩荣,刘小松,吴自玉.“代数迭代重建算法在折射衬度CT中的应用”,物理学报,57 (6): 3410 (2008).
七、联系方式
线站负责人:袁清习,010-88235990,yuanqx@ihep.ac.cn
用户联系人:黄万霞,010-88235990,huangwx@ihep.ac.cn