Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Excitation of nonlinear second order betatron sidebands
for Knock Out slow extraction at the third-integer resonance

Philipp Niedermayer \orcidlink0000-0002-4722-6619 p.niedermayer@gsi.de    Rahul Singh \orcidlink0000-0002-7427-8000 GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
(August 1, 2024)
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]

H𝐻\displaystyle Hitalic_H =3πd(X2+X2)+S4(3XX2X3)absent3𝜋𝑑superscript𝑋2superscript𝑋2𝑆43𝑋superscript𝑋2superscript𝑋3\displaystyle=3\pi d\left(X^{2}+X^{\prime 2}\right)+\frac{S}{4}\left(3XX^{% \prime 2}-X^{3}\right)= 3 italic_π italic_d ( italic_X start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_X start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT ) + divide start_ARG italic_S end_ARG start_ARG 4 end_ARG ( 3 italic_X italic_X start_POSTSUPERSCRIPT ′ 2 end_POSTSUPERSCRIPT - italic_X start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT )
=6πdJxS2Jx3/2cos(3Θx)absent6𝜋𝑑subscript𝐽x𝑆2superscriptsubscript𝐽x323subscriptΘx\displaystyle=6\pi dJ_{\text{x}}-\frac{S}{\sqrt{2}}J_{\text{x}}^{3/2}\cos(3% \Theta_{\text{x}})= 6 italic_π italic_d italic_J start_POSTSUBSCRIPT x end_POSTSUBSCRIPT - divide start_ARG italic_S end_ARG start_ARG square-root start_ARG 2 end_ARG end_ARG italic_J start_POSTSUBSCRIPT x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT roman_cos ( start_ARG 3 roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG )

where S=βx(s)3/2k2l/2𝑆subscript𝛽𝑥superscript𝑠32subscript𝑘2𝑙2S=-\beta_{x}(s)^{3/2}k_{2}l/2italic_S = - italic_β start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ( italic_s ) start_POSTSUPERSCRIPT 3 / 2 end_POSTSUPERSCRIPT italic_k start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_l / 2 is the normalized sextupole strength and

X𝑋\displaystyle Xitalic_X =2Jxcos(Θx)=1βx(s)x(s)absent2subscript𝐽xsubscriptΘx1subscript𝛽x𝑠𝑥𝑠\displaystyle=\sqrt{2J_{\text{x}}}\cos(\Theta_{\text{x}})=\frac{1}{\sqrt{\beta% _{\text{x}}(s)}}\,x(s)= square-root start_ARG 2 italic_J start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG roman_cos ( start_ARG roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG ) = divide start_ARG 1 end_ARG start_ARG square-root start_ARG italic_β start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ( italic_s ) end_ARG end_ARG italic_x ( italic_s ) (1)
Xsuperscript𝑋\displaystyle X^{\prime}italic_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT =2Jxsin(Θx)=βx(s)2βx(s)x(s)+βx(s)x(s)absent2subscript𝐽xsubscriptΘxsubscriptsuperscript𝛽x𝑠2subscript𝛽x𝑠𝑥𝑠subscript𝛽x𝑠superscript𝑥𝑠\displaystyle=-\sqrt{2J_{\text{x}}}\sin(\Theta_{\text{x}})=-\frac{\beta^{% \prime}_{\text{x}}(s)}{2\sqrt{\beta_{\text{x}}(s)}}\,x(s)+\sqrt{\beta_{\text{x% }}(s)}\,x^{\prime}(s)= - square-root start_ARG 2 italic_J start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG roman_sin ( start_ARG roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG ) = - divide start_ARG italic_β start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ( italic_s ) end_ARG start_ARG 2 square-root start_ARG italic_β start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ( italic_s ) end_ARG end_ARG italic_x ( italic_s ) + square-root start_ARG italic_β start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ( italic_s ) end_ARG italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ( italic_s )

