Proceedings of the 11th Annual Meeting of Particle Accelerator Society of Japan
August 9-11, 2014, Aomori, Japan

PASJ2014-SAP028

二種類の周波数領域の高周波を用いる新たな遅いビーム取り出し法
AN ADVANCED METHOD FOR SLOW EXTRACTION
WITH DUAL BAND RF VOLTAGE
西内秀晶#, A ) , 平本和夫 B ) , えび名風太郎 A )Hideaki Nishiuchi#, A), Kazuo Hiramoto B), Futaro Ebina A)
A)
Hitachi Research Laboratory, Hitachi, Ltd.
B)
Research & Development Group, Hitachi, Ltd.

Abstract

We propose “Dual-Step RF Extraction method” (DS-RF method) as an advanced slow extraction method of the ion
synchrotron with non-linear resonance of the betatron oscillations. The DS-RF method applies dual band RF voltage, andeach RF voltage controls it independently. One of the dual band RF voltages, an extraction RF voltage near the resonancemake the beam exceed the separatrix by augmentation of the betatron oscillation amplitude of a circulating beam andextracts it from a synchrotron. The other RF voltage, a diffusion RF voltage far from the resonance is applied to diffuse
the beam. A circulating beam is diffused by the RF voltage for diffusion in separatrix, and the beam spread in theseparatrix neighborhood is extracted outside a synchrotron with the RF voltage for extraction. The control range of theextraction beam current is enlarged by controlling these RF voltages independently. In this article, we show the principleof the DS-RF method and the characteristic of this extraction method as a result of multi particle beam simulation.
RF-KO 法においては、セパラトリックス近傍の粒子
1. はじめにのベータトロン振動数に合わせた単一周波数の高周高エネルギービームの利用分野、例えば、核物理 波電圧を従来の FM 変調電圧に加えることで、照射研究や粒子線によるがん治療では、イオンシンクロ 中にセパラトリックス近傍の粒子を選択的に取り出[4]トロンからの遅いビーム取り出しが用いられてきた。 す制御法が示されている 。この制御法は、単一周これらのビーム利用では、多くの場合、ビーム電流 波数の高周波電圧を重畳して印加することでセパラの一定化が望まれる他、特に粒子線治療では、目標 トリックス近傍の粒子密度分布を低減し、シンクロを高精度で照射するため、ビーム位置の変化を抑え トロン振動による周回ビームの運動量変化や電磁石た取り出し法が望まれている。遅いビーム取り出し の電源リップルなどにより、高周波電圧を印加しな法では、共鳴現象によるベータトロン振動の増加を くても取り出されてしまうビームの発生を抑制でき利用してビームをシンクロトロンからビーム輸送系 る。
本研究では、ベータトロン振動の非線形共鳴を用へ取り出す。ベータトロン振動の共鳴は、ベータトいるイオンシンクロトロンの新たな遅いビーム取りロン振動の変位がシンクロトロンに設置した多極電磁石を励磁して形成するセパラトリックスを超えた 出し法として、Dual-Step RF Extraction 法（DS-RF 法）時に発生する。このような遅いビーム取り出し法の を提案する。DS-RF 法では、二種類の周波数領域のひとつに、セパラトリックスを一定に保持し、高周 高周波電圧を用い、それぞれ独立に制御する。二種波電圧でセパラトリックス内を周回するビームの 類の周波数領域のうち、共鳴に近い側の周波数を持ベータトロン振動を増大することでセパラトリック つ高周波電圧（取り出し用高周波電圧）により、スを越えさせて共鳴を発生する、高周波ビーム取り ベータトロン振動の振幅を増加させることでセパラ出し法 [1]がある。高周波ビーム取り出し法では、セ トリックスを越えさせ、シンクロトロンからビームパラトリックスを一定にしたままでベータトロン振 を取り出す。この際、共鳴から離れた側の周波数領動に共鳴を発生させるため、取り出しビームの位 域の高周波電圧（拡散用高周波電圧）を合わせて置・勾配を一定にできる特徴がある。また、取り出 ビームに印加する。この拡散用高周波電圧により周しビーム電流は、周回ビームに印加する高周波電圧 回ビームがセパラトリックス内で拡散し、この拡散されたビームが取り出し用高周波でシンクロトロンの振幅値にて容易に制御可能である。
高周波ビーム取り出し法の高周波電圧印加方法と 外に取り出される。このような二種類の周波数領域して、以下の二つが挙げられる。(1) RF-DE 法[2]は、 の高周波電圧各々を制御することで取り出しビーム周回ビームに印加する高周波電圧として帯域制限高 電流の制御範囲を広げるとともに、従来の一種類の周波電圧を採用し、高周波電圧の振幅値制御にて取 高周波電圧で取り出した際に生じる取り出しビームり出しビーム電流を制御する。帯域制限高周波電圧 電流の時間構造が抑制可能となる。本論では、DSは、複数の線スペクトルで構成される。(2) RF-KO RF 法の原理と、多粒子ビームシミュレーション結果
法[3]は、正弦波高周波電圧を FM 変調することで広 より本取り出し法の特徴を示す。
帯域の高周波電圧を周回ビームに印加する。更に、___________________________________________# hideaki.nishiuchi.un@hitachi.com- 524 -

