上海光源二期工程的垂直发射度研究

李昌亮 田顺强 张文志



上海光源二期工程的垂直发射度研究

李昌亮1,2田顺强1张文志1

1（中国科学院上海应用物理研究所 嘉定园区 上海 201800）2（中国科学院大学 北京 100049）

上海光源二期工程正在设计一个新的磁聚焦结构，它包含双斜插入件、超导二极铁以及超导扭摆器等新的元件。光谱亮度是表征同步辐射光源性能的一个重要参数，通过降低垂直发射度的方式来提高亮度是一种简便可行的方法。首先分别模拟和计算了各类磁铁准直误差对垂直发射度的影响，然后采用奇值分解法校正了误差引起的闭轨畸变，最后使用60个斜四极铁来校正垂直发射度。结果表明，在上海光源二期工程的磁铁准直公差与目前正在运行的上海光源相同的条件下，六极铁垂直准直误差对垂直发射度的影响最大；用60个斜四极铁来校正垂直发射度，能使耦合度保持在0.5%以下。

上海光源二期工程，磁铁准直误差，垂直发射度，耦合校正

上海同步辐射光源(Shanghai Synchrotron Radiation Facility, SSRF)是一台标准能量为3.5 GeV的第三代同步辐射储存环光源。上海光源储存环由20个双二极铁消色差单元(Double Bend Achromatic, DBA)组成，它包含4个超周期结构，每个单元包含两个二极铁、10个四极铁和7个六极铁。储存环总长为432 m，其自然发射度为3.9 nm·rad[1‒2]。全环共有140个束流位置监测器(Beam Position Monitor, BPM)，水平与垂直各80个校正铁。近年来同步辐射装置吸引了越来越多的实验用户，这极大地促进了上海光源二期线站工程的升级改造。上海光源二期工程目前正处于设计阶段，它的主要任务就是建造超过16个线站供用户使用，其中包含4个低温永磁波荡器、一个超导扭摆器和两个超导二极铁单

元[3‒6]。升级改造的过程中，大部分二极铁、四极铁、六极铁、校正铁和BPM的位置都保留不变，改造后储存环的主要参数见表1。

表1 上海光源二期工程主要参数 Table 1 Main parameters in SSRF Phase II project.

光谱亮度是表征同步辐射光源性能的一个重要参数，高亮度的同步辐射光在照射样品时能保持极高的辐射能量密度，以很高的信噪比在很短时间内获得满意的实验结果，因而同步辐射实验用户对同步辐射光的亮度极其地关心。亮度的表达式可以简化为：

式中：beam为束流流强；ε、ε表示水平和垂直发射度；为耦合度，它表示垂直发射度与水平发射度的比值，其中水平发射度由磁聚焦结构决定，束流流强受限于束流集体效应。因此，通过降低垂直发射度的方式来提高光源亮度，进而改善光源的性能是一种简便可行的方法。尽管在降低垂直发射度的同时，会引起束流寿命的降低，但采用恒流(Top-up)的注入方式可以有效地解决这一不利影响。此外，在降低垂直发射度的同时，也获得了其他方面的优点，例增加了横向相干性、减小了垂直束斑尺寸。

目前国际上同类光源，如瑞士光源(Swiss Light Source, SLS)[7]、澳大利亚光源(Australian Synchrotron Project, ASP)[8]等，为提高各自光源的性能，都积极地开展了降低束流垂直发射度的理论和实验研究。上海光源已高性能、高稳定地运行了6年多，在调束和运行阶段，光源物理组在降低束流垂直发射度上做了大量详细深入的研究[9]。

本文旨在研究上海光源二期工程的垂直发射度，主要分析了各磁铁准直误差对垂直发射度的影响，并且使用斜四极铁来校正了这些误差所引起的垂直发射度，使上海光源二期在以后正常运行时能保持与目前上海光源运行时相同的耦合度。

1 垂直发射度的来源

垂直发射度来源于三个方面的贡献：1) 量子极限；2) 垂直色散，此为纵向与垂直方向的耦合；3) 横向耦合，此为水平方向与垂直方向的耦合。以上三种来源中，第一种来源是不可消除的，它是由于同步辐射光并不是准确沿切线方向辐射，而是沿着1/的张角辐射而产生的。对于某一储存环来讲，量子极限所决定的垂直发射度为一不变的常数，并且一般小于1×10‒3nm rad，其计算公式为[10]：