are the normalized phase space coordinates with the action Jxsubscript𝐽xJ_{\text{x}}italic_J start_POSTSUBSCRIPT x end_POSTSUBSCRIPT and angle ΘxsubscriptΘx\Theta_{\text{x}}roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT [3]. The relation to the physical coordinate111Corrected for closed orbit and dispersion x𝑥xitalic_x and divergence xsuperscript𝑥x^{\prime}italic_x start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT is thereby given by the beta function βx(s)subscript𝛽x𝑠\beta_{\text{x}}(s)italic_β start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ( italic_s ) which describes the optical properties of the lattice at the location s𝑠sitalic_s. The (small) quantity d=QxQres𝑑subscript𝑄xsubscript𝑄resd=Q_{\text{x}}-Q_{\text{res}}italic_d = italic_Q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT - italic_Q start_POSTSUBSCRIPT res end_POSTSUBSCRIPT is the distance of the tune222Including chromatic detuning: Qx=Qx,ref+ξxδsubscript𝑄xsubscript𝑄xrefsubscript𝜉x𝛿Q_{\text{x}}=Q_{\text{x},\text{ref}}+\xi_{\text{x}}\deltaitalic_Q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT = italic_Q start_POSTSUBSCRIPT x , ref end_POSTSUBSCRIPT + italic_ξ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT italic_δ Qxsubscript𝑄xQ_{\text{x}}italic_Q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT to the third-integer resonance Qressubscript𝑄resQ_{\text{res}}italic_Q start_POSTSUBSCRIPT res end_POSTSUBSCRIPT with 3Qres3subscript𝑄res3Q_{\text{res}}\in\mathbb{N}3 italic_Q start_POSTSUBSCRIPT res end_POSTSUBSCRIPT ∈ blackboard_N. In the following, uppercase symbols will be used to denote absolute tune values, while the lowercase q0.5𝑞0.5q\leq 0.5italic_q ≤ 0.5 refers to the fractional tune as distance to the nearest integer. At the separatrix, the Kobayashi Hamiltonian reaches a value of Hsep=(4πd)3/S2subscript𝐻sepsuperscript4𝜋𝑑3superscript𝑆2H_{\text{sep}}={\left(4\pi d\right)^{3}}/{S^{2}}italic_H start_POSTSUBSCRIPT sep end_POSTSUBSCRIPT = ( 4 italic_π italic_d ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT / italic_S start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT.

In the linear case (S=0𝑆0S=0italic_S = 0) particles perform harmonic betatron oscillations with a constant phase advance per turn μ1=2πQxsubscript𝜇12𝜋subscript𝑄x\mu_{1}=2\pi Q_{\text{x}}italic_μ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 2 italic_π italic_Q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT. 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 ΔX=SX2Δsuperscript𝑋𝑆superscript𝑋2\Delta X^{\prime}=SX^{2}roman_Δ italic_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = italic_S italic_X start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. For the angle after the sextupole one can thus write

tan(Θx,1)=X+SX2X=tan(Θx)S2Jxcos(Θx)subscriptΘx1superscript𝑋𝑆superscript𝑋2𝑋subscriptΘx𝑆2subscript𝐽xsubscriptΘx\tan(\Theta_{\text{x},1})=-\frac{X^{\prime}+SX^{2}}{X}=\tan(\Theta_{\text{x}})% -S\sqrt{2J_{\text{x}}}\cos(\Theta_{\text{x}})roman_tan ( start_ARG roman_Θ start_POSTSUBSCRIPT x , 1 end_POSTSUBSCRIPT end_ARG ) = - divide start_ARG italic_X start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + italic_S italic_X start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_X end_ARG = roman_tan ( start_ARG roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG ) - italic_S square-root start_ARG 2 italic_J start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG roman_cos ( start_ARG roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG )

Considering that the sextupole kick is small and using the small angle approximation

tan(Θx+ΔΘx)tan(Θx)+ΔΘxcos2(Θx)+(ΔΘx2)subscriptΘxΔsubscriptΘxsubscriptΘxΔsubscriptΘxsuperscript2subscriptΘxorderΔsuperscriptsubscriptΘx2\tan(\Theta_{\text{x}}+\Delta\Theta_{\text{x}})\approx\tan(\Theta_{\text{x}})+% \frac{\Delta\Theta_{\text{x}}}{\cos^{2}(\Theta_{\text{x}})}+\order{\Delta% \Theta_{\text{x}}^{2}}roman_tan ( start_ARG roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT + roman_Δ roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG ) ≈ roman_tan ( start_ARG roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG ) + divide start_ARG roman_Δ roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG start_ARG roman_cos start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ) end_ARG + ( start_ARG roman_Δ roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG )

