WO2013175600A1 - 荷電粒子ビーム輸送系及び粒子線治療装置 - Google Patents

荷電粒子ビーム輸送系及び粒子線治療装置 Download PDF

Info

Publication number
WO2013175600A1
WO2013175600A1 PCT/JP2012/063268 JP2012063268W WO2013175600A1 WO 2013175600 A1 WO2013175600 A1 WO 2013175600A1 JP 2012063268 W JP2012063268 W JP 2012063268W WO 2013175600 A1 WO2013175600 A1 WO 2013175600A1
Authority
WO
WIPO (PCT)
Prior art keywords
particle beam
charged particle
phase
deflection unit
gantry
Prior art date
Application number
PCT/JP2012/063268
Other languages
English (en)
French (fr)
Inventor
菅原 賢悟
周平 小田原
克久 吉田
Original Assignee
三菱電機株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to CN201280073432.4A priority Critical patent/CN104335686B/zh
Priority to EP12877384.3A priority patent/EP2858463A4/en
Priority to US14/382,637 priority patent/US9078335B2/en
Priority to JP2014516579A priority patent/JP5748912B2/ja
Priority to PCT/JP2012/063268 priority patent/WO2013175600A1/ja
Priority to TW101137237A priority patent/TWI486188B/zh
Publication of WO2013175600A1 publication Critical patent/WO2013175600A1/ja

Links

Images

Classifications

    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05HPLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H7/00Details of devices of the types covered by groups H05H9/00, H05H11/00, H05H13/00
    • H05H7/04Magnet systems, e.g. undulators, wigglers; Energisation thereof
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1042X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy with spatial modulation of the radiation beam within the treatment head
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1077Beam delivery systems
    • A61N5/1079Sharing a beam by multiple treatment stations
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N5/1077Beam delivery systems
    • A61N5/1081Rotating beam systems with a specific mechanical construction, e.g. gantries
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05HPLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H11/00Magnetic induction accelerators, e.g. betatrons
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05HPLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H13/00Magnetic resonance accelerators; Cyclotrons
    • H05H13/04Synchrotrons
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05HPLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H7/00Details of devices of the types covered by groups H05H9/00, H05H11/00, H05H13/00
    • H05H7/001Arrangements for beam delivery or irradiation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N5/00Radiation therapy
    • A61N5/10X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy
    • A61N2005/1085X-ray therapy; Gamma-ray therapy; Particle-irradiation therapy characterised by the type of particles applied to the patient
    • A61N2005/1087Ions; Protons
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05HPLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H7/00Details of devices of the types covered by groups H05H9/00, H05H11/00, H05H13/00
    • H05H7/001Arrangements for beam delivery or irradiation
    • H05H2007/002Arrangements for beam delivery or irradiation for modifying beam trajectory, e.g. gantries
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05HPLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H7/00Details of devices of the types covered by groups H05H9/00, H05H11/00, H05H13/00
    • H05H7/04Magnet systems, e.g. undulators, wigglers; Energisation thereof
    • H05H2007/046Magnet systems, e.g. undulators, wigglers; Energisation thereof for beam deflection
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05HPLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H7/00Details of devices of the types covered by groups H05H9/00, H05H11/00, H05H13/00
    • H05H7/04Magnet systems, e.g. undulators, wigglers; Energisation thereof
    • H05H2007/048Magnet systems, e.g. undulators, wigglers; Energisation thereof for modifying beam trajectory, e.g. gantry systems
    • HELECTRICITY
    • H05ELECTRIC TECHNIQUES NOT OTHERWISE PROVIDED FOR
    • H05HPLASMA TECHNIQUE; PRODUCTION OF ACCELERATED ELECTRICALLY-CHARGED PARTICLES OR OF NEUTRONS; PRODUCTION OR ACCELERATION OF NEUTRAL MOLECULAR OR ATOMIC BEAMS
    • H05H2277/00Applications of particle accelerators
    • H05H2277/10Medical devices
    • H05H2277/11Radiotherapy