Proceedings of the 11th Annual Meeting of Particle Accelerator Society of Japan
August 9-11, 2014, Aomori, Japan

PASJ2014-SAP028

2.

DS-RF 法の原理

Number of
circulating
particlesに印加する高周波電圧の周波数構造が取り出しビーム電流に重畳されにくくなり、高精度な取り出し
DS-RF 法の原理を従来法と比較しながら説明する。
ビーム電流の制御が実現できる。
図 1 に DS-RF 法の原理を示す。
図 2 に従来法と DS-RF 法でのセパラトリックス内のビーム分布を示す。従来法（図2(a)）ではセパラParticle distribution of RF voltage
トリックス内のビームは、図1(a)に点線で示したよRF voltage for
うに、セパラトリックス近傍まで広く分布している。
diffusion and extraction
一方、DS-RF 法（図 2(b)）では、ビームはセパラトリックス近傍まで広がっておらず、取り出し用高周波電圧を印加しないとビームが取り出されない分布を実現している。
Horizontal tune (νx)νres
Diffusion and extraction velocity ofthe extraction RF voltage(a) Conventional RF beam extraction method
0.01
Particle distribution ofdiffusion RF voltage

X’[a.u.]Number of
circulating
particles

RF voltage for extraction

RF voltage for
diffusion
νres

Horizontal tune (νx)

0-0.01-1.0

Diffusion velocity ofExtraction velocity ofthe extraction RF voltage the diffusion RF voltage

0X [a.u.]

1.0(a) Conventional method(b) DS-RF method

Figure 1: Principle of DS-RF method.

X’[a.u.]従来の高周波ビーム取り出し法 [1]（図 1(a)）では、0.01
高周波電圧の印加範囲によりビームの拡散速度を調整し、振幅値により取り出しビーム電流を調整している。しかし、これらの機能は完全に分離できない。
0
図 1(a)に点線で示したように、ビーム停止直後のシンクロトロン内のビーム分布は、セパラトリックス限界まで広く分布している。取り出しビーム電流の高精度かつ取り出しビーム電流制御範囲を広げる制-0.01
御を実現するには、セパラトリックス内のビーム分
-1.0
01.0
布に応じてビームの拡散速度と高周波電圧の振幅値
X [a.u.]を独立かつ高精度に制御する必要がある。
(b) DS-RF method
DS-RF 法（図 1(b)）では、周回ビームの拡散速度を制御する高周波電圧（拡散用高周波電圧）とビー
Figure 2: Comparison of the beam distribution of
ムを取り出す高周波電圧（取り出し用高周波電圧）
conventional method and DS-RF method.の二つを用いる。拡散用高周波電圧は、ビームが取図 3 に DS-RF 法での高周波電圧の印加周波数を示り出されない範囲で拡散速度を調整する。つまり、拡散用高周波電圧は、図 1(b)に点線で示したように、 す 。図 3 で は、高 周波電 圧の印加 周波数 を水平ビームを拡散するがセパラトリックスを越えさせな チューン νx に換算している。DS-RF 法では、拡散用いよう、高周波電圧の印加範囲と振幅値を調整する。 高周波電圧 Vdif をスペクトル間隔 νdiv の帯域制限高周
そのため、従来の高周波ビーム取り出し法と比較す 波電圧とし、帯域の中心 νcd は周回ビームのベータトると拡散用高周波電圧の印加範囲は狭く、振幅値も ロン振動周波数の重心値、帯域幅 νwd は周回ビーム低い。また、取り出し用高周波電圧は、拡散用高周 のベータトロン振動の周波数分布よりも狭く設定す
波電圧で拡散したビームを安定に取り出せるように、 る。また、取り出し用高周波電圧 Vext は、正弦波発
拡散用高周波電圧よりも共鳴側に印加し、振幅値で 振器から出力される高周波電圧とし、拡散用高周波取り出し速度を制御することでシンクロトロン外に 電圧 Vdif により拡散された周回ビームを効率良く取ビームを取り出す。さらに、取り出し用高周波電圧 り出し可能な周波数 νce を設定する。に単一周波数の高周波電圧を用いることで、ビーム- 525 -Proceedings of the 11th Annual Meeting of Particle Accelerator Society of Japan
August 9-11, 2014, Aomori, JapanPASJ2014-SAP028