the change in angle follows to first order as

ΔΘx=Θx,1Θx=S2Jxcos3(Θx)1ΔsubscriptΘxsubscriptΘx1subscriptΘx𝑆2subscript𝐽xsuperscript3subscriptΘxmuch-less-than1\Delta\Theta_{\text{x}}=\Theta_{\text{x},1}-\Theta_{\text{x}}=-S\sqrt{2J_{% \text{x}}}\cos^{3}(\Theta_{\text{x}})\ll 1roman_Δ roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT = roman_Θ start_POSTSUBSCRIPT x , 1 end_POSTSUBSCRIPT - roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT = - italic_S square-root start_ARG 2 italic_J start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG roman_cos start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ) ≪ 1

As a result, for the nonlinear betatron oscillation, the phase advance per turn

μ~1=2πQxS2Jxcos3(Θx)subscript~𝜇12𝜋subscript𝑄x𝑆2subscript𝐽xsuperscript3subscriptΘx\tilde{\mu}_{1}=2\pi Q_{\text{x}}-S\sqrt{2J_{\text{x}}}\cos^{3}(\Theta_{\text{% x}})\\ over~ start_ARG italic_μ end_ARG start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 2 italic_π italic_Q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT - italic_S square-root start_ARG 2 italic_J start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_ARG roman_cos start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT ( roman_Θ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ) (2)

shows an amplitude and phase dependent detuning.

As a particle with constant H𝐻Hitalic_H 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.

Refer to caption
Figure 1: Nonlinear detuning of the phase advance per turn as function of angle, calculated with equation 2 for constant H=0𝐻0H=0italic_H = 0 (blue) to H=Hsep𝐻subscript𝐻sepH=H_{\text{sep}}italic_H = italic_H start_POSTSUBSCRIPT sep end_POSTSUBSCRIPT (purple). The dashed line marks the resonance condition. The inset shows the corresponding phase space trajectories calculated with equation 1.

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 f=(n±qx)frev𝑓plus-or-minus𝑛subscript𝑞xsubscript𝑓revf=(n\pm q_{\text{x}})f_{\text{rev}}italic_f = ( italic_n ± italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ) italic_f start_POSTSUBSCRIPT rev end_POSTSUBSCRIPT below and above each harmonic n𝑛n\in\mathbb{N}italic_n ∈ blackboard_N of the revolution frequency frevsubscript𝑓revf_{\text{rev}}italic_f start_POSTSUBSCRIPT rev end_POSTSUBSCRIPT. 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 k𝑘k\in\mathbb{N}italic_k ∈ blackboard_N around each harmonic n𝑛nitalic_n:

f=(n±kqx)frev𝑓plus-or-minus𝑛𝑘subscript𝑞xsubscript𝑓rev{f}=\left(n\pm k\,q_{\text{x}}\right){f_{\text{rev}}}italic_f = ( italic_n ± italic_k italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ) italic_f start_POSTSUBSCRIPT rev end_POSTSUBSCRIPT

II.1 Measurement of transfer function

A beam transfer function (BTF) measurement is performed by exciting the beam with a fixed frequency fexsubscript𝑓exf_{\text{ex}}italic_f start_POSTSUBSCRIPT ex end_POSTSUBSCRIPT 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.

Refer to caption
Figure 2: Measured response to a sinusoidal excitation of varying frequency (BTF measurement) of an Au19765+{}^{197}\text{Au}{}^{65+}start_FLOATSUPERSCRIPT 197 end_FLOATSUPERSCRIPT Au start_FLOATSUPERSCRIPT 65 + end_FLOATSUPERSCRIPT beam at 800 MeV/nucleontimes800MeVnucleon800\text{\,}\mathrm{M}\mathrm{e}\mathrm{V}\mathrm{/}\mathrm{n}\mathrm{u}% \mathrm{c}\mathrm{l}\mathrm{e}\mathrm{o}\mathrm{n}start_ARG 800 end_ARG start_ARG times end_ARG start_ARG roman_MeV / roman_nucleon end_ARG. Vertical lines indicate the location of the betatron sidebands and third-integer resonance.
Refer to caption
Figure 3: Transverse Schottky spectra for the linear case (sextupoles off) and three nonlinear cases (sextupoles on) with and without excitation as obtained from a simulation. The zoom parts show beam oscillations at first and second order sidebands (vertical black lines) and at the excitation bands (yellow and blue shaded regions).

