Geant4-based Monte Carlo simulations of a transport beam line for - - PowerPoint PPT Presentation

geant4 based monte carlo simulations of a transport beam
SMART_READER_LITE
LIVE PREVIEW

Geant4-based Monte Carlo simulations of a transport beam line for - - PowerPoint PPT Presentation

MCMA 2017 Napoli, Italy 17 th October 2017 Geant4-based Monte Carlo simulations of a transport beam line for multidisciplinary applications of laser-driven proton beams Giuliana Milluzzo On behalf of the ELIMED collaboration INFN -


slide-1
SLIDE 1
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Geant4-based Monte Carlo simulations of a transport beam line for multidisciplinary applications of laser-driven proton beams

Giuliana Milluzzo

1

On behalf of the ELIMED collaboration

INFN - Laboratori Nazionali del Sud, Catania, Italy Università degli Studi di Catania

MCMA 2017

Napoli, Italy
 17th October 2017

slide-2
SLIDE 2
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Conventional RF accelerators vs laser-driven for hadrontherapy

2

In a few microns mul0-Mev protons and ions are accelerated

  • An intense laser field

(> 1018 W/cm2) blows

  • ff electrons from a

target surface

  • Fast electrons

penetrate the foil and ionize atoms along their path à a strong electric field is created (≈ TV/m)

slide-3
SLIDE 3
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Laser-driven beams peculiarities

3

  • Large energy spread
  • Wide angular distribu2ons
  • Extremely high dose rate per pulse

107 -109 Gy/min (vs 10-50 Gy/min) New paradigms for beam dosimetry

slide-4
SLIDE 4
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

4

Collection and diagnostics Selection, transport and diagnostics Dosimetry and sample irradiation

Graphics by L. Allegra

ELIMED beam line layout

4

  • Front. Oncol. (in press)

NIM A, Vol. 829 (2016) JINST, Vol. 11 (2016) JINST, Vol. 10 (2015)

slide-5
SLIDE 5
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ELIMED WPs

5

slide-6
SLIDE 6
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ELIMED WPs

5

slide-7
SLIDE 7
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

The ELIMED application

6

Source PMQs Energy selector system Conventional transport elements Dosimetry devices

100 20 40 60 80 70 60 50 40 30 20 10

Application structure

Component realistic model Magnetic and electric field implementation Realistic laser-driven particle source Information scoring along the beamline

Requirements from ELI

Easily modify geometrical configurations Accurate transport in magnetic fields User friendly

PIC simulation

  • J. Pipek, F. Romano, G. Milluzzo et al.,, Journal of Instrumentation, Volume 12, March 2017

OUTPUT

slide-8
SLIDE 8
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

The ELIMED application as a tool for beam optimization and feasibility studies for multi-disciplinary applications

slide-9
SLIDE 9
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ΔE/E 11% 8% 7%

  • Tr. Eff.

9.9% 8.0% 4.9%

60 MeV

slide-10
SLIDE 10
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ΔE/E 11% 8% 7%

  • Tr. Eff.

9.9% 8.0% 4.9%

60 MeV

slide-11
SLIDE 11
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

2 m air + 100 um Ta 2 m air + 300 um Ta

+ 10 mm collimator + 20 mm collimator

10 % 10 mm 10 % 20 mm

slide-12
SLIDE 12
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Towards dose delivery of clinical relevance: feasibility studies for multi-disciplinary applications

slide-13
SLIDE 13
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active energy modulator?

11

PMQs configuration fixed for the maximum energy transport

  • ptimization (60 MeV)

10 20 30 40 50 60 70 10

1

10

2

10

3

10

4

Energy [MeV] Counts

slide-14
SLIDE 14
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active energy modulator?

11

PMQs configuration fixed for the maximum energy transport

  • ptimization (60 MeV)

10 20 30 40 50 60 70 10

1

10

2

10

3

10

4

Energy [MeV] Counts

range in water between 1 and 2 cm

slide-15
SLIDE 15
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active energy modulator?

11

PMQs configuration fixed for the maximum energy transport

  • ptimization (60 MeV)

10 20 30 40 50 60 70 10

1

10

2

10

3

10

4

Energy [MeV] Counts

range in water between 1 and 2 cm

Protons in the energy selector

slide-16
SLIDE 16
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active energy modulator?

11

PMQs configuration fixed for the maximum energy transport

  • ptimization (60 MeV)

10 20 30 40 50 60 70 10

1

10

2

10

3

10

4

Energy [MeV] Counts

range in water between 1 and 2 cm

Protons in the energy selector

slide-17
SLIDE 17
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active energy modulator?

11

PMQs configuration fixed for the maximum energy transport

  • ptimization (60 MeV)

10 20 30 40 50 60 70 10

1

10

2

10

3

10

4

Energy [MeV] Counts

range in water between 1 and 2 cm

Protons in the energy selector

slide-18
SLIDE 18
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active modulator: results

12

slide-19
SLIDE 19
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active modulator: results

12

3% of dose homogeneity in 2 cm of modula5on region!

slide-20
SLIDE 20
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active modulator: results

12

Normalization to typical doses per session delivered in ocular melanoma proton treatments (15 Gy) Re-calculation of the absolute weights (in Gy) of any single peak Computation of number of shots per single peak to achieve the required dose

3% of dose homogeneity in 2 cm of modula5on region!

slide-21
SLIDE 21
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ESS as an active modulator: results

12

Normalization to typical doses per session delivered in ocular melanoma proton treatments (15 Gy) Re-calculation of the absolute weights (in Gy) of any single peak Computation of number of shots per single peak to achieve the required dose

̴ 1000 shots in total

Considering both laser repetition rate and and ESS field frequency oh 1 Hz ̴ 16 minutes Considering laser repetition rate of 10 Hz and and ESS field frequency oh 1 Hz ̴ 2 minutes

3% of dose homogeneity in 2 cm of modula5on region!

slide-22
SLIDE 22
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Conclusions and future developments

13

Conclusions

  • ELIMED application has been designed to be a Geant4 User-oriented application
  • Geant4 fits well all the requirements from ELI management:
  • Simulate complex geometrical elements with the possibility to switch on/off components

modularity

  • Easy methods for changing geometrical configurations interactive commands
  • Reliability for particle transport in magnetic fields robust tracking
  • Accurate energy/dose depositions predictions and secondary particles production -> well

tested and validated physics models

  • User-friendly and easy-to-use for non experienced Users -> Qt interface with the

possibility

  • f changing all the key parameters/elements of the beamline

Future developments

  • Voxel phantom import (from DICOM)
  • Calculation of LET distributions (and RBE?) in the perspective of future in-vitro/in-vivo

experiments to be performed at ELI-Beamline

  • Implementation of Geant4-DNA models? spatio-temporal track correlations due to the extremely

high dose rate per pulse

slide-23
SLIDE 23
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Thank you

slide-24
SLIDE 24
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Validation with reference code: magnetic fields& tracking

15

Magnetic field gradient within the 160 mm quad Magnetic field intensity along the radial axis of the dipole

Comparison with COMSOL and SIMION software

slide-25
SLIDE 25
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Validation with reference code: magnetic fields& tracking

15

Magnetic field gradient within the 160 mm quad Magnetic field intensity along the radial axis of the dipole

Comparison with COMSOL and SIMION software

60 MeV p track within ESS system in the x-z plane 60 MeV p track within PMQ system in the x-z plane

slide-26
SLIDE 26
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Monoenerge'c

Adaptive source

16

slide-27
SLIDE 27
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Monoenerge'c

Analytic

Adaptive source

16

slide-28
SLIDE 28
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Monoenerge'c

Analytic

Adaptive source

16

Detailed laser-driven proton and ion beam source to use as input

PIC 3D PIC 2D

slide-29
SLIDE 29
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Monoenerge'c

Analytic

Adaptive source

16

Detailed laser-driven proton and ion beam source to use as input

PIC 3D PIC 2D

slide-30
SLIDE 30
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Monoenerge'c

Analytic

Adaptive source

16

Detailed laser-driven proton and ion beam source to use as input

PIC 3D PIC 2D

slide-31
SLIDE 31
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Weighting low energy conventionally accelerated proton beams for modulators

17

slide-32
SLIDE 32
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Weighting low energy conventionally accelerated proton beams for modulators

17

Tumor thickness

Spread Out Bragg Peak

slide-33
SLIDE 33
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Weighting low energy conventionally accelerated proton beams for modulators

17

Tumor thickness

Spread Out Bragg Peak

  • B. Jia, F. Romano, G.A.P. Cirrone, G. CuEone, M.H. Hadizadeh, A.A. Mowlavi, L.
  • Raffaele. NIM A (Jan 2016)
slide-34
SLIDE 34
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Weighting low energy conventionally accelerated proton beams for modulators

17

Tumor thickness