Definitions

  • the present invention relates to a charged particle beam transport system that transports a charged particle beam composed of charged particles such as protons and heavy particles, and a particle beam therapy system that irradiates an object to be irradiated such as an object or a human body with the transported charged particle beam. Is.
  • a particle beam therapy system is connected to a beam generator for generating a charged particle beam, an accelerator for accelerating the generated charged particle beam, and a charge emitted after being accelerated to energy set by the accelerator.
  • a charged particle beam transport system that transports a particle beam, and a particle beam irradiation device that is installed downstream of the beam transport system and that irradiates the irradiation target with the charged particle beam.
  • the motion of charged particles in a charged particle beam in a charged particle beam transport system is such that the traveling direction of the beam is s direction, perpendicular to the s direction, and the deflection direction of the accelerator is the x direction, the s direction, and the x direction. It is described based on a coordinate system in the y direction perpendicular to the direction.
  • the distribution of charged particles in the beam is not uniform, and in the case of the slow extraction method of the synchrotron as an accelerator, a Gaussian distribution is generally shown in the y direction, but the distribution in the x direction is not.
  • the slow extraction method is a method in which a charged particle beam is extracted little by little over a long period of time.
  • the traveling direction thereof differs individually and changes with time, and the entire beam vibrates at a constant period called betatron vibration.
  • the angle and width of the beam ellipse of the charged particles in the phase space can be changed by a quadrupole electromagnet or the like, the emittance (the area of the ellipse) is kept constant even in this case.
  • emittance is asymmetric in the x and y directions, and this difference cannot be canceled out by a quadrupole electromagnet or the like.
  • the emittance is asymmetric, it is difficult to ensure the uniformity of the dose irradiated to the affected area when planning a treatment plan for particle beam therapy. That is, the fact that the charged particles in the x direction in the affected area have a non-Gaussian distribution complicates treatment planning and causes non-uniformity in the actually irradiated dose, although it is within an allowable range.
  • the spot shape at the affected area of the beam changes depending on the rotation angle of the rotating gantry, which also makes it difficult to make a treatment plan.
  • Patent Document 1 discloses a charged particle beam transport system including an emittance adjusting means such as a skew quadrupole electromagnet or a solenoid electromagnet for the purpose of symmetrizing the emittance.
  • Patent Document 2 discloses a charged particle beam irradiation apparatus provided with a downstream electromagnet including a scatterer and a plurality of quadrupole electromagnets provided downstream of the scatterer in a charged particle beam transport system. By this scatterer and the downstream electromagnet, the emittance is symmetrized, and the distribution of charged particles in the x and y directions is made Gaussian.
  • JP-A-9-265000 (0005th to 0009th stage, FIG. 1) Japanese Patent Laying-Open No. 2006-351339 (stages 0011 to 0013, FIG. 3)
  • emittance adjusting means for symmetrizing emittance before being input to the deflecting electromagnet of the rotating gantry A scatterer and a downstream electromagnet were provided in the charged particle beam transport system. Since the conventional charged particle beam transport system has added emittance adjusting means, a scatterer, and a downstream electromagnet, there is a problem that the charged particle beam transport system becomes large.
  • the present invention has been made to solve the above-described problems.
  • the charged particle beam transport system absorbs the difference in emittance generated when the light is emitted from the accelerator with a slow take-off, and is less rotationally dependent at the isocenter.
  • the purpose is to realize the beam size.
  • the charged particle beam emitted from the accelerator is in a plane perpendicular to the traveling direction of the charged particle beam at the emission position of the accelerator, and the direction in the orbital plane of the accelerator is the x direction.
  • the direction perpendicular to the x direction in the plane perpendicular to the traveling direction of the charged particle beam is the y direction
  • the charged particle beam has a different emittance in the x direction and emittance in the y direction, and can rotate around the isocenter.
  • a charged particle beam transport system for transporting a charged particle beam to a particle beam irradiation apparatus mounted on a rotating gantry, the rotating deflector being mounted on the rotating gantry and rotating around the gantry rotation axis of the rotating gantry;
  • a fixed transport unit from the accelerator to the rotation deflection unit, and the rotation deflection unit has a plurality of deflection electromagnets.
  • the gantry angle is the gantry reference angle
  • the change in the phase space distribution of the charged particle beam that advances from the entrance of the rotating deflection unit to the isocenter is defined as the first phase advance, and the gantry angle is rotated by 90 ° from the gantry reference angle.
  • the phase in the phase space distribution of the charged particle beam is isolated from the entrance of the rotary deflection unit.
  • the phase difference of the phase space distribution of the charged particle beam at the entrance of the rotating deflection unit is set to the average value of the first phase advance and the second phase advance.
  • the phase is determined based on the calculation based on the above.
  • the fixed transport unit calculates the phase space distribution of the charged particle beam at the entrance of the rotation deflection unit, and the phase is an average value of the first phase advance and the second phase advance.
  • X phase emittance which is a feature of slow extraction, at the entrance of the rotating deflection unit of the charged particle beam transport system. Even for a beam with a very small beam size, it is possible to realize a beam size with little rotation dependence at the isocenter.
  • FIG. 1 is a schematic configuration diagram of a particle beam therapy system according to the present invention. It is a figure which shows the structure of the particle beam irradiation apparatus of FIG. It is a figure which shows the structure of the rotation deflection
  • FIG. 1 is a schematic configuration diagram of a particle beam therapy apparatus according to the present invention
  • FIG. 2 is a diagram illustrating a configuration of a particle beam irradiation apparatus according to the present invention
  • FIG. 3 is a diagram showing the configuration of the rotational deflection unit of the charged particle beam transport system according to the present invention
  • FIG. 4 is a side view of the rotational deflection unit of FIG. 3 viewed from the left side.
  • a particle beam therapy apparatus 51 includes a beam generation apparatus 52, a charged particle beam transport system 59, and particle beam irradiation apparatuses 58a and 58b.
  • the beam generator 52 includes an ion source (not shown), a pre-stage accelerator 53, and an accelerator 54.
  • the particle beam irradiation device 58b is installed in a rotating gantry (not shown).
  • the particle beam irradiation device 58a is installed in a treatment room having no rotating gantry.
  • the role of the charged particle beam transport system 59 is in communication between the accelerator 54 and the particle beam irradiation devices 58a and 58b.
  • the accelerator 54 and the particle beam irradiation devices 58 a and 58 b are connected via a vacuum duct 57, and the charged particle beam passes through the vacuum duct 57.
  • a part of the charged particle beam transport system 59 is installed in a rotating gantry (not shown), and a rotating deflection unit 60 installed in the rotating gantry includes a plurality of deflecting electromagnets 10a, 10b, and 10c, and a plurality of quadrupole electromagnets 11a, 11b, 11c, and 11d.
  • the rotation deflection unit 60 rotates with the rotation of the rotating gantry.
  • the fixed transport unit 61 which is a part other than the rotating deflection unit 60 of the charged particle beam transport system 59, includes a plurality of deflection electromagnets 48a, 48b, 48c, 48d, 48e, 48f, 48g, 48h, 48i, and a plurality of quadrupole electromagnets 49a. 49b, 49c, 49d, 49e, 49f, 49g, 49h, 49i, 49j.
  • the reference numeral 48 for the deflection electromagnet of the fixed transport unit 61 is generally used, and 48a to 48i are used in the case of distinction.
  • the reference sign of the quadrupole electromagnet of the fixed transport unit 61 is 49 in general, and 49a to 49j are used in the case of distinction.
  • the reference sign of the deflection electromagnet of the rotating deflection unit 60 is generally 10, and 10a, 10b, and 10c are used in the case of distinction.
  • the charged particle beam which is a particle beam such as a proton beam generated in the ion source, is accelerated by the pre-accelerator 53 and is incident on the accelerator 54 from the incident device 46.
  • the accelerator 54 will be described using a synchrotron as an example.
  • the charged particle beam is accelerated to a predetermined energy.
  • the charged particle beam emitted from the emission device 47 of the accelerator 54 is transported to the particle beam irradiation devices 58a and 58b through the charged particle beam transport system 59.
  • the particle beam irradiation devices 58a and 58b irradiate the irradiation target 45 (see FIG. 2) with a charged particle beam.
  • the reference numeral 58 of the particle beam irradiation apparatus is used as a whole, and 58a and 58b are used in the case of distinction.
  • the charged particle beam 31 emitted from the accelerator 54 is charged in a direction perpendicular to the traveling direction (s direction) of the charged particle beam 31 at the emission position of the accelerator 54 and the direction in the orbital surface of the accelerator 54 is the x direction.
  • the direction perpendicular to the x direction in the plane perpendicular to the traveling direction of the particle beam 31 is the y direction
  • the charged particle beam has a small emittance in the x direction and a large emittance in the y direction.
  • the beam coordinate system of the charged particle beam 31 at the emission position of the accelerator 54 is a beam coordinate system 73.
  • the beam coordinate system of the charged particle beam 31 at the entrance of the rotating deflection unit 60 of the charged particle beam transport system 59 is a beam coordinate system 74.
  • the s-axis of the beam coordinate system is referred to as the s-axis
  • the x-axis of the beam coordinate system is referred to as the x-axis
  • the y-axis of the beam coordinate system is referred to as the y-axis.
  • the particle beam irradiation device 58 includes an X-direction scanning electromagnet 32 and a Y-direction scanning electromagnet 33 that scan the charged particle beam 31 in the X direction and the Y direction that are perpendicular to the charged particle beam 31, and a position monitor 34.
  • the irradiation management device 38 includes an irradiation control computer 39 and an irradiation control device 40.
  • the dose data converter 36 includes a trigger generation unit 42, a spot counter 43, and an inter-spot counter 44.
  • the traveling direction of the charged particle beam 31 is the ⁇ Z direction.
  • the X-direction scanning electromagnet 32 is a scanning electromagnet that scans the charged particle beam 31 in the X direction
  • the Y-direction scanning electromagnet 33 is a scanning electromagnet that scans the charged particle beam 31 in the Y direction.
  • the position monitor 34 detects beam information for calculating a passing position (center of gravity position) and a size of a beam through which the charged particle beam 31 scanned by the X direction scanning electromagnet 32 and the Y direction scanning electromagnet 33 passes.
  • the beam data processing device 41 calculates the passing position (center of gravity position) and size of the charged particle beam 31 based on beam information made up of a plurality of analog signals (beam information) detected by the position monitor 34. Further, the beam data processing device 41 generates an abnormality detection signal indicating an abnormal position or size abnormality of the charged particle beam 31 and outputs this abnormality detection signal to the irradiation management device 38.
  • the dose monitor 35 detects the dose of the charged particle beam 31.
  • the irradiation management device 38 controls the irradiation position of the charged particle beam 31 on the irradiation object 45 based on treatment plan data created by a treatment planning device (not shown), is measured by the dose monitor 35, and is measured by the dose data converter 36.
  • the dose converted into digital data reaches the target dose, the charged particle beam 31 is stopped.
  • the scanning electromagnet power source 37 sets the set currents of the X direction scanning electromagnet 32 and the Y direction scanning electromagnet 33 based on control inputs (commands) to the X direction scanning electromagnet 32 and the Y direction scanning electromagnet 33 output from the irradiation management device 38. Change.
  • the scanning irradiation method of the particle beam irradiation apparatus 58 is changed to a hybrid scanning irradiation method (a method in which the beam is not stopped when the beam irradiation position (spot) is changed), specifically, the irradiation position of the charged particle beam 31 is set.
  • a hybrid scanning irradiation method a method in which the beam is not stopped when the beam irradiation position (spot) is changed
  • the irradiation position of the charged particle beam 31 is set.
  • a description will be given of a method in which the charged particle beam 31 is not stopped when changing, as in a raster scanning irradiation method, and a beam irradiation position is successively moved between spot positions as in a spot scanning irradiation method.
  • the spot counter 43 measures the irradiation dose while the beam irradiation position of the charged particle beam 31 is stopped.
  • the spot-to-spot counter 44 measures the irradiation dose while the beam irradiation position of the charged particle beam 31 is moving.
  • the trigger generation unit 42 generates a dose expiration signal when the dose of the charged particle beam 31 at the beam irradiation position reaches the target irradiation dose.
  • FIG. 3 shows the rotation deflecting unit 60 of the charged particle beam transport system 59 when the rotation angles (hereinafter also referred to as gantry angles as appropriate) rotating around the gantry rotation axis 15 are 0 ° and 180 °.
  • the solid line shows the case where the rotation angle is 0 °
  • the broken line shows the case where the rotation angle is 180 °.
  • the reference numerals for the broken lines are different from those for the solid lines.
  • the charged particle beam 31 passes through each electromagnet in the following order and reaches the particle beam irradiation device 58.
  • the order in which the charged particle beam 31 passes through the electromagnet is the deflection electromagnet 10a, the quadrupole electromagnets 11a and 11b, the deflection electromagnet 10b, the quadrupole electromagnets 11c and 11d, and the deflection electromagnet 10c.
  • the charged particle beam 31 passes through each electromagnet in the following order and reaches the particle beam irradiation device 68 (FIGS. 3 and 4).
  • the reference numerals of the particle beam irradiation apparatus in FIGS. 3 and 4 are different from those of 58 in order to distinguish the rotation angle.
  • the reference numerals of the particle beam irradiation apparatus at the rotation angles of 90 °, 180 °, and 270 ° are 67, 68, and 69, respectively.
  • the order in which the charged particle beam 31 passes through the electromagnet is the deflection electromagnet 22a, the quadrupole electromagnets 23a and 23b, the deflection electromagnet 22b, the quadrupole electromagnets 23c and 23d, and the deflection electromagnet 22c.
  • the isocenter IC is a point where the beam center lines 14 of the charged particle beam 31 intersect when the rotation deflection unit 60 rotates.
  • the beam coordinate system of the charged particle beam that has passed through the particle beam irradiation apparatus at rotation angles of 0 °, 90 °, 180 °, and 270 ° is a beam coordinate system 75, a beam coordinate system 77, a beam coordinate system 76, and a beam coordinate, respectively.
  • the gantry angle which is the rotation angle that rotates around the gantry rotation axis 15, is the beam center line 14 of the charged particle beam 31 in the particle beam irradiation device 58 and the beam coordinate system 74 of the charged particle beam 31 at the entrance of the rotation deflection unit 60. Is defined as an angle with the axis in the y direction.
  • the gantry reference angle is defined as the gantry angle when the beam center line 14 and the y-direction axis (y-axis) of the charged particle beam 31 at the entrance of the rotating deflection unit 60 overlap.
  • the charged particle beam 31 is charged by the charged particle beam transport system 59 so that the emittance in the x direction and the emittance in the y direction are separated from each other and the emittance is maintained. Will be transported by.
  • the gantry reference angle may be any gantry angle, but is 0 ° here.
  • FIG. 4 is a side view of the rotation deflection unit shown in FIG. 3 as viewed from the left side, but is also added when the rotation angles are 90 ° and 270 °.
  • the beam center line 16 is obtained by rotating the beam center line 14 passing through the particle beam irradiation device 58 having a rotation angle of 0 ° and the particle beam irradiation device 68 having a rotation angle of 180 ° by 90 °, and a particle beam having a rotation angle of 90 °.
  • This is a beam center line of the charged particle beam 31 passing through the irradiation device 67 and the particle beam irradiation device 69 having a rotation angle of 270 °.
  • FIG. 4 is a side view of the rotation deflection unit shown in FIG. 3 as viewed from the left side, but is also added when the rotation angles are 90 ° and 270 °.
  • the beam center line 16 is obtained by rotating the beam center line 14 passing through the particle beam irradiation device 58 having
  • the deflection electromagnets 10 a and 10 c having a rotation angle of 0 ° are changed in sign when the rotation angle is 90 ° and 270 °.
  • the deflection electromagnets 20a and 20c are used.
  • the deflection electromagnets 24a and 24c are used.
  • the beam design of the charged particle beam transport system 59 is to adjust the magnet arrangement of the deflecting electromagnets and the quadrupole electromagnets and the strength of each quadrupole electromagnet, and to determine the transport matrix shown in the following formula (1). .
  • the beam design of the rotating deflection unit 60 and the fixed transport unit 61 of the charged particle beam transport system 59 is the same as the beam design of the charged particle beam transport system 59.
  • the rotation deflection unit 60 of the charged particle beam transport system 59 will be described as an example.
  • the transport matrix indicates the relationship between the three-dimensional coordinates of the entrance and the momentum thereof, and the three-dimensional coordinates of the exit and the momentum.
  • (X 1 , y 1 , x ′ 1 , y ′ 1 , t 1 , p t1 ) on the right side are the position (x 1 , y 1 ) of the charged particle at the entrance of the rotating gantry and the inclination (x ′ 1 ) in the s direction.
  • the beam design of the rotation deflection unit 60 of the charged particle beam transport system 59 has a problem of determining the transport matrix with a small number of magnets and parameters.
  • the transport matrix is modified to produce a betatron function ⁇
  • the concept of momentum dispersion function ⁇ has been used.
  • the momentum dispersion function ⁇ is ⁇ pt / pt.
  • Equation (2) becomes Equation (3), and there is no coupling in the x and y directions.
  • the conditions of the entrance and exit of the rotation deflection unit 60 are as follows.
  • Adjustment parameters in equation (3) are m 11 , m 12 , m 21 , m 22 , m 33 , m 34 , m 43 , and m 44 .
  • FIG. 5 and FIG. 6 show the betatron function ⁇ and the momentum dispersion function ⁇ of the rotational deflection unit 60 in FIG.
  • FIG. 5 is a diagram showing a betatron function of charged particles in the isocenter IC from the rotation deflection unit according to the first embodiment of the present invention.
  • the horizontal axis is the distance in the s direction, and the vertical axis is the betatron function ⁇ .
  • the betatron function ⁇ is not shown.
  • FIG. 5A shows the betatron function ⁇ in the x direction
  • FIG. 5B shows the betatron function ⁇ in the y direction.
  • FIG. 5A shows the betatron function ⁇ in the x direction
  • FIG. 5B shows the betatron function ⁇ in the y direction.
  • FIG. 6 is a diagram showing a momentum dispersion function of charged particles in the isocenter IC from the rotation deflection unit according to the first embodiment of the present invention.
  • the horizontal axis is the distance in the s direction
  • the vertical axis is the momentum dispersion function ⁇ .
  • FIG. 6A shows the momentum dispersion function ⁇ in the x direction
  • FIG. 5B shows the momentum dispersion function ⁇ in the y direction.
  • the exit of the rotational deflection unit 60 is the left end.
  • S 14 m.
  • the four quadrupole electromagnets 11a, 11b, 11c, and 11d realize a transport matrix in which the charged particle beam 31 passes through the rotation deflection unit 60 and the beam characteristics of the isocenter IC are good. ing.
  • the beam size in the x direction and the y direction in the isocenter IC is different, in actual irradiation, the beam size is expanded by a scatterer or wobbler electromagnet, and the exit of the rotation deflection unit In this case, the charged particles are scattered and irradiated to such an extent that it does not depend on the original beam size.
  • the scanning irradiation method requires a beam size as small as possible, it is desired that the scattering is as small as possible.
  • the beam sizes of the accelerator 54 and the charged particle beam transport system 59 have the same beam sizes in the x and y directions, and the beam sizes in the x and y directions are those of the rotating gantry. It is desirable that there is no angle dependency.
  • the x-direction beam size ⁇ x and the y-direction beam size ⁇ y are generally expressed as in Expression (4).
  • ⁇ x and ⁇ y are the emittance in the x direction and the emittance in the y direction, respectively.
  • Emittances ⁇ x and ⁇ y are areas on the phase space.
  • the phase space in the x direction is represented as shown in FIG. 7A, where the horizontal axis is the distance x in the x direction, and the vertical axis is the slope x ′ in the x direction.
  • the phase space in the y direction is represented as shown in FIG. 7B, where the horizontal axis is the distance y in the y direction and the vertical axis is the slope y ′ in the y direction.
  • the betatron functions ⁇ x and ⁇ y are quantities (quantities depending on the transport matrix) that can be changed by the beam design of the rotating gantry.
  • the emittances ⁇ x and ⁇ y are amounts determined by the emission of the accelerator 54, and are constant when there is no scatterer or the like during transport of charged particles in the charged particle beam transport system 59. . Further, the emittances ⁇ x and ⁇ y depend on the take-out mechanism from the accelerator 54, and the emittances ⁇ x and ⁇ y have different values when emitted from the accelerator 54 with slow take-out.
  • emittance ⁇ y is much larger, and emittance ⁇ y is 2 to 10 times or more of emittance ⁇ x .
  • the difference in beam size due to the difference in emittance ⁇ x and ⁇ y is minimized by the beam design of the charged particle beam transport system 59.
  • the rotation deflecting unit 60 has a configuration in which a scatterer and a quadrupole electromagnet are combined, or the intensity of each quadrupole electromagnet is changed for each rotation angle of the rotating gantry even if the number of the quadrupole electromagnets of the rotating deflection unit 60 is small, It is possible to eliminate the difference in size as much as possible. However, it is important that the rotational deflection unit 60 can be realized at low cost by making no difference in beam size for each rotation angle without changing the strength of the quadrupole electromagnet.
  • Such a beam design method for the charged particle beam transport system 59 can be applied to an existing facility, and it is desirable that the rotating deflection unit 60 of the charged particle beam transport system 59 can be realized also in the existing facility. From the equation (3), the beam design of the charged particle beam transport system 59 is divided into a 2 ⁇ 2 transport matrix for x and x ′ and a 2 ⁇ 2 transport matrix for y and y ′. You can see that it can be done.
  • FIG. 7 is a diagram showing a phase space distribution of charged particles at the entrance of the rotational deflection unit according to Embodiment 1 of the present invention.
  • FIG. 7A shows the phase space distribution in the x direction, and FIG.
  • phase space distribution in the y direction shows the phase space distribution in the y direction.
  • the phase space distribution in the x direction at the entrance of the rotating deflection unit 60 is limited to a distribution in the x direction due to the characteristics specific to the emission from the accelerator 54, and is a distribution closer to a line than an ellipse.
  • the phase space distribution in the y direction is an ellipse because it is not limited to the y direction. Since the area in the phase space is emittance, there is a large difference between the x direction and the y direction.
  • the beam design of the charged particle beam transport system 59 includes the beam design of the rotation deflection unit 60 in which the charged particle beam 31 passes through the rotation deflection unit 60 and the beam characteristic of the isocenter IC is good. It is divided into the beam design of the fixed transport unit 61 that realizes the phase space distribution of the charged particles that are input to the rotation deflection unit 60 for realizing a beam size that is less dependent on rotation by the isocenter IC.
  • the phase space distribution of the charged particles to be input to the rotation deflection unit 60 for realizing a beam size with little rotation dependency by the isocenter IC is determined.
  • the phase space distribution in the x direction approximates a line
  • the phase space distribution in the y direction approximates an ellipse. Since the phase space distribution in the x direction is approximated by a line, the phase is important, not the area.
  • the phase of the phase space distribution in the x direction at the entrance of the rotation deflection unit 60 is determined by paying attention to the phase relationship of the rotation gantry. Specifically, the phase ⁇ ix of the phase space distribution in the x direction at the entrance of the rotational deflection unit 60 is the distance from the entrance to the exit of each rotational deflection unit 60 at two gantry angles ⁇ 1 and ⁇ 2 that differ by 90 °.
  • the phase is determined based on the calculation based on the average value of the phase advance (phase advance). For example, each of the fractional portions of the phase advance (phase advance) from the entrance to the exit at the gantry angles of 0 ° and 90 ° is averaged. The average value is a minus (-).
  • FIG. 3 is an electromagnet arrangement of the rotary deflection unit 60 shown in FIG. 3, and the transport matrix of the rotary deflection unit 60 is set so as to have the strength of each electromagnet that becomes the betatron function and the momentum dispersion function shown in FIGS. 5 and 6.
  • deviation part 60 is designed using is considered.
  • the phase advance ⁇ 1 is 1.382 (497.52 °), and 497 at the exit of the rotation deflection unit 60. .52 ° rotation.
  • the phase advance ⁇ 2 is 0.522. It is rotated 93.96 ° at the outlet.
  • the phase space distribution of the charged particle beam 31 at a gantry angle of 180 ° theoretically matches the phase space distribution of the charged particle beam 31 at a gantry angle of 0 °.
  • the beam size in the isocenter IC is substantially equal to the beam size at the exit of the rotary deflection unit 60.
  • the beam size in the isocenter IC slightly changes compared to the case where the charged particle beam 31 is not scanned. Even so, the beam size in the isocenter IC is almost equal to the beam size at the exit of the rotary deflection unit 60.
  • FIG. 8 to 15 show isocenter ICs according to Embodiment 1 of the present invention when the angles of the rotation deflection units are 0 °, 30 °, 60 °, 90 °, 120 °, 150 °, 180 °, and 210 °, respectively.
  • the horizontal axis is the distance x in the x direction
  • the vertical axis is the inclination x ′ in the x direction.
  • the horizontal axis is the distance y in the y direction
  • the vertical axis is the inclination y ′ in the y direction.
  • rms means root mean square.
  • the beam size (rms) shown in the figure is a size corresponding to the standard deviation ⁇ .
  • the beam size (rms) will be referred to as the beam size 1 ⁇ .
  • the beam sizes 1 ⁇ in the x and y directions are 1.29 mm and 1.29 mm, respectively, as shown in FIG.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotating deflection unit is 30 ° are 1.30 mm and 1.30 mm, respectively.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotation deflection unit is 60 ° are 1.30 mm and 1.30 mm, respectively.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotating deflection unit is 90 ° are 1.30 mm and 1.30 mm, respectively.
  • the beam size 1 ⁇ in the x direction and the y direction is 1.29 mm and 1.30 mm, respectively, as shown in FIG.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotating deflection unit is 150 ° are 1.29 mm and 1.29 mm, respectively.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotation deflection unit is 180 ° are 1.29 mm and 1.30 mm, respectively.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotation deflection unit is 210 ° are 1.30 mm and 1.30 mm, respectively.
  • the shape and phase of the phase space distribution in the x direction and the y direction are different, but the beam sizes 1 ⁇ in the x direction and the y direction are substantially equal.
  • the reason for adding a minus ( ⁇ ) to the average value of the phase advance ⁇ 1 and ⁇ 2 will be described.
  • the selection of the phase at the entrance of the rotary deflection unit 60 is important. If the phase at the entrance is poorly selected, the phase space distribution of the beam in the phase space at the exit will be in the vertical direction, so that the beam size becomes too small to achieve the desired beam size, which is not desirable.
  • phase advance ⁇ 1 and ⁇ 2 cancels the phase advance of the gantry at the exit (isocenter) of the rotation deflection unit 60. Since the phase advances from the entrance to the exit in the rotating gantry, the state before the phase advance is desirable at the entrance. For example, if the rotation gantry advances 60 °, setting the phase to ⁇ 60 ° cancels the phase advance in the rotation deflection unit 60 at the exit (isocenter) of the rotation deflection unit 60.
  • FIG. 16 is a diagram showing the relationship between the angle of the rotation deflection unit and the beam size in the isocenter IC according to the first embodiment of the present invention.
  • the horizontal axis is the gantry angle ⁇ of the rotating deflection unit 60
  • the vertical axis is the beam size 3 ⁇ which is a size corresponding to three times the standard deviation ⁇ .
  • a characteristic 71 indicated by a solid line is an x-direction characteristic
  • a characteristic 72 indicated by a broken line is a y-direction characteristic.
  • the phase ⁇ ix of the phase space distribution in the x direction at the entrance of the rotation deflector 60 is the phase advance (phase advance) from the entrance to the exit of the respective rotation deflectors 60 at two gantry angles ⁇ 1 and ⁇ 2 that differ by 90 °.
  • phase ⁇ ix of the phase space distribution in the x direction at the entrance of the rotating deflection unit 60 is added to the fixed transport unit 61 of the charged particle beam transport system 59 or the quadrupole Shift the position of the electromagnet or change the setting value of the excitation amount.
  • the fixed transport unit 61 of the charged particle beam transport system 59 is designed to satisfy the following three conditions.
  • Condition 1 is that the charged particle beam 31 passes through the fixed transport unit 61.
  • Condition 2 is that the momentum dispersion function ⁇ and the phase space distribution of the charged particles in the y direction satisfy the conventional conditions.
  • Condition 3 is that the phase space distribution of the charged particles in the x direction at the entrance of the fixed transport unit 61 is a predetermined phase space distribution input to the rotation deflection unit 60 for realizing a beam size with little rotation dependency in the isocenter IC. It is to become.
  • the predetermined phase space distribution in the rotation deflection unit 60 of the charged particle beam transport system 59 of the first embodiment is the phase ⁇ ix of the phase space distribution in the x direction at the entrance of the rotation deflection unit 60 described above, and the gantry angle is 90.
  • ° Minus (-) to the average value obtained by averaging the fractional parts obtained by subtracting the half-integer part from the phase advance (phase advance) from the entrance to the exit of each of the rotational deflection parts 60 at two different gantry angles ⁇ 1 and ⁇ 2. It is attached.
  • FIG. 40 shows a case where the phase at the entrance of the rotating gantry is
  • the phase line 101 is the phase at the entrance of the rotating gantry
  • the phase line 102 is the phase of the isocenter IC when the rotating gantry is 0 °
  • the phase line 103 is the phase when the rotating gantry is 90 °. This is the phase of the isocenter IC.
  • the part excluding the half-integer part of 2.9 is 0.4.
  • the part excluding the half integer part of 0.8 is 0.3.
  • the following equation is obtained and the result is the same.
  • the expansion of the isocenter IC in the x direction is almost the same.
  • the angle between 0 ° and 90 ° is not simple because there is coupling in the x direction and the y direction.
  • the beam size can be made constant regardless of the gantry angle because it is approximately between these angles.
  • FIG. 17 is a diagram for explaining the phase space distribution of the charged particles in the x direction in the charged particle beam transport system of the present invention.
  • FIG. 17A is an example in which the x-direction phase space distribution becomes the x-direction phase distribution 98 at the exit of the fixed transport section 61
  • FIG. 17B is the x-direction phase space at the exit of the fixed transport section 61.
  • the distribution is an x-direction phase distribution 99.
  • the x-direction phase distribution 97 is a phase space distribution at the exit of the accelerator 54.
  • the x-direction phase distributions 98 and 99 are phase space distributions at the exit of the fixed transport unit 61, that is, the entrance of the rotational deflection unit 60.
  • the x-direction phase distributions 97, 98, and 99 are shown as linear distributions.
  • the coordinates of the end point p1 of the x-direction phase distribution 97 are (x 1 , x ′ 1 ), and the coordinates of the end point p2 of the x-direction phase distribution 97 are ( ⁇ x 1 , ⁇ x ′ 1 ).
  • the phase ⁇ 1 is a phase in the x-direction phase distribution 97
  • the phase ⁇ 2 is a phase in the x-direction phase distribution 98 or the x-direction phase distribution 99.
  • the end point p1 at the entrance of the rotation deflection unit 60 is the end point p11 (x 2 , x ′ 2 ), and the end point p2 at the entrance of the rotation deflection unit 60 is the end point p12 ( ⁇ x 2). , ⁇ x ′ 2 ).
  • the difference between the x-direction phase distribution 98 and the x-direction phase distribution 99 is that the quadrant in the phase space where the end point p11 exists is reversed.
  • Condition 2 is that ⁇ y , ⁇ y , ⁇ x , ⁇ ′ x , ⁇ y , and ⁇ ′ y satisfy the conventional conditions at the entrance of the rotating deflection unit 60, that is, ⁇ x , ⁇ ′ x , ⁇ y and ⁇ ′ y are all 0, and ⁇ y and ⁇ y are values satisfying the beam condition in the isocenter IC.
  • This condition 2 will be referred to as constraint condition 1.
  • Condition 3 is that an x-direction phase distribution 98 having an end point p11 (x 2 , x ′ 2 ) and an end point p12 ( ⁇ x 2 , ⁇ x ′ 2 ) at the exit of the fixed transport unit 61 (inlet of the rotation deflection unit 60). And determining the transport matrix M so that the x-direction phase distribution 99 is obtained.
  • This condition 3 will be referred to as constraint condition 2. Since the transport matrix M is a function of the entrance and exit twist parameters of the fixed transport unit 61 of the charged particle beam transport system 59 and the phase advance, a condition can be expressed by using the twist parameters as shown in Equation (5). become. Note that the twist parameter is a term used when the betatron functions ⁇ and ⁇ are used as parameters.
  • ⁇ 1 and ⁇ 1 are twist parameters at the entrance of the fixed transport section 61
  • ⁇ 2 and ⁇ 2 are twist parameters at the exit of the fixed transport section 61
  • is a phase from the entrance to the exit of the fixed transport section 61.
  • Advance phase advance
  • x 1 and x ′ 1 are determined by the output design of the accelerator 54
  • x 2 and x ′ 2 are determined by the design of the rotating gantry (rotational deflecting unit 60) when the beam size is not angle-dependent.
  • the phase ⁇ ix of the phase space distribution in the x direction is different from each of two gantry angles ⁇ 1 and ⁇ 2 that differ by 90 ° in gantry angle.
  • the average value obtained by averaging the fractional parts of the phase advance from the entrance to the exit of the rotation deflection unit 60 is minus ( ⁇ ). That is, the phase of the x-direction phase distribution 98 (or 99) having the end point p11 (x 2 , x ′ 2 ) and the end point p12 ( ⁇ x 2 , ⁇ x ′ 2 ) is two gantry angles that are 90 ° different from each other.
  • the phase ⁇ ix is obtained by adding a minus ( ⁇ ) to the average value obtained by averaging the fractional parts of the phase advance (phase advance) from the entrance to the exit of each of the rotation deflectors 60 of ⁇ 1 and ⁇ 2.
  • the phase ⁇ 2 of the x-direction phase distribution 98 (or 99) is the phase ⁇ ix .
  • Step S1 As in the past, the arrangement and polarity of the quadrupole electromagnet 49 are determined using the twist parameter as a condition. Here, if the arrangement of the quadrupole electromagnet 49 is added to the degree of freedom, the subsequent process becomes too long, and thus the base of the quadrupole electromagnet 49 is made. By doing so, the subsequent steps can be shortened.
  • Step S2 An arbitrary ellipse is written so as to be in contact with the phase space distribution in the x direction at the exit of the fixed transport section 61 (inlet of the rotation deflection section 60), and the twist parameter is calculated. The phase ⁇ 2 of the ellipse is also calculated. Since the beam distribution at the exit of the fixed transport section 61 (the entrance of the rotation deflection section 60) is important, the exit is used as a reference.
  • Step S3 The twist parameter at the entrance of the fixed transport unit 61 is calculated using the equation (5).
  • is assumed to be an arbitrary value
  • the inlet phase ⁇ 1 is calculated from ⁇ 2 and ⁇ .
  • the steps S2 and S3 are repeated, and the ellipse (emitter is naturally the value of the exit) corresponding to the entrance twist parameter is the beam distribution at the entrance. Is included. This is not the case when tracking is performed every time and the beam size is confirmed.
  • Step S4 The fixed transport unit 61 is designed so as to satisfy the twist parameter and the constraint condition 2 at the entrance and the exit of the fixed transport unit 61.
  • the conventional design method corresponds to the execution of step S4.
  • Step S5 Steps S2 to S4 are repeated until an appropriate solution is found in step S4. When there is no solution, it returns to step S1 and arrangement
  • the rotating deflection unit 60 of the charged particle beam transport system 59 of the first embodiment is a conventional one without changing the strength of the electromagnets of the three deflecting electromagnets 10a, 10b, and 10c and the four quadrupole electromagnets 11a, 11b, 11c, and 11d.
  • the emittance is reduced at the entrance of the rotary deflection unit 60 of the charged particle beam transport system 59. Even if there is an asymmetry, a beam size with little rotation dependency can be realized with an isocenter IC.
  • the charged particle beam 31 emitted from the accelerator 54 is in a plane perpendicular to the traveling direction of the charged particle beam 31 at the emission position of the accelerator 54.
  • the direction in the orbital plane of the accelerator is the x direction and the direction perpendicular to the x direction in the plane perpendicular to the traveling direction of the charged particle beam 31 is the y direction
  • the emittance in the x direction is different from the emittance in the y direction.
  • a charged particle beam transport system 59 that transports a charged particle beam 31 to a particle beam irradiation device 58 mounted on a rotating gantry that is a charged particle beam and is rotatable about an isocenter IC, and is mounted on the rotating gantry.
  • a rotating deflection unit 60 that rotates around the gantry rotating shaft 15 of the rotating gantry, and a fixed transport unit 61 from the accelerator 54 to the rotating deflection unit 60.
  • the polygonal portion 60 has a plurality of deflection electromagnets 10.
  • the charged particle beam transport system 59 when viewed from the gantry rotation axis direction, the beam center line 14 of the particle beam irradiation device 58 and the y of the charged particle beam 31 at the entrance of the rotation deflection unit 60 are displayed.
  • the rotating gantry through which the charged particle beam 31 is transported so that the emittance in the x direction and the emittance in the y direction are separated from each other and the emittance in the y direction is maintained at an exit position of the accelerator 54 is an angle with the direction axis.
  • the phase phase distribution of the charged particle beam 31 advances from the entrance of the rotary deflection unit 60 to the isocenter IC as the first phase advance.
  • the position of the charged particle beam 31 is A change in the phase of the spatial distribution from the entrance of the rotary deflection unit 60 to the isocenter IC is defined as a second phase advance, and the fixed transport unit 61 has a phase of the phase space distribution of the charged particle beam 31 at the entrance of the rotary deflection unit 60.
  • the fixed transport unit 61 is provided at the entrance of the rotation deflection unit 60.
  • the phase ⁇ ix of the phase space distribution in the x-direction at x can be set to a predetermined phase, that is, the phase space distribution of the charged particle beam 31 at the entrance of the rotary deflecting unit 60 has the first phase advance and the second phase.
  • a predetermined phase space distribution such as a phase determined based on an arithmetic operation based on an average value of the advance can be obtained, which is a feature of slow extraction at the entrance of the rotating deflection unit 60 of the charged particle beam transport system 59. Also with respect to the direction emittance is very small such beams, it is possible to realize a small beam sizes rotation dependencies isocenter IC.
  • the charged particle beam 31 is generated, the beam generation apparatus 52 that accelerates the charged particle beam 31 by the accelerator 54, and the charge accelerated by the accelerator 54.
  • the particle beam 31 is in a plane perpendicular to the traveling direction of the charged particle beam 31 at the exit position of the accelerator 54, the direction in the orbital plane of the accelerator 54 is the x direction, and in the plane perpendicular to the traveling direction of the charged particle beam 31. If the direction perpendicular to the x direction is the y direction, the charged particle beam transport system 59 for transporting the charged particle beam 31 and the charged particle beam 31 are charged particle beams having different emittance in the x direction and emittance in the y direction.
  • a particle beam irradiation device 58 for irradiating the irradiation target 45 with the charged particle beam 31 transported by the transport system 59 and a particle beam irradiation device are mounted, and an isocenter is mounted.
  • the charged particle beam transport system 59 is mounted on the rotating gantry and rotated from the rotating deflector 60 rotating around the gantry rotating shaft 15 of the rotating gantry and the accelerator 54.
  • a rotating transport unit 61 having a plurality of deflecting electromagnets 10 when viewed from the direction of the gantry rotation axis, and the beam center line 14 of the particle beam irradiation device 58;
  • An angle between the charged particle beam 31 and the y-direction axis at the entrance of the rotating deflection unit 60 is a gantry angle, and the emittance in the x direction and the emittance in the y direction at the emission position of the accelerator 54 are separated to maintain each emittance.
  • the angle of the rotating gantry through which the charged particle beam 31 is transported is the gantry reference angle, and the gantry angle is the gantry reference angle.
  • the gantry angle is an angle rotated by 90 ° from the gantry reference angle.
  • a change in the phase of the charged particle beam 31 in the phase space distribution that advances from the entrance of the rotary deflection unit 60 to the isocenter IC is defined as a second phase advance, and the fixed transport unit 61 uses the charged particle beam at the entrance of the rotary deflection unit 60.
  • the fixed transport unit 61 but the phase phi ix phase space distribution in the x direction at the inlet of the rotating deflecting unit 60 can be a predetermined phase, i.e., the charged particle beam at the entrance of the rotating deflecting unit 60 1 phase space distribution can be made into a predetermined phase space distribution such that the phase is a phase determined based on an operation based on an average value of the first phase advance and the second phase advance, and the charged particle beam At the entrance of the rotation deflection unit 60 of the transport system 59, a beam size with little rotation dependency can be realized with the isocenter IC even for a beam having a very small x-direction emittance, which is a characteristic of slow extraction. .
  • FIG. A charged particle beam transport system 59 according to the second embodiment of the present invention will be described.
  • the phase space distribution in the x direction is approximated to a linear distribution having end points, and the linear distribution having end points calculated using a matrix obtained by extracting a part of the transport matrix of Equation (3) is obtained.
  • the fixed transport unit 61 of the charged particle beam transport system 59 is configured to be realized at the entrance of the rotation deflection unit 60 of the charged particle beam transport system 59. In this way, by controlling the phase space distribution in the x direction of the charged particle beam 31 at the entrance of the rotating deflection unit 60, the difference in emittance generated when the charged particle beam 31 is extracted from the accelerator 54 by slow extraction is charged. And an isocenter IC can realize a beam size with little rotation dependency.
  • the deflection electromagnet of the rotary deflection unit 60 is composed only of the deflection electromagnets that deflect in the parallel plane and the vertical plane of the deflection electromagnet of the fixed transport unit 61, the x direction and the y direction are the same as in the first embodiment.
  • a transport matrix as shown in Equation (3) without coupling is obtained.
  • Adjustment parameters in equation (3) are m 11 , m 12 , m 21 , m 22 , m 33 , m 34 , m 43 , and m 44 .
  • FIG. 18 and 19 show the betatron function ⁇ and the momentum dispersion function ⁇ of the rotational deflection unit 60 (see FIG. 3) according to the second embodiment of the present invention.
  • FIG. 18 is a diagram illustrating a betatron function of charged particles in the isocenter IC from the rotation deflection unit according to the second embodiment of the present invention. The horizontal axis is the distance in the s direction, and the vertical axis is the betatron function ⁇ . In FIG. 18, the betatron function ⁇ is not shown.
  • FIG. 18A shows the betatron function ⁇ in the x direction
  • FIG. 18B shows the betatron function ⁇ in the y direction.
  • FIG. 19 is a diagram showing a momentum dispersion function of charged particles in the isocenter IC from the rotation deflection unit according to the second embodiment of the present invention.
  • the horizontal axis is the distance in the s direction, and the vertical axis is the momentum dispersion function ⁇ .
  • FIG. 19A shows the momentum dispersion function ⁇ in the x direction
  • FIG. 19B shows the momentum dispersion function ⁇ in the y direction.
  • the exit of the rotation deflection unit 60 is the left end.
  • S 14 m.
  • a transport matrix with good beam characteristics of the isocenter IC is realized by the four quadrupole electromagnets 11a, 11b, 11c, and 11d.
  • the electromagnet arrangement of the rotary deflection unit 60 according to the second embodiment of the present invention is the electromagnet arrangement of the rotary deflection unit 60 shown in FIG. 3, so that the betatron function and the momentum dispersion function shown in FIGS.
  • a part of the transport matrix of the rotating gantry is selected and extracted.
  • a matrix obtained by extracting a part of the transport matrix of Expression (3) is Expression (7).
  • a matrix obtained by selecting and extracting a part of the transport matrix will be referred to as a selection matrix.
  • Each of m 11 , m 12 , m 33 , and m 34 in the matrix of Expression (7) represents the beam size in the x direction and the beam size in the y direction in the isocenter IC, the rotational gantry entrance, that is, the entrance of the rotational deflection unit 60. Is a quantity relating the distribution in the phase space. Based on this amount, an ideal linear distribution on the phase space in the x direction at the entrance of the rotating gantry is determined as follows.
  • Equation (8) The quantity characterizing the linear beam on the phase space is the coordinates (x, x ') of the end point on the phase space. Therefore, the simultaneous equations of Equation (8) are derived using the transport matrix of Equation (7) with the coordinates (x, x ′) of the end points as unknowns.
  • equation (8) The meaning of equation (8) will be described.
  • the equation is for realizing the same beam size when the gantry angle is 0 ° or 180 ° and when the gantry angle is 90 ° or 270 °. If an optical system (charged particle beam transport system 59) capable of realizing substantially the same beam size at the two angles is realized, the difference between the two angles does not cause a problem.
  • an optical system charged particle beam transport system 59
  • the coordinates (x, x ′) of the end point obtained by explaining the equation (8) are taken as one end, and coordinates that are point-symmetric about the origin of the coordinates (x, x ′) of the end point, that is, ( ⁇ x, If a linear distribution in the phase space having ⁇ x ′) at the other end is realized at the gantry entrance (inlet of the rotation deflection unit 60), it can pass through the rotation deflection unit 60 designed by the conventional setting method. In the isocenter IC, the rotation dependency of the beam size is eliminated.
  • phase space distribution of the charged particles at the entrance of the rotary deflection unit 60 according to the second embodiment obtained by the above procedure is as shown in FIG. Further, the phase space distribution of the charged particles at the isocenter IC of the rotary deflection unit 60 according to the second embodiment is as shown in FIGS.
  • FIGS. 21 to 28 show isocenter ICs in the case where the angles of the rotation deflection units according to the second embodiment of the present invention are 0 °, 30 °, 60 °, 90 °, 120 °, 150 °, 180 °, and 210 °, respectively.
  • the horizontal axis is the distance x in the x direction
  • the vertical axis is the inclination x ′ in the x direction.
  • the horizontal axis is the distance y in the y direction
  • the vertical axis is the inclination y ′ in the y direction.
  • the beam size (rms) shown in the figure is the beam size 1 ⁇ described in the first embodiment.
  • the beam sizes 1 ⁇ in the x direction and the y direction are 1.29 mm and 1.29 mm, respectively, as shown in FIG.
  • the beam sizes 1 ⁇ in the x and y directions are 1.30 mm and 1.29 mm, respectively, as shown in FIG.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotating deflection unit is 60 ° are 1.30 mm and 1.30 mm, respectively.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotating deflection unit is 90 ° are 1.30 mm and 1.30 mm, respectively.
  • the beam size 1 ⁇ in the x direction and y direction is 1.29 mm and 1.30 mm, respectively, as shown in FIG.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotation deflection unit is 150 ° are 1.29 mm and 1.30 mm, respectively.
  • the beam sizes 1 ⁇ in the x direction and the y direction when the angle of the rotation deflection unit is 180 ° are 1.29 mm and 1.29 mm, respectively.
  • the angle of the rotation deflection unit is 210 °
  • the beam sizes 1 ⁇ in the x direction and the y direction are 1.30 mm and 1.29 mm, respectively, as shown in FIG.
  • FIG. 29 is a diagram showing the relationship between the angle of the rotation deflection unit and the beam size in the isocenter IC according to the second embodiment of the present invention.
  • the horizontal axis is the gantry angle ⁇ of the rotating deflection unit 60
  • the vertical axis is the beam size 3 ⁇ which is a size corresponding to three times the standard deviation ⁇ .
  • a characteristic 81 indicated by a solid line is an x-direction characteristic
  • a characteristic 82 indicated by a broken line is a y-direction characteristic.
  • FIGS. 33 to 37 show the entrance of the rotation deflection unit when the angles of the rotation deflection unit according to the second embodiment of the present invention are 0 °, 30 °, 60 °, 90 °, 120 °, and 150 °, respectively. It is a figure which shows the relationship of the beam size of the charged particle in an isocenter. The horizontal axis is the distance in the s direction, and the vertical axis is the beam size 3 ⁇ .
  • sections S1, S3, and S5 are sections that pass through the deflection electromagnets 10a, 10b, and 10c, respectively, a section S2 is a section where the quadrupole electromagnets 11a and 11b are arranged, and a section S4 is the quadrupole electromagnet 11c, This is a section in which 11d is arranged.
  • the section S6 is a section from the outlet of the deflection electromagnet 10c, that is, the outlet of the rotary deflection unit 60 to the isocenter IC.
  • FIG. 31 is a diagram showing the x-direction beam size when the angle of the rotating deflection unit according to the second embodiment of the present invention is 0 °
  • FIG. 32 shows the angle of the rotating deflection unit according to the second embodiment of the present invention being zero. It is a figure which shows the y direction beam size in the case of (degree).
  • the x direction characteristic 83 is a result obtained by simulating the phase space distribution in the x direction as the linear distribution of FIG.
  • the y-direction characteristic 84 is a result obtained by simulating the phase space distribution in the y-direction as the elliptic distribution in FIG.
  • the x-direction characteristic 85 is a result obtained by simulation using a twist parameter in the x-direction.
  • the y-direction characteristic 86 is a result obtained by simulation using the y-direction twist parameter.
  • the x-direction characteristic 83 and the x-direction characteristic 85 do not match very much. This indicates that when the phase space distribution in the x direction is approximated to an elliptic distribution, an accurate beam size cannot be obtained.
  • the y-direction characteristic 84 and the y-direction characteristic 86 are in good agreement. Although there is a shift in the section S1, the shift is small and negligible.
  • the rotating deflection unit 60 of the charged particle beam transport system 59 of the second embodiment is a conventional one without changing the strength of the electromagnets of the three deflection electromagnets 10a, 10b, and 10c and the four quadrupole electromagnets 11a, 11b, 11c, and 11d.
  • the predetermined phase space distribution in the x-direction is a point centered on the origin of the coordinates (x, x ′) of the end point with the end point coordinates (x, x ′) obtained by explaining Equation (8) as one end. It is a linear distribution on the phase space having symmetrical coordinates, that is, ( ⁇ x, ⁇ x ′) at the other end.
  • the phase space distribution in the predetermined x-direction is obtained by adding a quadrupole electromagnet in the fixed transport unit 61 of the charged particle beam transport system 59 or shifting the position of the quadrupole electromagnet and re-optimizing the excitation amount,
  • the required particle distribution is achieved at 60 inlets.
  • the charged particle beam transport system 59 of the second embodiment has fewer constraint conditions than the first embodiment and can achieve the object. Therefore, the number of adjustment elements can be reduced, and further downsizing can be achieved.
  • the phase space distribution in the x direction at the entrance of the rotating deflection unit 60 is changed to the phase space distribution in the predetermined x direction described above. Even if emittance asymmetry is present at the entrance of the rotation deflection unit 60, the isocenter IC can realize a beam size with little rotation dependency.

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Plasma & Fusion (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Optics & Photonics (AREA)
  • Radiation-Therapy Devices (AREA)
  • Particle Accelerators (AREA)

Abstract

 加速器から遅い取り出しで出射したときに発生するエミッタンスの違いを荷電粒子ビーム輸送系で吸収し、アイソセンタで回転依存性の少ないビームサイズを実現することを目的とする。 本発明の荷電粒子ビーム輸送系(59)は、その固定輸送部(61)が、回転ガントリのガントリ回転軸(15)の周りを回転する回転偏向部(60)の入口における荷電粒子ビーム(31)の位相空間分布の位相が、第1の位相進み及び第2の位相進みの平均値に基づいた演算を基に決定された位相となるようにしたことを特徴とする。第1の位相進みは、ガントリ角度がガントリ基準角度である場合に、位相空間分布における位相が回転偏向部(60)の入口からアイソセンタ(IC)まで進む変化分であり、第2の位相進みは、ガントリ角度がガントリ基準角度から90°回転した角度である場合の上記変化分である。

Description

荷電粒子ビーム輸送系及び粒子線治療装置
 本発明は、陽子や重粒子などの荷電粒子からなる荷電粒子ビームを輸送する荷電粒子ビーム輸送系及び、輸送された荷電粒子ビームを物体、人体等の被照射体に照射する粒子線治療装置に関するものである。
 一般に粒子線治療装置は、荷電粒子ビームを発生するビーム発生装置と、ビーム発生装置につながれ、発生した荷電粒子ビームを加速する加速器と、加速器で設定されたエネルギーまで加速された後に出射される荷電粒子ビームを輸送する荷電粒子ビーム輸送系と、ビーム輸送系の下流に設置され、荷電粒子ビームを照射対象に照射するための粒子線照射装置とを備える。
 一般に、荷電粒子ビーム輸送系における荷電粒子ビーム中の荷電粒子の運動は、ビームの進行方向をs方向、s方向に対して垂直であり、かつ加速器の偏向方向であるx方向、s方向及びx方向に垂直なy方向の座標系に基づいて記述される。ここで、荷電粒子ビーム輸送系におけるビーム中の荷電粒子の分布は均一ではなく、加速器であるシンクロトロンの遅い取り出し方式の場合は、一般にy方向にはガウス分布を示すが、x方向については非ガウス分布となっており、ビームの位相空間上の面積に相当するエミッタンスは、x方向、y方向で非対称になる。なお、遅い取り出し方式は、長期間に渡って少しずつ荷電粒子ビームを取り出す方式である。
 また、個々の荷電粒子の運動に注目すると、その進行方向は個々に異なり、かつ、経時的に変化しており、ビーム全体ではベータトロン振動と呼ばれる一定周期の振動をしている。位相空間上における荷電粒子のビーム楕円の角度や幅は四極電磁石等によって変えることはできるが、この場合でもエミッタンス(楕円の面積)は一定に保たれる。一般に、荷電粒子ビーム輸送系おいては、x方向とy方向ではエミッタンスが非対称になっており、四極電磁石等によってもこの差を打ち消すことはできない。
 上記のようにエミッタンスが非対称であることにより、粒子線治療の治療計画を立案するとき、患部に照射される線量の均一性を確保することが困難になる。すなわち、患部でx方向の荷電粒子が非ガウス分布であることは治療計画立案を複雑にし、また許容範囲内ではあるが、実際に照射される線量の不均一を生じる原因になる。また、患部に集中的に粒子線を照射しつつ、正常組織の被爆量を最小限にするには、患部の周囲に粒子線照射装置を旋回させる回転ガントリを用いることが好ましい。しかし、エミッタンスの非対称性があると、回転ガントリの回転角度によりビームの患部におけるスポット形状が変化することになり、これも治療計画の立案を困難にする原因となる。
 特許文献1には、エミッタンスを対称化することを目的として、スキュー4極電磁石やソレノイド電磁石であるエミッタンス調整手段を備えた荷電粒子ビーム輸送系が示されている。特許文献2には、荷電粒子ビームの輸送系に散乱体とこの散乱体の下流側に設けられた複数の四極電磁石からなる下流側電磁石を備えた荷電粒子線照射装置が示されている。この散乱体及び下流側電磁石により、エミッタンスの対称化を図り、かつ、x方向及びy方向の荷電粒子の分布をガウス分布化することを図っている。
特開平9-265000号公報(0005段乃至0009段、図1) 特開2006-351339号公報(0011段乃至0013段、図3)
 従来は、患者の患部におけるビームのスポット形状が、回転ガントリの回転角度により変化することを防止するために、回転ガントリの偏向電磁石に入力される前に、エミッタンスの対称化を図るエミッタンス調整手段や散乱体及び下流側電磁石を荷電粒子ビーム輸送系に設けていた。従来の荷電粒子ビーム輸送系はエミッタンス調整手段や、散乱体及び下流側電磁石を追加していたので、荷電粒子ビーム輸送系が大型化する問題があった。
 本発明は上記のような課題を解決するためになされたものであり、加速器から遅い取り出しで出射したときに発生するエミッタンスの違いを荷電粒子ビーム輸送系で吸収し、アイソセンタで回転依存性の少ないビームサイズを実現することを目的とする。
 本発明に係る荷電粒子ビーム輸送系は、加速器から出射された荷電粒子ビームが、加速器の出射位置における荷電粒子ビームの進行方向に垂直な面内で、加速器の周回軌道面内の方向をx方向、荷電粒子ビームの進行方向に垂直な面内でx方向と直交する方向をy方向とした場合、x方向のエミッタンスとy方向のエミッタンスが異なる荷電粒子ビームであり、アイソセンタを中心に回転可能な回転ガントリに搭載された粒子線照射装置に荷電粒子ビームを輸送する荷電粒子ビーム輸送系であって、回転ガントリに搭載されると共に、回転ガントリのガントリ回転軸の周りを回転する回転偏向部と、加速器から回転偏向部までの固定輸送部と、を備え、回転偏向部は複数の偏向電磁石を有する。本発明に係る荷電粒子ビーム輸送系によれば、ガントリ回転軸方向から見た場合において、粒子線照射装置のビーム中心線と、回転偏向部の入口における荷電粒子ビームのy方向の軸との角度をガントリ角度とし、加速器の出射位置におけるx方向のエミッタンスとy方向のエミッタンスとが分離されてそれぞれのエミッタンスが保たれるように荷電粒子ビームが輸送される回転ガントリの角度をガントリ基準角度とし、ガントリ角度がガントリ基準角度である場合に、荷電粒子ビームの位相空間分布における位相が回転偏向部の入口からアイソセンタまで進む変化分を第1の位相進みとし、ガントリ角度がガントリ基準角度から90°回転した角度である場合に、荷電粒子ビームの位相空間分布における位相が回転偏向部の入口からアイソセンタまで進む変化分を第2の位相進みとし、固定輸送部は、回転偏向部の入口における荷電粒子ビームの位相空間分布の位相が、第1の位相進み及び第2の位相進みの平均値に基づいた演算を基に決定された位相となるようにしたことを特徴とする。
 本発明に係る荷電粒子ビーム輸送系によれば、固定輸送部が、回転偏向部の入口における荷電粒子ビームの位相空間分布を、その位相が第1の位相進み及び第2の位相進みの平均値に基づいた演算を基に決定された位相である等の所定の位相空間分布となるようにしたので、荷電粒子ビーム輸送系の回転偏向部の入り口にて、遅い取り出しの特徴であるx方向エミッタンスが非常に小さいようなビームに対しても、アイソセンタで回転依存性の少ないビームサイズを実現することができる。
本発明による粒子線治療装置の概略構成図である。 図1の粒子線照射装置の構成を示す図である。 本発明による荷電粒子ビーム輸送系の回転偏向部の構成を示す図である。 図3の回転偏向部を左側から見た側面図である。 本発明の実施の形態1による回転偏向部からアイソセンタにおける荷電粒子のベータトロン関数を示す図である。 本発明の実施の形態1による回転偏向部からアイソセンタにおける荷電粒子の運動量分散関数を示す図である。 本発明の実施の形態1による回転偏向部の入口における荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度が0°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度が30°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度が60°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度が90°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度が120°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度が150°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度が180°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度が210°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態1による回転偏向部の角度とアイソセンタにおけるビームサイズの関係を示す図である。 本発明による荷電粒子ビーム輸送系の固定輸送部のx方向の荷電粒子の位相空間分布を説明する図である。 本発明の実施の形態2による回転偏向部からアイソセンタにおける荷電粒子のベータトロン関数を示す図である。 本発明の実施の形態2による回転偏向部からアイソセンタにおける荷電粒子の運動量分散関数を示す図である。 本発明の実施の形態2による回転偏向部の入口における荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度が0°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度が30°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度が60°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度が90°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度が120°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度が150°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度が180°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度が210°の場合におけるアイソセンタでの荷電粒子の位相空間分布を示す図である。 本発明の実施の形態2による回転偏向部の角度とアイソセンタにおけるビームサイズの関係を示す図である。 本発明の実施の形態2による回転偏向部の角度が0°の場合におけるビームサイズを示す図である。 本発明の実施の形態2による回転偏向部の角度が0°の場合におけるx方向ビームサイズを示す図である。 本発明の実施の形態2による回転偏向部の角度が0°の場合におけるy方向ビームサイズを示す図である。 本発明の実施の形態2による回転偏向部の角度が30°の場合におけるビームサイズを示す図である。 本発明の実施の形態2による回転偏向部の角度が60°の場合におけるビームサイズを示す図である。 本発明の実施の形態2による回転偏向部の角度が90°の場合におけるビームサイズを示す図である。 本発明の実施の形態2による回転偏向部の角度が120°の場合におけるビームサイズを示す図である。 本発明の実施の形態2による回転偏向部の角度が150°の場合におけるビームサイズを示す図である。 位相演算の概念を説明する規格化位相空間の図である。 位相演算の概念を説明する規格化位相空間の図である。 位相演算の概念を説明する規格化位相空間の図である。
実施の形態1.
 図1は本発明による粒子線治療装置の概略構成図であり、図2は本発明による粒子線照射装置の構成を示す図である。図3は本発明による荷電粒子ビーム輸送系の回転偏向部の構成を示す図であり、図4は図3の回転偏向部を左側から見た側面図である。図1において、粒子線治療装置51は、ビーム発生装置52と、荷電粒子ビーム輸送系59と、粒子線照射装置58a、58bとを備える。ビーム発生装置52は、イオン源(図示せず)と、前段加速器53と、加速器54とを有する。粒子線照射装置58bは回転ガントリ(図示せず)に設置される。粒子線照射装置58aは回転ガントリを有しない治療室に設置される。
 荷電粒子ビーム輸送系59の役割は加速器54と粒子線照射装置58a、58bの連絡にある。加速器54と粒子線照射装置58a、58bとは真空ダクト57を介して接続され、荷電粒子ビームは真空ダクト57を通過する。荷電粒子ビーム輸送系59の一部は回転ガントリ(図示せず)に設置され、回転ガントリに設置される回転偏向部60は、複数の偏向電磁石10a、10b、10c、及び複数の四極電磁石11a、11b、11c、11dを有する。回転偏向部60は、回転ガントリの回転に伴い回転する。荷電粒子ビーム輸送系59の回転偏向部60以外の部分である固定輸送部61は、複数の偏向電磁石48a、48b、48c、48d、48e、48f、48g、48h、48i、及び複数の四極電磁石49a、49b、49c、49d、49e、49f、49g、49h、49i、49jを有する。固定輸送部61の偏向電磁石の符号は、総括的に48を用い、区別して説明する場合に48a乃至48iを用いる。固定輸送部61の四極電磁石の符号は、総括的に49を用い、区別して説明する場合に49a乃至49jを用いる。回転偏向部60の偏向電磁石の符号は、総括的に10を用い、区別して説明する場合に10a、10b、10cを用いる。
 イオン源で発生した陽子線等の粒子線である荷電粒子ビームは、前段加速器53で加速され、入射装置46から加速器54に入射される。ここでは、加速器54は、シンクロトロンを例に説明する。荷電粒子ビームは、所定のエネルギーまで加速される。加速器54の出射装置47から出射された荷電粒子ビームは、荷電粒子ビーム輸送系59を経て粒子線照射装置58a、58bに輸送される。粒子線照射装置58a、58bは荷電粒子ビームを照射対象45(図2参照)に照射する。粒子線照射装置の符号は、総括的に58を用い、区別して説明する場合に58a、58bを用いる。加速器54から出射された荷電粒子ビーム31は、加速器54の出射位置における荷電粒子ビーム31の進行方向(s方向)に垂直な面内で、加速器54の周回軌道面内の方向をx方向、荷電粒子ビーム31の進行方向に垂直な面内でx方向と直交する方向をy方向とした場合、x方向のエミッタンスが小さく、y方向のエミッタンスが大きな荷電粒子ビームとなっている。加速器54の出射位置における荷電粒子ビーム31のビーム座標系は、ビーム座標系73である。荷電粒子ビーム輸送系59の回転偏向部60の入口における荷電粒子ビーム31のビーム座標系は、ビーム座標系74である。なお、ビーム座標系のs方向の軸をs軸、ビーム座標系のx方向の軸をx軸、ビーム座標系のy方向の軸をy軸と呼ぶことにする。
 ビーム発生装置52で発生され、所定のエネルギーまで加速された荷電粒子ビーム31は、荷電粒子ビーム輸送系59を経由し、粒子線照射装置58へと導かれる。図2において、粒子線照射装置58は、荷電粒子ビーム31に垂直な方向であるX向及びY方向に荷電粒子ビーム31を走査するX方向走査電磁石32及びY方向走査電磁石33と、位置モニタ34と、線量モニタ35と、線量データ変換器36と、ビームデータ処理装置41と、走査電磁石電源37と、粒子線照射装置58を制御する照射管理装置38とを備える。照射管理装置38は、照射制御計算機39と照射制御装置40とを備える。線量データ変換器36は、トリガ生成部42と、スポットカウンタ43と、スポット間カウンタ44とを備える。なお、荷電粒子ビーム31の進行方向は-Z方向である。
 X方向走査電磁石32は荷電粒子ビーム31をX方向に走査する走査電磁石であり、Y方向走査電磁石33は荷電粒子ビーム31をY方向に走査する走査電磁石である。位置モニタ34は、X方向走査電磁石32及びY方向走査電磁石33で走査された荷電粒子ビーム31が通過するビームにおける通過位置(重心位置)やサイズを演算するためのビーム情報を検出する。ビームデータ処理装置41は、位置モニタ34が検出した複数のアナログ信号(ビーム情報)からなるビーム情報に基づいて荷電粒子ビーム31の通過位置(重心位置)やサイズを演算する。また、ビームデータ処理装置41は、荷電粒子ビーム31の位置異常やサイズ異常を示す異常検出信号を生成し、この異常検出信号を照射管理装置38に出力する。
 線量モニタ35は、荷電粒子ビーム31の線量を検出する。照射管理装置38は、図示しない治療計画装置で作成された治療計画データに基づいて、照射対象45における荷電粒子ビーム31の照射位置を制御し、線量モニタ35で測定され、線量データ変換器36でデジタルデータに変換された線量が目標線量に達すると荷電粒子ビーム31を停止する。走査電磁石電源37は、照射管理装置38から出力されたX方向走査電磁石32及びY方向走査電磁石33への制御入力(指令)に基づいてX方向走査電磁石32及びY方向走査電磁石33の設定電流を変化させる。
 ここでは、粒子線照射装置58のスキャニング照射方式を、ハイブリッドスキャニング照射方式(ビーム照射位置(スポット)を変更する場合にビームを停止しない方式)、具体的には、荷電粒子ビーム31の照射位置を変えるときに荷電粒子ビーム31を停止させないラスタースキャニング照射方式のように行い、スポットスキャニング照射方式のようにビーム照射位置がスポット位置間を次々と移動していく方式として説明する。スポットカウンタ43は、荷電粒子ビーム31のビーム照射位置が停留している間の照射線量を計測するものである。スポット間カウンタ44は、荷電粒子ビーム31のビーム照射位置が移動している間の照射線量を計測するものである。トリガ生成部42は、ビーム照射位置における荷電粒子ビーム31の線量が目標照射線量に達した場合に、線量満了信号を生成するものである。
 図3は、ガントリ回転軸15の周りを回転する回転角度(以降、適宜、ガントリ角度とも称する。)が0°と180°における荷電粒子ビーム輸送系59の回転偏向部60を示している。実線で示したものは回転角度が0°の場合であり、破線で示したものは回転角度が180°の場合である。破線で示したものに対する符号は、実線のものと区別するために異なるものを付した。回転偏向部60において、荷電粒子ビーム31は次に示す順に各電磁石を通過して粒子線照射装置58に至る。荷電粒子ビーム31の電磁石通過順は、偏向電磁石10a、四極電磁石11a、11b、偏向電磁石10b、四極電磁石11c、11d、偏向電磁石10cである。回転角度が180°における回転偏向部60において、荷電粒子ビーム31は次に示す順に各電磁石を通過して粒子線照射装置68(図3及び図4)に至る。図3及び図4における粒子線照射装置の符号は、回転角度を区別するために、58とは異なる符号も用いた。回転角度が90°、180°、270°における粒子線照射装置の符号は、それぞれ67、68、69である。荷電粒子ビーム31の電磁石通過順は、偏向電磁石22a、四極電磁石23a、23b、偏向電磁石22b、四極電磁石23c、23d、偏向電磁石22cである。アイソセンタICは、回転偏向部60が回転したときの荷電粒子ビーム31のビーム中心線14がそれぞれ交わる点である。回転角度が0°、90°、180°、270°における粒子線照射装置を通過した荷電粒子ビームのビーム座標系は、それぞれ、ビーム座標系75、ビーム座標系77、ビーム座標系76、ビーム座標系78である。
 ここで、ガントリ角度の基準を考える。ガントリ回転軸15の周りを回転する回転角度であるガントリ角度は、粒子線照射装置58における荷電粒子ビーム31のビーム中心線14と、回転偏向部60の入口における荷電粒子ビーム31のビーム座標系74におけるy方向の軸との角度と定義する。ビーム中心線14と回転偏向部60の入口における荷電粒子ビーム31のy方向の軸(y軸)とが重なる場合のガントリ角度をガントリ基準角度とする。ガントリ角度がガントリ基準角度の場合は、加速器54の出射位置におけるx方向のエミッタンスとy方向のエミッタンスとが分離されてそれぞれのエミッタンスが保たれるように荷電粒子ビーム31が荷電粒子ビーム輸送系59で輸送されることになる。ガントリ基準角度は、任意のガントリ角度で構わないが、ここでは、0°とする。
 図4は、図3の回転偏向部を左側から見た側面図であるが、回転角度が90°と270°の場合も追加したものである。ビーム中心線16は、回転角度0°の粒子線照射装置58と回転角度180°の粒子線照射装置68を通過するビーム中心線14を90°回転したものであり、回転角度90°の粒子線照射装置67と回転角度270°の粒子線照射装置69を通過する荷電粒子ビーム31のビーム中心線である。図4において、回転角度0°の偏向電磁石10a、10cは、回転角度が90°、270°の場合に、その符号を変えている。回転角度が90°の場合は偏向電磁石20a、20cであり、回転角度が270°の場合は偏向電磁石24a、24cである。
 荷電粒子ビーム輸送系59のビーム設計とは、偏向電磁石や四極電磁石の磁石配列と、それぞれの四極電磁石の強度を調整することであり、下記式(1)に示す輸送行列を決定することにある。
Figure JPOXMLDOC01-appb-M000001
 荷電粒子ビーム輸送系59の回転偏向部60及び固定輸送部61のビーム設計も、荷電粒子ビーム輸送系59のビーム設計と同様である。荷電粒子ビーム輸送系59の回転偏向部60を例に説明する。輸送行列は、入口の3次元座標とその運動量と出口の3次元座標と運動量の関係を示すものである。右辺の(x1、y1、x’1、y’1、t1、pt1)は、回転ガントリ入口における荷電粒子の位置(x1、y1)と、s方向の傾き(x’1、y’1)と、時間t1と、荷電粒子の運動量pt1である。左辺の(x2、y2、x’2、y’2、t2、pt2)は、回転偏向部60の出口における荷電粒子の位置(x2、y2)と、s方向の傾き(x’2、y’2)と、時間t2と、荷電粒子の運動量pt2である。プライム‘は、sによる微分、すなわちd/dsを表す。
 荷電粒子ビーム輸送系59の回転偏向部60のビーム設計は、少ない磁石数やパラメータで上記輸送行列を決定することが課題であり、従来では、上記輸送行列を変形して、ベータトロン関数β、運動量分散関数ηという概念が用いられてきた。運動量分散関数ηは、Δpt/ptである。本発明では、運動量分散関数ηに関しては、従来法と同様に、回転偏向部60の入口において、ηx=0、η’x=0、ηy=0、η’y=0となるビームが来ることを想定し、回転偏向部60の出口でηx=0、η’x=0、ηy=0、η’y=0を満たすように設計する。この制約を課すことによって、回転偏向部60のビーム設計としては、下記式(2)のみを考えれば良い。
Figure JPOXMLDOC01-appb-M000002
 式(2)では、時間t1、t2、及び荷電粒子の運動量pt1、pt2を使用しないので、ガントリの回転に伴って回転偏向部60が回転しても、回転偏向部60の出口において、運動量の差によるビームサイズの変化や、ビームサイズの変化率の差が発生しない。
 次に、ガントリ角度が0°(または180°)のときと90°(または270°)の時に着目する。回転偏向部60の偏向電磁石が、垂直もしくは水平のみの偏向電磁石から構成される場合には、式(2)は式(3)のようになり、x方向とy方向のカップリングがなくなる。
Figure JPOXMLDOC01-appb-M000003
 さらに、回転偏向部60の入口において、ビームの位相空間上の特徴量であるベータトロン関数β、αをx方向とy方向で一致させることで、すなわち、βx=βy、αx=αyとすることで、ベータトロン関数βの差によるビームサイズの回転依存性はなくなる。なお、sに対するβの関数をベータトロン関数と呼ぶ。
 実施の形態1の回転偏向部60のビーム設計において、回転偏向部60の入口及び出口の条件を下記のようにする。
 回転偏向部60の入口の条件は、βx=βy、αx=αy、ηx=0、η’x=0、ηy=0、η’y=0である。
 回転偏向部60の出口の条件は、βx=βy、αx=αy、ηx=0、η’x=0、ηy=0、η’y=0である。
 回転偏向部60の入口及び出口の条件を上記のようにすることで、すなわち、できる限り簡易化することで、少ない調整パラメータで回転偏向部60のビーム設計を行うことができる。式(3)における調整パラメータは、m11、m12、m21、m22、m33、m34、m43、m44である。
 図3の回転偏向部60のベータトロン関数βと運動量分散関数ηを図5及び図6に示す。図5は、本発明の実施の形態1による回転偏向部からアイソセンタICにおける荷電粒子のベータトロン関数を示す図である。横軸はs方向の距離であり、縦軸はベータトロン関数βである。なお、図5においてベータトロン関数αは記載していない。図5(a)はx方向のベータトロン関数βであり、図5(b)はy方向のベータトロン関数βである。図6は、本発明の実施の形態1による回転偏向部からアイソセンタICにおける荷電粒子の運動量分散関数を示す図である。横軸はs方向の距離であり、縦軸は運動量分散関数ηである。図6(a)はx方向の運動量分散関数ηであり、図5(b)はy方向の運動量分散関数ηである。図5(a)、図5(b)、図6(a)、図6(b)のそれぞれにおいて、回転偏向部60の入口は、s=0mであり、回転偏向部60の出口は、左端のs=14m付近である。実施の形態1の例では、4個の四極電磁石11a、11b、11c、11dにより、荷電粒子ビーム31が回転偏向部60を通過して、アイソセンタICのビーム特性が良好となる輸送行列を実現している。
 スキャニング照射方式ではない従来照射法であれば、アイソセンタICにおけるx方向とy方向のビームサイズは異なっていても、実際の照射では、散乱体やワブラー電磁石によってビームサイズは広がり、回転偏向部の出口において元々のビームサイズの大きさによらない程度に荷電粒子が散乱させて照射するため、ビームサイズの違いは問題にならなかった。しかしながら、スキャニング照射方式では、出来る限り小さいビームサイズが要求されるため、散乱はなるべく少ないことが望まれる。散乱があまり期待できない場合には、加速器54や荷電粒子ビーム輸送系59のビーム設計として、x方向とy方向のビームサイズはそろっており、かつ、x方向とy方向のビームサイズは回転ガントリの角度依存性がないことが望まれる。
 x方向ビームサイズσx及びy方向ビームサイズσyは、一般に式(4)のように表される。εx、εyは、それぞれx方向のエミッタンス及びy方向のエミッタンスである。エミッタンスεx、εyは、位相空間上における面積である。x方向の位相空間は、図7(a)に示すように表され、横軸をx方向の距離xとし、縦軸をx方向の傾きx’として表された空間である。同様に、y向の位相空間は、図7(b)に示すように表され、横軸をy方向の距離yとし、縦軸をy方向の傾きy’として表された空間である。
Figure JPOXMLDOC01-appb-M000004
 ベータトロン関数βx、βyは、回転ガントリのビーム設計によって変えることのできる量(輸送行列に依存する量)である。これに対して、エミッタンスεx、εyは、加速器54の出射によって決まる量であり、荷電粒子ビーム輸送系59にて荷電粒子の輸送中に散乱体などが存在しない場合には、一定である。また、エミッタンスεx、εyは、加速器54からの取り出しメカニズムに依存し、加速器54から遅い取り出しで出射した場合は、エミッタンスεx、εyは異なる値になる。エミッタンスεxとエミッタンスεyを比べると、エミッタンスεyのほうが遥かに大きく、エミッタンスεyはエミッタンスεxの2乃至10倍以上になる。本発明では、このエミッタンスεx、εyの違いによるビームサイズの違いを、荷電粒子ビーム輸送系59のビーム設計で極力無くすようにした。
 回転偏向部60に散乱体と四極電磁石を組み合わせた構成が存在したり、回転偏向部60の四極電磁石が少なくても回転ガントリの回転角度ごとに各四極電磁石の強度を変更したりすれば、ビームサイズの違いを極力無くことは可能である。しかし、四極電磁石の強度を変えずとも回転角度ごとのビームサイズの違いが無いようにすることで、回転偏向部60を安価に実現できることは重要である。このような荷電粒子ビーム輸送系59のビーム設計方法は、既存の施設においても適用でき、既存の施設においても荷電粒子ビーム輸送系59の回転偏向部60を実現できることが望ましい。式(3)から、x、x’を対象にした2×2の輸送行列と、y、y’を対象にした2×2の輸送行列に分けて、荷電粒子ビーム輸送系59のビーム設計を実行できることが分かる。
 次に、エミッタンスが小さいことによる影響について議論する。荷電粒子ビーム31を加速器54から遅く取り出すビーム取り出し方式(遅いビームビーム取り出し方式)は、加速器54を周回する荷電粒子ビーム31を、加速器54の外周側(x方向)に移動させて、荷電粒子ビーム31の外周側のみを取り出すようにして、周回する荷電粒子ビーム31を長期間に渡って少しずつ加速器54から荷電粒子ビーム輸送系59に出射するものである。図7は、本発明の実施の形態1による回転偏向部の入口における荷電粒子の位相空間分布を示す図である。図7(a)はx方向の位相空間分布であり、図7(b)はy方向の位相空間分布である。回転偏向部60の入口でのx方向の位相空間分布は、加速器54からの出射固有の特性より、x方向の分布は制限されており、楕円というよりは線に近い分布となっている。これに対して、y方向の位相空間分布は、y方向に制限されることがないため楕円となっている。位相空間上の面積がエミッタンスであるから、x方向とy方向では大きな違いがある。
 本発明の実施の形態1による荷電粒子ビーム輸送系59のビーム設計は、荷電粒子ビーム31が回転偏向部60を通過し、アイソセンタICのビーム特性が良好となる回転偏向部60のビーム設計と、アイソセンタICで回転依存性の少ないビームサイズを実現するための回転偏向部60に入力する荷電粒子の位相空間分布を実現する固定輸送部61のビーム設計とに分かれている。固定輸送部61のビーム設計の際に、アイソセンタICで回転依存性の少ないビームサイズを実現するための回転偏向部60に入力する荷電粒子の位相空間分布を決定する。
 回転偏向部60に入力する荷電粒子の位相空間分布を決定する際に、x方向の位相空間分布は線と近似し、y方向の位相空間分布は楕円と近似する。x方向の位相空間分布は線で近似するので、面積ではなく位相が重要である。回転偏向部60の入口におけるx方向の位相空間分布の位相は、回転ガントリの位相関係に注目して決定する。具体的には、回転偏向部60の入口におけるx方向の位相空間分布の位相φixは、ガントリ角度が90°異なる2つのガントリ角度θ1、θ2のそれぞれの回転偏向部60の入口から出口までの位相進み(フェイズアドバンス)の平均値に基づいた演算を基に決定された位相であり、例えばガントリ角度0°及び90°の入口から出口までの位相進み(フェイズアドバンス)の小数部のそれぞれを平均した平均値にマイナス(-)を付けたものである。
 図3に示した回転偏向部60の電磁石配列であり、図5及び図6に示したベータトロン関数及び運動量分散関数となるような各電磁石の強度となるように回転偏向部60の輸送行列を用いて回転偏向部60が設計された場合を考える。ガントリ角度が0°の場合(図4の実線で記載した回転偏向部60の配置)の場合の位相進みφ1は、1.382(497.52°)であり、回転偏向部60の出口では497.52°回転している。ガントリ角度が90°の場合(図4の偏向電磁石20a、20c、粒子線照射装置67で記載した回転偏向部60の配置)の場合の位相進みφ2は、0.522であり、回転偏向部60の出口では93.96°回転している。ガントリ角度が180°における荷電粒子ビーム31の位相空間分布は、ガントリ角度が0°における荷電粒子ビーム31の位相空間分布と理論上一致する。ガントリ角度360°は一周であり、180°の2倍であるから、一周分は無視して、すなわち180°の位相進みで基準化して、位相進みφ1は、497.52°-360°=137.52°とする。
 ガントリ角度の変化により、x方向のビームサイズが変化しないようにするために、位相進みφ2である93.96°と、位相進みφ1である137.52°の中間の値(平均値)である114.465°を利用する。すなわち、回転偏向部60の入口におけるx方向の位相空間分布の位相を、-114.465°にしておけば、回転偏向部60の出口やアイソセンタICにおいて、図8及び図11のように、回転偏向部60のガントリ角度、すなわち回転ガントリの設置角度が0°の場合でも、90°の場合でもどちらでもほぼ同じビームサイズが実現できる。平均値にマイナス(-)を付ける理由は、後述する。
 X方向走査電磁石32やY方向走査電磁石33で荷電粒子ビーム31を走査しない場合は、アイソセンタICおけるビームサイズは、回転偏向部60の出口におけるビームサイズとほぼ等しい。また、X方向走査電磁石32やY方向走査電磁石33で荷電粒子ビーム31を走査する場合は、荷電粒子ビーム31を走査しない場合に比べて、アイソセンタICおけるビームサイズは、多少変化するが、この場合であっても、アイソセンタICおけるビームサイズは、回転偏向部60の出口におけるビームサイズとほぼ等しい。
 図8乃至図15は、本発明の実施の形態1による回転偏向部の角度がそれぞれ0°、30°、60°、90°、120°、150°、180°、210°の場合におけるアイソセンタICでの荷電粒子の位相空間分布を示す図である。各図の(a)はx方向の位相空間分布であり、各図の(b)はy方向の位相空間分布である。各図の(a)において、横軸はx方向の距離xであり、縦軸はx方向の傾きx’である。各図の(b)において、横軸はy方向の距離yであり、縦軸はy方向の傾きy’である。図中のrmsは、二乗平均平方根(root mean square)の意味である。図中に示したビームサイズ(rms)は、標準偏差σに相当する大きさである。ビームサイズ(rms)をビームサイズ1σと呼ぶことにする。
 回転偏向部60の角度が0°の場合のx方向及びy方向のビームサイズ1σは、図8に示すように、それぞれ1.29mm及び1.29mmである。回転偏向部の角度が30°の場合のx方向及びy方向のビームサイズ1σは、図9に示すように、それぞれ1.30mm及び1.30mmである。回転偏向部の角度が60°の場合のx方向及びy方向のビームサイズ1σは、図10に示すように、それぞれ1.30mm及び1.30mmである。回転偏向部の角度が90°の場合のx方向及びy方向のビームサイズ1σは、図11に示すように、それぞれ1.30mm及び1.30mmである。
 回転偏向部の角度が120°の場合のx方向及びy方向のビームサイズ1σは、図12に示すように、それぞれ1.29mm及び1.30mmである。回転偏向部の角度が150°の場合のx方向及びy方向のビームサイズ1σは、図13に示すように、それぞれ1.29mm及び1.29mmである。回転偏向部の角度が180°の場合のx方向及びy方向のビームサイズ1σは、図14に示すように、それぞれ1.29mm及び1.30mmである。回転偏向部の角度が210°の場合のx方向及びy方向のビームサイズ1σは、図15に示すように、それぞれ1.30mm及び1.30mmである。
 図8乃至図15に示すように、x方向及びy方向の位相空間分布の形状や位相は異なるが、x方向及びy方向のビームサイズ1σは、ほぼ等しくなっている。位相進みφ1及びφ2の平均値にマイナス(-)を付ける理由を説明する。回転偏向部60の入口における位相の選択は重要である。入口における位相の選び方が悪いと、出口における位相空間上でビームの位相空間分布が鉛直方向となってしまうため、ビームサイズが小さくなりすぎて所望のビームサイズにすることができず、望ましくない。位相進みφ1及びφ2の平均値にマイナス(-)を付けると、回転偏向部60の出口(アイソセンタ)でガントリの位相進み分がキャンセルされるからである。回転ガントリでは位相が入口から出口まで進むため、入口では位相が進む前の状態が望ましい。例えば、回転ガントリで60°進むとすると、-60°に設定しておけば、回転偏向部60の出口(アイソセンタ)で回転偏向部60中の位相進み分がキャンセルされることになる。
 図16は、本発明の実施の形態1による回転偏向部の角度とアイソセンタICにおけるビームサイズの関係を示す図である。横軸は回転偏向部60のガントリ角度θであり、縦軸は標準偏差σの3倍に相当する大きさであるビームサイズ3σである。実線で示した特性71はx方向特性であり、破線で示した特性72はy方向特性である。実施の形態1の回転偏向部60を上記のようにすることで、回転偏向部60の回転角度、すなわち回転ガントリの回転角度に依存せずに、ほぼ一定のビームサイズを実現できる。
 実施の形態1の荷電粒子ビーム輸送系59は、従来の設計における回転偏向部60の入口の条件である、βx=βy、αx=αy、ηx=0、η’x=0、ηy=0、η’y=0に加えて、入口におけるx方向の位相空間分布の位相φixを追加した入口の条件を用いて実現する。回転偏向部60の入口におけるx方向の位相空間分布の位相φixは、ガントリ角度が90°異なる2つのガントリ角度θ1、θ2のそれぞれの回転偏向部60の入口から出口までの位相進み(フェイズアドバンス)の小数部のそれぞれを平均した平均値にマイナス(-)を付けたものである。回転偏向部60の入口におけるx方向の位相空間分布の位相φixを設計した所定値にするために、例えば、荷電粒子ビーム輸送系59の固定輸送部61に四極電磁石を追加するか、もしくは四極電磁石の位置をずらす、もしくは励磁量の設定値を変更する。
 次に、本発明による荷電粒子ビーム輸送系59の固定輸送部61のビーム設計について説明する。荷電粒子ビーム輸送系59の固定輸送部61は、次の3つの条件を満たすように設計する。条件1は、荷電粒子ビーム31が固定輸送部61を通過することである。条件2は、運動量分散関数η及びy方向の荷電粒子の位相空間分布が従来の条件を満たすことである。条件3は、固定輸送部61の入口におけるx方向の荷電粒子の位相空間分布が、アイソセンタICで回転依存性の少ないビームサイズを実現するための回転偏向部60に入力する所定の位相空間分布になることである。実施の形態1の荷電粒子ビーム輸送系59の回転偏向部60における所定の位相空間分布は、前述した回転偏向部60の入口におけるx方向の位相空間分布の位相φixであり、ガントリ角度が90°異なる2つのガントリ角度θ1、θ2のそれぞれの回転偏向部60の入口から出口までの位相進み(フェイズアドバンス)から半整数部を引き算した小数部のそれぞれを平均した平均値にマイナス(-)を付けたものである。
 ビーム物理では一般に、横軸をx、縦軸をx’とした位相空間上の粒子の運動を考える。ビームの進行とともに位相空間上の粒子分布はベータトロン関数(β(s))にあわせて伸縮しつつ、位相進み(フェイズアドバンス)にあわせて回転する。βが同じ場合、ビームサイズを左右するのは粒子分布の位相である。位相のみに注目するために、ベータトロン関数による変形を規格化した規格化位相空間を考える。この規格化位相空間では、xの取りうる範囲は±1であり、x’の取りうる範囲は±1である。
 この規格化位相空間では大きさは変化せず、ビームは輸送されるごとに反時計周りに回転を行う。例えば、回転ガントリが0°のときのx方向位相進み(フェイズアドバンス)をν0=2.9とし、回転ガントリが90°のときのx方向の位相進み(フェイズアドバンス)をν90=0.8とする。図38~図40を用いて、位相演算の概念を説明する。
 図38は、回転ガントリ入り口の位相を0.15として、回転ガントリが0°のときのx方向位相進みが、ν0=2.9へと進んだ場合の規格化位相空間である。図39は、回転ガントリ入り口の位相を0.15として、回転ガントリが90°のときのx方向位相進みが、ν90=0.8へと進んだ場合の規格化位相空間である。図40は、回転ガントリ入り口の位相を0.15として、x方向位相進みがν0=2.9である場合とx方向位相進みがν90=0.8である場合とを比較するための図である。図38~図40において、位相線101は回転ガントリ入口の位相であり、位相線102は回転ガントリが0°のときのアイソセンタICの位相であり、位相線103は回転ガントリが90°のときのアイソセンタICの位相である。
 回転ガントリ入り口の位相は、-(2.9+0.8)/2=-1.85として計算した。位相は、半周回しても同じため、回転ガントリ入り口の位相は、-1.85±n/2=0.15±m/2とすることができる。なお、n、mは自然数である。したがって、以下の式に変形できる。
 -(2.9+0.8)/2=-0.85=0.15-2/2
これは、2.9の小数部0.9と0.8の小数部0.8の平均値を用いた場合と結果が同じということを示している。
 次に、半整数部について考える。2.9の半整数部を除いた部分は、0.4である。0.8の半整数部を除いた部分は、0.3である。この場合も上記と同様に、下記の式がなりたち、結果は同じである。
 -(0.4+0.3)/2=-0.35=0.15-1/2
 角度で表現すると、54°±180°を反時計周りにした角度である。
 ガントリ角度が0°のときと90°のときで、ほぼアイソセンタICのx方向の広がり(図38~図40では投影長さ)がほぼ同一にできている。0°と90°の間の角度では、x方向とy方向の結合があるため単純ではないが、ほぼこの間にはいるため、ビームサイズはガントリ角度に依存せず一定化可能になる。
 上述した条件1は当然なので、条件2及び3を具体的に説明する。図17は、本発明の荷電粒子ビーム輸送系のx方向の荷電粒子の位相空間分布を説明する図である。図17(a)は固定輸送部61の出口においてx方向の位相空間分布がx方向位相分布98になった例であり、図17(b)は固定輸送部61の出口においてx方向の位相空間分布がx方向位相分布99になった例である。x方向位相分布97は、加速器54の出口における位相空間分布である。x方向位相分布98、99は、固定輸送部61の出口、すなわち回転偏向部60の入口における位相空間分布である。x方向位相分布97、98、99は、線状分布として示す。x方向位相分布97の端点p1の座標は(x1、x’1)であり、x方向位相分布97の端点p2の座標は(-x1、-x’1)である。位相Ψ1はx方向位相分布97における位相であり、位相Ψ2はx方向位相分布98やx方向位相分布99における位相である。
 x方向位相分布98、99は、いずれも回転偏向部60の入口における端点p1が端点p11(x2、x’2)になり、回転偏向部60の入口における端点p2が端点p12(-x2、-x’2)になったものである。x方向位相分布98とx方向位相分布99との違いは、端点p11が存在する位相空間上の象限が逆になっていることである。
 条件2は、回転偏向部60の入口でβy、αy、ηx、η’x、ηy、η’yが従来どおりの条件を満たすことであり、すなわち、ηx、η’x、ηy、η’yが全て0にあり、βy、αyがアイソセンタICにおけるビーム条件を満たす値にすることである。この条件2を拘束条件1と呼ぶことにする。
 条件3は、固定輸送部61の出口(回転偏向部60の入口)において、端点p11(x2、x’2)及び端点p12(-x2、-x’2)を有するx方向位相分布98やx方向位相分布99となるように輸送行列Mを決定することである。この条件3を拘束条件2と呼ぶことにする。輸送行列Mは、荷電粒子ビーム輸送系59の固定輸送部61の入口及び出口のツイスパラメータと、位相進み(フェイズアドバンス)の関数なので、ツイスパラメータを用いて条件を表現すると式(5)のようになる。なお、ツイス(twiss)パラメータは、ベータトロン関数β、αをパラメータとする場合の呼び方である。
Figure JPOXMLDOC01-appb-M000005
 ここで、式(5)の計数A1、B1、C1は式(6)のように表される。
Figure JPOXMLDOC01-appb-M000006
 ただし、β1、α1は固定輸送部61の入口のツイスパラメータであり、β2、α2は固定輸送部61の出口のツイスパラメータであり、ΔΨは固定輸送部61の入口から出口の位相進み(フェイズアドバンス)であり、Ψ2-Ψ1に相当する。x1、x’1は加速器54の出射設計で決まるものであり、x2、x’2はビームサイズが角度依存しないようにする場合の回転ガントリ(回転偏向部60)の設計で決まる。
 本発明の実施の形態1による荷電粒子ビーム輸送系59では、回転偏向部60の入口において、x方向の位相空間分布の位相φixをガントリ角度が90°異なる2つのガントリ角度θ1、θ2のそれぞれの回転偏向部60の入口から出口までの位相進み(フェイズアドバンス)の小数部のそれぞれを平均した平均値にマイナス(-)を付けたものとなるようにする。すなわち、端点p11(x2、x’2)及び端点p12(-x2、-x’2)を有するx方向位相分布98(または99)の位相は、ガントリ角度が90°異なる2つのガントリ角度θ1、θ2のそれぞれの回転偏向部60の入口から出口までの位相進み(フェイズアドバンス)の小数部のそれぞれを平均した平均値にマイナス(-)を付けた位相φixとなるようにする。x方向位相分布98(または99)の位相Ψ2は、上記位相φixとなる。
 条件1乃至3を満たす荷電粒子ビーム輸送系59の固定輸送部61の設計フローは以下の通りである。
 ステップS1:従来どおり、ツイスパラメータを条件に用いて四極電磁石49の配置及び極性を決める。ここでは、四極電磁石49の配置を自由度に加えると,後の工程が長くなりすぎるので、四極電磁石49のベースを作る。こうすることで後の工程を短くすることができる。
 ステップS2:固定輸送部61の出口(回転偏向部60の入口)におけるx方向の位相空間分布に接するように任意の楕円を書き、ツイスパラメータを計算する。楕円の位相Ψ2も計算しておく。固定輸送部61の出口(回転偏向部60の入口)のビーム分布が重要なため、出口を基準にする。
 ステップS3:式(5)を用いて固定輸送部61の入口のツイスパラメータを計算する。このとき、ΔΨは任意の値を仮定し、入口の位相Ψ1はΨ2とΔΨから計算する。ツイスパラメータから荷電粒子ビーム31が固定輸送部61を通過できることを知るためには、ステップS2及び3を繰り返して,入口のツイスパラメータに対応する楕円(エミッタンスは当然出口の値)が入口のビーム分布を包含するようにする。毎回トラッキングを行い、ビームサイズを確認する場合はこの限りではない。
 ステップS4:固定輸送部61の入口と出口におけるツイスパラメータおよび拘束条件2を満たすように固定輸送部61を設計する。なお、従来の設計方法は、このステップS4の実行に相当する。
 ステップS5:ステップS4で適切な解が見つかるまでステップS2乃至ステップS4を繰り返す。解がない場合は、ステップS1に戻り、四極電磁石49の配置や極性変更を行う。
 なお、このように設計した設計結果が、従来の設計条件を満たしていることを確認する。この確認をすることにより、荷電粒子ビーム輸送系59の設置後の調整作業をスムーズに行うことができる。また、ステップS1で決めた四極電磁石49の配置及び極性を遵守する場合も、荷電粒子ビーム輸送系59の設置後の調整作業をスムーズに行うことができる。
 実施の形態1の荷電粒子ビーム輸送系59の回転偏向部60は、3つの偏向電磁石10a、10b、10cと4つの四極電磁石11a、11b、11c、11dの電磁石の強度を変更することなく、従来の設計における入口の条件に入口におけるx方向の位相空間分布の位相φixを追加した新たな入口の条件を適用することで、荷電粒子ビーム輸送系59の回転偏向部60の入口にてエミッタンスの非対称が存在しても、アイソセンタICで回転依存性の少ないビームサイズを実現することができる。
 以上のように実施の形態1による荷電粒子ビーム輸送系59によれば、加速器54から出射された荷電粒子ビーム31が、加速器54の出射位置における荷電粒子ビーム31の進行方向に垂直な面内で、加速器の周回軌道面内の方向をx方向、荷電粒子ビーム31の進行方向に垂直な面内でx方向と直交する方向をy方向とした場合、x方向のエミッタンスとy方向のエミッタンスが異なる荷電粒子ビームであり、アイソセンタICを中心に回転可能な回転ガントリに搭載された粒子線照射装置58に荷電粒子ビーム31を輸送する荷電粒子ビーム輸送系59であって、回転ガントリに搭載されると共に、回転ガントリのガントリ回転軸15の周りを回転する回転偏向部60と、加速器54から回転偏向部60までの固定輸送部61と、を備え、回転偏向部60は複数の偏向電磁石10を有する。実施の形態1による荷電粒子ビーム輸送系59によれば、ガントリ回転軸方向から見た場合において、粒子線照射装置58のビーム中心線14と、回転偏向部60の入口における荷電粒子ビーム31のy方向の軸との角度をガントリ角度とし、加速器54の出射位置におけるx方向のエミッタンスとy方向のエミッタンスとが分離されてそれぞれのエミッタンスが保たれるように荷電粒子ビーム31が輸送される回転ガントリの角度をガントリ基準角度とし、ガントリ角度がガントリ基準角度である場合に、荷電粒子ビーム31の位相空間分布における位相が回転偏向部60の入口からアイソセンタICまで進む変化分を第1の位相進みとし、ガントリ角度がガントリ基準角度から90°回転した角度である場合に、荷電粒子ビーム31の位相空間分布における位相が回転偏向部60の入口からアイソセンタICまで進む変化分を第2の位相進みとし、固定輸送部61は、回転偏向部60の入口における荷電粒子ビーム31の位相空間分布の位相が、第1の位相進み及び第2の位相進みの平均値に基づいた演算を基に決定された位相となるようにしたことを特徴とするので、固定輸送部61が、回転偏向部60の入口におけるx方向の位相空間分布の位相φixを所定の位相にでき、すなわち、回転偏向部60の入口における荷電粒子ビーム31の位相空間分布を、その位相が第1の位相進み及び第2の位相進みの平均値に基づいた演算を基に決定された位相である等の所定の位相空間分布にでき、荷電粒子ビーム輸送系59の回転偏向部60の入り口にて、遅い取り出しの特徴であるx方向エミッタンスが非常に小さいようなビームに対しても、アイソセンタICで回転依存性の少ないビームサイズを実現することができる。
 以上のように実施の形態1による粒子線治療装置51によれば、荷電粒子ビーム31を発生させ、この荷電粒子ビーム31を加速器54で加速させるビーム発生装置52と、加速器54により加速された荷電粒子ビーム31が、加速器54の出射位置における荷電粒子ビーム31の進行方向に垂直な面内で、加速器54の周回軌道面内の方向をx方向、荷電粒子ビーム31の進行方向に垂直な面内でx方向と直交する方向をy方向とした場合、x方向のエミッタンスとy方向のエミッタンスが異なる荷電粒子ビームであり、この荷電粒子ビーム31を輸送する荷電粒子ビーム輸送系59と、荷電粒子ビーム輸送系59で輸送された荷電粒子ビーム31を照射対象45に照射する粒子線照射装置58と、粒子線照射装置を搭載し、アイソセンタを中心に回転可能な回転ガントリと、を備え、荷電粒子ビーム輸送系59は、回転ガントリに搭載されると共に、回転ガントリのガントリ回転軸15の周りを回転する回転偏向部60と、加速器54から回転偏向部60までの固定輸送部61と、を備え、回転偏向部60は複数の偏向電磁石10を有し、ガントリ回転軸方向から見た場合において、粒子線照射装置58のビーム中心線14と、回転偏向部60の入口における荷電粒子ビーム31のy方向の軸との角度をガントリ角度とし、加速器54の出射位置におけるx方向のエミッタンスとy方向のエミッタンスとが分離されてそれぞれのエミッタンスが保たれるように荷電粒子ビーム31が輸送される回転ガントリの角度をガントリ基準角度とし、ガントリ角度がガントリ基準角度である場合に、荷電粒子ビーム31の位相空間分布における位相が回転偏向部60の入口からアイソセンタICまで進む変化分を第1の位相進みとし、ガントリ角度がガントリ基準角度から90°回転した角度である場合に、荷電粒子ビーム31の位相空間分布における位相が回転偏向部60の入口からアイソセンタICまで進む変化分を第2の位相進みとし、固定輸送部61は、回転偏向部60の入口における荷電粒子ビーム31の位相空間分布の位相が、第1の位相進み及び第2の位相進みの平均値に基づいた演算を基に決定された位相となるようにしたことを特徴とするので、固定輸送部61が、回転偏向部60の入口におけるx方向の位相空間分布の位相φixを所定の位相にでき、すなわち、回転偏向部60の入口における荷電粒子ビーム31の位相空間分布を、その位相が第1の位相進み及び第2の位相進みの平均値に基づいた演算を基に決定された位相である等の所定の位相空間分布にでき、荷電粒子ビーム輸送系59の回転偏向部60の入り口にて、遅い取り出しの特徴であるx方向エミッタンスが非常に小さいようなビームに対しても、アイソセンタICで回転依存性の少ないビームサイズを実現することができる。
実施の形態2.
 本発明の実施の形態2による荷電粒子ビーム輸送系59について説明する。実施の形態2では、x方向の位相空間分布を、端点を有する線状分布と近似し、式(3)の輸送行列の一部を取り出したマトリクスを用いて計算した端点を有する線状分布を、荷電粒子ビーム輸送系59の回転偏向部60の入口で実現するように、荷電粒子ビーム輸送系59の固定輸送部61を構成する。このように、回転偏向部60の入口における荷電粒子ビーム31のx方向の位相空間分布を制御することで、加速器54から遅い取り出しで出射したときに発生するエミッタンスの違いを荷電粒子ビーム輸送系59で吸収し、アイソセンタICで回転依存性の少ないビームサイズを実現することができる。
 本発明の実施の形態2による荷電粒子ビーム輸送系59の回転偏向部60のビーム設計について説明する。実施の形態1で説明したように、運動量分散関数ηに関しては、従来法と同様に、回転偏向部60の入口において、ηx=0、η’x=0、ηy=0、η’y=0となるビームが来ることを想定し、回転偏向部60の出口でηx=0、η’x=0、ηy=0、η’y=0を満たすように設計する。さらに、ガントリ角度が0°(または180°)のときと90°(または270°)の時に着目する。回転偏向部60の偏向電磁石が、固定輸送部61の偏向電磁石の平行面および垂直面内で偏向する偏向電磁石のみで構成される場合には、実施の形態1と同様に、x方向とy方向のカップリングがない、式(3)のような輸送行列が得られる。
 回転偏向部60の入口及び出口の条件において、実施の形態1とは異なり、αx=αyの制約は設けない。すなわち、実施の形態2の回転偏向部60のビーム設計において、回転偏向部60の入口及び出口の条件を下記のようにする。
 回転偏向部60の入口の条件は、βx=βy、ηx=0、η’x=0、ηy=0、η’y=0である。
 回転偏向部60の出口の条件は、ηx=0、η’x=0、ηy=0、η’y=0である。
 回転偏向部60の入口及び出口の条件を上記のようにすることで、すなわち、できる限り簡易化することで、少ない調整パラメータで回転偏向部60のビーム設計を行うことができる。式(3)における調整パラメータは、m11、m12、m21、m22、m33、m34、m43、m44である。
 本発明の実施の形態2による回転偏向部60(図3参照)のベータトロン関数βと運動量分散関数ηを図18及び図19に示す。図18は、本発明の実施の形態2による回転偏向部からアイソセンタICにおける荷電粒子のベータトロン関数を示す図である。横軸はs方向の距離であり、縦軸はベータトロン関数βである。なお、図18においてベータトロン関数αは記載していない。図18(a)はx方向のベータトロン関数βであり、図18(b)はy方向のベータトロン関数βである。図19は、本発明の実施の形態2による回転偏向部からアイソセンタICにおける荷電粒子の運動量分散関数を示す図である。横軸はs方向の距離であり、縦軸は運動量分散関数ηである。図19(a)はx方向の運動量分散関数ηであり、図19(b)はy方向の運動量分散関数ηである。図18(a)、図18(b)、図19(a)、図19(b)のそれぞれにおいて、回転偏向部60の入口は、s=0mであり、回転偏向部60の出口は、左端のs=14m付近である。実施の形態2の例では、4個の四極電磁石11a、11b、11c、11dにより、アイソセンタICのビーム特性が良好となる輸送行列を実現している。
 本発明の実施の形態2による回転偏向部60の電磁石配列が、図3に示した回転偏向部60の電磁石配列であり、図18及び図19に示したベータトロン関数及び運動量分散関数となるような各電磁石の強度の場合に、回転偏向部60の輸送行列を用いて回転偏向部60のビーム設計を行う例を考える。この場合の回転ガントリの輸送行列の一部を選択して取り出す。式(3)の輸送行列の一部を取り出したマトリクスは、式(7)である。輸送行列の一部を選択して取り出したマトリクスを、選択マトリクスと呼ぶことにする。
Figure JPOXMLDOC01-appb-M000007
 式(7)のマトリクスおけるm11、m12、m33、m34のそれぞれは、アイソセンタICにおけるx方向のビームサイズ、とy方向のビームサイズと、回転ガントリ入口、すなわち回転偏向部60の入口の位相空間上の分布を関連つける量である。この量を基に、回転ガントリ入口におけるx方向の理想的な位相空間上の線状分布を、以下のように決定する。
 位相空間上の線状ビームを特徴つける量は、位相空間上の端点の座標(x、x’)である。そこで、端点の座標(x、x’)を未知数として、式(7)の輸送行列を用いて、式(8)の連立方程式を導出する。
Figure JPOXMLDOC01-appb-M000008
 式(8)の意味を説明する。ガントリ角度を0°又は180°にした場合と、ガントリ角度を90°又は270°にした場合とで、同じビームサイズが実現するための式になっている。上記2つの角度でほぼ同じビームサイズを実現できる光学系(荷電粒子ビーム輸送系59)を実現すると、上記2つの角度の間の角度では問題にならない程度の差しか発生しない。
 式(8)を説いて得られた端点の座標(x、x’)を一端とし、この端点の座標(x、x’)の原点を中心にした点対称である座標、すなわち(-x、-x’)を他端に持つ位相空間上の線状分布を、ガントリ入口(回転偏向部60の入口)で実現すれば、従来の設定方法で設計された回転偏向部60を通過すことで、アイソセンタICにおいて、ビームサイズの回転依存性がなくなる。
 上記の手順で得られた実施の形態2による回転偏向部60の入口における荷電粒子の位相空間分布は、図20のようになる。また、実施の形態2による回転偏向部60のアイソセンタICでの荷電粒子の位相空間分布は、図21乃至図28のようになる。
 図21乃至図28は、本発明の実施の形態2による回転偏向部の角度がそれぞれ0°、30°、60°、90°、120°、150°、180°、210°の場合におけるアイソセンタICでの荷電粒子の位相空間分布を示す図である。各図の(a)はx方向の位相空間分布であり、各図の(b)はy方向の位相空間分布である。各図の(a)において、横軸はx方向の距離xであり、縦軸はx方向の傾きx’である。各図の(b)において、横軸はy方向の距離yであり、縦軸はy方向の傾きy’である。図中に示したビームサイズ(rms)は、実施の形態1で説明したビームサイズ1σである。
 回転偏向部60の角度が0°の場合のx方向及びy方向のビームサイズ1σは、図21に示すように、それぞれ1.29mm及び1.29mmである。回転偏向部の角度が30°の場合のx方向及びy方向のビームサイズ1σは、図22に示すように、それぞれ1.30mm及び1.29mmである。回転偏向部の角度が60°の場合のx方向及びy方向のビームサイズ1σは、図23に示すように、それぞれ1.30mm及び1.30mmである。回転偏向部の角度が90°の場合のx方向及びy方向のビームサイズ1σは、図24に示すように、それぞれ1.30mm及び1.30mmである。
 回転偏向部の角度が120°の場合のx方向及びy方向のビームサイズ1σは、図25に示すように、それぞれ1.29mm及び1.30mmである。回転偏向部の角度が150°の場合のx方向及びy方向のビームサイズ1σは、図26に示すように、それぞれ1.29mm及び1.30mmである。回転偏向部の角度が180°の場合のx方向及びy方向のビームサイズ1σは、図27に示すように、それぞれ1.29mm及び1.29mmである。回転偏向部の角度が210°の場合のx方向及びy方向のビームサイズ1σは、図28に示すように、それぞれ1.30mm及び1.29mmである。
 図21乃至図28に示すように、x方向及びy方向の位相空間分布の形状や位相は異なるが、x方向及びy方向のビームサイズ1σは、ほぼ等しくなっている。図29は、本発明の実施の形態2による回転偏向部の角度とアイソセンタICにおけるビームサイズの関係を示す図である。横軸は回転偏向部60のガントリ角度θであり、縦軸は標準偏差σの3倍に相当する大きさであるビームサイズ3σである。実線で示した特性81はx方向特性であり、破線で示した特性82はy方向特性である。実施の形態2の回転偏向部60を上記のようにすることで、回転偏向部60の回転角度、すなわち回転ガントリの回転角度に依存せずに、ほぼ一定のビームサイズを実現できる。
 実施の形態2の荷電粒子ビーム輸送系59の固定輸送部61は、固定輸送部61の出口を回転偏向部60の入口の条件である、βx=βy、ηx=0、η’x=0、ηy=0、η’y=0に加えて、上述した式(8)を説いて得られた端点の座標(x、x’)を一端とし、この端点の座標(x、x’)の原点を中心にした点対称である座標、すなわち(-x、-x’)を他端に持つ位相空間上の線状分布を用いて実現する。
 固定輸送部61の設計フローは、実施の形態1で説明した設計フローと同じである。ただし、実施の形態2では、実施の形態1とは異なり、拘束条件においてβx=βy、αx=αyの制約条件は設定されない。
 図30と図33乃至図37は、本発明の実施の形態2による回転偏向部の角度がそれぞれ0°、30°、60°、90°、120°、150°の場合における回転偏向部の入口からアイソセンタにおける荷電粒子のビームサイズの関係を示す図である。横軸はs方向の距離であり、縦軸はビームサイズ3σである。各図において、区間S1、S3、S5は、それぞれ偏向電磁石10a、10b、10cを通過する区間であり、区間S2は四極電磁石11a、11bが配置された区間であり、区間S4は四極電磁石11c、11dが配置された区間である。区間S6は、偏向電磁石10cの出口、すなわち回転偏向部60の出口からアイソセンタICまでの区間である。
 図31は本発明の実施の形態2による回転偏向部の角度が0°の場合におけるx方向ビームサイズを示す図であり、図32は本発明の実施の形態2による回転偏向部の角度が0°の場合におけるy方向ビームサイズを示す図である。x方向特性83は、x方向の位相空間分布を図20(a)の線状分布としてシミュレーションして得られた結果である。y方向特性84は、y方向の位相空間分布を図20(b)の楕円分布としてシミュレーションして得られた結果である。x方向特性85は、x方向のツイスパラメータを用いてシミュレーションして得られた結果である。y方向特性86は、y方向のツイスパラメータを用いてシミュレーションして得られた結果である。図31では、x方向特性83とx方向特性85はあまり一致しない。これは、x方向の位相空間分布を楕円分布と近似すると、正確なビームサイズが得られないことを示している。図32では、y方向特性84とy方向特性86はよく一致している。区間S1ではずれているが、そのずれは小さく無視でる程度である。
 実施の形態2の荷電粒子ビーム輸送系59の回転偏向部60は、3つの偏向電磁石10a、10b、10cと4つの四極電磁石11a、11b、11c、11dの電磁石の強度を変更することなく、従来の設計における入口の条件から、βx=βy、αx=αyの制約条件を外し、所定のx方向の位相空間分布を追加することで、アイソセンタICにおいて、ビームサイズの回転依存性がなくなる。所定のx方向の位相空間分布は、式(8)を説いて得られた端点の座標(x、x’)を一端とし、この端点の座標(x、x’)の原点を中心にした点対称である座標、すなわち(-x、-x’)を他端に持つ位相空間上の線状分布である。所定のx方向の位相空間分布は、荷電粒子ビーム輸送系59の固定輸送部61における四極電磁石を追加するか、もしくは四極電磁石の位置をずらした上で、励磁量を再最適化し、回転偏向部60の入口で必要な粒子分布を実現する。実施の形態2の荷電粒子ビーム輸送系59は、実施の形態1に比べて、拘束条件が少なく、課題を達成できるため、調整要素を少なくし、さらなるコンパクト化が図れる。
 実施の形態2の荷電粒子ビーム輸送系59は、回転偏向部60の入口におけるx方向の位相空間分布を、上記に示した所定のx方向の位相空間分布にしたので、荷電粒子ビーム輸送系59の回転偏向部60の入り口にてエミッタンスの非対称が存在しても、アイソセンタICで回転依存性の少ないビームサイズを実現することができる。
10、10a、10b、10c…偏向電磁石、14…ビーム中心線、15…ガントリ回転軸、31…荷電粒子ビーム、45…照射対象、51…粒子線治療装置、54…加速器、58、58a、58b…粒子線照射装置、59…荷電粒子ビーム輸送系、60…回転偏向部、61…固定輸送部、βx、βy、αx、αy…ツイスパラメータ、IC…アイソセンタ。

Claims (9)

  1.  加速器から出射された荷電粒子ビームが、前記加速器の出射位置における前記荷電粒子ビームの進行方向に垂直な面内で、前記加速器の周回軌道面内の方向をx方向、前記荷電粒子ビームの進行方向に垂直な面内で前記x方向と直交する方向をy方向とした場合、x方向のエミッタンスとy方向のエミッタンスが異なる荷電粒子ビームであり、アイソセンタを中心に回転可能な回転ガントリに搭載された粒子線照射装置に前記荷電粒子ビームを輸送する荷電粒子ビーム輸送系であって、
    前記回転ガントリに搭載されると共に、前記回転ガントリのガントリ回転軸の周りを回転する回転偏向部と、前記加速器から前記回転偏向部までの固定輸送部と、を備え、
    前記回転偏向部は複数の偏向電磁石を有し、
    前記ガントリ回転軸方向から見た場合において、前記粒子線照射装置のビーム中心線と、前記回転偏向部の入口における前記荷電粒子ビームの前記y方向の軸との角度をガントリ角度とし、
    前記加速器の出射位置における前記x方向のエミッタンスと前記y方向のエミッタンスとが分離されてそれぞれのエミッタンスが保たれるように前記荷電粒子ビームが輸送される前記回転ガントリの角度をガントリ基準角度とし、
    前記ガントリ角度が前記ガントリ基準角度である場合に、前記荷電粒子ビームの位相空間分布における位相が前記回転偏向部の入口から前記アイソセンタまで進む変化分を第1の位相進みとし、
    前記ガントリ角度が前記ガントリ基準角度から90°回転した角度である場合に、前記荷電粒子ビームの前記位相空間分布における位相が前記回転偏向部の入口から前記アイソセンタまで進む変化分を第2の位相進みとし、
    前記固定輸送部は、前記回転偏向部の入口における前記荷電粒子ビームの前記位相空間分布の位相が、前記第1の位相進み及び前記第2の位相進みの平均値に基づいた演算を基に決定された位相となるようにしたことを特徴とする荷電粒子ビーム輸送系。
  2.  前記固定輸送部は、前記回転偏向部の入口における前記荷電粒子ビームの前記位相空間分布の位相が、前記第1の位相進みを180°の位相進みで基準化した位相進みと前記第2の位相進みを180°の位相進みで基準化した位相進みとの平均値に基づいた位相となるようにしたことを特徴とする請求項1記載の荷電粒子ビーム輸送系。
  3.  前記固定輸送部は、前記回転偏向部の入口における前記荷電粒子ビームの前記位相空間分布の位相が、前記第1の位相進みを180°の位相進みで基準化した位相進みと前記第2の位相進みを180°の位相進みで基準化した位相進みとの平均値を負にした位相となるようにしたことを特徴とする請求項1記載の荷電粒子ビーム輸送系。
  4.  前記回転偏向部は、当該回転偏向部における前記x方向及び前記y方向のツイスパラメータが、入口及び出口において合うようにされたことを特徴とする請求項1乃至3のいずれか1項に記載の荷電粒子ビーム輸送系。
  5.  前記回転偏向部は、当該回転偏向部の入口及び出口において、前記x方向のツイスパラメータβxと前記y方向のツイスパラメータβyとが合うようにされ、かつ前記x方向のツイスパラメータαxと前記y方向のツイスパラメータαyとが合うようにされたことを特徴とする請求項4項記載の荷電粒子ビーム輸送系。
  6.  加速器から出射された荷電粒子ビームが、前記加速器の出射位置における前記荷電粒子ビームの進行方向に垂直な面内で、前記加速器の周回軌道面内の方向をx方向、前記荷電粒子ビームの進行方向に垂直な面内で前記x方向と直交する方向をy方向とした場合、x方向のエミッタンスとy方向のエミッタンスが異なる荷電粒子ビームであり、アイソセンタを中心に回転可能な回転ガントリに搭載された粒子線照射装置に前記荷電粒子ビームを輸送する荷電粒子ビーム輸送系であって、
    前記回転ガントリに搭載されると共に、前記回転ガントリのガントリ回転軸の周りを回転する回転偏向部と、前記加速器から前記回転偏向部までの固定輸送部と、を備え、
    前記回転偏向部は複数の偏向電磁石を有し、
    前記固定輸送部は、前記回転偏向部の入口において、前記荷電粒子ビームの前記x方向における位相空間分布が、所定の選択マトリクスを用いて算出された線状分布で近似できる分布であり、
    前記選択マトリクスは、
    前記回転偏向部の入口及び出口において、前記y方向及び前記x方向の運動量分散関数が0であり、かつ前記運動量分散関数における前記x方向及び前記y方向に垂直な前記荷電粒子ビームの進行方向の傾きが0である場合における前記回転偏向部の入口から前記アイソセンタまでの輸送行列の一部であって、
    前記回転偏向部の入口における前記荷電粒子ビームを構成する荷電粒子の位置(x1、y1)及びこの位置における前記荷電粒子ビームの進行方向の傾き(x’1、y’1)を、前記アイソセンタにおける前記荷電粒子の位置(x2、y2)及びこの位置における前記荷電粒子ビームの進行方向の傾き(x’2、y’2)に結び付ける輸送行列の成分m11、m12、m33、m34を有し、
    前記アイソセンタにおける前記y方向のビームサイズをσyとし、
    110+m12x’0=σyを第1式とし、
    330+m34x’0=σyを第2式とし、
    前記線状分布の一方の端点(x0、x’0)と他方の端点(-x0、-x’0)は、前記第1式及び前記第2式から算出されたことを特徴とする荷電粒子ビーム輸送系。
  7.  前記回転偏向部は、当該回転偏向部における前記x方向及び前記y方向のツイスパラメータが、入口及び出口において合うようにされたことを特徴とする請求項6項記載の荷電粒子ビーム輸送系。
  8.  前記回転偏向部は、当該回転偏向部の入口及び出口において、前記x方向のツイスパラメータβxと前記y方向のツイスパラメータβyとが合うようにされたことを特徴とする請求項7項記載の荷電粒子ビーム輸送系。
  9.  荷電粒子ビームを発生させ、この荷電粒子ビームを加速器で加速させるビーム発生装置と、
    前記加速器により加速された荷電粒子ビームが、前記加速器の出射位置における前記荷電粒子ビームの進行方向に垂直な面内で、前記加速器の周回軌道面内の方向をx方向、前記荷電粒子ビームの進行方向に垂直な面内で前記x方向と直交する方向をy方向とした場合、x方向のエミッタンスとy方向のエミッタンスが異なる荷電粒子ビームであり、この荷電粒子ビームを輸送する荷電粒子ビーム輸送系と、
    前記荷電粒子ビーム輸送系で輸送された荷電粒子ビームを照射対象に照射する粒子線照射装置と、
    前記粒子線照射装置を搭載し、アイソセンタを中心に回転可能な回転ガントリと、を備え、
    前記荷電粒子ビーム輸送系は、請求項1乃至8のいずれか1項に記載の荷電粒子ビーム輸送系であることを特徴とする粒子線治療装置。
PCT/JP2012/063268 2012-05-24 2012-05-24 荷電粒子ビーム輸送系及び粒子線治療装置 WO2013175600A1 (ja)

Priority Applications (6)

Application Number Priority Date Filing Date Title
CN201280073432.4A CN104335686B (zh) 2012-05-24 2012-05-24 带电粒子束输送系统及粒子射线治疗装置
EP12877384.3A EP2858463A4 (en) 2012-05-24 2012-05-24 LOAD PARTICLE BEAM TRANSPORT SYSTEM AND PARTICLE BEAM PROCESSING DEVICE
US14/382,637 US9078335B2 (en) 2012-05-24 2012-05-24 Charged particle beam transport system and particle beam therapy system
JP2014516579A JP5748912B2 (ja) 2012-05-24 2012-05-24 荷電粒子ビーム輸送系及び粒子線治療装置
PCT/JP2012/063268 WO2013175600A1 (ja) 2012-05-24 2012-05-24 荷電粒子ビーム輸送系及び粒子線治療装置
TW101137237A TWI486188B (zh) 2012-05-24 2012-10-09 帶電粒子射束輸送系統及粒子束治療裝置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2012/063268 WO2013175600A1 (ja) 2012-05-24 2012-05-24 荷電粒子ビーム輸送系及び粒子線治療装置

Publications (1)

Publication Number Publication Date
WO2013175600A1 true WO2013175600A1 (ja) 2013-11-28

Family

ID=49623331

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/063268 WO2013175600A1 (ja) 2012-05-24 2012-05-24 荷電粒子ビーム輸送系及び粒子線治療装置

Country Status (6)

Country Link
US (1) US9078335B2 (ja)
EP (1) EP2858463A4 (ja)
JP (1) JP5748912B2 (ja)
CN (1) CN104335686B (ja)
TW (1) TWI486188B (ja)
WO (1) WO2013175600A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015226672A (ja) * 2014-06-02 2015-12-17 株式会社日立製作所 粒子線治療システム及び装置、粒子線治療システムの制御方法
WO2017056087A1 (en) * 2015-09-30 2017-04-06 HIL Applied Medical Ltd. Beam delivery system for proton therapy for laser-accelerated protons
JP7259144B1 (ja) 2022-02-04 2023-04-17 株式会社東芝 粒子線治療システムの製造方法および粒子線治療システム

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6341655B2 (ja) * 2013-12-09 2018-06-13 株式会社東芝 円形加速器及び重粒子線治療装置
US9558915B2 (en) * 2014-01-13 2017-01-31 Frederick A. Flitsch Method and apparatus for a high resolution imaging system
EP3193567A4 (en) 2014-09-12 2018-05-09 Mitsubishi Electric Corporation Beam delivery system and particle beam therapy device
US10434337B2 (en) * 2014-12-04 2019-10-08 Kabushiki Kaisha Toshiba Particle beam adjustment device, particle beam adjustment method, and particle beam therapeutic device
US20180114667A1 (en) * 2015-09-11 2018-04-26 Mitsubishi Electric Corporation Electromagnet mounting frame, electromagnet device, and particle beam therapy system
JP6794302B2 (ja) * 2017-03-14 2020-12-02 株式会社東芝 回転照射装置、回転照射方法、及び回転照射治療装置
CN107596579B (zh) * 2017-10-12 2018-06-05 合肥中科离子医学技术装备有限公司 基于紧凑型超导回旋加速器的质子治疗系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09265000A (ja) 1996-03-29 1997-10-07 Mitsubishi Electric Corp 荷電粒子ビーム輸送系
JP2002541464A (ja) * 1999-04-01 2002-12-03 ジー エス アイ ゲゼルシャフト フュア シュベールイオーネンフォルシュンク エム ベー ハー イオン光学システムを有するガントリ
JP2003282300A (ja) * 2002-03-26 2003-10-03 Hitachi Ltd 粒子線治療システム
JP2005329252A (ja) * 2005-06-06 2005-12-02 Hitachi Ltd 粒子線治療システム
JP2006351339A (ja) 2005-06-15 2006-12-28 Natl Inst Of Radiological Sciences 荷電粒子線照射装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4870287A (en) * 1988-03-03 1989-09-26 Loma Linda University Medical Center Multi-station proton beam therapy system
JP3178381B2 (ja) * 1997-02-07 2001-06-18 株式会社日立製作所 荷電粒子照射装置
JPH11174198A (ja) 1997-12-12 1999-07-02 Mitsubishi Electric Corp 荷電粒子照射装置及びその運転方法
JP4691574B2 (ja) * 2008-05-14 2011-06-01 株式会社日立製作所 荷電粒子ビーム出射装置及び荷電粒子ビーム出射方法
US8063381B2 (en) * 2009-03-13 2011-11-22 Brookhaven Science Associates, Llc Achromatic and uncoupled medical gantry

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09265000A (ja) 1996-03-29 1997-10-07 Mitsubishi Electric Corp 荷電粒子ビーム輸送系
JP2002541464A (ja) * 1999-04-01 2002-12-03 ジー エス アイ ゲゼルシャフト フュア シュベールイオーネンフォルシュンク エム ベー ハー イオン光学システムを有するガントリ
JP2003282300A (ja) * 2002-03-26 2003-10-03 Hitachi Ltd 粒子線治療システム
JP2005329252A (ja) * 2005-06-06 2005-12-02 Hitachi Ltd 粒子線治療システム
JP2006351339A (ja) 2005-06-15 2006-12-28 Natl Inst Of Radiological Sciences 荷電粒子線照射装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP2858463A4 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015226672A (ja) * 2014-06-02 2015-12-17 株式会社日立製作所 粒子線治療システム及び装置、粒子線治療システムの制御方法
WO2017056087A1 (en) * 2015-09-30 2017-04-06 HIL Applied Medical Ltd. Beam delivery system for proton therapy for laser-accelerated protons
JP7259144B1 (ja) 2022-02-04 2023-04-17 株式会社東芝 粒子線治療システムの製造方法および粒子線治療システム
WO2023149170A1 (ja) * 2022-02-04 2023-08-10 株式会社 東芝 粒子線治療システムの製造方法および粒子線治療システム
JP2023114436A (ja) * 2022-02-04 2023-08-17 株式会社東芝 粒子線治療システムの製造方法および粒子線治療システム

Also Published As

Publication number Publication date
EP2858463A4 (en) 2015-12-23
JPWO2013175600A1 (ja) 2016-01-12
US20150038764A1 (en) 2015-02-05
US9078335B2 (en) 2015-07-07
EP2858463A1 (en) 2015-04-08
JP5748912B2 (ja) 2015-07-15
CN104335686A (zh) 2015-02-04
TWI486188B (zh) 2015-06-01
TW201347803A (zh) 2013-12-01
CN104335686B (zh) 2016-08-24

Similar Documents

Publication Publication Date Title
JP5748912B2 (ja) 荷電粒子ビーム輸送系及び粒子線治療装置
JP7245352B2 (ja) 粒子を用いた回転式の放射線治療を提供する方法
TWI404549B (zh) 膠塊土、膠塊土之製造方法、粒子束治療裝置,及治療計劃裝置
US7839973B2 (en) Treatment planning using modulability and visibility factors
JP5722559B2 (ja) 治療計画装置
US11938342B2 (en) Time optimized radiation treatment
JP6207755B2 (ja) ビーム輸送系及び粒子線治療装置
CN107949423A (zh) 粒子射束传输系统及其分段
CN107789749A (zh) 带电粒子束偏转装置及治疗系统
JP4639401B2 (ja) 荷電粒子線照射装置
JP2011206237A (ja) 荷電粒子ビーム輸送装置及び粒子線治療システム
EP2889055A1 (en) Control device for scanning electromagnet and particle beam therapy apparatus
US9630027B2 (en) Beam transport system and particle beam therapy system
JP2019195408A (ja) スキャニング照射装置および粒子線治療システム、ならびにスキャニング照射装置の調整方法
JP6117675B2 (ja) 粒子線治療装置
JP2019082389A (ja) ビーム輸送系および粒子線治療装置
JP7482048B2 (ja) 治療計画装置、治療計画方法及びコンピュータプログラム
WO2023238643A1 (ja) 粒子線ビーム制御システムおよび粒子線ビーム制御方法
WO2024082293A1 (en) Multicriterial treatment plan optimization using let cost functions
WO2020188890A1 (ja) 粒子線治療システム、計測粒子線ct画像生成方法、およびct画像生成プログラム
JP2015110171A (ja) 粒子線治療システム
WO2016027356A1 (ja) 粒子線治療システム及び位相空間粒子分布の計測方法
JP2014138669A (ja) 荷電粒子線治療装置
JP2021040899A (ja) 治療計画装置、粒子線治療システム及びコンピュータプログラム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12877384

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014516579

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 14382637

Country of ref document: US

WWE Wipo information: entry into national phase

Ref document number: 2012877384

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: DE