Excitation of nonlinear second order betatron sidebands
for Knock Out slow extraction at the third-integer resonance
Abstract
Radio Frequency Knock Out resonant slow extraction is a standard method for extracting stored particle beams from synchrotrons by transverse excitation. Classically, the beam is excited with an RF field comprising a frequency band around one of the betatron sidebands. This article demonstrates that the third-integer resonance commonly used for the slow extraction induces nonlinear motion, resulting in the appearance of additional sidebands of higher order at multiples of the betatron tune. Measured and simulated beam spectra are presented, revealing these sidebands and the beam’s response to being excited at first and second order sidebands. The feasibility of using a second order sideband for the purpose of slow extraction is demonstrated. This results in a significant improvement in the temporal structure (spill quality) of the extracted beam, but at the cost of higher excitation power requirements. This is observed both experimentally and in tracking simulations. The mechanism behind the observed improvement is explained using beam dynamics simulations.
I Introduction
Radio Frequency Knock Out (RF-KO) resonant slow extraction is used to extract stored particle beams from synchrotron rings [1, 2]. Therefore, a beam optics near a third-integer resonance is chosen where sextupole fields create a transverse, betatron amplitude dependent instability (separatrix) [3]. For storing and accelerating a beam such instabilities are undesired and investigated to be corrected [4]; but in the context of slow extraction they are exploited and enhanced by dedicated sextupole magnets. The beam is driven into the instability in a controlled manner by increasing the particle amplitude through transverse excitation. The excitation system consists of a radio frequency (RF) signal generator, RF amplifiers, and a stripline kicker [5], inside of which electromagnetic RF fields deflect traversing particles on each turn. When particles reach the separatrix, their motion becomes unbound and the betatron amplitude increases rapidly, such that septa can be used to deflect the extracted particles into an extraction beam line. This spill of extracted particles is then delivered to experiments or used for medical therapy.
I.1 Nonlinear betatron oscillation
With the working point (betatron tune) of the circular accelerator close to a third-integer resonance driven by nonlinear sextupole fields, the three-turn particle dynamics is described by the Kobayashi Hamiltonian [6]
where is the normalized sextupole strength and
(1) | ||||
are the normalized phase space coordinates with the action and angle [3].
The relation to the physical coordinate111Corrected for closed orbit and dispersion and divergence is thereby given by the beta function which describes the optical properties of the lattice at the location .
The (small) quantity is the distance of the tune222Including chromatic detuning: to the third-integer resonance with .
In the following, uppercase symbols will be used to denote absolute tune values, while the lowercase refers to the fractional tune as distance to the nearest integer. At the separatrix, the Kobayashi Hamiltonian reaches a value of .
In the linear case () particles perform harmonic betatron oscillations with a constant phase advance per turn . Near the resonance, however, the betatron oscillations becomes nonlinear and the sextupole field causes an amplitude and phase dependent detuning. In the thin lens and flat beam approximations, the kick induced by the sextupole is . For the angle after the sextupole one can thus write
Considering that the sextupole kick is small and using the small angle approximation
the change in angle follows to first order as
As a result, for the nonlinear betatron oscillation, the phase advance per turn
(2) |
shows an amplitude and phase dependent detuning.
As a particle with constant revolves in phase space, the detuning term effectively modulates the betatron frequency as depicted in figure 1. This frequency modulation generates anharmonic motion and sidebands of higher order.
II Nonlinear beam spectra
under excitation
Spectra of transverse beam motion are obtained by observing the position of the circulating beam as a function of time at a fixed location in the synchrotron. For the linear case, the resulting spectra are well understood and dominated by the betatron sidebands at below and above each harmonic of the revolution frequency . These harmonics are related to the transverse position of the beam centroid [7]. The sidebands are the result of the amplitude modulation of this position due to betatron oscillations and are used to determine linear quantities like tune and chromaticity [8].
In the nonlinear case the betatron oscillation are frequency modulated. Therefore, one expects multiple sidebands of order around each harmonic :
II.1 Measurement of transfer function
A beam transfer function (BTF) measurement is performed by exciting the beam with a fixed frequency and observing the magnitude and phase of the induced beam oscillations filtered at that specific excitation frequency. The measurement is repeated within a single machine cycle for a range of frequencies, yielding the beam response as function of frequency. For a linear system the transfer function is well defined as the complex quotient of the observed and induced oscillation [9]. However, in the nonlinear case the amplitude detuning affects the measurement by altering the oscillation frequencies present in the excited system, as the excitation causes particles to be depleted from the beam core to higher amplitudes, where they are detuned. Since the BTF is solely sensitive to oscillations at the excitation frequency, it can thus only represent a subset of the frequencies and nonlinear dynamics involved and becomes sensitive to the particle distribution and the speed and amplitude of the excitation frequency sweep [10].
Nevertheless, a BTF measurement can be used to probe the beam response at the frequency where the nonlinear sidebands of second order are expected. Such a measurement is presented in figure 2 for a typical RF-KO extraction setup from the Heavy Ion Synchrotron (SIS18) at GSI using a vector network analyzer connected to a stipline kicker and pickup. Alongside the nonlinear case with sextupoles, also the linear case with sextupole magnets switched off is shown. As the closed orbit was not corrected, the sextupole magnets cause a tune shift. Therefore the location of the first and second order sidebands is indicated for both cases as obtained by fitting of the respective phase response.
A response of the beam to an excitation near the second order betatron sideband at can only be observed for the nonlinear case with sextupoles switched on. The response is measured at the harmonic because of the higher transfer impedance of the terminated pickup. The measurement demonstrates the existence of the second order sideband and the feasibility of exciting it. While the response is very weak compared to the first order sideband, one has to keep in mind that the BTF does not give insight into the beam oscillations at any but the excited frequency. For example, spectral components near the second order sideband induced by excitation at the first order sideband are not visible. However, the dip in the response at the first order sideband suggest, that in this case the primary oscillation frequencies are shifted away from the excitation frequency where the BTF is not sensitive [10]. Therefore, in the following, the beam motion is investigated by means of Schottky spectra in simulations.
II.2 Simulation of Schottky spectra
Schottky diagnostics is a vital tool to study the incoherent particle motion by sampling the beam position at frequencies higher than the revolution frequency [11]. Unlike the BTF, it provides a complete picture of the intrinsic beam oscillations independent of an external excitation. A particle tracking simulation is performed with Xsuite [12] for a typical machine setup of a Knock Out extraction from SIS18 at GSI. A beam of ions with a rigidity of and revolution frequency of is prepared for extraction using a machine tune of and a corrected orbit.
Figure 3 shows the spectra of the horizontal beam centroid as obtained from the simulation using a sampling frequency of .
In the linear case (sextupoles off), the regular betatron sidebands of first order are visible, which are relatively broad due to the large chromaticity of .
When the sextupoles are switched on for the purpose of slow extraction, the chromaticity reduces to and the first order sidebands become narrower. In addition, the sidebands of higher order () become visible, most prominently the second order sidebands at .
But also the weaker third order sideband at and faint coupling bands at can be observed.
The latter result from the fact that the sextupole fields also depend on the vertical coordinate, thus coupling the betatron motion in both planes.
The observation that the emergence of second and higher order sidebands is related to the nonlinear dynamics discussed in section I.1 is further emphasized in figure 4.
Here, Schottky spectra of beam slices for discrete ranges of the Hamiltonian from the core to the separatrix are depicted (compare also figure 1).
The plot shows not only how the frequency of the sidebands shifts due to the amplitude detuning,
but also reveals a clear correlation of the particles’ Hamiltonian with the strength of the nonlinear second and third order sidebands.
An external excitation at the second order sideband will thus couple strongly to those particles close to the separatrix undergoing strong nonlinear motion, while only weakly affecting particles in the core.
In RF-KO extraction, the beam is excited transversely with a signal of a certain bandwidth to account for its intrinsic tune spread. While many different signal types can be used depending on the application, a band-limited noise signal is used here to avoid the bias of any specific non-uniform excitation spectrum. Under the excitation, the betatron oscillation amplitude is increased, and so does the magnitude in the respective Schottky spectra in figure 3. When the amplitude grows beyond the separatrix where particles are extracted, a significant DC leakage into the low frequencies is observed in the Schottky spectra. This is due to the slow variation in centroid position compared to the case without extraction, where the mean beam position remains at zero.
Traditionally, the excitation band is placed in the vicinity of one of the first order betatron sidebands. This is the consequence of modelling the excitation as a forced oscillation of individual particles undergoing simple harmonic motion under linear restoring force. Motivated by the nonlinear dynamics discussed above, a second case is considered where the excitation band is placed near the second order betatron sideband. The choice of two nearby sidebands around avoids a systematic bias on the spill quality stemming from a large difference of excitation frequencies. It also allows to distinguish the behaviour of the excited bands from the corresponding basebands around . As can be seen in figure 3, in both cases the excitation induces coherent beam oscillations only at the absolute excitation frequency, whereas the spectral response at any other betatron sideband at the other harmonics remains one order of magnitude weaker. This becomes especially obvious for the excitation band placed near , for which the oscillation magnitude of the excited second order sideband is increased even beyond the magnitude of the nearby first order sideband . Such a behaviour is only observed at the excited sideband, and not at any of the other sidebands, which are not directly excited. In general, the overall magnitudes observed in the Schottky spectrum for the case of excitation near the second order sideband are smaller compared to the case of excitation at the first order sideband , even though in the former case the amplitude of the excitation signal used in the simulation is larger to achieve an equal extraction rate. This suggests that the coupling to the beam is weaker at the higher order sidebands, but it is still sufficient for driving an RF-KO extraction, as demonstrated in the upcoming section.
III Knock Out extraction with second order sidebands
To study the effect of an excitation at a second order sideband on the RF-KO slow extraction process, the required excitation power and achievable spill quality is investigated. The spill quality is given by the amount of unwanted intensity fluctuations in the time structure of the extracted beam. To quantify the fluctuations, the spill is divided into time intervals of length and the number of extracted particles is counted in each interval. The standard deviation and mean , evaluated over a larger time span , is then used to determine the coefficient of variation or the equivalent spill duty factor [13]:
In addition to the spill quality, also the extraction efficiency is of importance.
To extract all particles from the storage ring, a certain excitation power is required which depends strongly on the beam rigidity and the excitation signal used.
In practice the available power of the RF amplification system can impose limits on the excitation signals which are feasible.
III.1 Machine measurement
The effects of excitation frequency and bandwidth on the RF-KO extraction are studied experimentally at the SIS18 of GSI. An beam with a kinetic energy of (, ) is prepared and excited with a band-limited noise signal. The extracted spill is observed with a particle detector based on a plastic scintillator [14], which serves two purposes: First, the signal is used with a feedback system to dynamically adjust the excitation signal amplitude over the duration of the spill in order to maintain a constant spill rate of [15]. Secondly, the intensity fluctuations measured by the detector are evaluated to determine the spill quality. Thereby was chosen to provide sufficient counting statistics for the quality analysis.
Figure 5 shows the results of a large scan of the excitation frequency and bandwidth, exceeding the traditionally used vicinity of the first order betatron sideband. The measurement proves, that RF-KO slow extraction by exciting the second order betatron sideband is possible, and a comparable number of particles can be extracted, namely . The excitation level required is with about twice as large as for the first order sideband. The optimal choice of excitation frequency and bandwidth from an efficiency point of view hence lies in the close vicinity of the first order betatron sideband.
However, the optimum in terms of spill quality was found on the opposite side of the resonance half way towards the second order betatron sideband (figure 5, left). Here, the measurement shows a reduction of the intensity fluctuations to compared to for the efficiency optimum at the first order sideband. This improvement in spill quality comes with an increased power requirement. While the transverse RF excitation system is capable of delivering these demands for the presented conditions, in general for higher rigidity beams a trade-off between spill quality and required power has to be found.
III.2 Simulation of spill quality
The previously introduced simulation code is used to model the machine experiment. Figure 6 shows the spill quality and required excitation power as a function of the excitation frequency as obtained from the simulation. Thereby, the bandwidth of the band-limited noise signal remains fixed. In agreement with the measurement, the global optimum in terms of spill quality is located on the opposite side of the resonance towards the second order betatron sideband. While the spill fluctuations are significantly reduced compared to an excitation band near the first order betatron sideband, the required excitation power is about higher, corresponding to a factor three larger signal amplitude. The increase in the required excitation power is consistent with the experimental findings and the observations from the Schottky spectra discussed in the previous section.
Another local optimum can be found towards the fourth order betatron sideband , where the excitation power required is even higher. Once the power limit is reached, only a small fraction of the beam can be extracted and the value of increases again.
III.3 Simulation of beam dynamics
To understand the mechanism behind the observed spill quality improvement, the beam dynamics is analysed in the particle tracking simulations. For the two excitation bands near the first and second order betatron sidebands discussed above, the particle distribution in horizontal phase space is recorded as a function of time. The Kobayashi Hamiltonian normalized to its value at the separatrix is used to analyse the particle distribution during the extraction by means of a histogram (figure 7). Thereby, corresponds to a particle in the centre of phase space with zero betatron amplitude and (chromatic) tune ; and for the betatron amplitude and nonlinear detuning increases (compare figure 1).
For both cases, the initial distribution is the same before the excitation is switched on at . The excitation band near mainly couples to particles which perform linear betatron oscillations and are located in the core of the phase space. When these particles are excited by the band-limited signal, they diffuse outwards in phase space towards the separatrix, causing the distribution to become flatter (figure 7, top). This effectively increases the particle density in the vicinity of the separatrix, which makes the process more vulnerable to fluctuations of the separatrix size caused by power supply ripples [13]. It also intensifies fluctuations of the spill intensity caused by the diffusive nature of the noise excitation itself during the transition of particles across the separatrix. On the other hand, the higher population of the phase space relatively close to the separatrix means, that a lower excitation strength is required to continually feed the extraction process with the desired particle rate. This can be understood as an extraction driven from within the core of the phase space.
In contrast, the excitation band near can only couple to particles whose betatron oscillations are anharmonic and frequency modulated, since only then they have a spectral component at this frequency range. This applies to particles towards the separatrix which undergo amplitude and phase detuning (see section I.1). As a result, the beam is extracted from the outside and the phase space distribution becomes steeper towards the separatrix, while the particle density in the core at is almost unaffected (figure 7, bottom). In this case, the particle density in the vicinity of the separatrix does not increase during the slow extraction, such that spill intensity fluctuations are effectively reduced compared to the former case. For an equal average extraction rate, the reduced density also demands a faster flow of particles across the separatrix, which likewise reduces the sensitivity of the process to ripples. However, it also means that a higher excitation strength is required to make particles from the core reach and cross the separatrix in order to maintain the same extraction rate.
IV Conclusion
The nonlinear particle dynamics at the third-integer resonance used for resonant slow extraction induces a frequency modulation of betatron oscillations, which results in the appearance of additional betatron sidebands of higher order.
The second order betatron sidebands are observed prominently in BTF measurements and Schottky spectra and can be used for transverse excitation of the beam.
For RF-KO slow extraction, excitation of these second order sidebands requires a larger excitation power, but improves the spill quality significantly.
As the excitation at the nonlinear second order band couples to particles near the separatrix, it effectively
extracts the beam by pushing out particles from the outside rather than from the inside of the phase space distribution.
Thereby, the accumulation of particles in the vicinity of the separatrix is reduced and at the same time the flow of particles across the separatrix is sped up, which reduces the vulnerability to processes diminishing the spill quality.
To improve the efficiency of the excitation process when sufficient power is not available, hybrid excitation signals can be used. For example, a weak noise-band at the sideband feeding the extraction can be combined with a stronger excitation signal placed on the other side of the resonance towards the sideband to benefit from the described spill quality improvement mechanism. This applies likewise to the components of more advanced excitation signals like Noise ++, a combination of a band-limited noise with sinusoidal excitation signals, which is described in detail in [16]. In figure 8 the spill quality is shown as a function of the frequency of the noise and one of the sine components of the Noise ++ excitation signal used at SIS18 at GSI for a beam at . For the noise component, frequencies between the first order sideband and the resonance are detrimental to the spill quality, and it is preferable to place the noise band at a sufficient distance to the band or on the other side of the resonance. For the frequency of the sine component a local optimum exists not only at towards the first order sideband, but also at towards the second order sideband on the other side of the resonance. With and the observed behaviour on the left side of the resonance is essentially a mirror image of the right side, scaled by a factor of two.
In general, finding the optimal excitation frequencies is an optimization problem.
By setting appropriate boundaries, one can allow the optimization algorithm to explore also the region of the second order sideband and then chose the local optimum according to the available excitation power.
Acknowledgements.
The beam instrumentation department at GSI is acknowledged for the support in carrying out the measurements. We thank C. Cortés for valuable discussions on nonlinear BTF measurements.References
- Hiramoto and Nishi [1992] K. Hiramoto and M. Nishi, Resonant beam extraction scheme with constant separatrix, Nucl. Instrum. Methods Phys. Res., Sect. A 322, 154 (1992).
- Tomizawa et al. [1993] M. Tomizawa, M. Yoshizawa, K. Chida, J. Yoshizawa, Y. Arakaki, R. Nagai, A. Mizobuchi, A. Noda, K. Noda, M. Kanazawa, A. Ando, H. Muto, and T. Hattori, Slow beam extraction at tarn ii, Nucl. Instrum. Methods Phys. Res., Sect. A 326, 399 (1993).
- Benedikt et al. [1999] M. Benedikt, P. J. Bryant, L. Badano, M. Crescenti, P. Holy, A. T. Maier, M. Pullia, S. Rossi, and P. Knaus, Proton-Ion Medical Machine Study (PIMMS): Part I, Tech. Rep. CERN/PS 99-010 (DI) (Geneva, 1999).
- Franchi et al. [2014] A. Franchi, L. Farvacque, F. Ewald, G. Le Bec, and K. B. Scheidt, First simultaneous measurement of sextupolar and octupolar resonance driving terms in a circular accelerator from turn-by-turn beam position monitor data, Phys. Rev. ST Accel. Beams 17, 074001 (2014).
- Belver-Aguilar et al. [2014] C. Belver-Aguilar, A. Faus-Golfe, F. Toral, and M. J. Barnes, Stripline design for the extraction kicker of compact linear collider damping rings, Phys. Rev. ST Accel. Beams 17, 071003 (2014).
- Kobayashi [1970] Y. Kobayashi, Theory of the resonant beam ejection from synchrotrons, Nuclear Instruments and Methods 83, 77 (1970).
- Singh et al. [2018a] R. Singh, P. Forck, P. Kowina, and A. Reiter, Non-destructive position measurement of coasting beam using Schottky signals, Review of Scientific Instruments 89, 113305 (2018a).
- Betz et al. [2017] M. Betz, O. R. Jones, T. Lefevre, and M. Wendt, Bunched-beam schottky monitoring in the lhc, Nucl. Instrum. Methods Phys. Res., Sect. A 874, 113 (2017).
- Borer et al. [1979] J. Borer, G. Guignard, A. Hofmann, E. Peschardt, F. Sacherer, and B. Zotter, Information from beam response to longitudinal and transverse excitation, IEEE Trans. Nucl. Sci. 26, 3405 (1979).
- [10] E. C. Cortés García, P. Niedermayer, R. Singh, R. Taylor, E. Benedetto, E. Feldmeier, M. Hun, and T. Haberer, Interpretation of the horizontal beam response near the third integer resonance, preprint.
- Boussard [1995] D. Boussard, Schottky noise and beam transfer function diagnostics, in Advanced Accelerator Physics Course, CAS - CERN Accelerator School (CERN, Rhodes, Greece, 1995) pp. 749–782.
- Iadarola et al. [2023] G. Iadarola, R. D. Maria, S. Łopaciuk, A. Abramov, X. Buffat, D. Demetriadou, L. Deniau, P. D. Hermes, P. Kicsiny, P. Kruyt, A. Latina, L. Mether, K. Paraschou, G. Sterbini, F. F. Van der Veken, P. Belanger, P. Niedermayer, D. Di Croce, T. Pieloni, L. Van Riesen-Haupt, and M. Seidel, Xsuite: An integrated beam physics simulation framework, in Proc. HB’23 (JACoW Publishing, Geneva, Switzerland, 2023) pp. 73–80.
- Singh et al. [2018b] R. Singh, P. Forck, P. Boutachkov, S. Sorge, and H. Welker, Slow extraction spill characterization from micro to milli-second scale, J. Phys. Conf. Ser. 1067, 072002 (2018b).
- Forck et al. [1997] P. Forck, T. Hoffmann, and A. Peters, Detectors for slowly extracted heavy ions at the gsi facility, in Proc. DIPAC’97, European Workshop on Diagnostics and Instrumentation for Particle Accelerators (INFN-LNF, Frascati, Italy, 1997) pp. 165–167.
- Niedermayer et al. [2023] P. J. Niedermayer, R. Singh, and R. Geißler, Software-defined radio based feedback system for beam spill control in particle accelerators, in Proc. GRCon23, GNU Radio Conference (GNU Radio Journal, Tempe, AZ, USA, 2023).
- Niedermayer and Singh [2024] P. Niedermayer and R. Singh, Excitation signal optimization for minimizing fluctuations in knock out slow extraction, Sci. Rep. 14, 10310 (2024).