式中：q=3.84×10‒13m；β为垂直方向的beta函数；J为垂直方向的阻尼因子；()为二极铁的曲率半径。利用式(2)计算得到上海光源二期磁聚焦结构中量子极限所决定的垂直发射度为1.4×10‒4nm·rad。

垂直发射度主要来源于第二和第三种贡献，即垂直色散和横向耦合，而垂直色散和横向耦合又是由磁铁的准直误差所产生的。垂直色散来源于：二极铁的转角误差、四极铁的转角误差、四极铁中的垂直准直误差及六极铁的垂直准直误差。横向耦合来源于四极铁转角误差和六极铁的垂直准直误差。

二极铁有转角时，此时磁场不再是单一的垂直场，而是产生了水平方向的磁场分量，进而产生了闭轨畸变(Closed Orbit Distortion, COD)，形成垂直色散。四极铁垂直方向上的非准直，在参考轨道上（理想的闭轨）上将产生水平磁场，进而也产生闭轨畸变。垂直方向具有色散之后，当辐射发出光子时，粒子纵向的能量发生变化，引起了闭轨的变化，进而引起betatron振荡的振幅与相位发生变化，粒子在一个新的闭轨上振荡，辐射引起了betatron振荡振幅均方根值的增大，从而增加垂直发射度。横向耦合有两种效应：一种是把水平色散耦合到垂直平面上去，增加垂直色散；另一种是耦合水平平面上的betatron振荡到垂直平面上去，两种效应都增加了垂直发射度。

2 磁铁准直误差对垂直发射度的影响

分别模拟二极铁转角误差ΔΘ、四极铁转角误差ΔΘ、四极铁中的垂直准直误差ΔY以及六极铁的垂直准直误差ΔY对垂直方向闭轨畸变的影响，在模拟某一类准直误差对垂直方向闭轨畸变的影响时，确保其他磁铁没有准直误差，这样做的目的在于单独分析每种误差的影响。使用的模拟软件是加速器工具箱(Accelerator Toolbox, AT)[11]，AT是基于MATLAB运行环境下的模拟软件，其本质为MATLAB的一个工具箱，是由斯坦福直线加速器中心(Stanford Linear Accelerator Center, SLAC)为斯坦福大学正负电子非对称环(Stanford Positron Electron Asymmetric Ring, SPEAR3)的设计和使用研发的。

首先模拟二极铁转角对垂直方向闭轨畸变的影响，分别设置二极铁转角误差为0.04 mrad、0.08mrad、0.12 mrad、0.16 mrad、0.2 mrad，每一个误差值处都选择了200个随机种子，该误差均为高斯分布的均方根(Root Mean Square, RMS)误差，并且在2.5处截断。模拟结果见图1(a)所示，图1(a)中线性拟合的直线的斜率表示准直误差对闭轨畸变的放大倍数。误差棒为平均值上下的一个标准偏差，采用同样的方法来模拟四极铁的转角误差、四极铁的垂直准直误差以及六极铁的垂直准直误差对垂直闭轨畸变的影响，得到图1(b)、(c)、(d)。

图1 由二极铁转角(a)、四极铁转角(b)、四极铁垂直非准直(c)、六极铁垂直非准直(d)所引起的垂直方向 闭轨畸变的均方根值 Fig.1 RMS of the vertical COD resulting from tilts of dipole (a) and quadrupoles (b), vertical misalignments on the quadrupole magnets (c) and the sextupole magnets (d).

定义所有BPM处的闭轨畸变均方根值与四极磁铁垂直准直误差的均方根值的比值为四极磁铁垂直非准直对闭轨畸变影响的轨道放大因子，其计算公式可以表示为[12]：

式中：ν为垂直工作点。在上海光源二期磁聚焦结构中，将Twiss参数和磁铁的强度代入到式(3)中，得到其放大因子为42.70，这与图1(c)的模拟结果40.28近似相符合。

单独分析每种磁铁准直误差对垂直闭轨畸变的影响表明：四极铁的垂直准直误差对垂直方向闭轨畸变的影响很大，这将引起束流在穿越六极铁时，距六极铁中心有很大的偏移，进而产生垂直色散和横向耦合，从而增加垂直发射度，这与六极铁本身存在垂直准直误差时具有相同的影响，因此不能单独模拟分析四极铁的垂直准直误差对垂直发射度的影响。其他三类准直误差对垂直闭轨畸变几乎没有影响，所以才可以单独模拟这三种准直误差对垂直发射度的影响。各磁铁准直误差引起的垂直色散和横向耦合所带来的垂直发射度的近似计算公式分别如下[10,12]：