A response of the beam to an excitation near the second order betatron sideband at 102qx102subscript𝑞x10-2q_{\text{x}}10 - 2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT can only be observed for the nonlinear case with sextupoles switched on. The response is measured at the 10thsuperscript10th10^{\text{th}}10 start_POSTSUPERSCRIPT th end_POSTSUPERSCRIPT harmonic because of the higher transfer impedance of the 50 times50ohm50\text{\,}\mathrm{\SIUnitSymbolOhm}start_ARG 50 end_ARG start_ARG times end_ARG start_ARG roman_Ω end_ARG 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 U23873+{}^{238}\text{U}{}^{73+}start_FLOATSUPERSCRIPT 238 end_FLOATSUPERSCRIPT U start_FLOATSUPERSCRIPT 73 + end_FLOATSUPERSCRIPT ions with a rigidity of Bρ=7 T m𝐵𝜌times7timesteslameterB\rho=$7\text{\,}\mathrm{T}\text{\,}\mathrm{m}$italic_B italic_ρ = start_ARG 7 end_ARG start_ARG times end_ARG start_ARG start_ARG roman_T end_ARG start_ARG times end_ARG start_ARG roman_m end_ARG end_ARG and revolution frequency of frev=785 kHzsubscript𝑓revtimes785kHzf_{\text{rev}}=$785\text{\,}\mathrm{k}\mathrm{H}\mathrm{z}$italic_f start_POSTSUBSCRIPT rev end_POSTSUBSCRIPT = start_ARG 785 end_ARG start_ARG times end_ARG start_ARG roman_kHz end_ARG is prepared for extraction using a machine tune of Qx,ref=4.327subscript𝑄xref4.327Q_{\text{x},\text{ref}}=4.327italic_Q start_POSTSUBSCRIPT x , ref end_POSTSUBSCRIPT = 4.327 and a corrected orbit.

Figure 3 shows the spectra of the horizontal beam centroid xdelimited-⟨⟩𝑥\left\langle{x}\right\rangle⟨ italic_x ⟩ as obtained from the simulation using a sampling frequency of 3.92 MHz=5frevtimes3.92MHz5subscript𝑓rev$3.92\text{\,}\mathrm{M}\mathrm{H}\mathrm{z}$=5f_{\text{rev}}start_ARG 3.92 end_ARG start_ARG times end_ARG start_ARG roman_MHz end_ARG = 5 italic_f start_POSTSUBSCRIPT rev end_POSTSUBSCRIPT. 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 ξx=ΔQx/δ=6.6subscript𝜉xΔsubscript𝑄x𝛿6.6\xi_{\text{x}}=\Delta Q_{\text{x}}/\delta=-6.6italic_ξ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT = roman_Δ italic_Q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT / italic_δ = - 6.6. When the sextupoles are switched on for the purpose of slow extraction, the chromaticity reduces to ξx=2.0subscript𝜉x2.0\xi_{\text{x}}=-2.0italic_ξ start_POSTSUBSCRIPT x end_POSTSUBSCRIPT = - 2.0 and the first order sidebands become narrower. In addition, the sidebands of higher order (k2𝑘2k\geq 2italic_k ≥ 2) become visible, most prominently the second order sidebands at n±2qxplus-or-minus𝑛2subscript𝑞xn\pm 2q_{\text{x}}italic_n ± 2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT. But also the weaker third order sideband at 13qx=3|d|013subscript𝑞x3𝑑01-3q_{\text{x}}=3\absolutevalue{d}\approx 01 - 3 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT = 3 | start_ARG italic_d end_ARG | ≈ 0 and faint coupling bands at n±2qyplus-or-minus𝑛2subscript𝑞yn\pm 2q_{\text{y}}italic_n ± 2 italic_q start_POSTSUBSCRIPT y end_POSTSUBSCRIPT 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.

