beam dynamics and linac simulation petr ostroumov fermilab accelerator advisory committee may 10 th...
TRANSCRIPT
![Page 1: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/1.jpg)
Beam Dynamics and Linac Simulation
Petr Ostroumov
Fermilab Accelerator Advisory Committee
May 10th – 12th , 2006
![Page 2: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/2.jpg)
2Fermilab
Outline
• Main specifications for the Linac• Basic concepts for the Linac design
– RFQ– MEBT– Front end, 325 MHz– High energy section, 1300 MHz
• Choice of lattice parameters• High-intensity beam physics• Detailed design and simulations• Issues to be solved near future• Conclusion
![Page 3: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/3.jpg)
3Fermilab
Main Linac Specifications
• Provide 8-GeV 1.561014 protons per cycle in the MI
• Beam time structure– Extraction kicker -0.7 msec
– Fit into MI 52.8 MHz rf structure without losses
• Repetition rate & pulse length– Initial configuration: 2.5 Hz at 3 msec, 0.5 MW at 8 GeV
– Ultimate configuration: 10 Hz, 1 msec, 2 MW at 8 GeV
• Consequences– Peak current for beam dynamics design is 40 mA
– Average current over the pulse is 25 mA
– Fast chopper in the MEBT (rise/fall time 2 nsec)
– Debuncher upstream of the MI
![Page 4: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/4.jpg)
4Fermilab
8-GeV Linac conceptual design
• RF power fan-out from one klystron to multiple cavity results in application of SC technology for the whole linac but an initial 10 MeV section
• Use two-frequency Linac option to produce multi-GeV hadron beams:– Apply 1300 MHz ILC cavities above ~1.2 GeV– Develop and use S-ILC cavities (beta=0.81) in the energy range ~400
MeV-1.2 GeV – Spoke loaded SC cavities operating at ILC sub-harmonic frequency in
the front end• Select sub-harmonic frequency for the front end: 1/4.
Motivation: spoke loaded SC cavities are developed at ~345-350 MHz. Requires 30% less number of cavities compared to 433 MHz option. Klystrons are available from JHF developments.
• Below 10 MeV: use the RFQ and 16 RT-CH.
![Page 5: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/5.jpg)
5Fermilab
Linac conceptual design (cont’d)
• 325 MHz SSR-1, SSR-2 and TSR from 10 MeV to ~418 MeV
• Apply SC solenoid focusing to obtain compact lattice in the front end including MEBT
• RFQ delivers axial-symmetric 2.5 MeV H-minus beam • MEBT consists of 2 re-bunchers and a chopper. Smooth
axial-symmetric focusing mitigates beam halo formation• Beam matching between the cryostats: adjust parameters
of outermost elements (solenoid fields, rf phase)• In the frequency transition at ~418 MeV, matching in (, W)-plane is provided by 90 “bunch rotation” • Avoid beam losses due to halo formation, machine errors
and H-minus stripping
![Page 6: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/6.jpg)
6Fermilab
Linac Structure
325 MHz 1300 MHz 1300 MHz
Major Linac Sections
Front end Squeezed ILC-style ILS-style
0.065 2.5 10 33 110 4100.065 2.5 10 33 110 4100.05 2.5 10 32 123 418
Will be installed in the Meson Lab
SSR-2
![Page 7: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/7.jpg)
7Fermilab
Radio Frequency Quadrupole
• Well established accelerator (SNS, J-PARC,….)
• Basic PD requirements:– Cost-effective
– Produce axially-symmetric beam
– Small longitudinal emittance
Average radius R0, cm 0.340
Inter-vane voltage U0, kV 90.45
Vane length, cm 302.428
Peak surface field, kV/cm 330
Output energy, MeV/u 2.498
Transverse emittance, rms, in/out, mm mrad 0.10/0.10
Transverse emittance, 99.5%, in/out, mm mrad 0.14/0.17
Long. emittance, rms, keV/u deg 133 Long. emittance, 99.5%, keV/u deg 1870
Transmission efficiency, % 97.8 Acceleration efficiency, % 95.9
0 50 100 150 200 250
1.0
1.2
1.4
1.6
1.8
2.0
-90
-80
-70
-60
-50
-40
-30
Cell number
Mod
ulat
ion
Syn
. ph
ase
(deg
)
0 50 100 150 200 250
1.0
1.2
1.4
1.6
1.8
2.0
-90
-80
-70
-60
-50
-40
-30
Cell number
Mod
ulat
ion
Syn
. ph
ase
(deg
)
- V.N. Aseev (ANL-PHY)
- A.A. Kolomiets (ITEP, Moscow)
![Page 8: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/8.jpg)
8Fermilab
RFQ Beam Parameters (2.5 MeV, 43 mA)
Phase (deg)
Phase (deg)
X (cm)
X (cm)
Y (cm)
Y (cm)
dX/d
Z(m
rad)
dY/d
Z(m
rad)
dW/W
(%
)C
ount
s
Cou
nts
Cou
nts
Phase (deg)
Phase (deg)
X (cm)
X (cm)
Y (cm)
Y (cm)
dX/d
Z(m
rad)
dY/d
Z(m
rad)
dW/W
(%
)C
ount
s
Cou
nts
Cou
nts
Emittance ( keV/u-nsec)Emittance ( keV/u-nsec)
-W/W XX YY -W/W
![Page 9: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/9.jpg)
9Fermilab
MEBT
S=solenoid; B=buncher; C=cavity
19 25.5 26.5 96 25.5 54.833
Distance between the center of elements (cm)
External surface of the RFQ end wall Space for the chopper First accelerating CH-type cavity
This solenoid belongs to the
first accelerating period
S B CBS S S
Space for Beam Diagnostics Box (BPM, Profile, Toroid, Beam Stop,
Steer, Bunch Length)
S=solenoid; B=buncher; C=cavity
19 25.5 26.5 96 25.5 54.833
Distance between the center of elements (cm)
External surface of the RFQ end wall Space for the chopper First accelerating CH-type cavity
This solenoid belongs to the
first accelerating period
S B CBS S S
Space for Beam Diagnostics Box (BPM, Profile, Toroid, Beam Stop,
Steer, Bunch Length)
![Page 10: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/10.jpg)
10Fermilab
Chopper
-1.5
-1
-0.5
0
0.5
1
1.5
0 0.2 0.4 0.6 0.8 1
Length (m)
Ap
ertu
re h
alf-
size
(cm
)
Chopper plate
Pulser voltage ± 1.9 kVRep. rate 53 MHzRise/fall time 2 nsec (at 10% of the voltage level)Beam target power: 37 kW pulsed, 370 W average
11 msec 0.7 msec
1/52.8 sec
![Page 11: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/11.jpg)
11Fermilab
High intensity beam physics
• Phase advances of transverse oscillations for zero current beams must be below 90
• Wave numbers of oscillations must change adiabatically along the linac despite of many lattice transitions with different types of focusing and inter-cryostat spaces, cavity TTF.
• Avoid strong space charge resonances (Hoffman’s Chart) • Provide equipartitioning of betatron and synchrotron
oscillation temperatures along the linac, primarily in the front end
• Beam matching in the lattice transitions is very important to avoid emittance growth and beam halo formation
• Short focusing periods in the Front End• Analyze HOM and avoid excessive power losses on cavity
walls
![Page 12: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/12.jpg)
12Fermilab
Properties of an ion SC linac
• The acceleration is provided with several types of cavities designed for fixed beam velocity. For the same SC cavity voltage performance there is a significant variation of real-estate accelerating gradient as a function of the beam velocity.
• The length of the focusing period for a given type of cavity is fixed.
• There is a sharp change in the focusing period length in the transitions between the linac sections with different types of cavities
• The cavities and focusing elements are combined into relatively long cryostats with an inevitable drift space between them. There are several focusing periods within a cryostat.
![Page 13: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/13.jpg)
13Fermilab
Iterative procedure of the lattice design
• Select the type and geometric beta of the cavities using a simplified formula for the cavity TTF. Optimize the electrodynamics and the mechanical design of the cavities. By numerical simulation, design the cavities to reduce the ratio of peak surface fields to the accelerating field.
• Assume experimentally proven peak surface fields in SC cavities.• Select the focusing lattice. Select the cryostat length and inter-
cryostat spaces working with cryogenic and mechanical engineers.• Develop lattice tuning for the beam without space charge. • Using rms envelope equations check the lattice tune to verify and
avoid strong space charge resonances.• Provide matching of the beam for the design peak current in all
lattice transitions. • Simulate beam dynamics using multi-particle codes. Study beam
losses using a large number of multi-particles, ~106. • Iterate this procedure to obtain a linac design which satisfies the
engineering requirements and provides high quality beams.
![Page 14: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/14.jpg)
14Fermilab
Accelerating cavities
![Page 15: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/15.jpg)
15Fermilab
Cavity parameters and focusing lattice
CH
SSR-1
SSR-2
TSR
S-ILC
ILC-1
ILC-2
Section CH SSR-1 SSR-2 TSR S-ILC ILC-1 ILC-2b
G - 0.2 0.4 0.6 0.83
# of res. 16 18 33 42 56 63 224# of cryost. - 2 3 7 7 9 28Epeak (MV/m) - 30 28 30 52
Focusing SR SR SRR FRDR FR2DR2 * FR4DR3 FR8DR8
LFocsuing, m 0.515-0.75 0.75 1.6 3.81 6.1 12.2 24.4
1
52
![Page 16: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/16.jpg)
16Fermilab
Voltage gain per cavity
0
5
10
15
20
25
30
0 100 200 300 400 500
Cavity number
Vo
lta
ge
(M
V)
![Page 17: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/17.jpg)
17Fermilab
Stability diagram (betatron oscillations)
2
3 20
sin1
2 ( )
f m ss
S eE
m c
b
-1
-0.8
-0.6
-0.4
-0.2
0
0.2
0.4
0.6
0.8
1
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8
Defocusing factor
cos(
s T)
PRSPhCHSSR-1SSR-2TSRS-ILCILC
![Page 18: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/18.jpg)
18Fermilab
Wave numbers of T- and L- oscillations
0
0.5
1
1.5
2
2.5
3
3.5
0 10 20 30 40 50 60 70 80 90 100 110
Period number
Wav
e-n
umbe
r (1
/m)
SSR2-TSR
SSR1-SSR2
kT
kL
![Page 19: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/19.jpg)
19Fermilab
Linac parameters variation
0
30
60
90
120
0 20 40 60 80 100 120
Period number
Ph
ase
ad
van
ce (
de
g)
sT
sL
-80
-70
-60
-50
-40
-30
-20
-10
0
0 100 200 300 400 500
Cavity number
Sy
nc
hro
no
us
ph
ase
(d
eg
)
Phase advance Synchronous phase
![Page 20: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/20.jpg)
20Fermilab
Peak RF power per cavity
0
100
200
300
400
500
600
700
1 51 101 151 201 251 301 351 401 451
Cavity Number
Ca
vit
y P
ow
er
(kW
)
0
10
20
30
40
50
60
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Cavity number
RF
po
we
r (k
W)
BeamCopper
Due to high shunt impedance, the RT-CH cavities dissipate less rf power than a DTL cavity by a factor of 2
![Page 21: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/21.jpg)
21Fermilab
End-to-end simulations
• End-to-end simulations: the TRACK code. In the RFQ - 108, LINAC- 106 macro-particles
• All fields are 3-D, resonator’s fields - MWS, solenoid fields – EMS
• Lattice is tuned for 45 mA - (RFQ), 43.25 mA - linac
• Some earlier designs have been simulated including machine errors, 100 seeds, 40K particles in each seed
• Linac simulations cross-check by using several codes:– TRACK, ANL (main workhorse)
– ASTRA, DESY code, J.-P. Carneiro (FNAL-AD)
– IMPACT, LBNL/LANL code, B. Mustapha (ANL-Physics)
![Page 22: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/22.jpg)
22Fermilab
Beam envelopes, 43 mA
0 100 200 300 400 500 600 700-4.0
-3.5
-3.0
-2.5
-2.0
-1.5
-1.0
-0.5
0.0
0.5
1.0
1.5
2.0
2.5
3.0
3.5
4.0
Be
am
siz
e (
cm)
Distance (m)
x_rms[cm] y_rms[cm] Xmax[cm] Ymax[cm]
0 100 200 300 400 500 600 700
1
10
100
Pha
se w
idth
(de
g)
Distance (m)
rms maximum
325 MHz 1300 MHz 325 MHz 1300 MHz
![Page 23: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/23.jpg)
23Fermilab
Beam losses (MEBT and RT)
0 5 10 15 20 25 30 35 401E-7
1E-6
1E-5
1E-4
1E-3
0.01
Re
lativ
e b
ea
m lo
sse
s
Distance (m)
0 2000 4000 6000 80000.1
1
10
100
1000
10000
100000
Co
unts
Beam energy (keV)
~1.8% are lost in the MEBT, particle energy <1.4 MeV 8 particles out of million are lost in the RT section, particle energy 4-10 MeV
![Page 24: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/24.jpg)
24Fermilab
RMS emittance growth (43 mA)
0 100 200 300 400 500 600 7000.8
1.0
1.2
1.4
1.6
1.8
2.0R
MS
em
ittan
ce g
row
th fa
ctor
Distance (m)
X Y Z
![Page 25: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/25.jpg)
25Fermilab
Phase space plots (8-GeV, 43 mA)
Ene
rgy
spre
ad,
MeV
X ,
mra
d
Y ,
mra
d
X, cm Y, cm Bunch width, nsec
Image of 1 million particles at the end of Linac without errors
Total emittance
RMS emittance= 50 95 85
![Page 26: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/26.jpg)
26Fermilab
High statistics simulations with rf errors
Option I is more tolerant to jitter errors Option II: kind of beam halo formation for (1.0%, 1 deg),
may require more careful design optimization.
Image of 40 million particles at the end of Linac with rf errors
![Page 27: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/27.jpg)
27Fermilab
Energy Jitter Correction (1.0%, 1.0)
It is more efficient to place the debuncher at the end of the 972.5 m drift between the Linac and the Main Injector.
The required voltage is:
27.5 MV same as for
(0.5%, 0.5 deg)
After correction:
the energy width is
+/- 3 MeV , the phase width is +/- 65 deg of 1300 MHz (130 deg total = 0.28 ns).
![Page 28: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/28.jpg)
28Fermilab
Recent ILC Lattice
![Page 29: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/29.jpg)
29Fermilab
250 300 350 400 450 500 550 600 650 7000.0
0.5
1.0
1.5
2.0
2.5
3.0
Be
am
S
ize
(cm
)
Distance (m)
x_rms[cm] y_rms[cm] Xmax[cm] Ymax[cm]
Beam dynamics in the new ILC lattice
3.2 GeV, 8 ILC RF Units
Simulation of 20K particles, 43.25 mA – NO beam lossesSimulation of 1M part., 43.25 mA – 210-5 relative losses (~13 Watts local beam losses)
![Page 30: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/30.jpg)
30Fermilab
Future work• Integrated lattice design should be continued: beam physics;
mechanical, cryogenic design; implement ILC lattice directly for H-minus acceleration, work with current version of the ILC cryostat design
• Finalize beam diagnostics specifications– Develop beam tuning & commissioning procedures
• High-statistic machine error studies on parallel computer – Beam correction
• Study of different options of the linac to provide the most cost-effective design – Example: frequency transition energy (110 MeV vs 420 MeV)
• Detailed studies of HOM in ILC-style cavities and in the TSRs• Code developments
– Implement FVM feedback model (SCREAM-1D code) into the TRACK code (3D, parallel computing)
– fitting in realistic fields with space charge– Include all H-stripping mechanisms
![Page 31: Beam Dynamics and Linac Simulation Petr Ostroumov Fermilab Accelerator Advisory Committee May 10 th – 12 th, 2006](https://reader035.vdocument.in/reader035/viewer/2022062500/5697bff91a28abf838cbf95a/html5/thumbnails/31.jpg)
31Fermilab
CONCLUSIONS
• New approach in hadron Linacs - “Pulsed SC Front End”-provides high-quality beams
• The concept of “current-independent” tune works well for the SC Linac: the same “43-mA tune” is good for all beam currents in the range from 0 to 43 mA
• Baseline design of the 8-GeV Linac: no beam losses (except H-minus stripping) at present stage of the simulations
• Preliminary study shows that 5 modules (15 cryostats, 120 cavities) of the new ILC RF unit can be used in the high energy end of the linac