research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Investigation of the damping wiggler effect and application on the PAL fourth-generation storage ring

CROSSMARK_Color_square_no_text.svg

aDepartment of Physics, POSTECH, Pohang, Gyungbuk 37673, South Korea, bPohang Accelerator Laboratory, POSTECH, Pohang, Gyungbuk 37673, South Korea, and cSLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA
*Correspondence e-mail: tlssh@postech.ac.kr

Edited by A. F. Craievich, University of São Paulo, Brazil (Received 31 January 2020; accepted 22 August 2020; online 22 October 2020)

An investigation of the damping wiggler effect to reduce the emittance in the Pohang Accelerator Laboratory (PAL), a fourth-generation storage ring (4GSR) that uses a multi-bend achromat, is presented. A 4GSR lattice which has reduced emittance and increased dynamic aperture to amplify the synergy with two existing light sources (PLS-II and PAL-XFEL) at PAL is described.

1. Introduction

A fourth-generation storage ring (4GSR) that uses the multi-bend achromat (MBA) lattice concept may be able to surpass the brightness and coherence that are attained using present third-generation storage rings (3GSRs). 4GSRs are part of a world-wide effort to exceed the brightness and coherence reached by present 3GSRs. MAX IV (MAX IV, 2020[MAX IV (2010). MAX IV DDR - The MAX IV Detailed Design Report, https://www.maxiv.lu.se/accelerators-beamlines/accelerators/accelerator-documentation/max-iv-ddr/.]) was the first MBA machine; the next, Sirius (Liu et al., 2013[Liu, L., Milas, N., Mukai, A. H. C., Resende, X. R., Rodrigues, A. R. D. & de Sá, F. H. (2013). Proceedings of the Fourth International Particle Accelerator Conference (IPAC2013), 12-17 May 2013, Shanghai, China, pp. 1874-1876. TUPWO001.]), is currently under construction. Other projects are being conducted to convert existing 3GSRs such as ESRF, APS, SPring-8 and ALS (Revol, 2017[Revol, J.-L. (2017). ESRF-EBS: The Extremely Brilliant Source Project, https://indico.psi.ch/conferenceDisplay.py?confId=5589.]; APS, 2020[APS (2020). APS Upgrade, https://www.aps.anl.gov/APS-Upgrade.]; SPring-8, 2014[SPring-8 (2014). SPring-8-II Conceptual Design Report, https://rsc.riken.jp/eng/pdf/SPring-8-II.pdf.]; Steier et al., 2016[Steier, C., Byrd, J. M., De Santis, S., Nishimura, H., Robin, D., Sannibale, F., Sun, C., Venturini, M. & Wan, W. (2016). Proceedings of the 7th International Particle Accelerator Conference (IPAC2016), 8-13 May 2016, Busan, Korea, pp. 2956-2958. WEPOW049.]) to 4GSRs. These modifications will reduce beam emittance to a few hundred picometres, or even to ≤100 pm, but many efforts have been made to reduce the natural emittance further. The natural emittance in an electron storage ring is given by

[{\varepsilon_0} = {C_q}{\gamma^2}\, {{ \left\langle{\cal{H}}/\left|\rho^3\right|\right\rangle }\over{ {j_x}\left\langle1/\rho^2\right\rangle }}, \eqno(1)]

where Cq is a constant (for electrons [{C_q} \simeq 3.832 \times {10^{ - 13}}\ {\rm{m}}]), [{{\gamma}}] is the relativistic factor for a particle, [{\cal{H}}] ≃ 1/ρ2 (where ρ is the unit-cell bending angle) and jx is the horizontal damping partition number (Borland et al., 2014[Borland, M., Decker, G., Emery, L., Sajaev, V., Sun, Y. & Xiao, A. (2014). J. Synchrotron Rad. 21, 912-936.]). The natural emittance scales with the inverse cube of ρ, so MBA lattices enable emittance reduction by one to two orders of magnitude compared with third-generation light sources.

The advent of 4GSRs enables a ring free-electron laser (FEL) to perform high-average-power extreme ultraviolet lithography (Luong, 2018[Luong, V. (2018). EUV Lithography, https://eng.kuleuven.be/en/research/phd/arenberg-youngster-seminars/ays-pdf/20180516-luong-ays-final.pdf.]) due to its ability to generate extremely low electron beam emittance. However, 4GSRs still have limitations such as relatively large energy spread and low repetition rate in FELs in the bypass line from the storage ring. A method to compensate for the effects of beam energy spread by introducing a transverse variation of the undulator magnetic field has been proposed (Huang et al., 2012[Huang, Z., Ding, Y. & Schroeder, C. B. (2012). Phys. Rev. Lett. 109, 204801.]), but the average power of the photon beam is insufficient because of the low repetition rate in a ring FEL that uses a bypass line (Cai et al., 2013[Cai, Y., Ding, Y., Hettel, R., Huang, Z., Wang, L. & Xiao, L. (2013). An X-ray Free Electron Laser Driven by an Ultimate Storage Ring. Report SLAC-PUB-15380. SLAC National Accelerator Laboratory, Menlo Park, California, USA.]). To further reduce the emittance, the concepts of the longitudinal gradient bend and reverse bends have been considered for a lattice design to achieve sub-theoretical minimum emittance (Riemann & Streun, 2019[Riemann, B. & Streun, A. (2019). Phys. Rev. Accel. Beams, 22, 021601.]). Use of a damping wiggler is an option to reduce emittance further, and its effect should be explored. The damping wiggler increases the repetition rate of the FEL beam by reducing the damping time; this effect must be quantified.

Recently a 4GSR has been designed at the Pohang Accelerator Laboratory (PAL-4GSR) using concepts from ESRF-EBS and APS-U. As a result, a lattice with 3 GeV energy and 90 pm emittance with a 570 m ring was proposed. We investigate the effect of a damping wiggler on PAL-4GSR, particularly focusing on emittance reduction and repetition rate increase. In this paper, we describe the PAL-4GSR lattice, and the overall damping-wiggler effect on PAL-4GSR. Section 2[link] describes the lattice design and beam dynamics studies for PAL-4GSR. Section 3[link] presents derivation of analytic formulas for the damping wiggler effect on beam parameters. Section 4[link] applies the analytic formulas to PAL-4GSR and compares them with simulations conducted using Elegant software (Borland, 2000[Borland, M. (2000). Elegant: A Flexible SDDS-Compliant Code for Accelerator Simulation. Technical Report LS-287. Advanced Photon Source, Argonne, IL, USA.]). Section 5[link] presents a summary and concluding remarks.

2. PAL-4GSR

The PAL-4GSR storage ring (Table 1[link]) is a hybrid seven-bend achromat lattice with a horizontal emittance of 90 pm. The ring is 570 m in circumference and is composed of 20 symmetrical cells. The length of the straight section is set as 6.5 m to accommodate two superconducting RF modules, as a result of experience on PLS-II (Shin et al., 2013[Shin, S., Kwon, S., Kim, D.., Kim, D.., Kim, M., Kim, S.., Kim, S.., Kim, J., Kim, C., Park, B., Park, S.., Park, S.., Park, E., Son, Y., Yoon, J., Lee, B., Lee, E., Lee, J.., Lee, H.., Joo, Y., Choi, J., Ha, T., Hwang, W., Hwang, I., Lee, J.., Oh, B., Lee, C.., Lee, H.., Kim, J.., Hwang, J. Y., Nam, S. H. & Cho, M. (2013). J. Instrum. 8, P01019.]). The PAL-4GSR lattice contains a 2 T super-bend in the central dipole to produce radiation with a critical energy of 12 keV (Fig. 1[link]).

Table 1
Main parameters of the PAL-4GSR

Parameter Value Unit
Energy 3 GeV
Emittance (flat / round) 89.4 / 58 pm
Circumference 570 m
Tune (x / y) 48.380 / 17.376
Natural chromaticity (x / y) −95.97 / −57.91
Radiation loss per turn 0.468 (without ID) MeV
Momentum compaction 1.45 × 10−4  
Damping partition (H / V / L) 1.82 / 1.00 / 1.18  
Damping time (H / V / L) 13.37 / 24.35 / 20.66 ms
[Figure 1]
Figure 1
Lattice functions for PAL-4GSR.

The concepts of the ESRF-EBS and APS-U lattices were adopted in the PAL-4GSR lattice. The dispersion was deliberately enlarged between the first and second dipoles and between the sixth and seventh dipoles, and chromatic sextupoles have been placed in this dispersion-bump region to reduce the strength required to control the chromaticity. The betatron phase advances between the two dispersion bumps were set to [\Delta{\varphi}_{x}] ≃ 3π in the horizontal plane and [\Delta{\varphi}_{y}]π in the vertical plane. As a result, nonchromatic effects of the sextupoles are canceled naturally to first order. To minimize natural emittance, five-step longitudinal gradient dipoles and reverse bending magnets were considered.

The sextupoles are two-cell periodic in order to increase the number of control knobs. Therefore, a total of ten sextupoles in two cells, except two sextupoles for correcting the linear chromaticity (Fig. 2[link]), are used in the optimization with a multi-objective genetic algorithm (MOGA) (Deb, 2001[Deb, K. (2001). Multi-Objective Optimization Using Evolutionary Algorithms. Chichester: John Wiley & Sons.]). The objectives are an on-momentum dynamic aperture for off-axis beam injection and an off-momentum dynamic aperture (4%) to extend the Touschek lifetime. The result (Fig. 3[link]) of MOGA optimization seems to have converged after 150 generations, at which the strengths of ten harmonic sextupoles reach a compromise between on-momentum dynamic aperture and 4% off-momentum dynamic aperture. The candidate point for the studies (marked by the star in Fig. 3[link]) has a maximum on-momentum dynamic aperture area among the flat region for the off-momentum dynamics aperture area. As a result, large dynamic apertures were attained, enough for the off-axis injection with random errors (Fig. 4[link]). Here, 30 µm r.m.s. offset, 200–500 µrad r.m.s. rotation and 0.1–0.2% r.m.s. field errors are assumed to generate 100 random errors. Despite assuming this harmful error scenario, a horizontal dynamic aperture from −9 mm to 12 mm is promising for 6 mm off-axis beam injection.

[Figure 2]
Figure 2
Configuration of sextupoles and other magnetic elements in two cells.
[Figure 3]
Figure 3
MOGA optimization for large dynamic aperture area. The star on the Pareto front shows the candidate point for the studies; it has a maximum on-momentum dynamic aperture area among the flat region for the off-momentum dynamics aperture area.
[Figure 4]
Figure 4
Dynamic aperture with errors for beam injection.

To obtain maximum Touschek lifetime, the RF voltage is scanned in the particle tracking with the optimized sextupole configuration (Fig. 5[link]). At 3.1 MV RF voltage, which corresponds to 6% RF acceptance, Touschek lifetimes are 6 h for 1% coupling flat beam and 58 h for coupling round beam. Local momentum aperture tracking including error distributions is also compared with ideal local momentum aperture and RF momentum aperture (Fig. 6[link]). Evaluated Touschek lifetimes in realistic conditions are 1.3 h for 1% coupling flat beam and 9.7 h for round beam at 1.9 MV optimum RF voltage, which corresponds to 4.3% RF acceptance. These lifetime reductions invoke the need of a harmonic cavity. Note that bunch lengthening from 3.45 mm r.m.s. to 20 mm r.m.s. with a harmonic cavity enables Touschek lifetimes to be increased up to 7 h for 1% coupling flat beam and 51 h for round beam. Here the factor 5.8 of bunch lengthening from the higher harmonic cavity is assumed to obtain sufficient lifetime in hybrid fill pattern operation with high single bunch current. In addition to the Touschek lifetime, the elastic scattering lifetime is 61 h, assuming a minimum half physical aperture of 2 mm, and the bremsstrahlung lifetime is 79 h.

[Figure 5]
Figure 5
Touschek lifetime along RF voltage. At 3.1 MV RF voltage, Touschek lifetimes are 6 h for 1% coupling flat beam and 58 h for round beam.
[Figure 6]
Figure 6
Local momentum acceptance. Green line: lattice acceptance without error; blue line: lattice acceptance with error; red line: RF acceptance.

The results of Elegant simulation of intra-beam scattering (Fig. 7[link]) show a change in the zero-charge emittance along beam energy, which scales as ∼E2. In the simulation, the algorithm uses the Bjorken & Mtingwa formula (Bjorken & Mtingwa, 1983[Bjorken, J. D. & Mtingwa, S. K. (1983). Part. Accel. 13, 115-143.]) to simulate the intra-beam scattering (IBS) effect. For a 0.845 nC bunch, which is a single bunch current for a 400 mA multi-bunch, round beam operation must be used to avoid a large emittance increase due to the IBS effect. In addition, a harmonic cavity can be considered for bunch lengthening up to 20 mm. Here, the natural bunch length is 3.45 mm and the lengthening factor with harmonic cavity is 5.8.

[Figure 7]
Figure 7
Emittance along beam energy. (a) Round beam (100% coupling). (b) Flat beam (1% coupling).

3. Analytic formulas for the damping wiggler effect

In an electron storage ring, representative ring parameters can be expressed using synchrotron radiation integrals

[{\epsilon }_{x} = {C}_{q}{\gamma }^{2}\,{{{I}_{5}}\over{{I}_{2}-{I}_{4}}}, \eqno(2)]

[\sigma _\delta ^2 = {C_q}{\gamma ^2}\,{{{I_3}} \over {2{I_2} + {I_4}}}, \eqno(3)]

[{U_0} = {{{C_\gamma }} \over {2\pi }} \, {E^4}{I_2}, \eqno(4)]

[{\tau _x} = {{4\pi {T_0}} \over {{C_\gamma }{E^3}\left({{I_2} - {I_4}} \right)}}, \eqno(5)]

[{\tau _y} = {{2\pi {T_0}} \over {{C_\gamma }{E^3}{I_2}}}, \eqno(6)]

[{\tau _\delta } = {{4\pi {T_0}} \over {{C_\gamma }{E^3}\left({2{I_2} + {I_4}} \right)}}, \eqno(7)]

[{\alpha _c} = {{{I_1}} \over C}, \eqno(8)]

where Cq ≃ 3.84 × 10−13 m, [{C}_{\gamma }] ≃ 8.846 × 10−5 m GeV−3, [{\epsilon}_{x}] is the natural horizontal emittance, E is the electron energy and τ is the damping time. The definition of the synchrotron radiation integral is given by

[{I_1} = \mathop \oint \nolimits^ {\rm{d}}s\left\{ {{\eta \over \rho }} \right\},]

[{I_2} = \mathop \oint \nolimits^ {\rm{d}}s\left\{ {{1 \over {{\rho ^2}}}} \right\},]

[{I_3} = \mathop \oint \nolimits^ {\rm{d}}s\left\{ {{\eta \over {{{\left| \rho \right|}^3}}}} \right\}, \eqno(9)]

[{I_4} = \mathop \oint \nolimits^ {\rm{d}}s\left\{ {{\eta \over \rho }\left({{1 \over {{\rho ^2}}} + 2{k_1}} \right)} \right\},]

[{I_5} = \mathop \oint \nolimits^ {\rm{d}}s\left\{ {{H \over {{{\left| \rho \right|}^3}}}} \right\},]

where k1 = [({1/{B\rho }})({{\partial {B_y}}/{\partial x}})], H = [\gamma {\eta ^2}] + [2\alpha \eta \eta^\prime] + [\beta \eta^{\prime2}] (Wiedemann, 1993[Wiedemann, H. (1993). Particle Accelerator Physics I. Berlin: Springer.]), and ρ is the curvature of a reference particle. These integrals are changed by considering a damping wiggler that has a vertical magnetic field given by

[{B_y} = {B_{\rm{w}}}\sin{k_{\rm{w}}}s, \eqno(10)]

where Bw = [B\rho/{\rho_{\rm{w}}}] is the peak field, kw = [2\pi/{\lambda_{\rm{w}}}] is the wavenumber, and [{\lambda _{\rm{w}}}] is the wiggler period. For the 3 GeV electron storage ring, the magnet rigidity [B\rho] is ∼10 T m. The synchrotron radiation integrals for the wiggler in the straight section are (Appendix A[link])

[{I_{1,{\rm{wig}}}} = - {{\lambda _{\rm{w}}^2B_{\rm{w}}^2{L_{\rm{w}}}} \over {8{\pi ^2} \times {{10}^2}}},]

[{I_{2,{\rm{wig}}}} = {{B_{\rm{w}}^2{L_{\rm{w}}}} \over {2 \times {{10}^2}}},]

[{I_{3,{\rm{wig}}}} = {{4B_{\rm{w}}^3{L_{\rm{w}}}} \over {3\pi \times {{10}^3}}}, \eqno(11)]

[{I_{4,{\rm{wig}}}} = - {{3\lambda _{\rm{w}}^2B_{\rm{w}}^4{L_{\rm{w}}}} \over {32{\pi ^2} \times {{10}^4}}},]

[{I_{5,{\rm{wig}}}} = {{\langle{\beta _x}\rangle\lambda _{\rm{w}}^2B_{\rm{w}}^5{L_{\rm{w}}}} \over {15{\pi ^3} \times {{10}^5}}},]

where Lw is the length of the wiggler and [{\beta _x}] is a horizontal beta-function. Equation (11)[link] applies to a sinusoidal wiggler field versus s, but in reality the field of a long-period insertion device (ID) includes high harmonics of kw.

The emittance ratio between the storage ring emittances without and with wiggler is given by

[\eqalignno{ {{\epsilon_{\rm{tot}}}\over{\epsilon_{\rm{ring}}}} & = {{ C_q\gamma^2{{I_{5,{\rm{tot}}}}\over{I_{2,{\rm{tot}}}\,-\,I_{4,{\rm{ring}}}}} }\over{ C_q\gamma^2{{I_{5,{\rm{ring}}}}\over{I_{2,{\rm{ring}}}\,-\,I_{4,{\rm{ring}}}}} }} = {{ C_q\gamma^2{{I_{5,{\rm{ring}}}\,+\,I_{5,{\rm{wig}}}}\over{I_{2,{\rm{ring}}}\,+\,I_{2,{\rm{wig}}}\,-\,I_{4,{\rm{ring}}}\,-\,I_{4,{\rm{wig}}} }} }\over{ C_q\gamma^2{{I_{5,{\rm{ring}}}}\over{I_{2,{\rm{ring}}}\,-\,I_{4,{\rm{ring}}}}} }} \cr& = {{ 1+{{I_{5,{\rm{wig}}}}\over{I_{5,{\rm{ring}}}}} }\over{ 1+{{I_{2,{\rm{wig}}}\,-\,I_{4,{\rm{wig}}}}\over{I_{2,{\rm{ring}}}\,-\,I_{4,{\rm{ring}}}}} }} &(12)}]

where

[{I_{4,{\rm{wig}}}} = -{{\lambda _{\rm{w}}^2B_{\rm{w}}^2} \over {16{\pi ^2} \times {{10}^2}}}\,{I_{2,{\rm{wig}}}},]

which is very small compared with I2,wig, so I4,wig can be neglected for typical wiggler values in equation (12)[link]. The emittance ratio without I4,wig becomes

[{{ \epsilon_{\rm{tot}} }\over{ \epsilon_{\rm{ring}} }} = {{ 1+bL_{\rm{w}} }\over{ 1+aL_{\rm{w}} }}, \eqno(13)]

where

[a\equiv {{{B}_{\rm{w}}^{2}}\over{2\left({I}_{2,{\rm{ring}}}-{I}_{4,{\rm{ring}}}\right)\times {10}^{2}}}]

and

[b\equiv {{\langle {\beta }_{x}\rangle {\lambda }_{\rm{w}}^{2}{B}_{\rm{w}}^{5}}\over{15{\pi }^{3}{I}_{5,{\rm{ring}}}\times {10}^{5}}}.]

The form of the function (1 + bLw)/(1 + aLw) depends on the ratio of a and b, which are always positive. [{\epsilon}_{{\rm{tot}}}/{\epsilon}_{{\rm{ring}}}] = 1 for a = b, [{\epsilon}_{\rm{tot}}/{\epsilon}_{{\rm{ring}}}] = >1 for [a\,\lt\,b], and [{\epsilon}_{{\rm{tot}}}/{\epsilon}_{{\rm{ring}}}] < 1 for [a\,\gt\,b]. Here our main interest is the case of [{{\epsilon_{\rm{tot}}}/{\epsilon_{\rm{ring}}}}] = <1. When a = 1.389Bw2 and b = 0.486Bw5 the emittance ratio (Fig. 8[link]) along the length of a damping wiggler is limited on the asymptotic line given by b/a, despite the extremely long wiggler length. Therefore, short wiggler period, small horizontal beta-function, long wiggler length and low magnetic field are favorable conditions to reduce total emittance when the wiggler is infinitely long. Minimum emittance reduction ratio along natural storage ring emittance (Fig. 9[link]) implies that the emittance reduction by a damping wiggler is not so effective when natural storage-ring emittance is extremely small.

[Figure 8]
Figure 8
Emittance ratio with the function of wiggler length.
[Figure 9]
Figure 9
Minimum emittance reduction ratio versus natural storage ring emittance.

In the same manner, the ratios of energy spread, damping time and momentum compaction factor are also formulated as (Helm et al., 1973[Helm, R. H., Lee, M. J., Morton, P. L. & Sands, M. (1973). IEEE Trans. Nucl. Sci. 20, 900-901.])

[{{ \sigma_{\delta,{\rm{tot}}}^{2} }\over{ \sigma_{\delta,{\rm{ring}}}^{2} }} = {{1+dL_{\rm{w}}}\over{1+cL_{\rm{w}}}},\eqno(14)]

where

[c \equiv {{ B_{\rm{w}}^{2} }\over{ (2I_{2,{\rm{ring}}}+I_{4,{\rm{ring}}})\times10^2 }} \quad{\rm{and}}\quad d \equiv {{ 4B_{\rm{w}}^{3} }\over{ 3\pi I_{3,{\rm{ring}}}\times10^3 }}]

[{{\tau_{\rm{tot}}}\over{\tau_{\rm{ring}}}} = {{1}\over{1+eL_{\rm{w}}}},\eqno(15)]

where

[\eqalign{ & e_x \equiv {{ B_{\rm{w}}^{2} }\over{ 2(I_{2,{\rm{ring}}}-I_{4,{\rm{ring}}})\times10^2 }} \quad{\rm{and}}\quad e_y \equiv {{ B_{\rm{w}}^{2} }\over{ 2I_{2,{\rm{ring}}}\times10^2 }}\quad{\rm{and}}\cr& e_\delta \equiv {{ B_{\rm{w}}^{2} }\over{ (2I_{2,{\rm{ring}}}+I_{4,{\rm{ring}}})\times10^2 }}, }]

[{{\alpha_{c,{\rm{tot}}}}\over{\alpha_{c,{\rm{ring}}}}} =1-fL_{\rm{w}} \quad{\rm{where}}\quad f \equiv {{\lambda_{\rm{w}}^2B_{\rm{w}}^2}\over{8\pi^2\alpha_{c,{\rm{ring}}}\,c\times10^2}}.\eqno(16)]

The equation for energy spread ratio has a similar form to that of the emittance ratio. Energy spread reduction is also limited by an asymptotic line. Damping time is always reduced by the damping wiggler and the reduction of momentum compaction factor is negligible because f in equation (16)[link] is small.

4. Application on PAL-4GSR

The damping wiggler effects on the main beam parameters in the lattice of PAL-4GSR were investigated using the analytic formulae derived in Section 3[link] and the Elegant simulation. Four straight sections among a total of 20 straight sections are considered for each 4.5 m-long damping wiggler. In this case, four ID straight sections should be sacrificed for insertion of damping wigglers. The synchrotron radiation integrals (Table 2[link]) for PAL-4GSR were obtained using simulation in Elegant. The theoretical results were also compared with the results of the Elegant simulation.

Table 2
Synchrotron radiation integrals for the PAL-4GSR

Parameter Value Unit
I1,ring 9.251 ×10 - 2 m
I2,ring 4.108 ×10 - 1 m−1
I3,ring 4.254 ×10 - 1 m−2
I4,ring - 3.317 ×10 - 1 m−1
I5,ring 6.068 ×10 - 6 m−1

Assuming a damping wiggler that had a total length of 18 m, the emittance-reduction ratios were investigated by scanning the wiggler period and the wiggler peak field (Fig. 10[link]). The function plot from theoretical results agreed well with the simulation result obtained using Elegant. For a given wiggler period and wiggler length, an optimal wiggler peak field achieves the smallest emittance-reduction ratio (Fig. 10[link], stars); this peak field increases as the wiggler period decreases.

[Figure 10]
Figure 10
Emittance change by using a damping wiggler. Dots: simulation results; lines: function plot; stars: optimal combination.

The allowable peak field should be obtained for a given short wiggler period. The Halbach formula was used to determine the optimal wiggler field along the wiggler period and magnetic field versus wavelength (Fig. 11[link]). For the total 18 m wiggler length in PAL-4GSR, the minimum emittance reduction ratio is at most ∼0.83. As already shown in Fig. 9[link], emittance reduction by the damping wiggler is not so effective due to the small natural emittance.

[Figure 11]
Figure 11
Optimum magnetic field along wiggler period at a given total wiggler length from Fig. 10[link]. Black line: effect of wiggler period on wiggler magnetic field, calculated using the Halback formula at the given minimum gap of 5 mm. The region above the red line indicates a minimum gap < 5 mm.

For the FEL option through the bypass line, damping-time reduction was also investigated to estimate the increase of repetition rate for high average FEL power. Horizontal, vertical and longitudinal damping times have the same form; the only difference is in the coefficients [equation (15)[link]]. The coefficients are independent of the wiggler period. In PAL-4GSR, they are

[\eqalign{ e_x & = 6.73\times10^{-3}\,B_{\rm{w}}^{2}, \cr e_y & = 1.22\times10^{-2}\,B_{\rm{w}}^{2}, \cr e_\delta & = 2.04 \times10^{-2}\,B_{\rm{w}}^{2}.} \eqno(17)]

The damping-time reduction (Fig. 12[link]) by a wiggler was calculated. A half-damping time (or double repetition rate) can be expected in a magnetic field >2.5 T. Similarly, the energy spread reduction can also be estimated (Fig. 13[link]); an increase in energy spread is inevitable in a magnetic field >2 T. These results show that the damping wiggler does not have much effect on the increase of ring FEL power.

[Figure 12]
Figure 12
Reduction of damping times. Lines: damping time ratio plot; dots: simulation results. Red: horizontal damping time; blue: vertical damping time; green: longitudinal damping time.
[Figure 13]
Figure 13
Energy spread versus magnetic peak field. The critical field that makes energy spread >1 is ∼2 T.

5. Conclusion

Analytic analysis and Elegant simulation were performed to quantify the damping wiggler effects on PAL-4GSR which is a 3 GeV and 90 pm natural emittance storage ring with 570 m circumference. Analytic analysis agreed well with Elegant simulation, and described reductions in energy spread, damping time, and emittance. Emittance reduction by the damping wiggler was less effective than in third-generation storage rings due to the small natural emittance in the fouth-generation storage ring. Application of a damping wiggler on a ring FEL was also investigated in the aspect of increasing average brightness. Unfortunately the damping wiggler causes an increase in energy spread, and therefore cannot considerably increase the ring-FEL power.

APPENDIX A

Synchrotron radiation integrals

Helm et al. (1973[Helm, R. H., Lee, M. J., Morton, P. L. & Sands, M. (1973). IEEE Trans. Nucl. Sci. 20, 900-901.]) evaluated the synchrotron radiation integrals. Synchrotron integrals in the damping wiggler have been investigated by many authors (Katoh & Kamiya, 1987[Katoh, M. & Kamiya, Y. (1987). Proceedings of the 1987 IEEE Particle Accelerator Conference, 16-19 March 1987, Washington DC, USA, Vol. 1, pp. 437-439.]; Wiedemann, 1988[Wiedemann, H. (1988). Nucl. Instrum. Methods Phys. Res. A, 266, 24-31.]; Walker, 1995[Walker, R. P. (1995). CAS - CERN Accelerator School: 5th Advanced Accelerator Physics Course, 20 September-1 October 1993, Rhodes, Greece, pp. 807-835.]). This section presents specific details of the damping wiggler's contribution to synchrotron integrals.

Synchrotron integrals are described using η, k1 and [{1/\rho}]. Therefore, to obtain synchrotron integrals in the damping wiggler, the form of η, k1 and [{1/\rho}] in the damping wiggler should be realized. First, we obtained the dispersion function; assuming that the damping wiggler has a simple harmonic vertical field given by

[{B_y} = {B_{\rm{w}}}\sin{k_{\rm{w}}}z, \eqno(18)]

where Bw is the peak wiggler field and kw = [2\pi/{\lambda_{\rm{w}}}], with [{\lambda_{\rm{w}}}] being the wavelength of the damping wiggler, a Lorentz force [{\bf{F}}] is exerted on the electron in the damping wiggler as

[{\bf{F}} = q\,{\bf{v}} \times {\bf{B}}, \eqno(19)]

where [{\bf{F}}] is the force, q is the charge, [{\bf{v}}] is the velocity vector and [{\bf{B}}] is the magnetic field. The time variable is usually replaced by the position variable z along the longitudinal coordinate. With the constant-velocity approximation, the left-hand side of the equation (19)[link] becomes

[{\bf{F}} = {{{\rm{d}}{\bf{p}}}\over{{\rm{d}}t}} = {{{\rm{d}}\gamma{m}{\bf{v}}}\over{{\rm{d}}t}} \simeq {\gamma}mv_z^2\,{{{\rm{d}}^2{\bf{r}}}\over{{\rm{d}}{z^2}}} = {\gamma}mv_z^2x^{\prime\prime},]

where γ is the Lorentz factor, m is the mass, [{\bf{r}}] is a position vector, [x^{\prime}] is dx/dz and vz = dz/dt.

From now on, a prime denotes the derivative with respect to z. The right-hand side of equation (19)[link] becomes

[{F_x} = -e\,{\bf{v}}\times{\bf{B}} = e{v_z}{B_y},]

where e is the charge on an electron.

Assuming that the longitudinal momentum is much larger than the transverse momentum, then the horizontal equation of motion is

[\eqalignno{ x^{\prime\prime} & ={{e}\over{\gamma{m}v_z}}\,B_y \,\simeq\, {{e}\over{p}}\,B_y \,=\, {{e}\over{(1+\delta)\,p_0}}\,B_y \,\simeq\, (1-\delta)\,{{e}\over{p_0}}\,B_y \cr & = \,(1-\delta)\,{{B_y}\over{B_{\rm{w}}\,\rho_{\rm{w}}}} \quad \left(\,{{e}\over{p_0}}={{1}\over{B\rho}}\right), &(20)}]

where p is the magnitude of the momentum, p0 is the reference momentum, [B\rho] is the magnetic rigidity, and δ = p-p0)/p0.

Inserting equation (18)[link] into equation (20)[link] with the condition that δ is zero yields

[x_{\rm{ref}}^{\prime\prime} = {{1}\over{\rho_{\rm{w}}}}\sin{k_{\rm{w}}}z.\eqno(21)]

When δ = 0, x is a reference orbit. By solving the differential equation, x can be obtained as

[x_{\rm{ref}} = -{{1}\over{\rho_{\rm{w}}k_{\rm{w}}^2}} \,\sin{k_{\rm{w}}}z. \eqno(22)]

When δ ≠ 0, the orbit has an energy-dependent deviation from the reference orbit. The equation of motion is similar to equation (21)[link],

[x^{\prime\prime} = \left({1-\delta}\right)\,{1\over{{\rho_{\rm{w}}}}}\sin{k_{\rm{w}}}z. \eqno(23)]

Solving the differential equation (23)[link] yields

[x = {x_{\rm{ref}}} + \left(\,{{1\over{{\rho_{\rm{w}}}k_{\rm{w}}^2}} \,\sin{k_{\rm{w}}}z} \right)\,\delta. \eqno(24)]

The definition of the dispersion is given by

[x = {x_{\rm{ref}}} + \eta \delta.\eqno(25)]

Comparing equations (24)[link] and (25)[link], the dispersion function can be expressed as

[\eta = {1 \over {{\rho _{\rm{w}}}k_{\rm{w}}^2}}\,\sin {k_{\rm{w}}}z. \eqno(26)]

The normalized quadrupole strength k1 should be expressed with z. The definition of k1 is

[{k_1} = {1\over{{B_{\rm{w}}}B\rho}} \, {{{\rm{d}}{B_y}} \over {{\rm{d}}\xi}}. \eqno(27)]

The slope of ξ is perpendicular to the slope of the tangent line. The slope of the tangent line to the curve is dx/dz = [x^{\prime}]. The product of [x^{\prime}] and [{{{\rm{d}}\xi}/{{\rm{d}}z}}] is −1 ([{{{\rm{d}}\xi}/{{\rm{d}}z}} \times x^{\prime}] = −1), because they are perpendicular to each other. Then equation (27)[link] becomes

[{k_1} = {1\over{{B_{\rm{w}}}{\rho_{\rm{w}}}}} {{{{{\rm{d}}{B_y}}/{{\rm{d}}z}}}\over{{{{\rm{d}}\xi}/{{\rm{d}}z}}}} = {1\over{{B_{\rm{w}}}{\rho_{\rm{w}}}}} {{{{{\rm{d}}{B_y}}/{{\rm{d}}z}}}\over{-{1/{x^{\prime}}}}} = -{1\over {{B_{\rm{w}}}{\rho_{\rm{w}}}}}\,x^{\prime}\,{{{\rm{d}}{B_y}}\over{{\rm{d}}z}}. \eqno(28)]

Inserting equations (18)[link] and (22)[link] into (28)[link] yields k1 as

[{k_1} = {1\over{\rho_{\rm{w}}^2}} \,\cos^2k_{\rm{w}}z. \eqno(29)]

Lastly, the bending radius is rewritten with the variable z. It is related to the magnetic field as

[{1 \over \rho } = {{{B_y}} \over {{B_{\rm{w}}}{\rho _{\rm{w}}}}} = {{{B_{\rm{w}}}\sin {k_{\rm{w}}}z} \over {{B_{\rm{w}}}{\rho _{\rm{w}}}}} = {1 \over {{\rho _{\rm{w}}}}}\,\sin {k_{\rm{w}}}z. \eqno(30)]

The damping wiggler contribution to a synchrotron integral will be denoted by the sub-index wig. The integration is from the entrance of the wiggler to the exit of the wiggler, i.e. for wiggler Lw. By inserting equations (26)[link], (29)[link] and (30)[link] into the synchrotron integrals, the contributions of the damping wiggler to the synchrotron integrals can be obtained as

[{I}_{1,{\rm{wig}}} = \int_{0}^{{L}_{\rm{w}}}{\rm{d}}z \, {{\eta }\over{\rho}} = {{{\lambda}_{\rm{w}}^{2}{B}_{\rm{w}}^{2}{L}_{\rm{w}}}\over{8{\pi}^{2}{\left(B\rho\right)}^{2}}},]

[Z{I}_{2,{\rm{wig}}} = \int_{0}^{{L}_{\rm{w}}}\,{\rm{d}}z \, {{1}\over{{\rho}^{2}}} = {{{B}_{\rm{w}}^{2}{L}_{\rm{w}}}\over{2{\left(B\rho\right)}^{2}}},]

[{I}_{3,{\rm{wig}}} = \int_{0}^{{L}_{\rm{w}}}{\rm{d}}z \, {{1}\over{\left|{\rho}^{3}\right|}} = {{4{B}_{\rm{w}}^{3}{L}_{\rm{w}}}\over{3\pi{\left(B\rho\right)}^{3}}},]

[{I}_{4,{\rm{wig}}} = \int_{0}^{{L}_{\rm{w}}} {\rm{d}}z \, {{\eta}\over{\rho}} \left(\,{{1}\over{{\rho}^{2}}}+2{k}_{1}\right) = {{5{\lambda}_{\rm{w}}^{2}{B}_{\rm{w}}^{\,4}{L}_{\rm{w}}}\over{32{\pi}^{2}{\left(B\rho\right)}^{4}}},]

where [{\beta_x}] is the average value of the horizontal beta function from z = 0 to z = Lw and [{\alpha_x}], [{\beta_x}] and [{\gamma_x}] are the Twiss functions. The calculation of I5,wig uses an approximation

[{\gamma_x}{\eta^2} + 2{\alpha_x}\eta\eta^{\prime} + {\beta_x}\eta^{\prime2} \simeq {\beta_x}\eta^{\prime2}.]

Usually [{\eta/{\eta^{\prime}}}]1/kw [\ll] 1 m. For example, when the wiggler wavelength is 12.57 cm, 1/kw = 5 ×10 - 3 m. An assumption is also used: that [{\beta_x}] changes smoothly in the damping wiggler, so the integral of [{{{\beta_x}\eta^{\prime2}}/{|{{\rho^3}}|}}] is approximately the same as the product of [{\beta_x}] and the integrals of [{{\eta^{\prime2}}/{|{{\rho^3}}|}}]. Generally, the synchrotron radiation integrals are defined in a curve-linear coordinate. Here the synchrotron radiation integrals are integrated with z; they cannot represent the ring parameters related to path length. For example, I1,wig does not represent the contribution of the damping wiggler to a momentum compaction factor. Then the new synchrotron radiation integrals

[\eqalign{ {I}_{1\,{\rm{new,wig}}} & = \int_{0}^{{L}_{\rm{w}}} {\rm{d}}z \, {x_{\rm{ref}}^{\prime}} \eta^{\prime} = -{{ {\lambda}_{\rm{w}}^{2}{B}_{\rm{w}}^{\,2}{L}_{\rm{w}} }\over{ 8{\pi}^{2}{\left(B\rho\right)}^{2}}},\cr {I}_{4\,{\rm{new,wig}}} & = \int _{0}^{{L}_{\rm{w}}} {\rm{d}}z \, {{{x_{\rm{ref}}^{\prime}} \eta^{\prime}}\over{{\rho }^{2}}} = -{{3{\lambda}_{\rm{w}}^{2}{B}_{\rm{w}}^{\,4}{L}_{\rm{w}}}\over{32{\pi}^{2}{\left(B\rho\right)}^{4}}} }]

are defined to denote the ring parameters appropriately. With these new parameters, the damping wiggler contribution to the momentum compaction is expressed as [{\alpha_{c,{\rm{wig}}}} ] = I1,wig/C.

Acknowledgements

We thank H. Wiedemann for providing helpful information and many useful discussions.

Funding information

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF-2019R1C1C1003412) and the Basic Science Research Program through the National Research Foundation of Korea (NRF-2019R1A2C1004862).

References

First citationAPS (2020). APS Upgrade, https://www.aps.anl.gov/APS-UpgradeGoogle Scholar
First citationBjorken, J. D. & Mtingwa, S. K. (1983). Part. Accel. 13, 115–143.  CAS Google Scholar
First citationBorland, M. (2000). Elegant: A Flexible SDDS-Compliant Code for Accelerator Simulation. Technical Report LS-287. Advanced Photon Source, Argonne, IL, USA.  Google Scholar
First citationBorland, M., Decker, G., Emery, L., Sajaev, V., Sun, Y. & Xiao, A. (2014). J. Synchrotron Rad. 21, 912–936.  Web of Science CrossRef IUCr Journals Google Scholar
First citationCai, Y., Ding, Y., Hettel, R., Huang, Z., Wang, L. & Xiao, L. (2013). An X-ray Free Electron Laser Driven by an Ultimate Storage Ring. Report SLAC-PUB-15380. SLAC National Accelerator Laboratory, Menlo Park, California, USA.  Google Scholar
First citationDeb, K. (2001). Multi-Objective Optimization Using Evolutionary Algorithms. Chichester: John Wiley & Sons.  Google Scholar
First citationHelm, R. H., Lee, M. J., Morton, P. L. & Sands, M. (1973). IEEE Trans. Nucl. Sci. 20, 900–901.  CrossRef Web of Science Google Scholar
First citationHuang, Z., Ding, Y. & Schroeder, C. B. (2012). Phys. Rev. Lett. 109, 204801.  Web of Science CrossRef PubMed Google Scholar
First citationKatoh, M. & Kamiya, Y. (1987). Proceedings of the 1987 IEEE Particle Accelerator Conference, 16–19 March 1987, Washington DC, USA, Vol. 1, pp. 437–439.  Google Scholar
First citationLiu, L., Milas, N., Mukai, A. H. C., Resende, X. R., Rodrigues, A. R. D. & de Sá, F. H. (2013). Proceedings of the Fourth International Particle Accelerator Conference (IPAC2013), 12–17 May 2013, Shanghai, China, pp. 1874–1876. TUPWO001.  Google Scholar
First citationLuong, V. (2018). EUV Lithography, https://eng.kuleuven.be/en/research/phd/arenberg-youngster-seminars/ays-pdf/20180516-luong-ays-final.pdfGoogle Scholar
First citationMAX IV (2010). MAX IV DDR – The MAX IV Detailed Design Report, https://www.maxiv.lu.se/accelerators-beamlines/accelerators/accelerator-documentation/max-iv-ddr/Google Scholar
First citationRevol, J.-L. (2017). ESRF-EBS: The Extremely Brilliant Source Project, https://indico.psi.ch/conferenceDisplay.py?confId=5589Google Scholar
First citationRiemann, B. & Streun, A. (2019). Phys. Rev. Accel. Beams, 22, 021601.  CrossRef Google Scholar
First citationShin, S., Kwon, S., Kim, D.., Kim, D.., Kim, M., Kim, S.., Kim, S.., Kim, J., Kim, C., Park, B., Park, S.., Park, S.., Park, E., Son, Y., Yoon, J., Lee, B., Lee, E., Lee, J.., Lee, H.., Joo, Y., Choi, J., Ha, T., Hwang, W., Hwang, I., Lee, J.., Oh, B., Lee, C.., Lee, H.., Kim, J.., Hwang, J. Y., Nam, S. H. & Cho, M. (2013). J. Instrum. 8, P01019.  Google Scholar
First citationSPring-8 (2014). SPring-8-II Conceptual Design Report, https://rsc.riken.jp/eng/pdf/SPring-8-II.pdfGoogle Scholar
First citationSteier, C., Byrd, J. M., De Santis, S., Nishimura, H., Robin, D., Sannibale, F., Sun, C., Venturini, M. & Wan, W. (2016). Proceedings of the 7th International Particle Accelerator Conference (IPAC2016), 8–13 May 2016, Busan, Korea, pp. 2956–2958. WEPOW049.  Google Scholar
First citationWalker, R. P. (1995). CAS – CERN Accelerator School: 5th Advanced Accelerator Physics Course, 20 September–1 October 1993, Rhodes, Greece, pp. 807–835.  Google Scholar
First citationWiedemann, H. (1988). Nucl. Instrum. Methods Phys. Res. A, 266, 24–31.  CrossRef Web of Science Google Scholar
First citationWiedemann, H. (1993). Particle Accelerator Physics I. Berlin: Springer.  Google Scholar

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds