Errors and optics study of a permanent magnet quadrupole system F. - - PowerPoint PPT Presentation

errors and optics study of a permanent magnet quadrupole
SMART_READER_LITE
LIVE PREVIEW

Errors and optics study of a permanent magnet quadrupole system F. - - PowerPoint PPT Presentation

Errors and optics study of a permanent magnet quadrupole system F. Schillaci IoP-ASCR, ELI-Beamlines Prague, Czech Republic And MEDical application @ ELI-Beamlines INFN-LNS Catania, Italy 25 th International Conference on Magnet Technology,


slide-1
SLIDE 1

Errors and optics study of a permanent magnet quadrupole system

  • F. Schillaci

IoP-ASCR, ELI-Beamlines

Prague, Czech Republic And INFN-LNS Catania, Italy francesco.schillaci@eli-beams.eu MEDical application @ ELI-Beamlines

25th International Conference on Magnet Technology, Amsterdam 27 August – 01 September 2017

slide-2
SLIDE 2

Laser-driven ion beams

2

  • Large proton number: 1010 ÷ 1013
  • Short bunch duration: few psec
  • High Beam Current: kA
  • !Low Emittance!: 5x10-3 π mm mrad

(microscale spot size but...)

  • Wide Angular Aperture: 10 – 20°

(if we are lucky!)

  • High Energy Spread: ΔE/E >> 10%
  • Low shot-to-shot reproducibilty
slide-3
SLIDE 3

Laser-driven ion beams

3

  • Large proton number: 1010 ÷ 1013
  • Short bunch duration: few psec
  • High Beam Current: kA
  • !Low Emittance!: 5x10-3 π mm mrad
  • Wide Angular Aperture: 10 – 20° (if we are lucky)
  • High Energy Spread: ΔE/E >> 10%
  • Low shot-to-sho reproducibilty
  • High dose-rate per bunch: ~109 Gy/sec

PIC simulations by J. Psikal Expected @ ELI Beamlines

slide-4
SLIDE 4

Laser-driven hadrontherapy

4

slide-5
SLIDE 5

Laser-driven hadrontherapy

5

slide-6
SLIDE 6

ELIMAIA & ELIMED

6

slide-7
SLIDE 7

ELIMAIA & ELIMED

7 Beam line elements: 1) Collection system 2) Selection system 3) Standard transport elements (quadrupoles and steerers) 4) in air dosimetry and irradiation Beam line features: 1) Tunability (deliver ion beams from 5 up to 60 MeV/u) with a controllable energy spread (5% up to 20%) and 106- 1011 ions/pulse 2) Large acceptance 3) Flexibility to meet different User requirements

slide-8
SLIDE 8

8

Magnets for laser-driven particles

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)
  • 20 mm long dipole
  • 50 mm gap
  • C-shape
  • NdFeBo magnets + iron yoke
slide-9
SLIDE 9

9

Magnets for laser-driven particles

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)
  • 20 mm long dipole
  • 50 mm gap
  • C-shape
  • NdFeBo magnets + iron yoke
  • Electron spectrometer!

The general idea of laser-people is: “I need X Telsa, just put a random magnet there and it will work“

slide-10
SLIDE 10

10

Magnets for laser-driven particles

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)
  • 20 mm long dipole
  • 50 mm gap
  • C-shape
  • NdFeBo magnets + iron yoke

Sample: 2-10 Mev

40 mrad uniform divergence

MeV

shot #81

N 6.4 bar N 3.2 bar

shot #317

> 40 MeV electrons N 6.4 bar

shot #225

Laser axis Sample: 10-40 Mev 2 mrad uniform divergence 40 mrad pointing down

slide-11
SLIDE 11

11

Magnets for laser-driven particles

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)
  • 20 mm long dipole
  • 50 mm gap
  • C-shape
  • NdFeBo magnets + iron yoke

Sample: 2-10 Mev

40 mrad uniform divergence

MeV

shot #81

N 6.4 bar N 3.2 bar

shot #317

> 40 MeV electrons N 6.4 bar

shot #225

Laser axis Sample: 10-40 Mev 2 mrad uniform divergence 40 mrad pointing down

slide-12
SLIDE 12

12

Magnets for laser-driven particles

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)
  • 20 mm long dipole
  • 50 mm gap
  • C-shape
  • NdFeBo magnets + iron yoke

