Scheme for generating 1 nm X-ray beams carrying orbital angular momentum at the SXFEL

Scheme for generating 1 nm X-ray beams carrying orbital angular momentum at the SXFEL 1 2 first-author https://orcid.org/0000-0001-8663-7528 3 corresp Corresponding author, chenjianhui@zjlab.org.cn Corresponding author, chenjianhui@zjlab.org.cn chenjianhui@zjlab.org.cn 3 Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China University of Chinese Academy of Sciences, Beijing 100049, China Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201210, China Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201210, China Optical vortices have the main features of helical wavefronts and spiral phase structures, and carry orbital angular momentum (OAM). This special structure of visible light has been produced and studied for various applications. These notable characteristics of photons were also tested in the extreme-ultraviolet and X-ray regimes. In this article, we simulate the use of a simple afterburner configuration by directly adding helical undulators after the SASE undulators with the Shanghai Soft X-ray FEL to generate high intensity X-ray vortices with wavelengths ∼ 1 nm. Compared to other methods, this approach is easier to implement, cost-effective, and more efficient. X-ray Orbital angular momentum (OAM) Synchrotron light source Free-electron laser (FEL) 1 Introduction 1 Structured light, whose phase structures rotate when they propagate through space, has been a topic of intense research. With the development of tools and technologies for creating and detecting structured light, applications have begun to appear steadily [ 1 1 ]. Among these studies, the optics and laser communities are interested in phase singularity[ 2 2 ] and visible photons carrying orbital angular momentum (OAM)[ 3 3 ], to provide new explanations of optical effects and develop distinct applications ranging from manipulations of small particles, driving micro-machines, and quantum computation, to optical communications[ 4 4 - 14 14 4-14 ]. The noteworthy attributes of photons were also demonstrated in the extreme-ultraviolet and X-ray or ever γ ray[ 15 15 ] regions, pioneered by Peele and co-workers[ 16 16 , 17 17 ]. Triggered by this pioneering effort, various generation methods based on particle accelerators[ 18 18 - 25 25 18-25 ], high-power lasers[ 26 26 - 29 29 26-29 ], and plasma[ 30 30 , 31 31 ], as well as possible novel applications[ 32 32 - 38 38 32-38 ], are growing steadily. Undulators are the workhorses in third-generation synchrotron radiation light sources and free-electron lasers (FELs) for generating X-ray beams. Sasaki and McNulty[ 39 39 ] were the first to indicate that when electrons traveled along a spiral path in helical undulators, the harmonics produced by them carried a spiral phase. Therefore, it is proposed that the helical undulator could be used as a source for generating intense X-rays carrying OAM light (XOAM). Bahrdt and coworkers at BESSY-II[ 40 40 ] first measured the special helical structure in spontaneous undulator emission. This was later measured in the XUV (X-ray and extreme ultraviolet) region at UVSOR-III[ 41 41 ] by Tatsuo and his coworkers, and in coherent radiation emitted from energy and density modulated (bunched) electron beams[ 42 42 ] by the US group at SLAC. In 2014, Hemsing and coworkers described a simpler afterburner configuration based on the HGHG setup to generate OAM light[ 43 43 ]. Based on Hemsing’s HGHG scheme, Primoz Rebernik Ribic et al. in 2014 proposed the use of an optical phase mask to produce a transverse phase modulation in the incident seed laser, which could also produce an XUV or X-ray optical vortex[ 44 44 ]. Three years later, in 2017, they demonstrated two techniques to obtain intense femtosecond coherent optical vortices at the FERMI FEL in the extreme ultraviolet[ 45 45 ]. The first method used the helical undulator to produce spiral beams at the second harmonic, while the second method relied on a spiral zone plate by putting it directly into the propagation path of the transverse Gaussian optical beam. In this study, we simulate the use of a simple afterburner configuration by directly adding helical undulators after the SASE undulators at the Shanghai Soft X-ray FEL (SXFEL)[ 46 46 ] to generate high intensity X-ray vortices with wavelengths of ∼ 1 nm. This configuration has two distinct advantages. First, the setup is very simple, since it only requires the addition of a few more helical undulators to the existing polarization control configuration. Second, it does not need special devices to produce and preserve the spiral bunched e beam structure during e beam transport, compared to previous methods. 2 Theory and modeling of higher harmonics from helical undulators 1 As mentioned earlier, X-ray photons carrying OAM can be generated in particle accelerators by converting an X-ray beam[ 18 18 - 23 23 18-23 ] or using helical undulators[ 39 39 ]. This triggers a careful re-examination of the angular momentum properties of emission from charged particles in circular and spiral motion[ 47 47 - 51 51 47-51 ]. We know that there are different kinds of undulator magnetic fields and possibly various structures of external lasers in diverse undulator devices. When the relativistic electron beam ( e beam) is injected into the undulator, its energy and the eventual density distribution are affected by the magnetic and laser fields. In other words, the energy and final density distribution of the relativistic e beam are determined by the three-dimensional structure of the resonant interaction between the e beam, undulator field, and external injected laser. After the resonant interaction, the microbunched e beam is produced. These microbunched e beams can then be used to radiate electromagnetic waves in different downstream radiation facilities. Special light beams can be generated by resonant coupling with higher order harmonics. Therefore, it may be used to radiate fields with special 3D structures by properly tuning the injected microbunching and harmonic resonances in the downstream radiator. The helical light beams carrying OAM, generated in Ref.([ 42 42 ]), arise from the transmission of the structural information of the helical microbunched e beams to the optical field. We use the standard linearized FEL equation to describe the interaction between the injected e beams and input laser in the planar undulator. After some discussion, the density perturbation n 1 ( r ) can be expressed as the sum of the expanded eigenmodes. The density modulation expansion is[ 52 52 ]: <math display="block" id="M1"><mrow><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>k</mi><mi>w</mi></msub><msub><mi>ϵ</mi><mn>0</mn></msub></mrow><mi>e</mi></mfrac><mstyle displaystyle="true"><munder><mo>∑</mo><mi>j</mi></munder></mstyle><msub><mi>a</mi><mi>j</mi></msub><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><msub><mover accent="true"><mi>ε</mi><mo>˜</mo></mover><mrow><mo>⊥</mo><mi>j</mi></mrow></msub><mo stretchy="false">(</mo><msub><mstyle mathvariant="bold" mathsize="normal"><mi>r</mi></mstyle><mo>⊥</mo></msub><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where k w is the wave number (the undulator wavelength is λ w =2 π / k w ) and ϵ 0 is the vacuum dielectric constant. The evolution of the harmonic amplitudes can be expressed as[ 52 52 ]: <math display="block" id="M2"><mrow><mfrac><mrow><msup><mtext>d</mtext><mn>2</mn></msup></mrow><mrow><mtext>d</mtext><msup><mi>z</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>a</mi><mi>q</mi></msub><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><mo>+</mo><msubsup><mi>θ</mi><mi>p</mi><mn>2</mn></msubsup><mstyle displaystyle="true"><munder><mo>∑</mo><mi>j</mi></munder></mstyle><msub><mi>F</mi><mrow><mi>q</mi><mn>,</mn><mi>j</mi></mrow></msub><msub><mi>a</mi><mi>j</mi></msub><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><mo>=</mo><mo>−</mo><mstyle displaystyle="true"><munder><mo>∑</mo><msup><mi>q</mi><mo>′</mo></msup></munder></mstyle><msubsup><mi>D</mi><mrow><mi>q</mi><mn>,</mn><msup><mi>q</mi><mo>′</mo></msup></mrow><mrow><mo stretchy="false">(</mo><mi>h</mi><mo stretchy="false">)</mo></mrow></msubsup><msub><mi>C</mi><msup><mi>q</mi><mo>′</mo></msup></msub><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>i</mi><msubsup><mi>θ</mi><msup><mi>q</mi><mo>′</mo></msup><mrow><mo stretchy="false">(</mo><mi>h</mi><mo stretchy="false">)</mo></mrow></msubsup><mi>z</mi></mrow></msup><mn>,</mn></mrow></math> where θ p is the relativistic longitudinal plasma wave number and a q (0), d a q (0)/d z , and C q (0) indicate the bunching, velocity modulation, and amplitude, respectively, of the incident field in the initial case. The second term on the left-hand side of equation(2) comes from the longitudinal space charge. The coefficient <math id="M3"><mrow><msubsup><mi>D</mi><mrow><mi>q</mi><mn>,</mn><msup><mi>q</mi><mo>′</mo></msup></mrow><mrow><mo stretchy="false">(</mo><mi>h</mi><mo stretchy="false">)</mo></mrow></msubsup></mrow></math> [ 52 52 ]on the right-hand side gives the coupling condition between the electron beams and harmonic fields. <math display="block" id="M4"><mrow><msubsup><mi>D</mi><mrow><mi>q</mi><mn>,</mn><msup><mi>q</mi><mo>′</mo></msup></mrow><mrow><mo stretchy="false">(</mo><mi>h</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>=</mo><msub><mover accent="true"><mi>g</mi><mo>^</mo></mover><mo>⊥</mo></msub><mfrac><mrow><msubsup><mi>θ</mi><mi>p</mi><mn>2</mn></msubsup><mi>c</mi><mi>K</mi><msup><mrow><mrow><mo>[</mo><mrow><msub><mi>k</mi><mrow><mi>z</mi><msup><mi>q</mi><mo>′</mo></msup></mrow></msub><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo><mo>+</mo><msub><mi>k</mi><mi>w</mi></msub></mrow><mo>]</mo></mrow></mrow><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>γ</mi><mi>ω</mi><msub><mi>k</mi><mi>w</mi></msub></mrow></mfrac><msubsup><mi>F</mi><mrow><mi>q</mi><mn>,</mn><msup><mi>q</mi><mo>′</mo></msup></mrow><mrow><mo stretchy="false">(</mo><mi>h</mi><mo stretchy="false">)</mo></mrow></msubsup><mn>,</mn></mrow></math> where c is the speed of light, K is the undulator parameter, ω is the frequency, <math id="M5"><mrow><msub><mover accent="true"><mi>g</mi><mo>^</mo></mover><mo>⊥</mo></msub></mrow></math> denotes the electron motion, <math id="M6"><mrow><msubsup><mi>F</mi><mrow><mi>q</mi><mn>,</mn><msup><mi>q</mi><mo>′</mo></msup></mrow><mrow><mo stretchy="false">(</mo><mi>h</mi><mo stretchy="false">)</mo></mrow></msubsup></mrow></math> indicates the overlay effect of the e beam and incident laser lights, and f ( r ⊥ ) represents the transverse density distribution of electron beams. After interacting with a typical harmonic field incident mode, the e beam modulation in equation (1) can be expressed as follows[ 52 52 ]: <math display="block" id="M7"><mrow><msub><mover accent="true"><mi>ε</mi><mo>˜</mo></mover><mrow><mo>⊥</mo><mi>q</mi></mrow></msub><mo stretchy="false">(</mo><msub><mi>r</mi><mo>⊥</mo></msub><mo stretchy="false">)</mo><mo>=</mo><msubsup><mi>R</mi><mi>p</mi><mi>l</mi></msubsup><mo stretchy="false">(</mo><mi>r</mi><mo stretchy="false">)</mo><msup><mi>e</mi><mrow><mi>i</mi><mi>l</mi><mi>ϕ</mi></mrow></msup><mn>.</mn></mrow></math> In the above formula, the mode index q includes two coefficients, q =( p , l ), where p and l correspond to the radial and azimuthal modes, respectively. When the e beam distribution is axisymmetric f ( r ⊥ )= f ( r ), the overlay effect of the e beam and incident fields can be straightforward[ 52 52 ] <math display="block" id="M8"><mrow><msubsup><mi>F</mi><mrow><mo stretchy="false">(</mo><mi>p</mi><mn>,</mn><mi>l</mi><mo stretchy="false">)</mo><mn>,</mn><mo stretchy="false">(</mo><msup><mi>p</mi><mo>′</mo></msup><mn>,</mn><msup><mi>l</mi><mo>′</mo></msup><mo stretchy="false">)</mo></mrow><mrow><mo stretchy="false">(</mo><mi>h</mi><mo stretchy="false">)</mo></mrow></msubsup><mo>∝</mo><mn>2</mn><mi>π</mi><msub><mi>δ</mi><mrow><msup><mi>l</mi><mo>′</mo></msup><mn>,</mn><mi>l</mi><mo>±</mo><mo stretchy="false">(</mo><mi>h</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></msub><mn>.</mn></mrow></math> We can see the relationship between l , l’ , h and the polarization direction of the undulator, ±, where l and l’ represent the radial modes about the azimuth of the electron beams and electromagnetic fields. There are some simple cases. For example, when an incident electromagnetic wave in a random l’ mode is injected into the undulator in the fundamental mode ( h =1), the homologous azimuthal modes can be obtained in the electron beams l = l’ . When an axisymmetric electromagnetic wave l’ =0 is injected into the undulator working at the second harmonic mode ( h =2), a spiral bunched electron beam is formed, where the helical undulator can be either right-handed ( l =-1) or left-handed ( l =1). Therefore, corkscrewed micro-bunched electron beams at the second harmonic can be produced by the interaction between an initially unmodulated electron beam and incident Gaussian electromagnetic field (laser field). We now discuss the properties of emission from helical undulators. The time-dependent complex field amplitude A of radiation from a pure circular undulator is[ 39 39 ]: <math display="block" id="M9"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mi>A</mi></mtd><mtd columnalign="left"><mrow><mo>=</mo><mo stretchy="false">(</mo><msub><mi>A</mi><mi>x</mi></msub><mo>−</mo><mi>i</mi><msub><mi>A</mi><mi>y</mi></msub><mo stretchy="false">)</mo><mo>/</mo><msqrt><mn>2</mn></msqrt></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>=</mo><msqrt><mn>2</mn></msqrt><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>exp</mi></mrow></mstyle><mrow><mi>i</mi><mo stretchy="false">(</mo><mi>n</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo><mi>ϕ</mi></mrow></msup><mrow><mo>[</mo><mrow><mrow><mo>(</mo><mrow><mi>γ</mi><mi>θ</mi><mo>−</mo><mfrac><mrow><mi>n</mi><mi>K</mi></mrow><mi>X</mi></mfrac></mrow><mo>)</mo></mrow><msub><mi>J</mi><mi>n</mi></msub><mo stretchy="false">(</mo><mi>X</mi><mo stretchy="false">)</mo><mo>−</mo><mi>K</mi><msub><msup><mi>J</mi><mo>′</mo></msup><mi>n</mi></msub><mo stretchy="false">(</mo><mi>X</mi><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> Here, X =2 nξγθK , ξ=1/(1+ γ 2 θ 2 + K 2 ), K is the undulator deflecting parameter, θ is the viewing angle of the propagation axis, ϕ is the azimuth of the observer’s location in the transverse plane perpendicular to the direction of propagation, γ is the Lorentz factor, J n is the n th order Bessel function of the first kind, and <math id="M10"><mrow><msub><msup><mi>J</mi><mo>′</mo></msup><mi>n</mi></msub><mo stretchy="false">(</mo><mi>X</mi><mo stretchy="false">)</mo></mrow></math> is the first derivative of J n with respect to X . It is clear that the above equation is typical of Laguerre-Gaussian (LG) modes with a topological charge of l = n -1. LG beams are a group of solutions of the paraxial wave equation in the cylindrical coordinate system, and they are characterized by the OAM of a photon with lℏ . The LG modes are expressed by[ 3 3 ]: <math display="block" id="M11"><mrow><msub><mi>u</mi><mrow><mtext>LG</mtext></mrow></msub><msub><mrow><mo stretchy="false">(</mo><mi>r</mi><mn>,</mn><mi>θ</mi><mn>,</mn><mi>z</mi><mo stretchy="false">)</mo></mrow><mrow><mi>p</mi><mn>,</mn><mi>l</mi></mrow></msub><mo>∝</mo><msqrt><mrow><msub><mi>I</mi><mrow><mi>p</mi><mn>,</mn><mi>l</mi></mrow></msub><mo stretchy="false">(</mo><mi>r</mi><mn>,</mn><mi>z</mi><mo stretchy="false">)</mo></mrow></msqrt><mi>exp</mi><mrow><mo>[</mo><mrow><mi>i</mi><mo stretchy="false">(</mo><mn>2</mn><mi>p</mi><mo>+</mo><mi>l</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><mi>tan</mi><mi>z</mi><mo>/</mo><msub><mi>z</mi><mi>R</mi></msub><msup><mrow/><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>]</mo></mrow><mi>exp</mi><mo stretchy="false">(</mo><mo>−</mo><mi>i</mi><mi>l</mi><mi>θ</mi><mo stretchy="false">)</mo><mn>.</mn></mrow></math> Similar to the above example, the last term on the right in equation (7) shows that there is an azimuthal component of the phase, where p is the radial index which decides the radial intensity distributions and l is the azimuthal index. The wavefront (or the phase structure) is a spiral surface when l ≠0. As an example, we perform a simulation using the current SSRF (Shanghai Synchrotron Radiation Facility)[ 53 53 - 55 55 53-55 ] beam parameters (3.5 GeV beam energy, 3 nm-rad horizontal emittance, and 1% coupling) and an actual magnetic field ( K∼ 0.78) of the 32-mm-period helical undulator, whose total length is 4.5 m. The radiation properties of such helical undulators are simulated using the SPECTRA code[ 56 56 ]. The trajectory of an electron moving inside the undulator is a spiral along the axis in the direction of propagation. The radiation produced by the electron is emitted along the direction of propagation within a very small divergence angle. The spectral flux density of the light from the helical undulator is shown in Fig. 1 Fig. 1 , which indicates that the fundamental and second harmonic modes are at 2256 eV and about 4413 eV within two orders of magnitude, respectively. We then calculate the complex amplitude spatial profile of the radiation field on the cross section perpendicular to the direction of propagation, and estimate the phases of the first and second harmonic radiation using the real and imaginary parts of E x and E y , which are shown in Fig. 2 Fig. 2 . We observe that the transverse flux density profile of the fundamental radiation is Gaussian, while the transverse flux density distribution of the second harmonic is nearly annular. The typical transverse phase distribution of the second harmonic light is a spiral surface along the direction of propagation, which indicates that the photons carry an OAM (in this case, a topological charge l =-1). It turns out that a relativistic e beam will move spirally along the axis of propagation after entering the helical undulator. Owing to this special motion mode, the light produced by the moving electrons forms a spiral structure at higher harmonics indicating that the photons carry OAM. 3 Theory of Self-Amplified Spontaneous Emission in Free-Electron Lasers 1 A revolutionary development in X-ray production was the progress of XFEL[ 57 57 ]. When an electron beam enters a periodic magnetic structure, the coherent radiation is amplified and this is the basis of the FEL. The radiation generated in an XFEL interacts with the beam of electrons in return to get more excellent X-ray than light sources based on storage rings in terms of intensity and coherence. At present, the high-energy e beams produced by a linac can be used to drive high-gain XFEL amplifiers using a long undulator. The gain of the FEL amplifier can be high enough to allow the initial incoherent undulator radiation to evolve into a high intensity quasi-coherent electromagnetic field——SASE[ 58 58 ]. In principle, an FEL can work at any wavelength, limited only by the energy and quality of the e beams obtained by the accelerator, and thus, can be used to fill gaps in the spectrum that cannot be covered by other coherent sources. At present, FEL oscillators based on electrostatic accelerators, cyclotrons, and rf-linear accelerators have been built. FEL oscillators in the storage rings are suitable for producing coherent radiation from visible to ultraviolet wavelengths[ 58 58 ]. The main technique for producing coherent radiation at shorter wavelengths involves a high gain FEL that does not require a reflector. In the late 1990s and early 2000s, devices were first designed and built to produce high-intensity ultraviolet and vacuum ultraviolet radiation. In 2005, the FLASH installed in Hamburg, Germany, began offering soft XFEL to scientific users[ 59 59 ]. In 2009, the Linac Coherent Light Source(LCLS) at SLAC ushered in the era of hard XFEL[ 60 60 ]. Subsequently, the XFEL devices in Japan(SACLA)[ 61 61 ] and Italy(FERMI)[ 71 71 ] were successfully debugged and provided services to users. At present, the soft XFEL device in Shanghai, China has been built and successfully tested[ 46 46 ]. We also mention the European XFEL in Germany[ 62 62 ], PAL XFEL at Posco in South Korea[ 63 63 ], and SwissFEL at PSI in Switzerland[ 64 64 ]. An FEL can operate in several modes like any laser, such as the FEL oscillator, amplifier, and SASE. The oscillator uses reflectors to constrain the radiation so that the radiation field passes through the undulator many times. An FEL can also function as a linear amplifier to amplify radiation with a central frequency close to the resonant condition. Owing to the need for seed light of the corresponding wavelength, this mode of operation is most easily implemented in the spectral region already covered by other types of light sources. FEL amplifiers extending to shorter wavelengths often require advanced techniques such as harmonic generation, so it is challenging to apply them to X-ray wavelengths. It is for this reason that the first scale and built XFEL devices were based on the SASE mode[ 58 58 ]. SASE originates from the growth and saturation of spontaneous emission in a long undulator, without the reflector and oscillator. Although the SASE radiation is not completely coherent, the gain process increases the brightness of the source to about 10 orders of magnitude higher than that of traditional synchrotron sources. In order to learn more about SASE, we consider the radiation field of the FEL process given by[ 58 58 ] <math display="block" id="M12"><mrow><mi>a</mi><mo stretchy="false">(</mo><mover accent="true"><mi>z</mi><mo>^</mo></mover><mo stretchy="false">)</mo><mo>≈</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mrow><mo>[</mo><mrow><mi>a</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo><mo>−</mo><mi>i</mi><mfrac><mrow><mi>b</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow><mrow><msub><mi>μ</mi><mn>3</mn></msub></mrow></mfrac><mo>−</mo><mi>i</mi><msub><mi>μ</mi><mn>3</mn></msub><mi>P</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow><mo>]</mo></mrow><msup><mi>e</mi><mrow><mo>−</mo><mi>i</mi><msub><mi>μ</mi><mn>3</mn></msub><mover accent="true"><mi>z</mi><mo>^</mo></mover></mrow></msup><mn>.</mn></mrow></math> The equation(8) is approximately obtained considering a long propagation distance, where <math id="M13"><mrow><mover accent="true"><mi>z</mi><mo>^</mo></mover><mo>≫</mo><mn>1</mn></mrow></math> . The first term in square brackets of the formula(8) represents the coherent amplification of the external radiation signal, while the second and third terms in square brackets represent the FEL radiation caused by density and energy modulation, respectively, of the electron beams. When these modulations occur due to the shot noises of e beams, the exponential growth is known as SASE. In equation (8), where the normalized ordinate is <math id="M14"><mrow><mover accent="true"><mi>z</mi><mo>^</mo></mover><mo>≡</mo><mn>2</mn><msub><mi>k</mi><mi>u</mi></msub><mi>ρ</mi><mi>z</mi></mrow></math> , and ρ is the nondimensional Pierce parameter[ 65 65 ], <math display="block" id="M15"><mrow><mi>ρ</mi><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><mfrac><mrow><msup><mi>e</mi><mn>2</mn></msup><msup><mi>K</mi><mn>2</mn></msup><msup><mrow><mrow><mo>[</mo><mrow><mi>J</mi><mi>J</mi></mrow><mo>]</mo></mrow></mrow><mn>2</mn></msup><msub><mi>n</mi><mi>e</mi></msub></mrow><mrow><mn>32</mn><msub><mi>ϵ</mi><mn>0</mn></msub><msubsup><mi>γ</mi><mi>r</mi><mn>3</mn></msubsup><mi>m</mi><msup><mi>c</mi><mn>2</mn></msup><msubsup><mi>k</mi><mi>u</mi><mn>2</mn></msubsup></mrow></mfrac></mrow><mo>]</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mn>.</mn></mrow></math> Here, the factor [ JJ ] comes from the longitudinal oscillation of the electrons in the plane undulator, which alters the average coupling between the particle and radiation field[ 58 58 ]. <math display="block" id="M16"><mrow><mrow><mo>[</mo><mrow><mi>J</mi><mi>J</mi></mrow><mo>]</mo></mrow><mo>≡</mo><msub><mi>J</mi><mn>0</mn></msub><mo stretchy="false">(</mo><mfrac><mrow><msup><mi>K</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><mo>+</mo><mn>2</mn><msup><mi>K</mi><mn>2</mn></msup></mrow></mfrac><mo stretchy="false">)</mo><mo>−</mo><msub><mi>J</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mfrac><mrow><msup><mi>K</mi><mn>2</mn></msup></mrow><mrow><mn>4</mn><mo>+</mo><mn>2</mn><msup><mi>K</mi><mn>2</mn></msup></mrow></mfrac><mo stretchy="false">)</mo><mn>.</mn></mrow></math> The parameters μ 1 , μ 2 , and μ 3 [ 66 66 ]are <math display="block" id="M17"><mrow><msub><mi>μ</mi><mn>1</mn></msub><mo>=</mo><mn>1,</mn><msub><mi>μ</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mo>−</mo><mn>1</mn><mo>−</mo><msqrt><mn>3</mn></msqrt><mi>i</mi></mrow><mn>2</mn></mfrac><mn>,</mn><msub><mi>μ</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><mo>−</mo><mn>1</mn><mo>+</mo><msqrt><mn>3</mn></msqrt><mi>i</mi></mrow><mn>2</mn></mfrac><mn>.</mn></mrow></math> We then define three parameters to analyze the characteristics of SASE radiation[ 58 58 ]. <math display="block" id="M18"><mrow><mfrac><mrow><mtext>d</mtext><mi>a</mi></mrow><mrow><mtext>d</mtext><mover accent="true"><mi>z</mi><mo>^</mo></mover></mrow></mfrac><mo>=</mo><mo>−</mo><mi>b</mi></mrow></math> <math display="block" id="M19"><mrow><mfrac><mrow><mtext>d</mtext><mi>b</mi></mrow><mrow><mtext>d</mtext><mover accent="true"><mi>z</mi><mo>^</mo></mover></mrow></mfrac><mo>=</mo><mo>−</mo><mi>i</mi><mi>P</mi></mrow></math> <math display="block" id="M20"><mrow><mfrac><mrow><mtext>d</mtext><mi>P</mi></mrow><mrow><mtext>d</mtext><mover accent="true"><mi>z</mi><mo>^</mo></mover></mrow></mfrac><mo>=</mo><mi>a</mi></mrow></math> where a is the radiation field which represents the coherent radiation drives energy modulation, b is the cluster factor which reflects the cluster- generated coherent radiation, and P is the collective momentum which reflects the transformation from energy modulation to density modulation. The resonance condition of the undulator is <math display="block" id="M21"><mrow><mfrac><mrow><msub><mi>k</mi><mi>u</mi></msub></mrow><mi>k</mi></mfrac><mo>=</mo><mfrac><mi>λ</mi><mrow><msub><mi>λ</mi><mi>u</mi></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><mn>1</mn><mo>+</mo><msup><mi>K</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mrow><mrow><mn>2</mn><msup><mi>γ</mi><mn>2</mn></msup></mrow></mfrac><mn>.</mn></mrow></math> For convenience, we assume that there is no external field( a (0)=0) and the divergence of beam energy is 0 ( P (0)=0). Thus, in the exponential growth range, the radiation field is[ 58 58 ] <math display="block" id="M22"><mrow><mrow><mo>〈</mo><mrow><mo>|</mo><mi>a</mi><mo stretchy="false">(</mo><mover accent="true"><mi>z</mi><mo>^</mo></mover><msup><mrow><mo stretchy="false">)</mo><mo>|</mo></mrow><mn>2</mn></msup></mrow><mo>〉</mo></mrow><mo>≈</mo><mfrac><mn>1</mn><mn>9</mn></mfrac><mrow><mo>〈</mo><mrow><mo>|</mo><mi>b</mi><msup><mrow><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo><mo>|</mo></mrow><mn>2</mn></msup></mrow><mo>〉</mo></mrow><msup><mi>e</mi><mrow><msqrt><mn>3</mn></msqrt><mover accent="true"><mi>z</mi><mo>^</mo></mover></mrow></msup><mn>.</mn></mrow></math> Here, by normalizing the propagation distance <math id="M23"><mrow><msqrt><mn>3</mn></msqrt><mover accent="true"><mi>z</mi><mo>^</mo></mover><mo>=</mo><msqrt><mn>3</mn></msqrt><mo stretchy="false">(</mo><mn>2</mn><msub><mi>k</mi><mi>u</mi></msub><mi>z</mi><mi>ρ</mi><mo stretchy="false">)</mo><mo>=</mo><mi>z</mi><mo>/</mo><msub><mi>L</mi><mrow><mi>G</mi><mn>0</mn></mrow></msub></mrow></math> , the ideal one-dimensional power gain length is <math display="block" id="M24"><mrow><msub><mi>L</mi><mrow><mi>G</mi><mn>0</mn></mrow></msub><mo>≡</mo><mfrac><mrow><msub><mi>λ</mi><mi>u</mi></msub></mrow><mrow><mn>4</mn><mi>π</mi><msqrt><mn>3</mn></msqrt><mi>ρ</mi></mrow></mfrac><mn>.</mn></mrow></math> The cluster factor 〈| b (0)| 2 〉 at the entrance of the undulator comes from the initial shot noises of the e beams, which have been amplified during FEL. The level of shot noises is determined by the number of particles in the radiation coherence length. The cluster factor is[ 58 58 ], <math display="block" id="M25"><mrow><mrow><mo>〈</mo><mrow><mo>|</mo><mi>b</mi><msup><mrow><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo><mo>|</mo></mrow><mn>2</mn></msup></mrow><mo>〉</mo></mrow><mo>=</mo><mrow><mo>〈</mo><mrow><msup><mrow><mfrac><mn>1</mn><mrow><msub><mi>N</mi><mrow><msub><mi>l</mi><mrow><mtext>coh</mtext></mrow></msub></mrow></msub></mrow></mfrac></mrow><mn>2</mn></msup><mo>|</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>j</mi><mo>∈</mo><msub><mi>l</mi><mrow><mtext>coh</mtext></mrow></msub></mrow></munder></mstyle><msup><mi>e</mi><mrow><mo>−</mo><mi>i</mi><msub><mi>θ</mi><mi>j</mi></msub></mrow></msup><msup><mo>|</mo><mn>2</mn></msup></mrow><mo>〉</mo></mrow><mo>≈</mo><mfrac><mn>1</mn><mrow><msub><mi>N</mi><mrow><msub><mi>l</mi><mrow><mtext>coh</mtext></mrow></msub></mrow></msub></mrow></mfrac><mn>,</mn></mrow></math> where N l coh is the number of electrons in the scale of the coherence length, normalized bandwidth of SASE is Δ ω / ω ∼ ρ , coherence time is t coh ∼ λ 1 / cρ , and coherence length is l coh ∼ λ 1 / ρ . Thus, the coherence length can be regarded as the shift of the radiation relative to the electron beam over several gain lengths. Therefore, the starting shot noise of the SASE FEL is represented by the following equation: <math display="block" id="M26"><mrow><msub><mi>N</mi><mrow><msub><mi>l</mi><mrow><mtext>coh</mtext></mrow></msub></mrow></msub><mo>∼</mo><mfrac><mn>1</mn><mrow><mi>e</mi><mi>c</mi></mrow></mfrac><mfrac><mrow><msub><mi>λ</mi><mn>1</mn></msub></mrow><mi>ρ</mi></mfrac><mn>.</mn></mrow></math> We note that any SASE pulse is essentially an amplified undulator radiation. The one-dimensional power spectrum density can be used to calculate the radiated energy of the undulator[ 58 58 ]: <math display="block" id="M27"><mrow><msub><mi>U</mi><mrow><mtext>und</mtext></mrow></msub><mo>=</mo><mi>T</mi><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext>d</mtext><mi>ω</mi><mtext>d</mtext><mi>ϕ</mi><mfrac><mrow><mtext>d</mtext><mi>P</mi></mrow><mrow><mtext>d</mtext><mi>ω</mi><mtext>d</mtext><mi>ϕ</mi></mrow></mfrac><mo>→</mo><mi>T</mi><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><mtext>d</mtext><mi>ω</mi><mfrac><mrow><msup><mi>λ</mi><mn>2</mn></msup></mrow><mrow><msub><mi mathvariant="script">A</mi><mrow><mtext>tr</mtext></mrow></msub></mrow></mfrac><mfrac><mrow><mtext>d</mtext><mi>P</mi></mrow><mrow><mtext>d</mtext><mi>ω</mi><mtext>d</mtext><mi>ϕ</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>ϕ</mi><mo>=</mo><mn>0</mn></mrow></msub></mrow></math> Here λ 2 / tr can be understood as the characteristic angular divergence of a source, whose area is tr . Therefore, <math display="block" id="M28"><mrow><mfrac><mrow><mtext>d</mtext><mi>P</mi></mrow><mrow><mtext>d</mtext><mi>ω</mi></mrow></mfrac><msub><mo>|</mo><mrow><mn>1</mn><mi>D</mi></mrow></msub><mo>=</mo><mfrac><mrow><msup><mi>λ</mi><mn>2</mn></msup></mrow><mrow><msub><mi mathvariant="script">A</mi><mrow><mtext>tr</mtext></mrow></msub></mrow></mfrac><mi>δ</mi><mo stretchy="false">(</mo><mi>ϕ</mi><mo stretchy="false">)</mo><mfrac><mrow><mtext>d</mtext><mi>P</mi></mrow><mrow><mtext>d</mtext><mi>ω</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msup><mi>λ</mi><mn>2</mn></msup></mrow><mrow><msub><mi mathvariant="script">A</mi><mrow><mtext>tr</mtext></mrow></msub></mrow></mfrac><mfrac><mrow><mtext>d</mtext><mi>P</mi></mrow><mrow><mtext>d</mtext><mi>ω</mi><mtext>d</mtext><mi>ϕ</mi></mrow></mfrac><msub><mo>|</mo><mrow><mi>ϕ</mi><mo>=</mo><mn>0</mn></mrow></msub><mn>.</mn></mrow></math> Then, at the saturation length of the FEL ( N u ≈1/ ρ ), the saturated FEL energy can be rewritten as[ 58 58 ]: <math display="block" id="M29"><mrow><msub><mi>U</mi><mrow><mtext>FEL</mtext></mrow></msub><mo>=</mo><msub><mi>N</mi><mtext>e</mtext></msub><mi>ρ</mi><msubsup><mi>γ</mi><mrow><mtext>rmc</mtext></mrow><mn>2</mn></msubsup><mo>=</mo><mfrac><mrow><msub><mi>N</mi><mtext>e</mtext></msub></mrow><mrow><mi>ρ</mi><msub><mi>ω</mi><mn>1</mn></msub><mi>T</mi></mrow></mfrac><mfrac><mrow><msub><mi>U</mi><mrow><mtext>und</mtext></mrow></msub></mrow><mrow><mn>8</mn><mi>π</mi></mrow></mfrac></mrow></math> <math display="block" id="M30"><mrow><mo>∼</mo><mfrac><mrow><msub><mi>t</mi><mrow><mtext>coh</mtext></mrow></msub><msub><mi>N</mi><mtext>e</mtext></msub></mrow><mi>T</mi></mfrac><msub><mi>U</mi><mrow><mtext>und</mtext></mrow></msub><mo>=</mo><msub><mi>N</mi><mrow><msub><mi>l</mi><mrow><mtext>coh</mtext></mrow></msub></mrow></msub><msub><mi>U</mi><mrow><mtext>und</mtext></mrow></msub><mo>=</mo><mfrac><mi>T</mi><mrow><msub><mi>t</mi><mrow><mtext>coh</mtext></mrow></msub></mrow></mfrac><msubsup><mi>N</mi><mrow><msub><mi>l</mi><mrow><mtext>coh</mtext></mrow></msub></mrow><mn>2</mn></msubsup><mfrac><mrow><msub><mi>U</mi><mrow><mtext>und</mtext></mrow></msub></mrow><mrow><msub><mi>N</mi><mtext>e</mtext></msub></mrow></mfrac><mn>.</mn></mrow></math> The equation (22) indicates that the output of the saturated FEL is much stronger than the undulator emissions. Their ratio is the number( <math id="M31"><mrow><msub><mi>N</mi><mrow><msub><mi>l</mi><mrow><mtext>coh</mtext></mrow></msub></mrow></msub><msup><mrow><mo>></mo><mn>10</mn></mrow><mn>5</mn></msup></mrow></math> ) of particles in the coherence time. The energy of the FEL is the product of the single electron radiation energy of the undulator, square of the number of electrons in the coherence length and number of coherent regions T/t coh . In conclusion, the SASE FEL emits radiation from shot noise and the electron beams in the undulator can be distributed in the scale of the radiation wavelength to form periodic density modulations and achieve coherent enhancement of radiation intensity. The density modulation originates from the resonant interaction between the electron beams and X-rays in the periodic undulator. Hence, when the undulator is long enough, current of the electron beams is sufficiently large, and electron beams are of good quality(low emissivity and low energy divergence), the electrons will generate periodic density modulation at the scale of the resonant X-ray wavelength after interacting with the radiation. This results in intensities higher than incoherent undulator radiation, even though the electron beam is much longer than the radiation wavelength. In this work, we use the SASE FEL to obtain periodic microbunching e beams for generating high-quality XOAM. 4 Scheme for generating coherent X-rays carrying OAM at the SXFEL 1 As indicated in the recent study of the interaction between atoms and XUV spiral light[ 67 67 ], it is necessary to observe the non-dipole transition optionally when the photon beam has a single clear phase singularity or the position of the sample atom is close to the phase singularity. It is suggested[ 67 67 ] that the diffraction-limited storage ring[ 68 68 ] generates such XUV or X-ray spiral light under special conditions and an FEL could produce coherent vortex beams taking advantage of diverse schemes. The phased project SXFEL[ 46 46 ] mainly includes the SXFEL test facility (SXFEL-TF) and SXFEL user facility (SXFEL-UF). The SXFEL-TF is mainly used to promote the research of FEL in China, such as investigating the feasibility of seed XFEL with two-stage cascade HGHG-HGHG, proposing a brand-new plan based on an EEHG-HGHG cascade, and promoting the research and development of critical XFEL technology. Finally, the Chinese Academy of Sciences dec ided to build the SXFEL-TF at the SSRF in Shanghai. The initial commissioning was successfully performed from 2014 until the beginning of 2020. The Shanghai Institute of Applied Physics (SINAP) cooperated with the Shanghai-Tech University to upgrade the SXFEL-TF to a water window SXFEL-UF, including the development of the scientific experimental line station, upgradation of the energy of the linear accelerator to 1.5 GeV, construction of the second SASE FEL facility, as well as the construction of other necessary public facilities. It is expected to open to the public user by the end of this year. In the newly added second SASE FEL branch, a typical polarization control setup (using helical undulators as afterburners) is being planned. We use a configuration with minor modifications to the baseline design to generate X-rays carrying OAM, as shown in Fig. 3 Fig. 3 . This special method is analyzed in greater detail in ref[ 69 69 ]. The first SASE undulators are used for generating microbunching from the electron beam, and the microbunched electron beams are then sent to the downstream helical undulators to generate coherent X-rays carrying OAM with wavelengths of ∼1 nm at the second harmonic of the helical undulators. The main parameters are summarized in Table 1 Table 1 . 1 1 Parameters 1 1 Value 1 1 e - beam energy (GeV) 1 1 1.5 1 1 e - energy spread 1 1 0.01% 1 1 e - bunch charge (nC) 1 1 0.5 1 1 e - normalized emittance(mm ⋅ mrad) 1 1 0.6 1 1 e - pulse length (FWHM) ( μ m) 1 1 50 1 1 e - repetition rate (Hz) 1 1 10 1 1 SASE undulator K 1 1 1.6 1 1 SASE undulator period (mm) 1 1 15 1 1 SASE undulator length/segment (m) 1 1 4 1 1 SASE undulator No.of segment 1 1 4 1 1 Helical undulator K 1 1 2.3 1 1 Helical undulator period (mm) 1 1 15 1 1 Helical undulator length/segment (m) 1 1 4 1 1 Helical undulator No.of segment 1 1 2 We use the SIMPLEX code[ 70 70 ] to simulate this process with time dependence. First, we simulated the SASE FEL part with six planar undulators, and the gain curve of the first harmonic of FEL radiation at λ 1 =2 nm is shown in Fig. 4(left) Fig. 4(left) . In order to achieve better FEL performance downstream, we selected four SASE undulator segments owing to the best total bunching factor at z =18.9 m (shown in Fig. 4(Right) Fig. 4(Right) ) and relatively lower energy spread, compared to the saturation point. In our simulation, the focusing magnets are periodically arranged to be compatible with the layout of the undulator. We select the FUDU(QF-u-QD-U) lattice type to satisfy our requirements. Since the gain of the FEL is related to the e beam size and angular spread, there is an optimum betatron function that maximizes the gain. Thus, the lattice function should be selected to realize the optimum betatron function over the whole undulator line. The field gradient of the focusing quadrupole magnets in the upstream planar undulator segments is 6.50744 T/m and the field gradient of the defocusing quadrupole magnets is -3.50995 T/m. The focusing and defocusing quadrupole magnets are both 0.2 m long. In addition to the lattice type, the initial Twiss parameters at the entrance of the undulator are listed as follows. The horizontal ( x ) and vertical ( y ) betatron functions are 11.8631 m and 4.61501 m, respectively. while the horizontal ( x ) and vertical ( y ) alpha functions(corresponding to the derivative of the betatron function) are 1.59554 and -0.68688, respectively. The horizontal( x ) and vertical( y ) betatron functions along all the four SASE undulator segments are shown in Fig. 5(left) Fig. 5(left) . Moreover, the lattice type in the afterburner helical undulator segments is FUDU(QF-U-QD-U), and the field gradients of the focusing and defocusing quadrupole magnets are 4.26216 T/m and -4.27089 T/m, respectively. The focusing and defocusing quadrupole magnets are both 0.2 m long. The horizontal ( x ) and vertical ( y ) betatron functions are 11.462 m and 5.04241 m, respectively, while the horizontal ( x ) and vertical ( y ) alpha functions are 0.966439 and -0.520785. The horizontal ( x ) and vertical ( y ) betatron functions along the afterburner are shown in Fig. 5(right) Fig. 5(right) . In the SASE FEL, there is no special distribution of e beams when they first enter the undulator. However, the unstable SASE FEL makes the electrons produce shot noises which act on the electrons after amplification and redistribute the e beams. The profile of the emission from the device is determined by the microbunching structure and geometric configuration of the FEL radiator. We know that the emission at the fundamental modes of both planar and helical undulators has a strong peak on the axis and the screw-type phase of the XOAM originates from the spiral microbunched e beams. Based on the above discussions, these fine-tuned microbunched e beams require precise control over the electron profile upstream of the circularly polarized undulator with a particular input laser field. The fine-tuned microbunched e beam is then injected into the downstream radiator to emit and amplify coherent XOAM light. Thus, the above method requires the careful preservation of the spiral structure of the microbunching e beams when propagating into the downstream undulator. While it is relatively difficult to maintain the spiral structures of the microbunched e beams, our simulated method does not require the preservation of the spiral structures. High-quality periodic micro-bunched electron beams are obtained with the FUDU lattice design. The micro-bunched e beams are then injected into the downstream helical undulator. The helical undulator has the same periodic structural parameters as the upstream planar undulator to obtain X-ray vortices at the second harmonic at the wavelength λ 2 =1 nm. The gain curve along z , transverse intensity, and horizontal phase at the exit z =8.9 m are shown in Fig. 6 Fig. 6 . The typical donut shaped spatial distribution of the intensity distribution and spiral phase indicate that the second harmonic carries OAM with an azimuthal mode index of l =-1. The X-ray OAM intensity is at the sub- μ J level. We use the SASE FEL to produce microbunched e beams with a specific periodic structure when the preservation of the spiral microbunched e beams is unnessesary. Hence, this is an easier method for generating high energy XOAM. We also discuss the stability of the pulse energy of the second harmonic radiation using the same parameters for the SASE and helical undulator but different initial macroparticle distributions, i.e., random numbers auto seeding the SASE FEL,under the time-dependent state simulation mode. In the first section of the undulator, spontaneous radiation and periodic micro-bunched e beams are generated, which then act as source e beams for the X-ray carrying OAM generation process. After ten simulations with various initial shot noise seeding in the SASE section, we find that the evolution of the bunching factor at the first harmonic along the direction of propagation in the planar undulator and the pulse energy of the second harmonic radiation in the helical undulator are similar, and stable growth within the error range. These are shown in Fig. 7 Fig. 7 . Therefore, using the SASE FEL mode as the upstream setup, it is possible to generate similar well-bunched microbunching e beams, despite fluctuations of the initial macroparticle distribution producing shot noise seeding in the SASE. The pulse energy of the second harmonic of the helical undulator is relatively stable, guaranteeing the stability of this method for generating XOAM. It is known that when microbunched relativistic electron beams pass through the helical undulator and produce resonant interactions, the axisymmetric incident laser field transmits the structural information to the electron beams. Thus, an electron beam propagating through a long helical undulator can generate many microbunches of electrons. Each microbunch of electron beams is a helix, which can produce OAM radiation in harmonics, whose radiation powers are relatively weak. To increase the intensity of the vortex light, the special helical microbunched electron beams could be regarded as sources of OAM emission in a radiator, while carefully maintaining the spiral density distribution during propagation in the downstream undulator. In our simulation, the well-defined microbunches of electron beams are generated in the upstream planar undulator, which are then injected directly into the downstream short helical undulator. This is tuned to the harmonic and thus, the bunched e beam radiates coherently in the afterburner. The resultant vortex radiation carrying OAM is generated on the harmonics. In the downstream afterburner, only the parameter K of the helical undulator is changed. This enables it to emit coherent radiation at the second harmonic of the helical undulator, but owing to the weaker interaction between the microbunched e beam and helical undulator radiation at the second harmonic, the intensity of the OAM radiation emitted at the second harmonic is smaller than that at the first harmonic. In order to compare with the FERMI-like but different all-helical-undulator configuration[ 71 71 ], which is not externally seeded possesses the same set of parameters for the accelerator beam and undulators, and only differing in the use of a helical undulator rather than linear one, we simulate its FEL output at the second harmonic. It is shown that in the afterburner configuration ( Fig. 3 Fig. 3 ), the FEL intensity quickly saturates within 10 m, while in the FERMI-like all-helical-undulator configuration, the FEL intensity grows slowly to saturation in 25 m, with similar FEL intensity output. However, if we use the exact FERMI configuration with external seeding, the afterburner type setup can be easily mimicked to obtain fully coherent OAM. To summarize, there are two advantages of this simple afterburner setup using the existing devices in SXFEL at Shanghai. First, the short helical undulators can be directly added after the existing SASE planar undulators, without the need to generate and preserve the complicated helical electron distribution. Second, this afterburner configuration is more cost-effective and efficient since it is based on existing devices and there is no need to build a new all-helical undulator. 5 Conclusion 1 We simulate a straightforward and easy-to-implement helical undulator afterburner configuration at the SXFEL to generate high power X-rays carrying OAM with wavelengths of ∼ 1 nm. We use four SASE segments from the SASE FEL line at the SXFEL to produce periodic microbunched e beams. These are a consequence of the interaction between the initial randomly distributed electron beams and shot noise in the SASE process. The microbunched e beams are then injected into the downstream helical undulator (two segments). The scheme does not need additional dispersion devices along the transport channel of the electron beams, which can then be used to obtain optical spiral radiation with topological charge l =-1. The minus sign of the topological charge is due to the right-handed helicity of the circularly polarized undulator. We select the FUDU(QF-u-QD-U) lattice type in the SASE planar undulator and downstream helical undulator to satisfy our requirements. Using different shot noises in the initial SASE undulator, the relatively stable results are obtained. After simulation, the short helical undulators can be added directly after the existing SASE planar undulator. This configuration is more cost-effective and efficient as it is based on existing devices. As far as we know, this is the first demonstration of the generation of intense X-ray vortices whose wavelength is ∼1 nm at the SXFEL in China. We expect that exciting experiments with OAM carrying X-rays at the SXFEL will be realized as soon as the SASE facility is completed and the structure of the helical undulator is well processed. These X-rays with OAM property in the "magnetic window" can be used in various fields of research, especially for exploring magnetic dichroism [ 34 34 ] and detecting microscopic magnetic vortices[ 35 35 ]. References H. Rubinsztein-Dunlop , A. Forbes , M. V. Berry et al ., Roadmap on structured light . Journal of Optics 19 , 013001 ( 2016 ). https://doi.org/10.1088/2040-8978/19/1/013001 https://doi.org/10.1088/2040-8978/19/1/013001 Roadmap on structured light J. F. Nye , M. V. Berry , Dislocations in wave trains . Proc. R. Soc. Lond. A 336 , 0012 ( 1974 ). https://doi.org/10.1098/rspa.1974.0012 https://doi.org/10.1098/rspa.1974.0012 Dislocations in wave trains L. Allen , M. W. Beijersbergen , R. J. C. Spreeuw , J. P. Woerdman , Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes . Phys.Rev. A 45 , 8185 ( 1992 ). https://doi.org/10.1103/PhysRevA.45.8185 https://doi.org/10.1103/PhysRevA.45.8185 Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes L. Allen , S. M. Barnett , M. J. Padgett (Eds.), Optical Angular Momentum ( Institute of Physics Pub. , Bristol, Philadelphia , 2003 ). Optical Angular Momentum G. Molina-Terriza , J. P. Torres , L. Torner , Twisted photons . Nat. Phys . 3 , 305 ( 2007 ). https://doi.org/10.1038/nphys607 https://doi.org/10.1038/nphys607 Twisted photons S. Franke-Arnold , L. Allen , M. Padgett , Advances in optical angular momentum . Laser Photonics Rev . 2 , 299 ( 2008 ). https://doi.org/10.1002/lpor.200810007 https://doi.org/10.1002/lpor.200810007 Advances in optical angular momentum A. M. Yao , M. J. Padgett , Orbital angular momentum: origins, behavior and applications . Adv. Opt. Photonics 3 , 161 ( 2011 ). https://doi.org/10.1364/AOP.3.000161 https://doi.org/10.1364/AOP.3.000161 Orbital angular momentum: origins, behavior and applications J. P. Torres , L. Torner (Eds.), Twisted photon: Applications of Light with Orbital Angular Momentum . 2011 WILEY-VCH Verlag GmbH & Co. KGaA , Weinheim . Twisted photon: Applications of Light with Orbital Angular Momentum D. L. Andrews , M. Babiker (Eds.), The Angular Momentum of Light , Cambridge University Press 2013 . The Angular Momentum of Light R. W. Bowman , M. J. Padgett , Optical trapping and binding . Rep. Prog. Phys . 76 , 026401 ( 2013 ). https://doi.org/10.1088/0034-4885/76/2/026401 https://doi.org/10.1088/0034-4885/76/2/026401 Optical trapping and binding M. J. Padgett , Orbital angular momentum 25 years on . Opt. Express 25 , 11265 ( 2017 ). https://doi.org/10.1364/OE.25.011265 https://doi.org/10.1364/OE.25.011265 Orbital angular momentum 25 years on Y.J. Shen , X.J. Wang , Z.W. Xie et al ., Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities . Light: Science & Applications 8 , 90 ( 2019 ). https://doi.org/10.1038/s41377-019-0194-2 https://doi.org/10.1038/s41377-019-0194-2 Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities D. S. Bykov , S.R. Xie , R. Zeltner et al ., Long-range optical trapping and binding of microparticles in hollow-core photonic crystal fibre . Light: Science & Applications 7 , 22 ( 2018 ). https://doi.org/10.1038/s41377-018-0015-z https://doi.org/10.1038/s41377-018-0015-z Long-range optical trapping and binding of microparticles in hollow-core photonic crystal fibre N.S. Huang , K. Li , H.X. Deng , BRIGHT: the three-dimensional X-ray crystal Bragg diffraction code . Nucl. Sci. Tech . 30 , 39 ( 2019 ). https://doi.org/10.1007/s41365-019-0559-5 https://doi.org/10.1007/s41365-019-0559-5 BRIGHT: the three-dimensional X-ray crystal Bragg diffraction code H.L. Wu , J.H. Chen , Z. T. Zhao et al ., Proposal for a high brightness γ-ray source at the SXFEL . Nucl. Sci. Tech . 26 , 050103 ( 2015 ). https://doi.org/10.13538/j.1001-8042/nst.26.050103 https://doi.org/10.13538/j.1001-8042/nst.26.050103 Proposal for a high brightness γ-ray source at the SXFEL A. G. Peele , P. J. McMahon , D. Paterson et al ., Observation of an x-ray vortex . Opt. Lett . 27 , 1752 ( 2002 ). https://doi.org/10.1364/OL.27.001752 https://doi.org/10.1364/OL.27.001752 Observation of an x-ray vortex A. G. Peele , K. A. Nugent , A. P. Mancuso et al ., X-ray phase vortices: theory and experiment . J. Opt. Soc. Am. A 21 , 1575 ( 2004 ). https://doi.org/10.1364/JOSAA.21.001575 https://doi.org/10.1364/JOSAA.21.001575 X-ray phase vortices: theory and experiment B. Terhalle , A. Langner , V. A. Guzenko et al ., Generation of extreme ultraviolet vortex beams using computer generated holograms . Optics Letters 36 , 4143 ( 2011 ). https://doi.org/10.1364/OL.36.004143 https://doi.org/10.1364/OL.36.004143 Generation of extreme ultraviolet vortex beams using computer generated holograms A. Sakdinawat , Y.W. Liu , Soft-x-ray microscopy using spiral zone plates . Optics Letters 32 , 2635 ( 2007 ). https://doi.org/10.1364/OL.32.002635 https://doi.org/10.1364/OL.32.002635 Soft-x-ray microscopy using spiral zone plates D. Cojoc , B. Kaulich , A. Carpentiero et al ., X-ray vortices with high topological charge . Microelectronic Engineering 83 , 1360 ( 2006 ). https://doi.org/10.1016/j.mee.2006.01.066 https://doi.org/10.1016/j.mee.2006.01.066 X-ray vortices with high topological charge Y. Kohmura , K. Sawada , M. Taguchi et al ., Formation of x-ray vortex dipoles using a single diffraction pattern and direct phase measurement using interferometry . Appl. Phys. Lett . 94 , 101112 ( 2009 ). https://doi.org/10.1063/1.3095828 https://doi.org/10.1063/1.3095828 Formation of x-ray vortex dipoles using a single diffraction pattern and direct phase measurement using interferometry J. C. T Lee , S. J. Alexander , S. D. Kevan et al ., Laguerre-Gauss and Hermite-Gauss soft X-ray states generated using diffractive optics . Nat. Photon . 13 , 205 ( 2019 ). https://doi.org/10.1038/s41566-018-0328-8 https://doi.org/10.1038/s41566-018-0328-8 Laguerre-Gauss and Hermite-Gauss soft X-ray states generated using diffractive optics L. Loetgering , M. Baluktsian , K. Keskinbora et al ., Generation and characterization of focused helical x-ray beams . Sci. Adv . 6 , eaax8836 ( 2020 ). https://dx.doi.org/10.1126/sciadv.aax8836 https://dx.doi.org/10.1126/sciadv.aax8836 Generation and characterization of focused helical x-ray beams D. Seipt , A. Surzhykov , S. Fritzsche , Structured x-ray beams from twisted electrons by inverse Compton scattering of laser light . Phys. Rev. A 90 , 012118 ( 2014 ). https://dx.doi.org/10.1103/PhysRevA.90.012118 https://dx.doi.org/10.1103/PhysRevA.90.012118 Structured x-ray beams from twisted electrons by inverse Compton scattering of laser light V. Petrillo , G. Dattoli , I. Drebot et al ., Compton Scattered X-Gamma Rays with Orbital Momentum . Phys. Rev. Lett . 117 , 123903 ( 2016 ). https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.117.123903 https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.117.123903 Compton Scattered X-Gamma Rays with Orbital Momentum M. Zurch , C. Kern , P. Hansinger et al ., Strong-field physics with singular light beams . Nature physics 8 , 743 ( 2012 ). https://doi.org/10.1038/nphys2397 https://doi.org/10.1038/nphys2397 Strong-field physics with singular light beams G. Gariepy , J. Leach , K. T. Kim et al ., Creating High-Harmonic Beams with Controlled Orbital Angular Momentum . Phys. Rev. Lett . 113 , 153901 ( 2014 ). https://dx.doi.org/10.1103/PhysRevLett.113.153901 https://dx.doi.org/10.1103/PhysRevLett.113.153901 Creating High-Harmonic Beams with Controlled Orbital Angular Momentum C. Hernandez-Garcia , J. Vieira , J. T. Mendonca et al ., Generation and Applications of Extreme-Ultraviolet Vortices . Photonics 4 , 28 ( 2017 ). https://doi.org/10.3390/photonics4020028 https://doi.org/10.3390/photonics4020028 Generation and Applications of Extreme-Ultraviolet Vortices L. Rego , K. M. Dorney , N. J. Brooks et al ., Generation of extreme-ultraviolet beams with time-varying orbital angular momentum . Science 364 , eaaw9486 ( 2019 ). https://doi.org/10.1126/science.aaw9486 https://doi.org/10.1126/science.aaw9486 Generation of extreme-ultraviolet beams with time-varying orbital angular momentum J. Vieira , R.M.G.M. Trines , E.P. Alves et al ., Amplification and generation of ultra-intense twisted laser pulses via stimulated Raman scattering . Nature Communications 7 , 10371 ( 2016 ). https://doi.org/10.1038/ncomms10371 https://doi.org/10.1038/ncomms10371 Amplification and generation of ultra-intense twisted laser pulses via stimulated Raman scattering J. Vieira , R. M. G. M. Trines , E. P. Alves et al ., High Orbital Angular Momentum Harmonic Generation . Phys. Rev. Lett . 117 , 265001 ( 2016 ). https://doi.org/10.1103/PhysRevLett.117.265001 https://doi.org/10.1103/PhysRevLett.117.265001 High Orbital Angular Momentum Harmonic Generation A. Picon , A. Benseny , J. Mompart et al ., Transferring orbital and spin angular momenta of light to atoms . New Journal of Physics 12 , 083053 ( 2010 ). https://doi.org/10.1088/1367-2630/12/8/083053 https://doi.org/10.1088/1367-2630/12/8/083053 Transferring orbital and spin angular momenta of light to atoms A. Picon , J. Mompart , J. R. V. de Aldana , et al ., Photoionization with orbital angular momentum beams . Opt. Express 18 , 3660 ( 2010 ). https://doi.org/10.1364/OE.18.003660 https://doi.org/10.1364/OE.18.003660 Photoionization with orbital angular momentum beams M. v. Veenendaal , I. McNulty , Prediction of Strong Dichroism Induced by X Rays Carrying Orbital Momentum . Phys. Rev. Lett . 98 , 157401 ( 2007 ). https://link.aps.org/doi/10.1103/PhysRevLett.98.157401 https://link.aps.org/doi/10.1103/PhysRevLett.98.157401 Prediction of Strong Dichroism Induced by X Rays Carrying Orbital Momentum M. v. Veenendaal , Interaction between x-ray and magnetic vortices . Phys. Rev. B 92 , 245116 ( 2015 ). https://doi.org/10.1103/PhysRevB.92.245116 https://doi.org/10.1103/PhysRevB.92.245116 Interaction between x-ray and magnetic vortices G. Friesecke , R. D. James , D. Juestel , Twisted X-Rays: Incoming Waveforms Yielding Discrete Diffraction Patterns for Helical Structures . SIAM J. Appl. Math . 76 , 1191 ( 2016 ). https://doi.org/10.1137/15M1043418 https://doi.org/10.1137/15M1043418 Twisted X-Rays: Incoming Waveforms Yielding Discrete Diffraction Patterns for Helical Structures D. Juestel , G. Friesecke , R. D. James , Bragg-von Laue diffraction generalized to twisted X-rays . Acta Cryst. A 72 , 190 ( 2016 ). https://dx.doi.org/10.1107/S2053273315024390 https://dx.doi.org/10.1107/S2053273315024390 Bragg-von Laue diffraction generalized to twisted X-rays L. Ye , J. R. Rouxel , S. Asban et al ., Probing Molecular Chirality by Orbital-Angular-Momentum-Carrying X-ray Pulses . J. Chem. Theory Comput . 15 , 4180 ( 2019 ). https://doi.org/10.1021/acs.jctc.9b00346 https://doi.org/10.1021/acs.jctc.9b00346 Probing Molecular Chirality by Orbital-Angular-Momentum-Carrying X-ray Pulses S. Sasaki , I. McNulty , Proposal for Generating Brilliant X-Ray Beams Carrying Orbital Angular Momentum . Phys. Rev. Lett . 100 , 124801 ( 2008 ). https://doi.org/10.1103/PhysRevLett.100.124801 https://doi.org/10.1103/PhysRevLett.100.124801 Proposal for Generating Brilliant X-Ray Beams Carrying Orbital Angular Momentum J. Bahrdt , K. Holldack , P. Kuske et al ., First Observation of Photons Carrying Orbital Angular Momentum in Undulator Radiation . Phys. Rev. Lett . 111 , 034801 ( 2013 ). https://doi.org/10.1103/PhysRevLett.111.034801 https://doi.org/10.1103/PhysRevLett.111.034801 First Observation of Photons Carrying Orbital Angular Momentum in Undulator Radiation T. Kaneyasu , Y. Hikosaka , M. Fujimoto et al ., Observation of an optical vortex beam from a helical undulator in the XUV region . J. Synchrotron Rad . 24 , 934 ( 2017 ). https://doi.org/10.1107/S1600577517009626 https://doi.org/10.1107/S1600577517009626 Observation of an optical vortex beam from a helical undulator in the XUV region E. Hemsing , A. Knyazik , M. Dunning et al ., Coherent optical vortices from relativistic electron beams . Nat. Phys . 9 , 549 ( 2013 ). https://doi.org/10.1038/nphys2712 https://doi.org/10.1038/nphys2712 Coherent optical vortices from relativistic electron beams E. Hemsing , M. Dunning , C. Hast et al ., First Characterization of Coherent Optical Vortices from Harmonic Undulator Radiation . Phys. Rev. Lett . 113 , 134803 ( 2014 ). https://doi.org/10.1103/PhysRevLett.113.134803 https://doi.org/10.1103/PhysRevLett.113.134803 First Characterization of Coherent Optical Vortices from Harmonic Undulator Radiation P. R. Ribic , D. Gauthier , G. D. Ninno , Generation of Coherent Extreme-Ultraviolet Radiation Carrying Orbital Angular Momentum . Phys. Rev. Lett . 112 , 203602 ( 2014 ). https://doi.org/10.1103/PhysRevLett.112.203602 https://doi.org/10.1103/PhysRevLett.112.203602 Generation of Coherent Extreme-Ultraviolet Radiation Carrying Orbital Angular Momentum P. R. Ribic , B. Rosner , D. Gauthier et al ., Extreme-Ultraviolet Vortices from a Free-Electron Laser . Phys. Rev. X 7 , 031036 ( 2017 ). https://doi.org/10.1103/PhysRevX.7.031036 https://doi.org/10.1103/PhysRevX.7.031036 Extreme-Ultraviolet Vortices from a Free-Electron Laser Z.T. Zhao , D. Wang , Q. Gu et al ., Status of the SXFEL Facility . Appl. Sci . 7 , 607 ( 2017 ). https://doi.org/10.3390/app7060607 https://doi.org/10.3390/app7060607 Status of the SXFEL Facility A. Afanasev , A. Mikhailichenko , On generation of photons carrying orbital angular momentum in the helical undulator . arXiv:1109.1603 ( 2011 ). On generation of photons carrying orbital angular momentum in the helical undulator M. Katoh , M. Fujimoto , H. Kawaguchi et al ., Angular Momentum of Twisted Radiation from an Electron in Spiral Motion . Phys. Rev. Lett . 118 , 094801 ( 2017 ). https://doi.org/10.1103/PhysRevLett.118.094801 https://doi.org/10.1103/PhysRevLett.118.094801 Angular Momentum of Twisted Radiation from an Electron in Spiral Motion M. Katoh , M. Fujimoto , N. S. Mirian et al ., Helical Phase Structure of Radiation from an Electron in Circular Motion . Sci. Rep . 7 , 6130 ( 2017 ). https://doi.org/10.1038/s41598-017-06442-2 https://doi.org/10.1038/s41598-017-06442-2 Helical Phase Structure of Radiation from an Electron in Circular Motion H. Kawaguchi , M. Katoh , Orbital angular momentum of Lienard-Wiechert fields . Prog. Theor. Exp. Phys . 2019 , 083A02 ( 2019 ). https://doi.org/10.1093/ptep/ptz059 https://doi.org/10.1093/ptep/ptz059 Orbital angular momentum of Lienard-Wiechert fields V. Epp , U. Guselnikova , Angular momentum of radiation from a charge in circular and spiral motion . Physics Letters A 383 , 2668 ( 2019 ). https://doi.org/10.1016/j.physleta.2019.05.038 https://doi.org/10.1016/j.physleta.2019.05.038 Angular momentum of radiation from a charge in circular and spiral motion E. Hemsing , P. Musumeci , S. Reiche et al ., Helical electron-beam microbunching by harmonic coupling in a helical undulator . Physical Review Letters 102 , 174801 ( 2009 ). https://doi.org/10.1103/PhysRevLett.102.174801 https://doi.org/10.1103/PhysRevLett.102.174801 Helical electron-beam microbunching by harmonic coupling in a helical undulator M.H.J , X. Yang , H.J. Xu et al ., Shanghai Synchrotron Radiation Facility . Chinese Science Bulletin 54 , 4171 ( 2009 ). https://doi.org/10.1007/s11434-009-0689-y https://doi.org/10.1007/s11434-009-0689-y Shanghai Synchrotron Radiation Facility Z.C. Chen , Y.B. Leng , R.X. Yuan et al ., Study of algorithms of phase advance measurement between BPMs and its application in SSRF . Nuclear Science and Techniques 24 , 010102 ( 2013 ). https://doi.org/10.13538/j.1001-8042/nst.2013.01.007 https://doi.org/10.13538/j.1001-8042/nst.2013.01.007 Study of algorithms of phase advance measurement between BPMs and its application in SSRF Z.C. Chen , Y.B. Leng , R.X. Yuan et al ., Beam position monitor troubleshooting by using principal component analysis in Shanghai Synchrotron Radiation Facility . Nuclear Science and Techniques 25 , 020102 ( 2014 ). https://doi.org/10.13538/j.1001-8042/nst.25.020102 https://doi.org/10.13538/j.1001-8042/nst.25.020102 Beam position monitor troubleshooting by using principal component analysis in Shanghai Synchrotron Radiation Facility T. Tanaka , H. Kitamura , SPECTRA: a synchrotron radiation calculation code . Journal of Synchrotron Radiation 8 , 1221 ( 2001 ). https://doi.org/10.1107/S090904950101425X https://doi.org/10.1107/S090904950101425X SPECTRA: a synchrotron radiation calculation code C. Feng , H.X. Deng , Review of fully coherent free-electron lasers . Nuclear Science and Techniques 29 , 160 ( 2018 ). https://doi.org/10.1007/s41365-018-0490-1 https://doi.org/10.1007/s41365-018-0490-1 Review of fully coherent free-electron lasers K.-J. Kim , Z.R. Huang , R. Lindberg , Synchrotron Radiation and Free-Electron Lasers: Principles of Coherent X-Ray Generation . Cambridge University Press 2017 . Synchrotron Radiation and Free-Electron Lasers: Principles of Coherent X-Ray Generation W. Ackermann , G. Asova , K. Zapfe et al ., Operation of a free-electron laser from the extreme ultraviolet to the water window . Nature Photonics 1 , 336 ( 2007 ). https://doi.org/10.1038/nphoton.2007.76 https://doi.org/10.1038/nphoton.2007.76 Operation of a free-electron laser from the extreme ultraviolet to the water window J. Galayda , J. Arthur , P. Anfinrud et al ., Linac Coherent Light Source (LCLS) conceptual design report. SLAC, Report SLAC-R-593 , 2002 Linac Coherent Light Source (LCLS) conceptual design report. SLAC, Report SLAC-R-593 T. Shintake et al ., sPring-8 Compact SASE Source Conceptual Design Report . 2005 . DOI: http://www-xfel.spring8.or.jp http://www-xfel.spring8.or.jp sPring-8 Compact SASE Source Conceptual Design Report R. Brinkmann , et al . TESLA XFEL: First stage of the X-ray laser laboratory (Technical Design Report, Supplement). DESY, Report TESLA FEL 2002-09 , 2002 . TESLA XFEL: First stage of the X-ray laser laboratory (Technical Design Report, Supplement). DESY, Report TESLA FEL 2002-09 H.-S. Kang , K.-W. Kim , I.-S. Ko , Current status of PAL-XFEL project . Proceedings of IPAC 2014 , p. 2897 , 2014 Current status of PAL-XFEL project R. Ganter , et al . SwissFEL-conceptual design report . Paul Scherrer Institute (PSI) , Report, 2010 SwissFEL-conceptual design report R. Bonifacio , C. Pellegrini , L. M. Narducci , Collective instabilities and high-gain regime in a free electron laser . Opt. Commun . 50 , 373 ( 1984 ). https://doi.org/10.1016/0030-4018(84)90105-6 https://doi.org/10.1016/0030-4018(84)90105-6 Collective instabilities and high-gain regime in a free electron laser N. Kroll , W. McMullin , Stimulated emission from relativistic electrons passing through a spatially periodic transverse magnetic field . Phys. Rev. A , 17 , 300 ( 1978 ). https://doi.org/10.1103/PhysRevA.17.300 https://doi.org/10.1103/PhysRevA.17.300 Stimulated emission from relativistic electrons passing through a spatially periodic transverse magnetic field T. Kaneyasu , Y. Hikosaka , M. Fujimoto et al ., Limitations in photoionization of helium by an extreme ultraviolet optical vortex . Phys. Rev. A 95 , 023413 ( 2017 ). https://doi.org/10.1103/PhysRevA.95.023413 https://doi.org/10.1103/PhysRevA.95.023413 Limitations in photoionization of helium by an extreme ultraviolet optical vortex M. Eriksson , J. F. van der Veen , C. Quitmann , Diffraction-limited storage rings - a window to the science of tomorrow . J. Synchrotron Radiat . 21 , 837 ( 2014 ). https://doi.org/10.1107/S1600577514019286 https://doi.org/10.1107/S1600577514019286 Diffraction-limited storage rings - a window to the science of tomorrow E. Hemsing , Coherent photons with angular momentum in a helical afterburner . Physical Review Accelerators and Beams 23 , 020703 ( 2020 ). https://doi.org/10.1103/PhysRevAccelBeams.23.020703 https://doi.org/10.1103/PhysRevAccelBeams.23.020703 Coherent photons with angular momentum in a helical afterburner T. Tanaka , SIMPLEX: simulator and postprocessor for free-electron laser experiments . Journal of Synchrotron Radiation 22 , 1319 ( 2015 ). https://doi.org/10.1107/S1600577515012850 https://doi.org/10.1107/S1600577515012850 SIMPLEX: simulator and postprocessor for free-electron laser experiments C. J. Bocchetta , et al . FERMI@Elettra CDR , http://www.elettra.trieste.it/FERMI http://www.elettra.trieste.it/FERMI , ( 2007 ). FERMI@Elettra CDR We appreciate the stimulating discussions with Nanshun Huang, Zhichu Chen, Hailong Wu, Bo Liu, and Dong Wang of Shanghai Advanced Research Institute, Chinese Academy of Sciences. 10.1007/s41365-020-00794-7 This work was supported by the National Development and Reform Commission ([2013]2347) and National Basic Research Program of China (No. 2015CB859700). 2020-05-29 2020-07-16 2020-07-17 2020-08-28 2020-09-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2020 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-020-00794-7 http://dx.doi.org/10.1007/s41365-020-00794-7 . 2020 Nuclear Science and Techniques 9 31