Discussion of 90° stopband in low-energy superconducting linear accelerators

Discussion of 90° stopband in low-energy superconducting linear accelerators 1 first-author https://orcid.org/0000-0002-0033-2107 1 1 1 1 1 corresp *lizhihui@scu.edu.cn *lizhihui@scu.edu.cn lizhihui@scu.edu.cn https://orcid.org/0000-0001-7892-8245 The Key Laboratory of Radiation Physics and Technology of Ministry of Education, Institute of Nuclear Science and Technology, Sichuan University, Chengdu 610065, China The Key Laboratory of Radiation Physics and Technology of Ministry of Education, Institute of Nuclear Science and Technology, Sichuan University, Chengdu 610065, China Superconducting linear accelerators (SCL) have a high acceleration gradient and are capable of operating in a high-duty factor mode. For high-power and high-intensity SCL, the design of beam dynamics generally follows the principle that the zero-current periodic phase advance (σ 0 ) of each degree of freedom is less than 90° to avoid envelope instability caused by space charge. However, this principle is obtained under the condition of a completely periodic focusing channel, and it is ambiguous for pseudo-periodic structures, such as linear accelerators. Although transverse beam dynamics without acceleration have been studied by other researchers, it appears that there are some connections between pure 2D and 3D beam dynamics. Based on these two points, five focusing schemes for the solenoid and quadrupole doublet channels were designed to simulate the beam behavior with non-constant σ 0 . Among them, the four schemes follow the characteristics of variation in the zero-current longitudinal phase advance (σ 0l ) under a constant acceleration gradient and synchronous phase. The zero-current transverse phase advance (σ 0t ) is consistent with σ 0l , based on the equipartition requirement. The initial σ 0t was set to 120°, 110°, 100°, and 90°, and was then gradually decreased to approximately 40° at the end of the channel. The last scheme maintains the maximum σ 0t of 88° by reducing the acceleration gradient of the corresponding cavities, until the point at which σ 0t equals 88° with a normal gradient. Using the stopbands obtained from the linearized envelope equations and multiparticle particle-in-cell (PIC) simulations, the transport properties of both continuous and 3D-bunched beams with the acceleration of the five focusing schemes were studied. It was found that for a CW beam, when tune depression > 0.7, σ 0t can break through 90° when the beams were transported in both solenoid and quadrupole doublet periodic focusing channels. When tune depression < 0.7, the conclusions were different. For the solenoid focusing system, σ 0t can partially break through 90°, and the beam quality is not significantly affected. For the quadrupole doublet focusing system, a partial breakthrough of 90° has a greater impact on the beam quality. The same conclusions were obtained for a bunched beam with acceleration. Proton beam Superconducting linear accelerators Envelope instability Periodic focusing structure Resonance. Introduction 1 Linear accelerators (linacs) have been investigated for several applications [ 1 1 ]. Since the successful first operation of a particle accelerator at the beginning of the last century, remarkable advances have been made in the physics and technology of linacs [ 2 2 ]. Theoretical and experimental studies on the stability of high-intensity beams in linacs began in the 1960s. Most of these theoretical studies have been based on two-dimensional transport models without acceleration [ 3 3 ]. The associated strong space-charge effects present challenges in the transport of high-intensity beams [ 4 4 ]. One of the conclusions of these studies is that when transverse phase advance per period σ 0t is greater than 90°, the space charge force may drive envelope instability and degrade the beam quality [ 5 5 , 6 6 ]. To avoid envelope instability in high-intensity linacs, one of the rules that most designers follow is to maintain σ 0t < 90° [ 7 7 ]. When the phase advance is near 90°, fourth-order resonance and envelope instability occur [ 8 8 ]. The first experimental study of the fourth-order resonance in an operational high-intensity accelerator was undertaken at GSI UNILAC, which confirmed the emittance degradation [ 9 9 ]. For Kapchinskij-Vladmirskij (K-V) beams, using perturbation analysis of the linearized Vlasov-Poisson equations, stopbands of instability can be obtained [ 5 5 ]. Stopbands are defined as regions in which any eigenvalue lies outside the complex unit cycle. In stopbands, the emittance increases exponentially in a nonzero state. Higher-order and lower-order stopbands may coincide, but instability typically decreases with increasing order [ 10 10 ]. Meanwhile, the special case of the "second-order even" mode is known as the envelope instability [ 11 11 ]. Envelope instability is characterized by parametric and confluent resonance. Parametric resonance represents the resonance between the focusing lattices and one of the envelope oscillatory modes that occurs only in the solenoid channel. When parametric resonance occurs, only one of the breathing and quadrupole modes excites resonance with the beam core. This is also called half-integer resonance because it represents the envelope eigenmode that oscillates at half the frequency of the periodically varying focusing (including space charge). Confluent resonance represents the resonance between both envelope oscillation modes and lattices, which occurs only in the quadrupole channel [ 12 12 ]. The space charge force not only causes envelope instability, but may also cause coupling of the particle motions and emittance exchange in different directions [ 13 13 ]. Therefore, Jameson introduced the equipartition concept into the beam dynamics design of linear accelerators [ 9 9 , 14 14 ]. By properly setting the longitudinal and transverse focusing strengths, the oscillation energies of different degrees of freedom were balanced, and there was no free energy for resonance to occur. For the equipartition design, the motions of the particles in the transverse and longitudinal directions are decoupled. Therefore, the transverse and longitudinal motions of the superconducting linear accelerator (SCL) based on equipartition design can be studied separately. However, σ 0t and σ 0l maintain the same trend because of the requirement for equipartition. SCL can provide high-power beams and have many applications in various fields. It has negligible wall power losses, which makes it possible to operate in the high-duty factor mode with a higher accelerating gradient. The SCL is typically composed of different acceleration sections, and each section is composed of identical cells [ 15 15 ]. The significant difference between SCL and periodic focusing channels is that SCL have radio-frequency (RF) cavities. The RF cavity can not only accelerate beams, but can also provide focusing in the longitudinal direction and defocusing in the transverse direction to the particles [ 16 16 ]. One of the most important outcomes for the SCL is the pseudo-periodic focusing structure, which may lead to different dynamics. When the acceleration gradient and synchronous phase are constant, the square of the focusing strength is inversely proportional to the 3D power of the βγ of particles [ 17 17 ], as follows: <math display="block" id="M1"><mrow><mtable><mtr><mtd><mrow><msubsup><mi>σ</mi><mrow><mn>0</mn><mi>l</mi></mrow><mn>2</mn></msubsup><mo>=</mo><mfrac><mrow><mn>2</mn><mi>π</mi><mi>q</mi><msub><mi>E</mi><mn>0</mn></msub><mi>T</mi><msup><mi>S</mi><mn>2</mn></msup><mtext>sin</mtext><mrow><mo>(</mo><mrow><mo>−</mo><msub><mi>ϕ</mi><mtext>s</mtext></msub></mrow><mo>)</mo></mrow></mrow><mrow><mi>m</mi><msup><mi>c</mi><mn>2</mn></msup><msubsup><mi>β</mi><mi>s</mi><mn>3</mn></msubsup><msubsup><mi>γ</mi><mi>s</mi><mn>3</mn></msubsup><mi>λ</mi></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> Thus, the maximum σ 0l occurs at the beginning of each section and decreases as the beam energy increases with the SCL. However, cavities at the beginning of each section must work under a reduced acceleration gradient because of the principle that σ 0 < 90°. This significantly reduces the acceleration efficiency. The remainder of this paper is as follows. In Sect. 2, five different focusing schemes are designed and presented. In Sect. 3, the types and characteristics of periodic focusing channels and an approach for calculating the stopbands are introduced in detail. In Sect. 4, the beam dynamics of five different focusing schemes are shown. Finally, the effects of σ 0l exceeding 90° in a real SCL section are studied and verified. The conclusions and discussion are presented in Sect. 6. Different focusing schemes 1 To study the transverse beam dynamics when σ 0 partially exceeds 90°, five different focusing schemes were designed to simulate σ 0t . The maximum σ 0t of the beam was designed to be 120°(①), 110°(②), 100°(③), and 90°(④). The purpose was to simulate the trend of σ 0t under different acceleration gradients in the presence of an acceleration field. A higher acceleration gradient corresponds to a larger σ 0t . The first scheme corresponds to the highest acceleration gradient and the fourth scheme corresponds to the lowest acceleration gradient. To avoid excessive sensitivity of the beam to the focusing parameters, for every focusing scheme, the minimum value of σ 0t was generally set to 40° after the beam passes through 80 focusing periods. Therefore, we designed focusing schemes based on a logistic function. This function is more suitable for the changes in σ 0t . In addition, for the case in which the maximal σ 0t was 120°, a control scheme was set to avoid envelope instability on the beam. The maximal σ 0t of this control scheme was 88°(⑤), where envelope instability does not occur. This can be controlled by reducing the acceleration gradient of the corresponding cavity. Once σ 0t of ① was less than 90°, ⑤ and ① remain within the same gradient. The common logistic function is given by <math display="block" id="M2"><mrow><mtable><mtr><mtd><mrow><mi>y</mi><mo>=</mo><msub><mi>A</mi><mn>2</mn></msub><mo>+</mo><mfrac><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>−</mo><msub><mi>A</mi><mn>2</mn></msub></mrow><mrow><mn>1</mn><mo>+</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mi>x</mi><mrow><msub><mi>x</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mi>p</mi></msup></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> In the design functions, x represents the period number, and y represents the value of σ 0t . The specific parameters of the functions are listed in Table. 1 Table. 1 . Additionally, σ 0t of the five schemes are shown in Fig. 2 Fig. 2 . 1 1 Scheme 1 1 Max  σ 0t   (°) 1 1 Min σ 0t (°) 1 1 A 1 1 1 A 2 1 1 x 0 1 1 p 1 1 ① 1 1 120 1 1 40 1 1 120 1 1 29.8312 1 1 26.3333 1 1 1.87755 1 1 ② 1 1 110 1 1 40 1 1 110 1 1 31.1023 1 1 26.3333 1 1 1.87755 1 1 ③ 1 1 100 1 1 40 1 1 100 1 1 32.3734 1 1 26.3333 1 1 1.87755 1 1 ④ 1 1 90 1 1 40 1 1 90 1 1 33.6445 1 1 26.3333 1 1 1.87755 The focusing channels and stopbands 1 Periodic focusing channels 2 According to the transverse motion equation, the RF electromagnetic field has an evident defocusing effect and axial symmetry at the low-energy end of the linear accelerator. Considering that the beam has a strong space-charge effect at the low-energy end, a relatively compact transverse focusing lattice is required to maintain an appropriate beam size within our design. Magnetic quadrupoles, solenoids, and axisymmetric electrostatic lattices were chosen by designers for beam transport systems at low-energy ends [ 18 18 ]. In this study, a solenoid, which has axial symmetry, and a quadrupole doublet, which is non-axisymmetric, were chosen as the focusing elements. The specific parameters of the focusing elements used for the following discussion are shown in Fig. 1 Fig. 1 . For the solenoid periodic focusing system, the <math id="M3"><mi>κ</mi></math> ( s ) is given by [ 12 12 ] <math display="block" id="M4"><mrow><mtable><mtr><mtd><mrow><mi>κ</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>κ</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>κ</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><mfrac><mrow><mi>q</mi><mi>B</mi></mrow><mrow><mn>2</mn><msub><mi>m</mi><mn>0</mn></msub><mi>c</mi><mi>β</mi><mi>γ</mi></mrow></mfrac></mrow><mo>]</mo></mrow></mrow><mn>2</mn></msup></mrow></mtd></mtr></mtable></mrow></math> where B is the solenoid lens strength. For a magnetic quadrupole channel, <math display="block" id="M5"><mrow><mtable><mtr><mtd><mrow><mi>κ</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>κ</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mo>−</mo><msub><mi>κ</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>G</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>c</mi><mi>β</mi><mi>γ</mi></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> where G is the quadrupole linear strength gradient. The distribution of beam particles in phase space plays a significant role because it is directly related to the space charge in high-intensity proton accelerators. The K-V distribution is an ellipse with uniform density in any projected phase plane, which is beneficial for calculating the space charge force. This is because the space charge force is a linear defocusing force that can be compensated for by an external field. Therefore, the K-V distribution is most commonly used for researching space charges in simulations. From the research of Sacher, the exact form of the distribution hardly works on the linear part of the self-field, while the RMS size of the distribution plays an important role. The envelope equations can be used directly in simulations when there is little variability in the distribution during beam transport. Simultaneously, the results of previous research have shown that beams with different distributions but with the same beam emittance terms can be equivalently considered [ 19 19 ]. Generally, it is justified to use the results from the K-V distribution for the other beam distributions. When particles are randomly generated in the K-V transverse ellipsoid, the envelopes where <math id="M6"><mrow><mi>X</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></math> and <math id="M7"><mrow><mi>Y</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></math> are in the two transverse directions of the periodic focusing channels are described by [ 12 12 ] <math display="block" id="M8"><mrow><mtable><mtr><mtd><mrow><msup><mi>X</mi><mo>″</mo></msup><mo>+</mo><msub><mi>κ</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>⋅</mo><mi>X</mi><mo>−</mo><mfrac><mrow><msup><mi>ε</mi><mn>2</mn></msup></mrow><mrow><msup><mi>X</mi><mn>3</mn></msup></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>2</mn><mi>K</mi></mrow><mrow><mi>X</mi><mo>+</mo><mi>Y</mi></mrow></mfrac><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mrow></math> <math display="block" id="M9"><mrow><mtable><mtr><mtd><mrow><msup><mi>Y</mi><mo>″</mo></msup><mo>+</mo><msub><mi>κ</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>⋅</mo><mi>Y</mi><mo>−</mo><mfrac><mrow><msup><mi>ε</mi><mn>2</mn></msup></mrow><mrow><msup><mi>Y</mi><mn>3</mn></msup></mrow></mfrac><mo>−</mo><mfrac><mrow><mn>2</mn><mi>K</mi></mrow><mrow><mi>X</mi><mo>+</mo><mi>Y</mi></mrow></mfrac><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable></mrow></math> where K=2eI/4πε 0 mc 3 β 3 γ 3 is the generalized provenance. Meanwhile, <math id="M10"><mrow><msub><mi>κ</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></math> and <math id="M11"><mrow><msub><mi>κ</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></math> represent the transverse focusing functions, and <math id="M12"><mrow><mi>ε</mi><mo>=</mo><msub><mi>ε</mi><mi>x</mi></msub><mo>=</mo><msub><mi>ε</mi><mi>y</mi></msub></mrow></math> during beam transport. Calculating the stopbands 2 The mismatch between the periodic focusing channels and particles can cause envelope resonance. This is an important factor in space-charge-dominated beams. At present, there are several methods for calculating envelope instability stopbands caused by resonance. First, using the theory established by Struckmeier and Reiser, it is solved by linearizing differential equations with disturbance terms [ 11 11 ]. Second, the perturbation theory based on the Chernin matrix can also obtain a complete set of second-order instabilities, including even and odd modes. The sum effects of the resonance on the beam were also analyzed [ 20 20 ]. Finally, perturbation theory based on the Vlasov-Poisson equations introduced by Hofmann provides an approach for calculating the stopbands of higher-order perturbations [ 5 5 ]. Because the perturbation theory established by Struckmeier and Reiser is relatively mature, widely used, and highly reliable, envelope instability stopbands are calculated based on this theory. With <math id="M13"><mrow><mi>z</mi><mo>=</mo><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><msup><mi>x</mi><mo>′</mo></msup><mo>,</mo><mi>y</mi><mo>,</mo><msup><mi>y</mi><mo>′</mo></msup></mrow><mo>)</mo></mrow></mrow></math> , an ensemble of phase space trajectories is represented which satisfies <math id="M14"><mrow><mi>Z</mi><mrow><mo>(</mo><mrow><mi>s</mi><mo>+</mo><mi>S</mi></mrow><mo>)</mo></mrow><mo>=</mo><mi>M</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mi>Z</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></math> in a periodic dynamics system. It is useful to define M ( s ) as the transfer matrix, where <math id="M15"><mrow><mi>M</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mi>M</mi><mrow><mo>(</mo><mrow><mi>s</mi><mo>+</mo><mi>S</mi></mrow><mo>)</mo></mrow></mrow></math> is also satisfied in a periodic dynamics system. Overall, <math id="M16"><mi>λ</mi></math> denotes the eigenvalues of the transfer matrix that can be greater than 1 or less than 1. If <math id="M17"><mrow><mrow><mo>|</mo><mi>λ</mi><mo>|</mo></mrow><mo>=</mo><mn>1</mn></mrow></math> , the system is stable; otherwise, the system is unstable. Meanwhile, <math id="M18"><mi>λ</mi></math> satisfies <math id="M19"><mrow><mi>λ</mi><mo>=</mo><mrow><mo>|</mo><mi>λ</mi><mo>|</mo></mrow><msup><mi>e</mi><mrow><mi>i</mi><mi>ϕ</mi></mrow></msup></mrow></math> , where <math id="M20"><mi>ϕ</mi></math> represents the eigenphase Arg [ M ( s )] of the transfer matrix. Additionally, <math id="M21"><mi>ϕ</mi></math> connects with the envelope modes [ 21 21 ]. Finally, by converting <math id="M22"><mi>ϕ</mi></math> to σ, a stopband with a specific σ 0t can be obtained under a certain focusing structure. Since σ 0t is different in each period, different values of σ 0t corresponds to different stopbands. Therefore, we need to separately calculate the stopbands corresponding to different σ 0t values and summarize the results. Finally, the stopbands that change with the period can be obtained using Mathematica, as shown in Fig. 2 Fig. 2 [ 22 22 ]. Because σ 0t of schemes ④ and ⑤ are designed to be less than or equal to 90°, there are no stopbands. Hence, this is not shown within this paper. Here, <math id="M23"><mrow><msub><mi>σ</mi><mtext>t</mtext></msub><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mi>x</mi></msub><mo>+</mo><msub><mi>σ</mi><mi>y</mi></msub></mrow><mn>2</mn></mfrac></mrow></math> , which represents an isotropic beam. Simulation results of beam evolution under different focusing schemes 1 TraceWin conducts simulations of any initial beam distribution and three-dimensional space charge force tracking simulation with a wonderful accelerator design capability; therefore, TraceWin was chosen for our simulation [ 23 23 ]. PARTRAN was used for multi-particle simulation and tracking. Subsequently, the K-V distribution was selected because it has many advantages in analyzing the effect of space charge. The initial kinetic energy of the proton beam was 10 MeV. The initial normalized RMS emittance of the CW beam was set to be 0.2 π mm mrad in both the x-x ’ and y-y ’ planes. The number of particles was set as 100000. σ 0t follows a focusing scheme. It is well established that a tune depression of 0.71 represents a dividing line. When the tune depression is greater than 0.71, the beam is dominated by the emittance. When the tune depression is less than 0.71, the beam is dominated by space charge [ 24 24 ]. Accordingly, the range of tune depression values was chosen from 1.0-0.6 to ensure a simulation that adapted well to various situations. By setting the initial conditions, the beam dynamics were simulated using TraceWin code. Here, we mainly discuss the parameters related to the beam quality, such as the RMS emittance and the particle distribution in the phase space. The normalized RMS emittance <math id="M25"><mrow><mi>ε</mi><mo>=</mo><mfrac><mrow><msub><mi>ε</mi><mi>x</mi></msub><mo>+</mo><msub><mi>ε</mi><mi>y</mi></msub></mrow><mn>2</mn></mfrac></mrow></math> is defined because there is not much difference between <math id="M26"><mrow><msub><mi>ε</mi><mi>x</mi></msub></mrow></math> and <math id="M27"><mrow><msub><mi>ε</mi><mi>y</mi></msub></mrow></math> which can also effectively reduce the errors caused by the software calculations. Finally, the <math id="M28"><mi>ε</mi></math> of the output during beam transport was obtained, as shown in Fig. 3 Fig. 3 . It is clear that there are significant differences in <math id="M30"><mi>ε</mi></math> between the quadrupole doublet and solenoid focusing channels. For the quadrupole doublet channels, when tune depression 0.7 and 0.6, the <math id="M31"><mi>ε</mi></math> of focusing scheme ① increased significantly, while the beam designed by other focusing schemes exhibited no such behavior. For the solenoid channels, the tune depression and focusing schemes seem to have no significant effect on the beam quality. This shows that the beam quality can be improved using a periodic focusing solenoid channel. This can be explained as follows. It can be seen from Fig. 2 Fig. 2 that the reduction speed of σ is consistent with σ 0t ; beams transported in the solenoid channel take less time to cross the stopbands under the same focusing scheme and tune depression. With the enhancement of the space charge effect and the increase in the maximum value of σ 0t , this phenomenon becomes more obvious. In stopbands, the beams are affected by envelope instability, which deteriorates the beam quality. Moreover, most of the beams transported in the solenoid channels do not enter the stopbands at the beginning. This implies that the beam is less affected by the initial state. Second, from previous studies by other researchers, envelope instability is caused by linearized perturbed oscillations of the breathing mode and quadrupole mode. The perturbations of the beam envelope are in phase when the breathing mode occurs, whereas out-of-phase perturbations of the beam exist when the quadrupole mode occurs. Parametric resonance can occur only in the breathing mode or quadrupole mode. Certainly, there is another type of resonance, called confluent resonance. Confluent resonance occurs when the two modes occur simultaneously. Confluent resonance means that the resonance effect is between the two envelope modes and focusing lattice. In a solenoid focusing channel, only one breathing or quadrupole mode can occur. This also means that, for a beam propagating in a solenoid, only parametric resonances cause instability. Meanwhile, the decoupling of the breathing and quadrupole modes means that the stopbands of parametric resonance for the two modes do not overlap. In contrast, in a quadrupole doublet focusing channel, only confluent resonance occurs, which involves locking of the breathing and quadrupole modes. The breathing and quadrupole modes are closely related and cannot be strictly separated [ 25 25 ]. If one mode is activated, the other is activated, which can result in a larger emittance. Third, compared to a solenoid channel, in a quadrupole doublet channel, the high-order stopbands overlap more with the envelope instability stopband, which can lead to an increase in emittance. From Fig. 3(a) Fig. 3(a) , it can be seen that when tune depression equals 0.7 or 0.6, the <math id="M32"><mi>ε</mi></math> corresponding to focus scheme ① reaches a maximum, and when tune depression equals 0.6, <math id="M33"><mi>ε</mi></math> has exceeded 10% of the initial emittance. To analyze this, particle distribution diagrams when beams were transported in the quadrupole doublet focusing channel with different tuning depressions are shown in Fig. 4 Fig. 4 . Simultaneously, relevant information about the beam transported in the solenoid focusing channel when the tune depression is equal to 0.6 is shown in Fig. 5 Fig. 5 , corresponding to focusing scheme ①. From Fig. 4 Fig. 4 , it is clear that there is a four-fold structure and envelope instability structure in the phase space distribution of the particles, and there is an obvious growth in <math id="M35"><mi>ε</mi></math> when the tune depression is equal to 0.6. The appearance of these structures is mainly due to the fourth-order resonance and envelope instability. Along with the beam envelope transport, envelope instability is excited from the mismatch generated by the fourth-order resonance [ 26 26 ]. Moreover, the fourth-order resonance appears earlier than the envelope instability. This also means that the fourth-order resonance stopband has some overlapping regions with an envelope instability stopband. Furthermore, as shown in Fig. 2(a) Fig. 2(a) , the beam envelope enters the envelope instability stopband for a relatively long time in the beginning when the tune depression is equal to 0.7 and 0.6, respectively. Hence, it is easy to observe an unstable structure in the phase-space distribution of the particles. In addition, as shown in Fig. 2 Fig. 2 , with a tune depression equal to 0.6, the beam takes more time to cross through the stopband. It seems that the stopband has a longer time to attract the beam owing to the space-charge effect. This further increased the emittance. As for the case in which the tune depression was equal to 0.7, in the first 20 periods, the distribution of particles in the phase space shows a four-fold structure. Hence, there is also a small increase in the emittance. The fourth-order resonance of the particle does not cause instability in the beam envelope, and the stopband has no obvious attraction to the beam. Finally, we analyzed the dynamics of the beam in the solenoid channel, as shown in Fig. 5 Fig. 5 . It is clear that the emittance does not exhibit a significant growth, regardless of the value of the tune depression and the scheme of focusing. Only the distribution of the particles in the phase space was elongated. Simulation results of beam evolution with acceleration 1 The above simulation results are all obtained using pure 2D transverse dynamics, which do not consider acceleration. The beam dynamics when σ 0t locally exceeds 90° under acceleration are ambiguous. In this section, we consider 3D beam dynamics and simulate the evolution of the beam under acceleration. To explore this, we assumed two periodic acceleration channels, where each period element contains, for example, focusing lattices and acceleration structures, as shown schematically in Fig. 6 Fig. 6 . The length of each periodic element was 1236.0 mm, which is similar to that of the periodic element shown in Fig. 1 Fig. 1 . To avoid coupling resonance between different degrees of freedom, it is necessary to optimize the selection of dynamic parameters using equipartition theory, so that the beam oscillation energy is balanced [ 27 27 ]. Under the condition of equipartition, σ 0l and σ 0t not only maintained the same changing trend but also satisfied <math id="M36"><mrow><mfrac><mrow><msub><mi>ε</mi><mrow><mi>l</mi><mi>n</mi></mrow></msub></mrow><mrow><msub><mi>ε</mi><mrow><mi>t</mi><mi>n</mi></mrow></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>σ</mi><mtext>t</mtext></msub></mrow><mrow><msub><mi>σ</mi><mtext>l</mtext></msub></mrow></mfrac></mrow></math> . Therefore, we assume that the maximum value of σ 0l is 2/3 of the maximum value of σ 0t . It is worth noting that σ 0t is consistent with the scheme shown in Fig. 7 Fig. 7 . Based on this, the longitudinal normalized RMS emittance is set to 0.3 [π mm mrad], which falls in the non-resonance region of the Hofmann chart. Subsequently, the number of proton particles was set to 100000, the transverse particle distribution was based on K-V, and the longitudinal particle distribution was based on three Gaussian standard deviations. The simulation results obtained using TraceWin are shown in Fig. 8 Fig. 8 . From Fig. 8 Fig. 8 , we can see that ε increased as the tune depression decreased (the current intensity increased), regardless of which periodic element was used. For the element that provides focusing strength by the quadrupole lattice, ε exhibits a large increase when tune depression was equal to 0.6. There are two reasons for the increase in ε: one is the influence of the space charge, and the other is the nonlinear influence of the accelerating field. Under the action of a non-linear force, particles with larger amplitude will feel a much stronger force, which makes them deviate from the law of harmonic oscillatory motion that leads to an increase in beam emittance. As the beams are focused by the solenoid lattices, <math id="M40"><mi>ε</mi></math> increases slightly. Compared to the results obtained by the 2D simulation, <math id="M41"><mi>ε</mi></math> is mostly consistent. This means that the results obtained by the 2D simulation can be extended to 3D simulations. However, there are some differences. When the tune depression was equal to 0.6, <math id="M42"><mi>ε</mi></math> was larger under focusing scheme ①. Based on this phenomenon, we separately analyzed the beam evolution under focusing schemes ① and ②. Fig. 9 Fig. 9 shows the corresponding particle distribution diagrams when the beam passed through the 9th, 24th, and 80th periods. When σ 0t was changed according to focusing scheme ①, the distribution of particles in the phase space had a typical four-fold structure, and filamentation occurs, but the beam core maintained a four-fold structure within a certain number of periods. Finally, the bunch was transported in the acceleration channel in an unstable structure. Therefore, in this case, the increase in the beam emittance results from the fourth-order resonance and envelope instability. This result was consistent with our conclusions. When σ 0t was changed according to focusing scheme ②, the distribution of particles in the phase space demonstrated a slight filamentation, but this stabilizes as the beam is transported. Conclusion 1 In this study, we analyzed the beam dynamics when σ 0 partially exceeded 90°. Under the condition of equipartition, we used the logistic function to fit σ 0 under acceleration and then calculated the stopband with a non-constant σ. From the simulation of the beam's transverse motion, it is clear that the beam quality is affected by the focusing scheme, tune depression, and periodic focusing channel. By comparing the time at which the beam crosses the stopband and analyzing the particle distribution in the phase space, we concluded that when the tune depression was ≤ 0.7, the beam quality was significantly affected by the space charge. Furthermore, when σ 0t partially exceeded 90°, the emittance of the beam transported in the quadrupole doublet focusing channel increased significantly, which is inseparable from the superposition of the oscillation modes in this focusing structure. When the maximum σ 0t was 120°, the distribution of the particles in the phase space exhibited an obvious fourth-fold structure over a long period of time. In addition, the main reasons for the increase in emittance are the fourth-order resonance and envelope instability. Of course, higher-order instability also affects the beam quality to a certain extent. However, for the solenoid focusing system, there was no significant increase in the emittance under the same conditions. When Tune depression was greater than 0.7, the periodic focusing structure type and focusing scheme demonstrated little effect on the beam quality. We then simulated the beam dynamics while considering acceleration. Beams are focused on solenoid lattices, which seem to have result in high beam quality. However, beams focused on quadrupole doublet lattices show different results. When the tune depression was 0.6, <math id="M44"><mi>ε</mi></math> reaches its maximum value. This is mainly due to the envelope instability and nonlinear forces in the longitudinal phase space. By comparing the results obtained from the 2D simulation with those obtained in 3D, we find that they are highly consistent. Therefore, the results obtained from 2D simulations can be extended to 3D simulations. References J.Y. Liu , J.R. Shi , H. Zha et al. , Analytic RF design of a linear accelerator with a SLED-I type RF pulse compressor . Nucl. Sci. Tech. 31 , 107 ( 2020 ). doi: 10.1007/s41365-020-00815-5 http://doi.org/10.1007/s41365-020-00815-5 Analytic RF design of a linear accelerator with a SLED-I type RF pulse compressor X.L. Lu , Y. Zhang , J.R. Wang et al. , Transport characteristics of space charge-dominated multi-species deuterium beam in electrostatic low-energy beam line . Nucl. Sci. Tech. 29 , 51 ( 2018 ). doi: 10.1007/s41365-018-0384 http://doi.org/10.1007/s41365-018-0384 Transport characteristics of space charge-dominated multi-species deuterium beam in electrostatic low-energy beam line J. Qiang , Bunched Beam Envelope Instability in a Periodic Focusing Channel . J. Phys: Conf. Ser , 1067 , 062015 ( 2018 ). doi: 10.1088/1742-6596/1067/6/062015 http://doi.org/10.1088/1742-6596/1067/6/062015 . Bunched Beam Envelope Instability in a Periodic Focusing Channel H.Y. Li , X.M. Wan , W. Chen et al . Optimization of the S-band side-coupled cavities for proton acceleration . Nucl. Sci. Tech. 31 , 23 ( 2020 ). doi: 10.1007/s41365-020-0735-7 http://doi.org/10.1007/s41365-020-0735-7 Optimization of the S-band side-coupled cavities for proton acceleration I. Hofmann , L. Laslett , L. Smith et al. , Stability of the kapchinskij-vladimirskij (k-v) distribution in long periodic transport systems . Part. Accel , 13 . 145 - 178 ( 1982 ). Stability of the kapchinskij-vladimirskij (k-v) distribution in long periodic transport systems T. Wangler , F. Merrill , L. Rybarcyk et al. , Space charge in proton linacs . AIP Conf. Pro , 448 , 3 - 14 ( 1998 ). doi: 10.1063/1.56764 http://doi.org/10.1063/1.56764 . Space charge in proton linacs Z.H. Li , P. Cheng , H.P. Geng et al. , Physics design of an accelerator for an accelerator-driven subcritical system . Phys. Rev. ST Accel. Beams , 16 . 080101 ( 2013 ). doi: 10.1103/PhysRevSTAB.16.080101 http://doi.org/10.1103/PhysRevSTAB.16.080101 . Physics design of an accelerator for an accelerator-driven subcritical system D.O. Jeon , Classification of space-charge resonances and instabilities in high-intensity linear accelerators . J. Korean Phys. Soc , 72 . 1523 - 1530 ( 2018 ). doi: 10.3938/jkps.72.1523 http://doi.org/10.3938/jkps.72.1523 . Classification of space-charge resonances and instabilities in high-intensity linear accelerators R.A. Jameson , Beam-intensity limitations in linear accelerators . IEEE T. Nucl. Sci , 28 . 2408 - 2412 ( 1981 ). doi: 10.1109/TNS.1981.4331708 http://doi.org/10.1109/TNS.1981.4331708 Beam-intensity limitations in linear accelerators C. Li , Q. Qin , R.A. Jameson et al. , Space Charge Induced Collective Modes and Beam Halo in Periodic Channels . Paper Presented at the 7th Int. Particle Accelerator Conf.(IPAC'16) , Busan, Korea , 8-13May, 2016 . Space Charge Induced Collective Modes and Beam Halo in Periodic Channels J. Struckmeier , M. Reiser , Theoretical studies of envelope oscillations and instabilities of mismatched intense charged particle beams in periodic focusing channels . Part. Accel , 14 . 227 - 260 ( 1983 ). Theoretical studies of envelope oscillations and instabilities of mismatched intense charged particle beams in periodic focusing channels M. Chung , Y.L. Cheon , H. Qin , Linear beam stability in periodic focusing systems: Krein signature and band structure . Nucl. Instrum. Meth. A , 962 . 163708 ( 2020 ). doi: 10.1016/j.nima.2020.163708 http://doi.org/10.1016/j.nima.2020.163708 . Linear beam stability in periodic focusing systems: Krein signature and band structure X.Q. Yan , R.A. Jameson , Y. Lu et al. , Matched and equipartitioned design method for modern high-intensity radio frequency quadrupole accelerators . Nucl. Instrum. Meth. A , 577 . 402 - 408 ( 2007 ). doi: 10.1016/j.nima.2007.03.031 http://doi.org/10.1016/j.nima.2007.03.031 . Matched and equipartitioned design method for modern high-intensity radio frequency quadrupole accelerators R.A. Jameson , An approach to fundamental study of beam loss minimization . Paper Presented at the AIP Conference Proceedings . American Institute of Physics , 480 , 21 - 30 ( 1999 ). doi: 10.1063/1.59501 http://doi.org/10.1063/1.59501 An approach to fundamental study of beam loss minimization P. Paul , G. Sprouse , Superconducting linear accelerators for heavy ions . Comments on Nuclear and Particle Physics , 11 . 217 - 229 ( 1983 ). Superconducting linear accelerators for heavy ions G. Geschonke , Superconducting structures for high intensity linac applications Proceedings of the XVIII International Linear Accelerator Conference , 8 October 1996 . pp. 910 - 914 . doi: 10.5170/CERN-1996-007.910 http://doi.org/10.5170/CERN-1996-007.910 . Superconducting structures for high intensity linac applications M. Weiss , Radio-frequency quadrupole . 5th Advanced Accelerator Physics Course . 959 - 991 ( 1995 ). doi: 10.5170/CERN-1995-006.959 http://doi.org/10.5170/CERN-1995-006.959 Radio-frequency quadrupole D.O. Jeon , J. Jang , H. Jin , Interplay of space-charge fourth order resonance and envelope instability . Nucl. Instrum. Meth. A , 832 . 43 - 50 ( 2016 ). doi: 10.1016/j.nima.2016.06.036 http://doi.org/10.1016/j.nima.2016.06.036 . Interplay of space-charge fourth order resonance and envelope instability F. J. Sacherer , Rms envelope equations with space charge . IEEE T. Nucl. Sci , 18 . 1105 - 1107 ( 1971 ). doi: 10.1109/TNS.1971.4326293 http://doi.org/10.1109/TNS.1971.4326293 Rms envelope equations with space charge D. Chernin , Evolution of rms beam envelopes in transport systems with linear x-y coupling . Part. Accel , 24 , 29 - 44 ( 1988 ). Evolution of rms beam envelopes in transport systems with linear x-y coupling C. Li , R.A. Jameson , Structure resonances due to space charge in periodic focusing channels . Phys. Rev. Accel. Beams , 21 , 024204 ( 2018 ). doi: 10.1103/PhysRevAccelBeams.21.024204 http://doi.org/10.1103/PhysRevAccelBeams.21.024204 . Structure resonances due to space charge in periodic focusing channels Wolfram mathematica , URL http://www.wolfram.com/mathematica/ http://www.wolfram.com/mathematica/ . Accessed 9 Jul 2022 . Wolfram mathematica TraceWin , URL http://irfu.cea.fr/dacm/en/logiciels http://irfu.cea.fr/dacm/en/logiciels . Accessed 9 Jul 2022 . TraceWin T. Wangler , K. Crandall , R. Ryne et al. , Particle-core model for transverse dynamics of beam halo . Physical review special topics-accelerators and beams , 1 , 084201 ( 1998 ). doi: 10.1103/PhysRevSTAB.1.084201 http://doi.org/10.1103/PhysRevSTAB.1.084201 . Particle-core model for transverse dynamics of beam halo Y.S. Yuan , O. Boine-Frankenheim , I. Hofmann , Modeling of second order space charge driven coherent sum and difference instabilities . Physical Review Accelerators and Beams , 20 , 104201 ( 2017 ). doi: 10.1103/PhysRevAccelBeams.20.104201 http://doi.org/10.1103/PhysRevAccelBeams.20.104201 . Modeling of second order space charge driven coherent sum and difference instabilities Y. Cheon , S. Moon , M. Chung et al. , Analysis on the stop band of fourth-order resonance in high-intensity linear accelerators . Physics of Plasmas. 27 , 063105 ( 2020 ). doi: 10.1063/5.0004651 http://doi.org/10.1063/5.0004651 . Analysis on the stop band of fourth-order resonance in high-intensity linear accelerators R.A. Jameson , Equipartitioning in linear accelerators . Los Alamos National Laboratory (Report) LA , Santa Fe, NM, USA , 19 October 1981 . Equipartitioning in linear accelerators Xin-Miao Wan and Zhi-Hui Li contributed to the study conception and design. Material preparation, data collection and analysis were performed by Xue-Jiang Pu, Wei-Hong Huang, Yu-Fei Yang Xin-Miao Wan and Zhi-Qiang Ren. The first draft of the manuscript was written by Xin-Miao Wan and Zhi-Hui Li. All authors commented on previous versions of the manuscript. 10.1007/s41365-022-01104-z This study was supported by the National Natural Science Foundation of China (Nos. 11375122 and 11875197). 2022-04-03 2022-07-11 2022-07-16 2022-09-24 2022-10-24 © 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. 2022 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-022-01104-z http://dx.doi.org/10.1007/s41365-022-01104-z . Springer Nature or its licensor holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law. 2022 Nuclear Science and Techniques 09 33