将上海光源二期工程磁聚焦结构的Twiss参数和磁铁强度代入式(4)、(5)，得到：

式中：ε_dispersion、ΔΘ、ΔY、ΔΘ和ΔY的单位分别为nm∙rad、mrad、mrad、mm和mm。可对式(6)、(7)分离出每种磁铁的准直误差对垂直发射度的影响，结果如下：

其中，二极铁的转角误差只对垂直色散有影响，故式(8)只有一项；四极铁转角误差对垂直色散和横向耦合都有影响，故式(9)实际上是对垂直发射度的两项贡献的和；六极铁的垂直准直误差对垂直色散和横向耦合也都有影响，故式(10)也是对垂直发射度的两项贡献的和。

在设计改造上海光源二期工程时，只是修改了磁聚焦结构中很少的一部分，大部分磁铁的位置没有改变，因而现在运行的上海光源磁铁的准直公差将会沿用到二期工程中。因此，在模拟上海光源二期工程中每类磁铁准直误差对垂直发射度的影响时，在理想的磁聚焦结构模型中引入的磁铁准直误差分别为：二极铁和四极铁的转角误差为0.2 mrad，四极铁和六极铁的垂直准直误差为0.15 mm，上述误差均为高斯分布的均方根误差，并且在2.5处截断。决定了各磁铁的准直误差后，选择1000个随机种子来分别模拟二极铁转角误差、四极铁转角误差和六极铁的垂直准直误差对垂直发射度的影响。模拟过程中，垂直发射度采用Ohmi介绍的方法来计算[13]，其原理为利用粒子追踪计算六维束流包络矩阵。模拟结果见图2，误差棒为平均值上下的一个标准偏差。

图2 二极铁转角(a)、四极铁转角(b)、六极铁垂直 准直误差(c)分别在1000个随机种子中引起的 垂直发射度的分布 Fig.2 Distribution of the vertical emittance from tilts of dipole (a) and quadrupoles (b), vertical misalignments of the sextupole magnets (c) for 1000 random seeds.

在图2中，每类磁铁误差引起的垂直发射度的分布中，有几个随机种子的垂直发射度远大于其他种子的垂直发射度，这是由于随机产生的较大误差值的磁铁处在具有较大色散函数或较大beta函数的位置上。模拟结果：二极铁转角、四极铁转角和六极铁垂直非准直所引起的垂直发射度分别在1000个随机种子中的平均值为4.0×10‒4nm·rad、1.1×10‒2nm·rad、0.18 nm·rad。由式(8)‒(10)可得到：

可以看出，近似公式计算结果与模拟结果近似相等，两者基本符合。

从模拟和近似公式的计算中得出结论，在上海光源二期工程中，保持各类磁铁的准直误差与目前正在运行的上海光源相同的情况下，由二极铁转角误差所带来的垂直发射度很小，垂直发射度主要来源于四极铁的转角和六极铁的垂直准直误差，其中六极铁的垂直准直误差对垂直发射度的贡献非常大。

3 上海光源二期工程的垂直发射度校正

虽然不能单独模拟分析四极铁的垂直准直误差对垂直发射度的影响，但是为使模拟更接近于真实的机器，在下面校正随机种子引起的垂直发射度时，将包含所有磁铁的准直误差，无论该误差对垂直发射度贡献的大小。首先设置各磁铁的准直误差，即二极铁和四极铁的转角误差为0.2 mrad，四极铁和六极铁的垂直准直误差为0.15 mm，上述误差均为高斯分布的均方根误差，并且在2.5处截断。为使模拟更具有合理性，共使用了200个随机种子。

由于随机种子误差中引入了四极铁的垂直准直误差，这很大程度地增加了垂直方向的闭轨畸变，导致束流在通过六极铁时与六极铁的中心具有很大的偏移，进而产生垂直色散和横向耦合，所以在采用斜四极铁来校正垂直发射度之前，先要使用校正铁来校正垂直方向的闭轨畸变。

3.1 垂直闭轨畸变校正

闭轨校正的目的是使束流轨道尽可能地接近理想轨道，所采用的方法为奇值分解法(Singular Values Decomposition, SVD)[14]。首先由BPM获取束流轨道信息，然后采用SVD方法计算校正铁的强度并改变校正铁的电流值使束流轨道尽可能接近理想轨道。响应矩阵(×)可以定义为：=，其中表示BPM的个数；表示校正铁的个数，表示轨道畸变向量；表示校正铁强度变化向量。响应矩阵为一个不可求逆的矩阵，但可以通过SVD来求响应矩阵的伪逆矩阵。根据SVD法，对于任意一个矩阵，都有可以写成=USV，其中：为×的对角阵，其对角线的值称为的奇值，它们都是大于或等于0的数；为×的正交矩阵；为×的正交矩阵。响应矩阵的伪逆矩阵为：‒1=‒1U，因此有：