Sinusoidal RF voltage
for extraction (Vext)

Band-limit RF voltage
for diffusion (Vdif)

ビームシミュレーションの初期粒子数は、FFT 解析時の統計量を考慮し 50 万点とした。

νdiv

Table 1: Simulation Conditions of the Medical Proton
SynchrotronAextRFamplitude Adifνresνceνcdνwd
Horizontal tune (νx)

Figure 3: Applied band of the RF voltage on the horizontaltune.
図 4 に DS-RF 法での高周波電圧の制御方法を示す。
拡散用高周波電圧 Vdif を構成する帯域制限高周波電圧は、逆フーリエ変換で生成した複数の線スペクトルで構成される[2]。この帯域幅 fwd、線スペクトル周波数間隔 fdiv の帯域制限高周波電圧は、帯域制限高周波電圧の中心周波数 fcd となる正弦波高周波電圧と乗算する。この乗算処理での中心周波数 fcd の調整により、周回ビームのエネルギーに応じた拡散用高周波電圧の印加範囲を容易に切り替えられる。

ADD

Amplitude signal
for extraction
Sinusoidal RF signal
for extraction

更に拡散用高周波電圧 Vdif の振幅値 Adif は、振幅変調回路 AMd で周回ビームのベータトロン振動振幅がセパラトリックスを越えない範囲に制御する。取り出し用高周波電圧 fce は、ビームの取り出し効率が高くなるように調整する。取り出し用高周波電圧
Vext の振幅値 Aext は、振幅変調回路 AMe で取り出しビーム電流を制御する。これら振幅値が制御された
拡 散 用高 周波 電圧 Vdif と 取 り出 し用 高周 波電圧
Vext を加算器で加算し、電力増幅器で増幅した後に高周波電極に印加する。我々は DS-RF 法を多粒子ビーム取り出しシミュレーションで評価し、拡散用高周波電圧と取り出し用高周波電圧の印加条件を整理した。詳細を以下に示す。

3.

ビームシミュレーション

DS-RF 法の原理は、六次元位相空間(x, x’, y, y’, φ,Δp/p)での多粒子トラッキングシミュレーション[5]で評価した。表 1 に示す医療用陽子シンクロトロンの運転条件を示し、表 2 に高周波電圧の印加条件を示す。これらのシミュレーション条件において、DSRF 法での取り出しビームの周波数成分を評価した。
高周波電圧の蹴り角は、周回ビームの 50%程度のビームを 500ms で取り出すことを想定し、拡散用高周波電圧による蹴り角を 0.3μrad、取り出し用高周波電圧による蹴り角を 0.7μrad でそれぞれ固定とした。

Circumference of synchrotron

18m

Extraction tune (νx, νy)(0.68, 1.46)

Chromaticity (ξx, ξy)(-2.0, 2.78)(5πmmmrad,2πmmmrad)

Momentum spread (Δp/p)

0.05%

Frequency of synchrotron oscillation

6.1kHzTable 2: Simulation Conditions of the DS-RF Method

To RF
AmplifierFigure 4. Block diagram of the RF signal generator for DSRF method.