Electron spectrometer! Radia Field uniformity ~30%!!!

slide-13
SLIDE 13

13

Magnets for laser-driven particles

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)
  • 20 mm long dipole
  • 50 mm gap
  • C-shape
  • NdFeBo magnets + iron yoke

Electron spectrometer! Radia Field uniformity ~30%!!!

slide-14
SLIDE 14

14

Permanent Magnet prototype test results @ LOA (Fr)

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)
slide-15
SLIDE 15

OUTLINE

15

  • Quadrupole features
  • Error source in magnets and modelling
  • Fixing the tolerances
  • Beam transport (simulations and experiment)
slide-16
SLIDE 16

OUTLINE

16

  • Quadrupole features
  • Error source in magnets and modelling
  • Fixing the tolerances
  • Beam transport (simulations and experiment)
slide-17
SLIDE 17

Quadrupole layout

17

4 PMQs features

(simulations)

  • 2 elements 40 mm long
  • 2 elements 80 mm long
  • 22 mm bore – 20 mm clearance
  • 100T/m field gradient
  • NdFeBo N50 permanent magnets
  • Gradient homogeneity: -6% @ R = 8mm
  • Integrated gradient homogeneity:
  • 1% @ R = 8mm
  • Harmonic content Bn/B2 < 2%
  • Cost-effective prototype

Iron

slide-18
SLIDE 18

Quadrupole layout

18

Iron

4 PMQs features

(simulations)

  • 2 elements 40 mm long
  • 2 elements 80 mm long
  • 22 mm bore – 20 mm clearance
  • 100T/m field gradient
  • NdFeBo N50 permanent magnets
  • Gradient homogeneity: -6% @ R = 8mm
  • Integrated gradient homogeneity:
  • 1% @ R = 8mm
  • Harmonic content Bn/B2 < 2%
  • Cost-effective prototype
slide-19
SLIDE 19

Quadrupole layout

19

Iron

Magnetic design and manufacturing Mechanics designed and manufactured at INFN

slide-20
SLIDE 20

2D Harmonic analysis

20

2D simulations:

  • r0 = 8 mm radius reference circle for B-field

post-processing and harmonic analysis

slide-21
SLIDE 21

2D Harmonic analysis

21

2D simulations:

  • r0 = 8 mm radius reference circle for B-field

post-processing and harmonic analysis

  • Modulus of induction |B| should be

constant

slide-22
SLIDE 22

2D Harmonic analysis

22

2D simulations:

  • r0 = 8 mm radius reference circle for B-field

post-processing and harmonic analysis

  • Modulus of induction |B| should be

constant

  • Radial component Brad = Bx (x/r0) + By (y/r0)

should be purely sinusolidal

slide-23
SLIDE 23

2D Harmonic analysis

23

2D simulations:

  • r0 = 8 mm radius reference circle for B-field

post-processing and harmonic analysis

  • Modulus of induction |B| should be

constant

  • Radial component Brad = Bx (x/r0) + By (y/r0)

should be purely sinusolidal

  • Fourier expansion of Brad gives the

magnitude of the harmonic components Cn:

  • Deviations from ideal behaviour affect the

field quality and the beam transport can show filamentation, emittance growth, steering

Cn= 1 N

N−1 k=1

Brad k r0 exp(ik(2π n N))

slide-24
SLIDE 24

OUTLINE

24

  • Quadrupole features
  • Error source in magnets and modelling
  • Fixing the tolerances
  • Beam transport (simulations and experiment)
slide-25
SLIDE 25

Error source in a magnet

25

  • Magnetization of permanent magnets (remanence, magnetization angle, ...)
  • Manufacturing errors (assembly, pole shimming, ...)
  • Alignment (skew components)
  • Eddy currents (see my talk Status and realization of an high efficiency transport beamline for laser-driven ion beamline [Wed-Mo-Or19])
  • ...

If one or more error sources are introduced symmetry is broken! In order to minimize the errors the tolerances have to be stated for each possible error source. The tighter are the tolerances the higher will be the cost!

slide-26
SLIDE 26

Error source in a magnet

26

  • Magnetization of permanent magnets (remanence, magnetization angle, ...)
  • Manufacturing errors (assembly, pole shimming, ...)
  • Alignment (skew components)
  • Eddy currents (see my talk Status and realization of an high efficiency transport beamline for laser-driven ion beamline [Wed-Mo-Or19])
  • ...