在闭轨校正时，使用了140个BPM和80个校正子，其位置与目前运行的上海光源中BPM与校正铁的位置相同。由于只关心垂直方法的闭轨，所以只校正了垂直方向的闭轨畸变，校正前后的闭轨畸变均方根值如图3所示，

图3 200个随机种子在校正前(a)和校正后(b)的垂直 闭轨畸变的均方根值 Fig.3 RMS of the vertical COD of 200 uncorrected (a) and corrected (b) random machine seeds.

校正之后，很大程度地降低了垂直方向的闭轨畸变，大部分随机种子的闭轨畸变均方根值都在0.06 mm以下，如图3(b)所示。

3.2 垂直发射度校正

在校正垂直方向的闭轨畸变之后，使用斜四极铁采用梯度下降法来校正垂直发射度。梯度下降法的基本原理为：一个有界可微的多变量函数在某一点处，函数值下降最快的方向为该点的负梯度的方向，即‒Ñ()的方向。这种方法需要先猜想一个初始点0为函数的局部最小值点，随后的点变化方式为：

则可以得到：(0)≥(1)≥(2)…，其中为一个非常小的数值，并且它可以在每次迭代中变化，的选取条件满足：(x‒γÑ(x))=min(x‒Ñ(x))。由于耦合度随斜四极铁强度变化的函数解析式未知，所以采用了非确定性线搜索并且满足Wolfe条件而定，即迭代寻找使耦合度在某一组斜四极铁强度处的值比上一组的值要小。

在模拟校正时，使用了60个斜四极铁来校正上海光源二期的垂直发射度，这些斜四极铁分别为S3（14个）、S8（4个）、S10（4个）、SFS（两个）、SD（36个）这60个六极铁的辅助线圈。

目前上海光源在正常运行时的耦合度为0.5%，所以在校正上海光源二期耦合度时，选定的目标值也为0.5%，其基本过程为扫描斜四极铁的强度值对耦合度的响应，得到其梯度值，然后采用梯度下降法来一步步循环迭代，当某一步计算的耦合度小于0.5%，则立即退出迭代循环。图4为校正前后的耦合度变化。

图4 200个随机种子在校正前(a)和校正后(b)的耦合度 Fig.4 Emittance coupling of 200 uncorrected (a) and corrected (b) random machine seeds.

从图4(b)中可以看出，200个随机种子的耦合度都被降到了0.5%以下。校正前与校正后的垂直色散的变化如图5所示。从图5(b)中可以看出，200个随机种子的垂直色散都得到了很大程度地减小。校正时所使用的60个斜四极磁铁的强度都处在目前上海光源的斜四极铁供电电源的供给范围内。

图5 200个随机种子在校正前(a)和校正后(b)的垂直 色散的均方根值 Fig.5 RMS of the vertical dispersion of 200 uncorrected (a) and corrected (b) random machine seeds.

4 结语

在上海光源二期工程的磁铁准直误差与目前正在运行的上海光源相同的条件下，垂直发射度主要来源于四极铁的转角和六极铁的垂直准直误差，其中六极铁准直误差对垂直发射度的贡献最大；对于由所有磁铁准直误差所引起的垂直发射度，研究结果表明使用60个斜四极铁来进行校正，能使耦合度保持在0.5%以下。

1 Xu H J, Zhao Z T. Current status and progressed of SSRF project[J]. Nuclear Science and Techniques, 2008, 19: 1‒6. DOI: 10.1016/S1001-8042(08)60013-5

2 Liu G M, Dai Z M, Li H H,. Lattice design for SSRF storage ring[J].and, 2006, 30(Suppl I): 144‒146. DOI: 10.3321/j.issn:0254- 3052.2006.z1.047

3 Zhao Z T, Yin L X, Leng Y B,Consideration on the future major upgrades of the SSRF storage ring[C]. Proceedings of IPAC2015, Richmond, VA, USA, 2015: 1672‒1674

4 Tian S Q, Jiang B C, Leng Y B,. Double-mini-βoptics design in the SSRF storage ring[J]. Nuclear Science and Techniques, 2014, 25: 030101. DOI: 10.13538/ j.1001-8042/nst.25.030101