Refer to caption
Figure 4: Transverse Schottky spectra (sextupoles on, excitation off) as function of phase space amplitude by means of the Kobayashi Hamiltonian H𝐻Hitalic_H as obtained from a simulation. The shading corresponds to four beam slices of particles in the range hi<H/Hsep<hi+1subscript𝑖𝐻subscript𝐻sepsubscript𝑖1h_{i}<H/H_{\text{sep}}<h_{i+1}italic_h start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT < italic_H / italic_H start_POSTSUBSCRIPT sep end_POSTSUBSCRIPT < italic_h start_POSTSUBSCRIPT italic_i + 1 end_POSTSUBSCRIPT with hi=(i/4)2subscript𝑖superscript𝑖42h_{i}=(i/4)^{2}italic_h start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = ( italic_i / 4 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and i𝑖i\in\mathbb{N}italic_i ∈ blackboard_N. The inset plot shows these ranges in normalized phase space.

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 f/frev1qx2qx2/3𝑓subscript𝑓rev1subscript𝑞x2subscript𝑞x23f/f_{\text{rev}}\approx 1-q_{\text{x}}\approx 2q_{\text{x}}\approx 2/3italic_f / italic_f start_POSTSUBSCRIPT rev end_POSTSUBSCRIPT ≈ 1 - italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ≈ 2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT ≈ 2 / 3 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 f/frev1/3𝑓subscript𝑓rev13f/f_{\text{rev}}\approx 1/3italic_f / italic_f start_POSTSUBSCRIPT rev end_POSTSUBSCRIPT ≈ 1 / 3. 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 2qx2subscript𝑞x2q_{\text{x}}2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT, for which the oscillation magnitude of the excited second order sideband is increased even beyond the magnitude of the nearby first order sideband 1qx1subscript𝑞x1-q_{\text{x}}1 - italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT. 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 2qx2subscript𝑞x2q_{\text{x}}2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT are smaller compared to the case of excitation at the first order sideband 1qx1subscript𝑞x1-q_{\text{x}}1 - italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT, 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 ΔtcountΔsubscript𝑡count{\Delta t_{\text{count}}}roman_Δ italic_t start_POSTSUBSCRIPT count end_POSTSUBSCRIPT and the number of extracted particles N𝑁Nitalic_N is counted in each interval. The standard deviation σ=(NN)2𝜎delimited-⟨⟩superscript𝑁delimited-⟨⟩𝑁2\sigma=\sqrt{\left\langle{(N-\left\langle{N}\right\rangle)^{2}}\right\rangle}italic_σ = square-root start_ARG ⟨ ( italic_N - ⟨ italic_N ⟩ ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⟩ end_ARG and mean μ=N𝜇delimited-⟨⟩𝑁\mu=\left\langle{N}\right\rangleitalic_μ = ⟨ italic_N ⟩, evaluated over a larger time span ΔtevalΔsubscript𝑡eval{\Delta t_{\text{eval}}}roman_Δ italic_t start_POSTSUBSCRIPT eval end_POSTSUBSCRIPT, is then used to determine the coefficient of variation cvsubscript𝑐vc_{\text{v}}italic_c start_POSTSUBSCRIPT v end_POSTSUBSCRIPT or the equivalent spill duty factor F𝐹Fitalic_F [13]:

cv=σμ=N2N21subscript𝑐v𝜎𝜇delimited-⟨⟩superscript𝑁2superscriptdelimited-⟨⟩𝑁21\displaystyle c_{\text{v}}=\frac{\sigma}{\mu}=\sqrt{\frac{\left\langle{N^{2}}% \right\rangle}{\left\langle{N}\right\rangle^{2}}-1}italic_c start_POSTSUBSCRIPT v end_POSTSUBSCRIPT = divide start_ARG italic_σ end_ARG start_ARG italic_μ end_ARG = square-root start_ARG divide start_ARG ⟨ italic_N start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⟩ end_ARG start_ARG ⟨ italic_N ⟩ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG - 1 end_ARG F=N2N2=11+cv2𝐹superscriptdelimited-⟨⟩𝑁2delimited-⟨⟩superscript𝑁211superscriptsubscript𝑐v2\displaystyle F=\frac{\left\langle{N}\right\rangle^{2}}{\left\langle{N^{2}}% \right\rangle}=\frac{1}{1+c_{\text{v}}^{2}}italic_F = divide start_ARG ⟨ italic_N ⟩ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG ⟨ italic_N start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⟩ end_ARG = divide start_ARG 1 end_ARG start_ARG 1 + italic_c start_POSTSUBSCRIPT v end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG

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

Refer to caption
Figure 5: Measured spill quality (left), number of particles extracted (middle) and required excitation strength (right) as a function of central frequency and bandwidth of a band-limited noise excitation signal. The coefficient of variation is calculated for Δtcount=100 µsΔsubscript𝑡counttimes100microsecond{\Delta t_{\text{count}}}=$100\text{\,}\mathrm{\SIUnitSymbolMicro s}$roman_Δ italic_t start_POSTSUBSCRIPT count end_POSTSUBSCRIPT = start_ARG 100 end_ARG start_ARG times end_ARG start_ARG roman_µ roman_s end_ARG and Δteval=1 sΔsubscript𝑡evaltimes1s{\Delta t_{\text{eval}}}=$1\text{\,}\mathrm{s}$roman_Δ italic_t start_POSTSUBSCRIPT eval end_POSTSUBSCRIPT = start_ARG 1 end_ARG start_ARG times end_ARG start_ARG roman_s end_ARG. The coloured contour plots are derived by triangulation from the sampled points (white dots). The red dot and span marks the respective global optimum and associated bandwidth. The vertical lines indicate the third-integer resonance and nearby betatron sidebands.

The effects of excitation frequency and bandwidth on the RF-KO extraction are studied experimentally at the SIS18 of GSI. An U23873+{}^{238}\text{U}{}^{73+}start_FLOATSUPERSCRIPT 238 end_FLOATSUPERSCRIPT U start_FLOATSUPERSCRIPT 73 + end_FLOATSUPERSCRIPT beam with a kinetic energy of 200 MeV/nucleontimes200MeVnucleon200\text{\,}\mathrm{M}\mathrm{e}\mathrm{V}\mathrm{/}\mathrm{n}\mathrm{u}% \mathrm{c}\mathrm{l}\mathrm{e}\mathrm{o}\mathrm{n}start_ARG 200 end_ARG start_ARG times end_ARG start_ARG roman_MeV / roman_nucleon end_ARG (Bρ=7 T m𝐵𝜌times7timesteslameterB\rho=$7\text{\,}\mathrm{T}\text{\,}\mathrm{m}$italic_B italic_ρ = start_ARG 7 end_ARG start_ARG times end_ARG start_ARG start_ARG roman_T end_ARG start_ARG times end_ARG start_ARG roman_m end_ARG end_ARG, frev=785 kHzsubscript𝑓revtimes785kHzf_{\text{rev}}=$785\text{\,}\mathrm{k}\mathrm{H}\mathrm{z}$italic_f start_POSTSUBSCRIPT rev end_POSTSUBSCRIPT = start_ARG 785 end_ARG start_ARG times end_ARG start_ARG roman_kHz end_ARG) 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 106 particles/stimesE6particless{10}^{6}\text{\,}\mathrm{p}\mathrm{a}\mathrm{r}\mathrm{t}\mathrm{i}\mathrm{c}% \mathrm{l}\mathrm{e}\mathrm{s}\mathrm{/}\mathrm{s}start_ARG start_ARG end_ARG start_ARG ⁢ end_ARG start_ARG power start_ARG 10 end_ARG start_ARG 6 end_ARG end_ARG end_ARG start_ARG times end_ARG start_ARG roman_particles / roman_s end_ARG [15]. Secondly, the intensity fluctuations measured by the detector are evaluated to determine the spill quality. Thereby Δtcount=100 µsfex1Δsubscript𝑡counttimes100microsecondmuch-greater-thansuperscriptsubscript𝑓ex1{\Delta t_{\text{count}}}=$100\text{\,}\mathrm{\SIUnitSymbolMicro s}$\gg f_{% \text{ex}}^{-1}roman_Δ italic_t start_POSTSUBSCRIPT count end_POSTSUBSCRIPT = start_ARG 100 end_ARG start_ARG times end_ARG start_ARG roman_µ roman_s end_ARG ≫ italic_f start_POSTSUBSCRIPT ex end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT 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 N2qx/N1qx=2.4(3)/2.7(3)=90 %subscriptdelimited-⟨⟩𝑁2subscript𝑞xsubscriptdelimited-⟨⟩𝑁1subscript𝑞xuncertain2.43uncertain2.73times90percent\left\langle{N}\right\rangle_{2q_{\text{x}}}/\left\langle{N}\right\rangle_{1-q% _{\text{x}}}=$2.4(3)$/$2.7(3)$=$90\text{\,}\mathrm{\char 37\relax}$⟨ italic_N ⟩ start_POSTSUBSCRIPT 2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_POSTSUBSCRIPT / ⟨ italic_N ⟩ start_POSTSUBSCRIPT 1 - italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_POSTSUBSCRIPT = start_ARG 2.4 end_ARG start_ARG ( 3 ) end_ARG / start_ARG 2.7 end_ARG start_ARG ( 3 ) end_ARG = start_ARG 90 end_ARG start_ARG times end_ARG start_ARG % end_ARG. The excitation level required is with Uex2qx/Uex1qx=0.21(4)/0.12(5)=1.7(6)subscriptdelimited-⟨⟩subscript𝑈ex2subscript𝑞xsubscriptdelimited-⟨⟩subscript𝑈ex1subscript𝑞xuncertain0.214uncertain0.125uncertain1.76\left\langle{U_{\text{ex}}}\right\rangle_{2q_{\text{x}}}/\left\langle{U_{\text% {ex}}}\right\rangle_{1-q_{\text{x}}}=$0.21(4)$/$0.12(5)$=$1.7(6)$⟨ italic_U start_POSTSUBSCRIPT ex end_POSTSUBSCRIPT ⟩ start_POSTSUBSCRIPT 2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_POSTSUBSCRIPT / ⟨ italic_U start_POSTSUBSCRIPT ex end_POSTSUBSCRIPT ⟩ start_POSTSUBSCRIPT 1 - italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT end_POSTSUBSCRIPT = start_ARG 0.21 end_ARG start_ARG ( 4 ) end_ARG / start_ARG 0.12 end_ARG start_ARG ( 5 ) end_ARG = start_ARG 1.7 end_ARG start_ARG ( 6 ) end_ARG 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 cv=0.76subscript𝑐v0.76c_{\text{v}}=$0.76$italic_c start_POSTSUBSCRIPT v end_POSTSUBSCRIPT = 0.76 compared to cv=1.21subscript𝑐v1.21c_{\text{v}}=$1.21$italic_c start_POSTSUBSCRIPT v end_POSTSUBSCRIPT = 1.21 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

Refer to caption
Figure 6: Spill quality (top) and excitation power (bottom) as function of excitation frequency for a band-limited noise signal (ΔQex=0.008Δsubscript𝑄ex0.008\Delta Q_{\text{ex}}=$0.008$roman_Δ italic_Q start_POSTSUBSCRIPT ex end_POSTSUBSCRIPT = 0.008) as obtained from simulations. The quality is calculated for Δtcount=100 µsΔsubscript𝑡counttimes100microsecond{\Delta t_{\text{count}}}=$100\text{\,}\mathrm{\SIUnitSymbolMicro s}$roman_Δ italic_t start_POSTSUBSCRIPT count end_POSTSUBSCRIPT = start_ARG 100 end_ARG start_ARG times end_ARG start_ARG roman_µ roman_s end_ARG and Δteval=100 msΔsubscript𝑡evaltimes100millisecond{\Delta t_{\text{eval}}}=$100\text{\,}\mathrm{ms}$roman_Δ italic_t start_POSTSUBSCRIPT eval end_POSTSUBSCRIPT = start_ARG 100 end_ARG start_ARG times end_ARG start_ARG roman_ms end_ARG. Vertical lines mark the resonance and nearby betatron sidebands, and the two excitation bands from figure 3 are indicated.
Refer to caption
Figure 7: Particle distributions before and during excitation near the first (top) and second order sideband (bottom) obtained from a simulation of 105superscript10510^{5}10 start_POSTSUPERSCRIPT 5 end_POSTSUPERSCRIPT particles. The band-limited noise excitation signals correspond to figures 3 and 6 (ΔQex=0.008Δsubscript𝑄ex0.008\Delta Q_{\text{ex}}=$0.008$roman_Δ italic_Q start_POSTSUBSCRIPT ex end_POSTSUBSCRIPT = 0.008). The location of the separatrix at H=Hsep𝐻subscript𝐻sepH=H_{\text{sep}}italic_H = italic_H start_POSTSUBSCRIPT sep end_POSTSUBSCRIPT is indicated by the solid red line.

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 9 dBtimes9dB9\text{\,}\mathrm{d}\mathrm{B}start_ARG 9 end_ARG start_ARG times end_ARG start_ARG roman_dB end_ARG 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 24qx24subscript𝑞x2-4q_{\text{x}}2 - 4 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT, 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 cv1/μproportional-tosubscript𝑐v1𝜇c_{\text{v}}\propto 1/\muitalic_c start_POSTSUBSCRIPT v end_POSTSUBSCRIPT ∝ 1 / italic_μ 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 H𝐻Hitalic_H normalized to its value at the separatrix Hsepsubscript𝐻sepH_{\text{sep}}italic_H start_POSTSUBSCRIPT sep end_POSTSUBSCRIPT is used to analyse the particle distribution during the extraction by means of a histogram (figure 7). Thereby, H=0𝐻0H=0italic_H = 0 corresponds to a particle in the centre of phase space with zero betatron amplitude and (chromatic) tune Qxsubscript𝑄xQ_{\text{x}}italic_Q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT; and for HHsep𝐻subscript𝐻sepH\to H_{\text{sep}}italic_H → italic_H start_POSTSUBSCRIPT sep end_POSTSUBSCRIPT 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 t=1.3 ms𝑡times1.3mst=$1.3\text{\,}\mathrm{m}\mathrm{s}$italic_t = start_ARG 1.3 end_ARG start_ARG times end_ARG start_ARG roman_ms end_ARG. The excitation band near 1qx1subscript𝑞x1-q_{\text{x}}1 - italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT 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 2qx2subscript𝑞x2q_{\text{x}}2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT 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 H=0𝐻0H=0italic_H = 0 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 k=1𝑘1k=1italic_k = 1 sideband feeding the extraction can be combined with a stronger excitation signal placed on the other side of the resonance towards the k=2𝑘2k=2italic_k = 2 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 Au19765+{}^{197}\text{Au}{}^{65+}start_FLOATSUPERSCRIPT 197 end_FLOATSUPERSCRIPT Au start_FLOATSUPERSCRIPT 65 + end_FLOATSUPERSCRIPT beam at 200 MeV/nucleontimes200MeVnucleon200\text{\,}\mathrm{M}\mathrm{e}\mathrm{V}\mathrm{/}\mathrm{n}\mathrm{u}% \mathrm{c}\mathrm{l}\mathrm{e}\mathrm{o}\mathrm{n}start_ARG 200 end_ARG start_ARG times end_ARG start_ARG roman_MeV / roman_nucleon end_ARG. 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 1qx1subscript𝑞x1-q_{\text{x}}1 - italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT band or on the other side of the resonance. For the frequency of the sine component a local optimum exists not only at Qex,10.677=10.323subscript𝑄ex10.67710.323Q_{\text{ex},1}\approx 0.677=1-0.323italic_Q start_POSTSUBSCRIPT ex , 1 end_POSTSUBSCRIPT ≈ 0.677 = 1 - 0.323 towards the first order sideband, but also at 0.648=2×0.3240.64820.3240.648=2\times 0.3240.648 = 2 × 0.324 towards the second order sideband on the other side of the resonance. With 1qx=2/3+|d|1subscript𝑞x23𝑑1-q_{\text{x}}=2/3+\absolutevalue{d}1 - italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT = 2 / 3 + | start_ARG italic_d end_ARG | and 2qx=2/32|d|2subscript𝑞x232𝑑2q_{\text{x}}=2/3-2\absolutevalue{d}2 italic_q start_POSTSUBSCRIPT x end_POSTSUBSCRIPT = 2 / 3 - 2 | start_ARG italic_d end_ARG | 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.

Refer to caption
Figure 8: Measured spill quality as a function of sine and noise frequency for an excitation with Noise ++. For the scan, the noise bandwidth ΔQex=0.01Δsubscript𝑄ex0.01\Delta Q_{\text{ex}}=0.01roman_Δ italic_Q start_POSTSUBSCRIPT ex end_POSTSUBSCRIPT = 0.01 and frequency of the second sine Qex,2=1.6715subscript𝑄ex21.6715Q_{\text{ex},2}=1.6715italic_Q start_POSTSUBSCRIPT ex , 2 end_POSTSUBSCRIPT = 1.6715 are kept constant.

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).