H+ 200MeVEmittance (εx, εy)Amplitude signal
for diffusion
Sinusoidal RF signal
for diffusion
Band-limit
DAC
RF signals

Beam

4.

fdif (convert to νx)

0.679-0.680

fext (convert to νx)

0.677

Adif

0.3μradAext

1.5μrad

fdiv

1kHz

Ninit

500kpts

Extraction time

500ms

DS-RF 法での取り出しビーム特性評価

図 5 にシミュレーション結果を示す。DS-RF 法での取り出し粒子波形は穏やかな変化を示しており（図 5(a)）、時間構造が密である。また、拡散用高周波電圧の線スペクトル間隔に対応する 1ms 周期の時間構造は生じていない（図 5(b)）。
図 6 に取り出し粒子波形の FFT 演算結果を示す。
FFT 演算結果でも、拡散用高周波電圧の線スペクトル周波数間隔に対応した 1kHz とその整数倍の周波数成分は見えず、一方で 20kHz から 30kHz の帯域に1kHz 間隔の周波数成分が僅かに確認できる。この20kHz から 30kHz の周波数成分は、取り出し用高周波電圧の周波数と拡散用高周波電圧の線スペクトル周波数によるビート成分である。このビート成分のスプリアスレベルは DC 成分に対して 42.3dB と十分低く抑えられており、取り出し粒子波形への影響も少ないことが分かる。このように、DS-RF 法では、従来の高周波ビーム取り出し法で生じていた、取り出し用高周波電圧の周波数に起因する低周波数のスプリアス成分が生じない。

- 526 -Proceedings of the 11th Annual Meeting of Particle Accelerator Society of Japan
August 9-11, 2014, Aomori, Japan

0100200300
400
Extraction time[ms]500

Number of Extracted Particles(a) Overview

35302520151050

Extracted particles
Moving average of extacted particles

60626466
68
Extraction time[ms]70(b) Detail view
Figure 5: Extracted particle waveform with DS-RF method.
(a) is overview of extraction waveform and (b) is detailview of extraction waveform from 60ms to 70ms.
FFT Magnitude [dB]806030302025201510500100200300
400
Extraction time [ms](a) Overview

500

Extracted particles
Moving average of extracted particles

251ms

2015105060626466
68
Extraction time [ms]70(b) Detail view

80

0-20-40Extracted particles
Moving average of extracted particles

Figure 7: Extract particle waveform when an extraction RF
signal set in 0.675 by DS-RF method. (a) is overview ofextraction particle waveform and (b) is detail view ofextraction particle waveform from 60ms to 70ms.

42.3dB
Beat frequency of fce and fwd

40

Number of Extracted Particles

した。図 7 に取り出し用高周波電圧を水平チューン換算で 0.6750 に印加した場合の取り出し粒子波形を示し、図 8 に図 7 に示した取り出し粒子波形の FFT演算結果を示す。

Extracted particles
Moving average of extacted particles

Number of Extracted Particles

353025201510500102030
Frequency [kHz]4050Figure 6: FFT operation result of the extraction particle
waveform of figure 5. Frequency division of FFT is 5Hz.

FFT Magnitude [dB]

Number of Extracted Particles

PASJ2014-SAP028

604027.0dB
Beat frequency of fcd and fwd

200-20-40ここで、取り出し粒子の時間構造に拡散用高周波
01020304050電圧の線スペクトル間隔に対応した 1kHz の整数倍
Frequency [kHz]の周波数成分が生じない原因を考察するため、取り出し用高周波電圧の周波数を拡散用高周波電圧の周 Figure 8: FFT operation result of the extraction particle
波数から離した条件で同様にシミュレーションした。 waveform of figure 7.
取り出し用高周波電圧の印加周波数は、拡散用高周取り出し用高周波電圧を拡散用高周波電圧から離波電圧の印加周波数から水平チューン換算で 0.677 したシミュレーション結果において、拡散用高周波（Δν=0.002）から 0.675（Δν=0.004）に離し、図 5 の 電圧の線スペクトル間隔に起因する 1kHz の整数倍シミュレーション結果と同程度の粒子数が取り出せ の周波数成分が確認できた。この周波数成分は、実るよう、取り出し用高周波電圧の蹴り角を 3.9μrad と- 527 -Proceedings of the 11th Annual Meeting of Particle Accelerator Society of Japan
August 9-11, 2014, Aomori, Japan