5 Tian S Q, Jiang B C, Zhou Q G,. Lattice design and optimization of the SSRF storage ring with super-bends[J]. Nuclear Science and Techniques, 2014, 25: 010102. DOI: 10.13538/j.1001-8042/nst.25.010102

6 Zhang Q L, Tian S Q, Jiang B C,. Beam dynamics of the superconducting wiggler on the SSRF storage ring[J]. Chinese Physics C, 2016, 40(3): 037001. DOI: 10.1088/ 1674-1137/40/3/037001

7 Aiba M, Böge M, Milas N,. Ultra low vertical emittance at SLS through systematic and random optimization[J].Instruments andinSection A, 2012, 694: 133‒139. DOI: 10.1016/j.nima.2012.08.012

8 Dowd R, Boland M, BlancLe G,. Achievement of ultralow emittance coupling in the Australian synchrotron storage ring[J]. Physical ReviewA & B, 2011, 14: 012804. DOI: 10.1103/PhysRevSTAB.14. 012804

9 Zhang M Z, Hou J, Li H H. Coupling measurement and correction at the SSRF storage ring[J]. Science- Physics Mechanics, 2011, 54(Suppl 2): s201‒s206. DOI: 10.1007/s11433-011-4513-8

10 Raubenheimer T O. Tolerances to limit the vertical emittance in future storage rings[J].Partners, 1991, 36(SLAC-PUB-4937): 75‒120

11 Terebilo A. Accelerator toolbox for MATLAB[R]. Stanford Linear Accelerator Center, Menlo Park, CA (US), 2001

12 Wolski A. Low-emittance storage rings[R]. CERN Accelerator School: Advanced Accelerator Physics Course, Trondheim, Norway, 2013. DOI: 10.5170/ CERN-2014-009.245

13 Ohmi K, Hirata K, Oide K. From the beam-envelope matrix to synchrotron-radiation integrals[J]. Physics Review E, 1994, 49: 751‒765. DOI: 10.1103/PhysRevE. 49.751

14 Hou J, Tian S Q, Zhang M Z,. Studies of closed orbit correction and slow orbit feedback for the SSRF storage ring[J]. Chinese Physics C, 2009, 33: 145‒150. DOI: 10.1088/1674-1137/33/2/013

上海光源维修改造经费成果奖励(No.1173000212)资助

Supported by the Maintenance and Upgrade Foundation of SSRF (No.1173000212)

Study on the vertical emittance for the SSRF Phase II project

LI Changliang1,2TIAN Shunqiang1ZHANG Wenzhi1

1(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China)2(University of Chinese Academy of Sciences, Beijing 100049, China)

Background: The spectral brightness is an important parameter for the synchrotron radiation light source. It is intimately associated with beam current, horizontal emittance and vertical emittance. The beam current is limited by collective effects and the horizontal emittance is decided by the lattice structure, which has almost been fixed in a machine. Reducing the vertical emittance to improve the brightness is a simple and feasible method. Other benefits such as the increase of transverse coherence and the reduction of vertical spot size can be obtained in reducing the vertical emittance. Purpose: The aim is to study the sensitivity of the vertical emittance to magnetic alignment errors and correct the vertical emittance for the Shanghai Synchrotron Radiation Facility (SSRF) Phase II project. Methods: Singular values decomposition (SVD) is used to correct the vertical closed orbit distortion (COD). Sixty skew quadrupole magnets are utilized to correct the emittance coupling with gradient descent algorithm. Results: An approximate mean of 0.06-mm vertical COD is obtained after orbit correction. The emittance coupling of each seed can be corrected below 0.5% by using 60 skew quadrupole magnets. Conclusion: If the alignment tolerances of all magnets in the SSRF Phase II project are the same as those in the currently operational SSRF storage ring, the vertical emittance is primarily determined by vertical misalignments of the sextupole magnets, and the emittance coupling can be corrected below 0.5% by using 60 skew quadrupole magnets.

SSRF Phase II project, Magnetic alignment errors, Vertical emittance, Coupling correction

LI Changliang, male, born in 1988, graduated from Hubei University of Science and Technology in 2013, master student, focusing on accelerator physics and beam dynamics

ZHANGWenzhi,E-mail:zhangwenzhi@sinap.ac.cn

TL594

10.11889/j.0253-3219.2016.hjs.39.050102

李昌亮，男，1988年出生，2013年毕业于湖北科技学院，现为硕士研究生，研究领域为加速器物理与束流动力学

张文志，E-mail: zhangwenzhi@sinap.ac.cn

2016-02-25，

2016-03-20