Spread Out Bragg Peak

  • B. Jia, F. Romano, G.A.P. Cirrone, G. CuEone, M.H. Hadizadeh, A.A. Mowlavi, L.
  • Raffaele. NIM A (Jan 2016)
slide-35
SLIDE 35
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Weighting low energy conventionally accelerated proton beams for modulators

17

Tumor thickness

Spread Out Bragg Peak

  • B. Jia, F. Romano, G.A.P. Cirrone, G. CuEone, M.H. Hadizadeh, A.A. Mowlavi, L.
  • Raffaele. NIM A (Jan 2016)
slide-36
SLIDE 36
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

In vacuum transport- focusing and energy selection

18

PMQs ESS

x [mm] z [mm] z [mm] y [mm] y [mm] x [mm] 35 MeV 45 MeV 60 MeV 30 MeV proton 60 MeV proton

slide-37
SLIDE 37
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

ELIMED Monte Carlo WP

19

ELIMED (LNS-INFN) @ ELI-Beamlines

WP2 (LNS-INFN) Monte Carlo Simulations

Resp.: F. Romano

  • J. Pipek
  • G. Milluzzo

WP1 (LNS-INFN) Ion Beam Transport WP3 (LNS-INFN) Ion Beam Diagnostics and Dosimetry

slide-38
SLIDE 38
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Optimization of transversal profiles for applications

slide-39
SLIDE 39
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Optimization of transversal profiles for applications

slide-40
SLIDE 40
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

…in short

21

Bucke&ul of par.cles! Controlled beams Conven.onal Laser-driven

slide-41
SLIDE 41
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Secondary radiation: prompt radiation

22 Fluence map of secondary prompt photons escaping from the beam line (top view) for 60 MeV p-configuration

p/cm2

Total number of prompt neutrons and photons escaped from the different elements for 60 MeV p-configuration

slide-42
SLIDE 42
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Secondary radiation study

23 Fluence map of secondary prompt photons escaping from the beam line (top view) for 60 MeV p-configuration Distribution of activated nuclei projected into the horizontal plane (60-MeV proton configuration)

slide-43
SLIDE 43
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

Activation

24 The total beta activity (including electron capture) of the beam-line elements as a function of time (60-MeV proton configuration). Distribution of activated nuclei projected into the horizontal plane (60-MeV proton configuration)

slide-44
SLIDE 44
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

User interface

25

Qt interface for visualization and interactive simulation

Macro commands

Remote simulation using network performances Batch mode

slide-45
SLIDE 45
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

The ELIMED application

26

Source PMQs Energy selector system Conventional transport elements Dosimetry devices

100 20 40 60 80 70 60 50 40 30 20 10

“Toolkit for the simulation of the passage of particles through matter.”

Application structure

Component realistic model Magnetic and electric field implementation Realistic laser-driven particle source Information scoring along the beamline

Requirements from ELI

Easily modify geometrical configurations Accurate transport in magnetic fields User friendly

http://www.geant4.org

PIC simulation

  • J. Pipek, F. Romano, G. Milluzzo et al.,, Journal of Instrumentation, Volume 12, March 2017

OUTPUT

slide-46
SLIDE 46
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

The ELIMED application

27

INPUT PMQs Energy selector system Conventional transport elements Dosimetry devices

“Toolkit for the simulation of the passage of particles through matter.”

Application structure Component realistic model Magnetic and electric field implementation Realistic laser-driven particle source Information scoring along the beamline

Requirements from ELI Easily modify geometrical configurations Accurate transport in magnetic fields User friendly

http://www.geant4.org

  • J. Pipek, F. Romano, G. Milluzzo et al.,, Journal of Instrumentation, Volume 12, March 2017

OUTPUT

slide-47
SLIDE 47
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

28

60 MeV

slide-48
SLIDE 48
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

28

ΔE/E 11% 8% 7%

  • Tr. Eff.

9.9% 8.0% 4.9%

60 MeV

slide-49
SLIDE 49
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

29

5 MeV 60 MeV

slide-50
SLIDE 50
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

29

ΔE/E 5.0% 4.3% 6.6%

  • Tr. Eff.

5.9% 4.5% 2.1% ΔE/E 11% 8% 7%

  • Tr. Eff.

9.9% 8.0% 4.9%

5 MeV 60 MeV

slide-51
SLIDE 51
  • G. Milluzzo- INFN-LNS (Italy) - gmilluzzo@lns.infn.it

5 MeV 60 MeV