If one or more error sources are introduced symmetry is broken! In order to minimize the errors the tolerances have to be stated for each possible error source. The tighter are the tolerances the higher will be the cost! The goal here is to have no more than 3% of total harmonic component

slide-27
SLIDE 27

2D Errors modelling

27

How to introduce errors in simulations:

Remanence: The remanence Mr of each rare-earth piece is

multiplied by a random number, rand1, with a fixed seed depending on the block identification number and on the ordinal number of themagnetic configuration produced (401 in total). rand1 is uniformly distributed around the mean value 1 with a range of ±0.03 and ±0.06, making the remanent magnetization increasing or decreasing up to 3% and 6%.

Assembly:

slide-28
SLIDE 28

2D Errors modelling

28

How to introduce errors in simulations:

Remanence: The remanence Mr of each rare-earth piece is

multiplied by a random number, rand1, with a fixed seed depending on the block identification number and on the ordinal number of themagnetic configuration produced (401 in total). rand1 is uniformly distributed around the mean value 1 with a range of ±0.03 and ±0.06, making the remanent magnetization increasing or decreasing up to 3% and 6%.

Assembly: The mechanical assembly errors is simulated introducing

a different displacement for each block controlled by a random number rand2 with fixed seed. The direction has been forced to avoid overlapping of the magnets (iron parts are considered fixed). The T-like pieces between two poles are treated as three independent blocks, even if they will be realized as a single one; this allow to take in account not only errors due to the assembly but also errors due to the machining of these parts. rand2 is been defined as uniformly distributed around the mean value 0 with a rangeof ±0.1 and ±0.2. In this way each block is shifted from the ideal position up to 100mm in the first case and up to 200mm in the second case.

slide-29
SLIDE 29

2D Errors modelling

29

How to introduce errors in simulations:

Remanence: The remanence Mr of each rare-earth piece is

multiplied by a random number, rand1, with a fixed seed depending on the block identification number and on the ordinal number of themagnetic configuration produced (401 in total). rand1 is uniformly distributed around the mean value 1 with a range of ±0.03 and ±0.06, making the remanent magnetization increasing or decreasing up to 3% and 6%.

Assembly: The mechanical assembly errors is simulated introducing

a different displacement for each block controlled by a random number rand2 with fixed seed. The direction has been forced to avoid overlapping of the magnets (iron parts are considered fixed). The T-like pieces between two poles are treated as three independent blocks, even if they will be realized as a single one; this allow to take in account not only errors due to the assembly but also errors due to the machining of these parts. rand2 is been defined as uniformly distributed around the mean value 0 with a rangeof ±0.1 and ±0.2. In this way each block is shifted from the ideal position up to 100μm in the first case and up to 200μm in the second case.

slide-30
SLIDE 30

Model validation I Ideal B1 = 0.092 units

30

1 mm shift Effect of the introduced dipole component

R A1 B1 Phase [°] 1 1.572 1.572 45 2 1.571 1.571 45 3 1.571 1.571 45 4 1.570 1.570 45 5 1.569 1.569 45 6 1.568 1.568 45 7 1.566 1.566 45 8 1.564 1.564 45 9 1.559 1.559 45

The radial displacement of the pole at 45° produces a small decrease in the peak of Brad at the same angle. The loss of symmetry produces a dipole contribution in the opposite direction of the pole shift. The real and imaginary parts of the coefficient C1 are equal to each other even if the field is analysed at different reference radii, which means that the phase of the dipole component is θ = arctan(B1)/(A1) = 45°, namely in the direction of the displaced pole.

slide-31
SLIDE 31

Model validation II Ideal B1 = 0.092 units

31

R A1 B1 Phase [°] 1 1.572 1.572 45 2 1.571 1.571 45 3 1.571 1.571 45 4 1.570 1.570 45 5 1.569 1.569 45 6 1.568 1.568 45 7 1.566 1.566 45 8 1.564 1.564 45 9 1.559 1.559 45

If the pole is in its ideal position but its remanence is increased by a factor of two there is a strong increase in the peak of Brad as the loss of symmetry produces a dipole contribution in the same direction of the pole magnetization direction.

2xMr Effect of the introduced dipole component

R A1 B1 Phase [°] 1

  • 12,644
  • 12,644

45 2

  • 12,644
  • 12,644

45 3

  • 12,643
  • 12,643