PASJ2014-SAP028

際にビームに印加している帯域制限高周波電圧の線 電圧で取り出した際に生じる線スペクトル周波数間スペクトル成分同士によるビートで生ずるものであ 隔に起因する 1kHz の整数倍の周波数の時間構造のる。この時のスプリアスレベルは DC 成分と 1kHz 成 抑制できることを確認した。引き続き、取り出し分で 27.0dB であり、取り出し用高周波電圧を 0.6770 ビーム電流制御性について評価を進める。に設定した場合より 14.3dB も強い。
以上のシミュレーション結果に基づき、DS-RF 法 参考文献
において 1kHz とその整数倍のスプリアス成分が抑 [1] K. Hiramoto, et al., Nucl. Instr. and Meth. A 322 (1992) 154160.
制される理由について考える。複数の線スペクトルで構成される拡散用高周波電圧は、それぞれの線ス [2] K. Hiramoto, et al., Proc. of the 16th Particle Accelerator
Conference, p113, Dallas, USA, 1995.
ペクトル周波数同士でビートを生じる。今回のシ[3] K. Noda, et al., Nucl. Instr. and Meth. A 374 (1996) 269-277.
ミュレーション条件では、線スペクトル間隔 fdiv は [4] K. Mizushima, et al., Nucl. Instr. and Meth. A 606 (2009) 3241kHz であるため、1kHz とその整数倍のビートが取329.
り出しビームにスプリアス成分として重畳する。ま [5] F. Ebina, et al., Proc. of the 5th Particle Accelerator SocietyMeeting, p382, Higashihiroshima, Japan, 2008.
た、DS-RF 法では、取り出し用高周波電圧に拡散用高周波電圧を重畳して印加するため、拡散用高周波電圧の周波数成分と取り出し用高周波電圧の周波数成分によりビートを生じる。今回のシミュレーションにおいて、二つの高周波電圧により生じるビートは 20kHz から 30kHz の間に 1kHz 間隔で生じる。
これら二つのビート周波数成分（1kHz 帯および、20kHz 帯）は、拡散用高周波電圧と取り出し用高周波電圧の拡散速度の差によって、取り出しビームに生じるスプリアス成分が変わると考える。図 5、図 6に示したシミュレーション結果では、取り出し用高周波電圧を拡散用高周波電圧に近付ける（Δν=0.002）
ことで、拡散用高周波電圧による拡散速度に対して取り出し用高周波電圧による移動速度が相対的に速くなる。そのため、拡散用高周波電圧と取り出し用高周波電圧によって生じる 20kHz 帯のビート成分が強く現れ、拡散用高周波電圧により生じる 1kHz 帯のビート成分が打ち消される。一方、図 7、図 8 に示したシミュレーション結果では、取り出し用高周波電圧の周波数が拡散用高周波電圧かの周波数から離れる（Δν=0.004）ことで、拡散用高周波電圧による拡散速度に対して取り出し用高周波電圧の移動速度が相対的に遅くなる。これにより、拡散用高周波電圧により生じる 1kHz 帯のビート成分が含まれたままビームが取り出される。
以上のことから、DS-RF 法において、拡散用高周波電圧に対して取り出し用高周波電圧の拡散速度を高めることで、拡散用高周波電圧の線スペクトル周波数成分を抑制が可能となり、取り出しビームに生じる低周波数の時間構造を抑制できる。

5.

まとめ

ベータトロン振動の非線形共鳴を用いるイオンシンクロトロンの新たな遅いビーム取り出し法として、
DS-RF 法を提案する。DS-RF 法では、ベータトロン振動振幅を増大させる高周波電圧で直接ビームを取り出さず、拡散用高周波電圧により周回ビームがセパラトリックス内で拡散し、セパラトリックス近傍に拡散されたビームが取り出し用高周波でシンクロトロン外に取り出す。DS-RF 法を多粒子シミュレーションにて評価した結果、二種類の周波数領域の高周波電圧各々を制御することで取り出しビーム電流制御範囲を広げるとともに、従来の一種類の高周波- 528 -