45 4

  • 12,643
  • 12,643

45 5

  • 12,642
  • 12,642

45 6

  • 12,642
  • 12,642

45 7

  • 12,639
  • 12,639

45 8

  • 12,637
  • 12,637

45 9

  • 12,634
  • 12,634

45

slide-32
SLIDE 32

OUTLINE

32

  • Quadrupole features
  • Error source in magnets and modelling
  • Fixing the tolerances
  • Beam transport (simulations and experiment)
slide-33
SLIDE 33

Random Errors

33

400 different simulations per range of variation of Mr and magnet position

slide-34
SLIDE 34

Random Errors

34

400 different simulations per range of variation of Mr and magnet position Ideal case results

slide-35
SLIDE 35

Random Errors

35

400 different simulations per range of variation of Mr and magnet position The normal content (Bn) does not increase significantly with the increasing of the errors The complex harmonics (Cn) are strongly affected by the errors and their contribution is about 3% of the main harmonic if the errors range in the wider interval.

slide-36
SLIDE 36

Random Errors

36

Combining errors on Mr and magnet position

each magnetic configuration is reproduced on all the different geometric configurations (400 x 400 simulations)

slide-37
SLIDE 37

Random Errors

37

Effects on the field quality Magnetic center nominal position Magnetic center shift

slide-38
SLIDE 38

Random Errors

38

Effects on the field quality

Magnetic measurement

slide-39
SLIDE 39

Random Errors

39

Effects on the field quality

Magnetic measurement

slide-40
SLIDE 40

OUTLINE

40

  • Quadrupole features
  • Error source in magnets and modelling
  • Fixing the tolerances
  • Beam transport (simulations and experiment)
slide-41
SLIDE 41

Beam Transport Simulations

41

slide-42
SLIDE 42

Beam Transport Simulations

42

Ideal case Perturbed case 1.4% Emittance growth

slide-43
SLIDE 43

Beam Transport Simulations

43

Ideal case Perturbed case 1.4% Emittance growth 0.5° bigger angular aperture and more filamentations

slide-44
SLIDE 44

Beam Transport Simulations

44

Ideal case Perturbed case 1.4% Emittance growth 0.5° bigger ang aperture and more filamentations Not negligible steering effect on the radial plane (as expected)

slide-45
SLIDE 45

Beam Transport Test

45

Wide big spot size for cell irradiation

slide-46
SLIDE 46

Beam Transport Test

@ LOA (Fr)

46

3.5 MeV 6.5 MeV

slide-47
SLIDE 47

47

3.5 MeV 6.5 MeV The centroid is shifted and the beam is rotated B1 + skew + missalignment

Beam Transport Test

@ LOA (Fr)

slide-48
SLIDE 48

48

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)

Conclusion

  • A model to study random errors in PMQs is proposed
  • Validated in simple cases
  • Effects of the harmonic contents on beam dynamics results in agreement with the

dipole component produced by the loss of symmetry due to the introduction of imperfection on magnets

  • The method results to be robust and reliable
  • This model is useful to state tolerances on magnet assembly
  • The model is completely general and can include any kind of error source... if you

have enough time to run and analyse thousands of simulations

slide-49
SLIDE 49

49

  • F. Schillaci et al., JINST 10 T05001 (2015)
  • F. Schillaci et al., JINST 11 T07005 (2016)

Thank you for your attention

http://www.eli-beams.eu/ INFN: F. Schillaci, M. Maggiore, G. A. P. Cirrone, G. Cuttone Thanks to W. Beeckman (SigmaPhi) for advice and discussion

slide-50
SLIDE 50

Quadrupole layout

50

4 PMQs features

(simulations)

  • 2 elements 40 mm long
  • 2 elements 80 mm long
  • 22 mm bore – 20 mm clearance
  • 100T/m field gradient
  • NdFeBo N50 permanent magnets
  • Gradient homogeneity: -6% @ R = 8mm
  • Integrated gradient homogeneity:
  • 1% @ R = 8mm
  • Harmonic content Bn/B2 < 2%
  • Cost-effective prototype

Iron

slide-51
SLIDE 51

Random Errors

51

400 different simulations per range of variation of Mr and magnet position The normal content (Bn) does not increase significantly with the increasing of the errors The complex harmonics (Cn) are strongly affected by the errors and their contribution is about 3% of the main harmonic if the errors range in the wider interval