Monte Carlo simulation for proton dosimetry in magnetic fields: Fano test and magnetic field correction factorskBfor Farmer-type ionization chambers

Objective. In this contribution we present a special Fano test for charged particles in presence of magnetic fields in the MC code TOol for PArticle Simulation (TOPAS), as well as the determination of magnetic field correction factorskBfor Farmer-type ionization chambers using proton beams.Approach. Customized C++ extensions for TOPAS were implemented to model the special Fano tests in presence of magnetic fields for electrons and protons. The Geant4-specific transport parameters,DRoverRandfinalRange,were investigated to optimize passing rate and computation time. ThekBwas determined for the Farmer-type PTW 30013 ionization chamber, and 5 custom built ionization chambers with same geometry but varying inner radius, testing magnetic flux density ranging from 0 to 1.0 T and two proton beam energies of 157.43 and 221.05 MeV.Main results. Using the investigated parameters, TOPAS passed the Fano test within 0.39 ± 0.15% and 0.82 ± 0.42%, respectively for electrons and protons. The chamber response (kB,M,Q) gives a maximum at different magnetic flux densities depending of the chamber size, 1.0043 at 1.0 T for the smallest chamber and 1.0051 at 0.2 T for the largest chamber. The local dose differencecBremained ≤ 0.1% for both tested energies. The magnetic field correction factorkB, for the chamber PTW 30013, varied from 0.9946 to 1.0036 for both tested energies.Significance. The developed extension for the special Fano test in TOPAS MC code with the adjusted transport parameters, can accurately transport electron and proton particles in magnetic field. This makes TOPAS a valuable tool for the determination ofkB. The ionization chambers we tested showed thatkBremains small (≤0.72%). To the best of our knowledge, this is the first calculations ofkBfor proton beams. This work represents a significant step forward in the development of MRgPT and protocols for proton dosimetry in presence of magnetic field.

Proton beam dosimetry in the presence of magnetic fields using Farmer-type ionization chambers of different radii

BACKGROUND

Magnetic resonance-guided proton therapy is promising, as it combines high-contrast imaging of soft tissue with highly conformal dose delivery. However, proton dosimetry in magnetic fields using ionization chambers is challenging since the dose distribution as well as the detector response are perturbed.

PURPOSE

This work investigates the effect of the magnetic field on the ionization chamber response, and on the polarity and ion recombination correction factors, which are essential for the implementation of a proton beam dosimetry protocol in the presence of magnetic fields.

METHODS

Three Farmer-type cylindrical ionization chambers, the 30013 with 3 mm inner radius (PTW, Freiburg, Germany) and two custom built chambers "R1" and "R6" with 1 and 6 mm inner radii respectively were placed at the center of an experimental electromagnet (Schwarzbeck Mess - Elektronik, Germany) 2 cm depth of an in-house developed 3D printed water phantom. The detector response was measured for a 3×10 cm² field of mono-energetic protons 221.05 MeV/u for the three chambers, and with an additional proton beam of 157.43 MeV/u for the chamber PTW 30013. The magnetic flux density was varied between 0.1 and 1.0 Tesla in steps of 0.1 Tesla.

RESULTS

At both energies, the ionization chamber PTW 30013 showed a non-linear response as a function of the magnetic field strength, with a decrease of the ionization chamber response of up to 0.27±0.06% (1 SD) at 0.2 Tesla, followed by a smaller effect at higher magnetic field strength. For the chamber R1, the response decreased slightly with the magnetic field strength up to 0.45±0.12% at 1 Tesla, and for the chamber R6, the response decreased up to 0.54±0.13% at 0.1 Tesla, followed by a plateau up to 0.3 Tesla, and a weaker effect at higher magnetic field strength. The dependence of the polarity and recombination correction factor on the magnetic field was ≤ $ \le $ 0.1% for the chamber PTW 30013.

CONCLUSIONS

The magnetic field has a small but significant effect on the chamber response in the low magnetic field region for the chamber PTW 30013 and for R6, and in the high magnetic field region for the chamber R1. Corrections may be necessary for ionization chamber measurements, depending on both the chamber volume and the magnetic flux density. No significant effect of the magnetic field on the polarity and recombination correction factor was detected in this work for the ionization chamber PTW 30013. This article is protected by copyright. All rights reserved.

End-to-end test for fractionated online adaptive MR-guided radiotherapy using a deformable anthropomorphic pelvis phantom

Objective.In MR-guided radiotherapy (MRgRT) for prostate cancer treatments inter-fractional anatomy changes such as bladder and rectum fillings may be corrected by an online adaption of the treatment plan. To clinically implement such complex treatment procedures, however, specific end-to-end tests are required that are able to validate the overall accuracy of all treatment steps from pre-treatment imaging to dose delivery.Approach.In this study, an end-to-end test of a fractionated and online adapted MRgRT prostate irradiation was performed using the so-called ADAM-PETer phantom. The phantom was adapted to perform 3D polymer gel (PG) dosimetry in the prostate and rectum. Furthermore, thermoluminescence detectors (TLDs) were placed at the center and on the surface of the prostate for additional dose measurements as well as for an external dose renormalization of the PG. For the end-to-end test, a total of five online adapted irradiations were applied in sequence with different bladder and rectum fillings, respectively.Main results.A good agreement of measured and planned dose was found represented by highγ-index passing rates (3%/3mmcriterion) of the PG evaluation of98.9%in the prostate and93.7%in the rectum. TLDs used for PG renormalization at the center of the prostate showed a deviation of-2.3%.Significance.The presented end-to-end test, which allows for 3D dose verification in the prostate and rectum, demonstrates the feasibility and accuracy of fractionated and online-adapted prostate irradiations in presence of inter-fractional anatomy changes. Such tests are of high clinical importance for the commissioning of new image-guided treatment procedures such as online adaptive MRgRT.

PAGAT gel dosimetry for everyone: gel production, measurement and evaluation

Polymer gel (PG) dosimetry is a valuable tool to measure complex dose distributions in 3D with a high spatial resolution. However, due to complex protocols that need to be followed for in-house produced PGs and the high costs of commercially available gels, PG gels are only rarely applied in quality assurance procedures worldwide. In this work, we provide an introduction to perform highly standardized dosimetric PG experiments using PAGAT (PolyAcrylamide Gelatine gel fabricated at ATmospheric conditions) dosimetry gel. PAGAT gel can be produced at atmospheric conditions, at low costs and is evaluated using magnetic resonance imaging (MRI). The conduction of PG experiments is described in great detail including the gel production, treatment planning, irradiation, MRI evaluation and post-processing procedure. Furthermore, a plugin in an open source image processing tool for post-processing is provided free of charge that allows a standardized and reproducible analysis of PG experiments.

Development of phantom materials with independently adjustable CT- and MR-contrast at 0.35, 1.5 and 3 T

Quality assurance in magnetic resonance (MR)-guided radiotherapy lacks anthropomorphic phantoms that represent tissue-equivalent imaging contrast in both computed tomography (CT) and MR imaging. In this study, we developed phantom materials with individually adjustable CT value as well as [Formula: see text]- and [Formula: see text]-relaxation times in MR imaging at three different magnetic field strengths. Additionally, their experimental stopping power ratio (SPR) for carbon ions was compared with predictions based on single- and dual-energy CT. Ni-DTPA doped agarose gels were used for individual adjustment of [Formula: see text] and [Formula: see text] at [Formula: see text] and 3.0 T. The CT value was varied by adding potassium chloride (KCl). By multiple linear regression, equations for the determination of agarose, Ni-DTPA and KCl concentrations for given [Formula: see text] [Formula: see text] and CT values were derived and employed to produce nine specific soft tissue samples. Experimental [Formula: see text] [Formula: see text] and CT values of these soft tissue samples were compared with predictions and additionally, carbon ion SPR obtained by range measurements were compared with predictions based on single- and dual-energy CT. The measured CT value, [Formula: see text] and [Formula: see text] of the produced soft tissue samples agreed very well with predictions based on the derived equations with mean deviations of less than [Formula: see text] While single-energy CT overestimates the measured SPR of the soft tissue samples, the dual-energy CT-based predictions showed a mean SPR deviation of only [Formula: see text] To conclude, anthropomorphic phantom materials with independently adjustable CT values as well as [Formula: see text] and [Formula: see text] relaxation times at three different magnetic field strengths were developed. The derived equations describe the material specific relaxation times and the CT value in dependence on agarose, Ni-DTPA and KCl concentrations as well as the chemical composition of the materials based on given [Formula: see text] and CT value. Dual-energy CT allows accurate prediction of the carbon ion range in these materials.

Feasibility of markerless fluoroscopic real-time tumor detection for adaptive radiotherapy: development and end-to-end testing

Respiratory-gated radiotherapy treatments of lung tumors reduce the irradiated normal tissue volume and potentially lower the risk of side effects. However, in clinical routine, the gating signal is usually derived from external markers or other surrogate signals and may not always correlate well with the actual tumor position. This study uses the kV-imaging system of a LINAC in combination with a multiple template matching algorithm for markerless real-time detection of the tumor position in a dynamic anthropomorphic porcine lung phantom. The tumor was realized by a small container filled with polymer dosimetry gel, the so-called gel tumor. A full end-to-end test for a gated treatment was performed and the geometric and dosimetric accuracy was validated. The accuracy of the tumor detection algorithm in SI- direction was found to be [Formula: see text] mm and the gel tumor was automatically detected in 98 out of 100 images. The measured 3D dose distribution showed a uniform coverage of the gel tumor and comparison with the treatment plan revealed a high 3D [Formula: see text]-passing rate of [Formula: see text] ([Formula: see text]). The simulated treatment confirmed the employed margin sizes for residual motion within the gating window and serves as an end-to-end test for a gated treatment based on a markerless fluoroscopic real-time tumor detection.

Absolute dosimetry with polymer gels-a TLD reference system

BACKGROUND AND PURPOSE

Absolute dosimetry in 3D with polymer gels (PG) is generally complicated and usually requires a second independent measurement with conventional detectors. This is why, PG are often used only for relative dosimetry. To overcome this drawback, we combine PG with a 1D thermoluminescence (TL) detector within the same measurement. The TL detector information is then used as additional information for calibration of the gel.

MATERIALS AND METHODS

The PAGAT dosimetry gel was used in combination with TLD600 (LiF:Mg,Ti). TL detectors were attached on the surface of the PG container placed inside a cylindrical phantom. To test the usability of this setup, two irradiation geometries were carried out: (a) homogeneous target coverage and (b) small-field irradiation. PG was evaluated with magnetic resonance imaging (MRI) and the TL detectors with a Harshaw 5500 hot gas reader.

RESULTS

PG dosimetry alone showed deviations of up to 4% as compared to calculations. Including additionally the dose information of the TL detectors for PG calibration, this deviation was decreased to less than 1% for both irradiation geometries. This is also reflected by the very high [Formula: see text]-passing rates of  >  96% (3%/3 mm) and  >93% (2%/2 mm), respectively.

CONCLUSION

This study presents a novel method combining 3D PG and TL dose measurements for the purpose of absolute 3D dose measurements that can also be applied in complex anthropomorphic phantoms using only a single measurement. The method was validated for two different irradiation geometries including a homogeneous large field as well as a small field irradiation with sharp dose gradients.

Feasibility of polymer gel-based measurements of radiation isocenter accuracy in magnetic fields

For conventional irradiation devices, the radiation isocenter accuracy is determined by star shot measurements on films. In magnetic resonance (MR)-guided radiotherapy devices, the results of this test may be altered by the magnetic field and the need to align the radiation and imaging isocenter may require a modification of measurement procedures. Polymer dosimetry gels (PG) may offer a way to perform both, the radiation and imaging isocenter test, however, first it has to be shown that PG reveal results comparable to the conventionally applied films. Therefore, star shot measurements were performed at a linear accelerator using PG as well as radiochromic films. PG were evaluated using MR imaging and the isocircle radius and the distance between the isocircle center and the room isocenter were determined. Two different types of experiments were performed: i) a standard star-shot isocenter test and (ii) a star shot, where the detectors were placed between the pole shoes of an experimental electro magnet operated either at 0 T or 1 T. For the standard star shot, PG evaluation was independent of the time delay after irradiation (1 h, 24 h, 48 h and 216 h) and the results were comparable to those of film measurements. Within the electro magnet, the isocircle radius increased from 0.39  ±  0.01 mm to 1.37  ±  0.01 mm for the film and from 0.44  ±  0.02 mm to 0.97  ±  0.02 mm for the PG-measurements, respectively. The isocenter distance was essentially dependent on the alignment of the magnet to the isocenter and was between 0.12  ±  0.02 mm and 0.82  ±  0.02 mm. The study demonstrates that evaluation of the PG directly after irradiation is feasible, if only geometrical parameters are of interest. This allows using PG for star shot measurements to evaluate the radiation isocenter accuracy with comparable accuracy as with radiochromic films.

Radiation dosimetry in magnetic fields with Farmer-type ionization chambers: determination of magnetic field correction factors for different magnetic field strengths and field orientations

The aim of this work was to determine magnetic field correction factors that are needed for dosimetry in hybrid devices for MR-guided radiotherapy for Farmer-type ionization chambers for different magnetic field strengths and field orientations. The response of six custom-built Farmer-type chambers irradiated at a 6 MV linac was measured in a water tank positioned in a magnet with magnetic field strengths between 0.0 T and 1.1 T. Chamber axis, beam and magnetic field were perpendicular to each other and both magnetic field directions were investigated. EGSnrc Monte Carlo simulations were compared to the measurements and simulations with different field orientations were performed. For all geometries, magnetic field correction factors, [Formula: see text], and perturbation factors were calculated. A maximum increase of 8.8% in chamber response was measured for the magnetic field perpendicular to chamber and beam axis. The measured chamber response could be reproduced by adjusting the dead volume layer near the chamber stem in the Monte Carlo simulations. For the magnetic field parallel to the chamber axis or parallel to the beam, the simulated response increased by 1.1% at maximum for field strengths up to 1.1 T. A complex dependence of the response was found on chamber radius, magnetic field strength and orientation of beam, chamber axis and magnetic field direction. Especially for magnetic fields perpendicular to beam and chamber axis, the exact sensitive volume has to be considered in the simulations. To minimize magnetic field correction factors and the influence of dead volumes on the response of Farmer chambers, a measurement set-up with the magnetic field parallel to the chamber axis or parallel to the beam is recommended for dosimetry.

Heavy Ion Therapy: Status and Perspectives

Starting with the pioneering work at the University of California in Berkeley in 1977, heavy ion radiotherapy has been of increasing interest especially in Japan and Europe in the last decade. There are currently 3 facilities treating patients with carbon ions, two of them in Japan within a clinical setting. In Germany, a research therapy facility is in operation and the construction of a new hospital based facility at the Heidelberg university will be started soon. An outline of the current status of heavy ion radiotherapy is given with emphasis to the technical aspects of the respective facilities. This includes a description of passive and active beam shaping systems, as well as their implications for treatment planning and dosimetry. The clinical trials and routine treatments performed at the German heavy ion facility are summarized.

An overview over the upcoming new facilities and their technical possibilities is given. It is discussed what the necessary improvements are to fully exploit the potential of these facilities. Especially the new Heidelberg facility with the possibility of active beam scanning in combination with the first isocentric gantry for ions and offering beams of protons, helium, oxygen and carbon ions has implications on treatment planning, dosimetry and quality assurance. The necessary and ongoing developments in these areas are summarized. The new facilities also offer the possibilities to perform more extensive clinical studies and to explore future indications for radiotherapy with heavy ions. An overview over the indications and treatment schemes is also given.

A voxel-based multiscale model to simulate the radiation response of hypoxic tumors

PURPOSE

In radiotherapy, it is important to predict the response of tumors to irradiation prior to the treatment. This is especially important for hypoxic tumors, which are known to be highly radioresistant. Mathematical modeling based on the dose distribution, biological parameters, and medical images may help to improve this prediction and to optimize the treatment plan.

METHODS

A voxel-based multiscale tumor response model for simulating the radiation response of hypoxic tumors was developed. It considers viable and dead tumor cells, capillary and normal cells, as well as the most relevant biological processes such as (i) proliferation of tumor cells, (ii) hypoxia-induced angiogenesis, (iii) spatial exchange of cells leading to tumor growth, (iv) oxygen-dependent cell survival after irradiation, (v) resorption of dead cells, and (vi) spatial exchange of cells leading to tumor shrinkage. Oxygenation is described on a microscopic scale using a previously published tumor oxygenation model, which calculates the oxygen distribution for each voxel using the vascular fraction as the most important input parameter. To demonstrate the capabilities of the model, the dependence of the oxygen distribution on tumor growth and radiation-induced shrinkage is investigated. In addition, the impact of three different reoxygenation processes is compared and tumor control probability (TCP) curves for a squamous cells carcinoma of the head and neck (HNSSC) are simulated under normoxic and hypoxic conditions.

RESULTS

The model describes the spatiotemporal behavior of the tumor on three different scales: (i) on the macroscopic scale, it describes tumor growth and shrinkage during radiation treatment, (ii) on a mesoscopic scale, it provides the cell density and vascular fraction for each voxel, and (iii) on the microscopic scale, the oxygen distribution may be obtained in terms of oxygen histograms. With increasing tumor size, the simulated tumors develop a hypoxic core. Within the model, tumor shrinkage was found to be significantly more important for reoxygenation than angiogenesis or decreased oxygen consumption due to an increased fraction of dead cells. In the studied HNSSC-case, the TCD50 values (dose at 50% TCP) decreased from 71.0 Gy under hypoxic to 53.6 Gy under the oxic condition.

CONCLUSIONS

The results obtained with the developed multiscale model are in accordance with expectations based on radiobiological principles and clinical experience. As the model is voxel-based, radiological imaging methods may help to provide the required 3D-characterization of the tumor prior to irradiation. For clinical application, the model has to be further validated with experimental and clinical data. If this is achieved, the model may be used to optimize fractionation schedules and dose distributions for the treatment of hypoxic tumors.

Estimation of the uncertainty of elastic image registration with the demons algorithm

The accuracy of elastic image registration is limited. We propose an approach to detect voxels where registration based on the demons algorithm is likely to perform inaccurately, compared to other locations of the same image. The approach is based on the assumption that the local reproducibility of the registration can be regarded as a measure of uncertainty of the image registration. The reproducibility is determined as the standard deviation of the displacement vector components obtained from multiple registrations. These registrations differ in predefined initial deformations. The proposed approach was tested with artificially deformed lung images, where the ground truth on the deformation is known. In voxels where the result of the registration was less reproducible, the registration turned out to have larger average registration errors as compared to locations of the same image, where the registration was more reproducible. The proposed method can show a clinician in which area of the image the elastic registration with the demons algorithm cannot be expected to be accurate.

Impact of enhancements in the local effect model (LEM) on the predicted RBE-weighted target dose distribution in carbon ion therapy

Biological optimization for treatment planning in carbon ion therapy is currently based on the first version of the local effect model (LEM I). Further developments implemented in the latest version (LEM IV) allowed to predict more accurately the Relative Biological Effectiveness (RBE) in-vitro. The main goal of this study is to compare the LEM IV against LEM I under treatment-like conditions for idealized target geometries. Therefore, physical dose distributions resulting from the biological optimization with LEM I were used to recalculate the RBE-weighted dose distribution based on LEM IV. Input parameters representing the clinical endpoints late toxicity in the central nervous system and the tumor control for chordoma were chosen to investigate the impact of changes on the predicted isoeffective dose levels. The recalculated RBE-weighted dose distributions show an increase within the target region, and the mean RBE-weighted dose values are dependent on the geometry and decrease with increasing target dimension. The differences between predictions of LEM IV and LEM I are less than 10% for typical tumor volumes treated in the pilot project at GSI. Median RBE-weighted doses predicted by LEM IV in the target region are consistent with clinically observed dose-response behavior as demonstrated by comparison to the 5-year local control curve for skull base chordoma.

Comparison of two respiration monitoring systems for 4D imaging with a Siemens CT using a new dynamic breathing phantom

Four-dimensional computed tomography (4D-CT) requires breathing information from the patient, and for this, several systems are available. Testing of these systems, under realistic conditions, requires a phantom with a moving target and an expandable outer contour. An anthropomorphic phantom was developed to simulate patient breathing as well as lung tumor motion. Using the phantom, an optical camera system (GateCT) and a pressure sensor (AZ-733V) were simultaneously operated, and 4D-CTs were reconstructed with a Siemens CT using the provided local-amplitude-based sorting algorithm. The comparison of the tumor trajectories of both systems revealed discrepancies up to 9.7 mm. Breathing signal differences, such as baseline drift, temporal resolution and noise level were shown not to be the reason for this. Instead, the variability of the sampling interval and the accuracy of the sampling rate value written on the header of the GateCT-signal file were identified as the cause. Interpolation to regular sampling intervals and correction of the sampling rate to the actual value removed the observed discrepancies. Consistently, the introduction of sampling interval variability and inaccurate sampling rate values into the header of the AZ-733V file distorted the tumor trajectory for this system. These results underline the importance of testing new equipment thoroughly, especially if components of different manufacturers are combined.

A stochastic approach to estimate the uncertainty involved in B-spline image registration

Uncertainties in image registration may be a significant source of errors in anatomy mapping as well as dose accumulation in radiotherapy. It is, therefore, essential to validate the accuracy of image registration. Here, we propose a method to detect areas where mono modal B-spline registration performs well and to distinguish those from areas of the same image, where the registration is likely to be less accurate. It is a stochastic approach to automatically estimate the uncertainty of the resulting displacement vector field. The coefficients resulting from the B-spline registration are subject to moderate and randomly performed variations. A quantity is proposed to characterize the local sensitivity of the similarity measure to these variations. We demonstrate the statistical dependence between the local image registration error and this quantity by calculating their mutual information. We show the significance of the statistical dependence with an approach based on random redistributions. The proposed method has the potential to divide an image into subregions which differ in the magnitude of their average registration error.

Treatment planning for carbon ion radiotherapy in Germany: Review of clinical trials and treatment planning studies

The GSI carbon ion radiotherapy facility established the first completely active beam shaping system for heavy ions, using energy variation on the synchrotron and pencil beam scanning. The introduction of an active beam shaping system for carbon ions has considerable impact on the design of the treatment planning system (TPS). The TPS has to account for the capability of the beam delivery and the biological modelling, which is needed to calculate the RBE for the resulting varying depth dose modulation. The TPS used in clinical routine with carbon ions is described and its use in treatment planning studies are outlined. A clinical trial with carbon ion therapy as primary therapy for chordoma and chondrosarcoma of the base of skull has been completed in 2001. Currently, carbon ion therapy as a boost treatment together with conventional conformal photon therapy or IMRT is under investigation in clinical trials for adenoid cystic carcinoma, chordoma and chondrosarcoma of the cervical spine and sacrococcygeal chordoma. Treatment planning studies comparing carbon ion therapy with IMRT, using optimization of combination therapy, and optimization of beam-line design have already been completed. Analysis of uncertainties in treatment planning has been started with the investigation of range uncertainties stemming from CT imaging. Uncertainties coming from the beam delivery play only a minor role. An attempt to asses the uncertainties introduced in treatment plans by the biological modelling, was done, using phantom verification of calculated cell survival levels. The clinical trials and planning studies are of special importance for the upcoming new clinical ion facility of the Heidelberg university hospital.

A calibration procedure for beam monitors in a scanned beam of heavy charged particles

An international code of practice (CoP) for dosimetry based on standards of absorbed dose to water has recently been published by the IAEA [Technical Report Series No. 398, 2000] (TRS-398). This new CoP includes procedures for proton and heavy ion beams as well as all other beam qualities. In particular it defines reference conditions to which dose measurements should refer to. For proton and ion beams these conditions include dose measurements in the center of all possible modulated Bragg peaks. The recommended reference conditions in general are used also for the calibration of beam monitors. For a dynamic beam delivery system using beam scanning in combination with energy variation, like, e.g., at the German carbon ion radiotherapy facility, this calibration procedure is not appropriate. We have independently developed a different calibration procedure. Similar to the IAEA CoP this procedure is based on the measurement of absorbed dose to water. This is translated in terms of fluence which finally results in an energy-dependent calibration of the beam monitor in units of particle number per monitor unit, which is unique for all treatment fields. In contrast to the IAEA CoP, the reference depth is chosen to be very small. The procedure enables an accurate and reliable determination of calibration factors. In a second step, the calibration is verified by measurements of absorbed dose in various modulated Bragg peaks by comparing measured against calculated doses. The agreement between measured and calculated doses is usually better than 1% for homogeneous fields and the mean deviation for more inhomogeneous treatment fields, as they are used for patient treatments, is within 3%. It is proposed that the CoP in general, and in particular the IAEA TRS-398 should include explicit recommendations for the beam monitor calibration. These recommendations should then distinguish between systems using static and dynamic beams.

Treatment planning intercomparison for spinal chordomas using intensity-modulated photon radiation therapy (IMRT) and carbon ions

Spinal chordomas cannot be treated with an effective dose using conventional radiation therapy (RT) without exceeding the tolerance dose of the spinal cord while ensuring sufficient target coverage at the same time. In this study we investigate the potential physical advantages of combined photon intensity-modulated radiation therapy (IMRT) and raster-scanned carbon ion RT over photon IMRT alone. For a representative patient we generated a carbon ion RT plan and a photon IMRT plan. Additionally, combined plans consisting of both carbon ions and photon IMRT were calculated using ratios of 20:40 GyE, 30:30 GyE and 40:20 GyE. The best target coverage was obtained using carbon ions alone. Using a combination of photon IMRT and carbon ions, the target coverage was better than with photon IMRT alone. Due to the applied dose constraints, the sparing of the spinal cord was comparable for all plans. Using carbon ions alone, the non-target tissue volume irradiated to at least 30 GyE/50.4 GyE was reduced by 72%/84% compared to photon IMRT alone. These advantages were evident even with combined techniques. The actually delivered dose distribution is expected to be more dependent on patient misalignment with carbon ions compared with photon IMRT. A combination of carbon ions and photon IMRT might be preferable in order to profit by the physical advantages of carbon ions while ensuring a safe treatment.

[Late radiation changes after small volume radiosurgery of the rat brain. Measuring local cerebral blood flow and histopathological studies]

BACKGROUND

The goal of this study was to investigate late effects following stereotactic single fraction and small volume irradiation on cerebral blood flow and histologic alterations in the rat brain parenchyma.

MATERIAL AND METHODS

66 Copenhagen rats, separated into eleven groups of six animals each received single doses of 20, 30, 40, 50 and 100 Gy using a 15 MV linear accelerator. Six rats served as controls. Two cylindrical collimators of 2 mm and 3 mm aperture were used. The diameters of the spherical 80% isodose were 3.7 and 4.7 mm, respectively (Table 1). Irradiation was applied to a predefined area in the right frontal lobe. 19 months after irradiation local cerebral blood flow (LCBF) was measured by the autoradiographic method in one animal of each dose group between 20 and 50 Gy. 9 and 19 months after irradiation, half of the animals of each group were sacrificed for brain histology. All animals irradiated with 100 Gy were sacrificed 7 months after irradiation.

RESULTS

An increase of local cerebral blood flow was measured in brain structures within the 80% isodose in animals irradiated with 50 Gy (Figure 3) compared to the contralateral hemisphere. Measurements close to necrotic areas showed a strong decrease of local cerebral blood flow (Figure 1). A volume increase of the irradiated hemisphere was seen after 19 months (Figure 2). The histologic examination after 19 months showed necrotic areas in the 30-50 Gy groups (Figure 4b) but not in the 20 Gy groups (Figure 4c). The animals who received 100 Gy demonstrated brain necrosis within 9 months after irradiation (Figure 4a). At both points in time the groups irradiated with the 3-mm collimator showed more pronounced histomorphologic and functional changes compared to the groups irradiated with the 2-mm collimator.

CONCLUSION

Alterations of the local cerebral blood flow were measured as a late effect after single dose irradiation. The alterations of the local cerebral blood flow could be explained by the histomorphologic changes of the blood vessels. Using a semi-quantitative classification a dose, time and volume dependence for the endpoint radionecrosis was seen.

Three-dimensional accuracy and interfractional reproducibility of patient fixation and positioning using a stereotactic head mask system

PURPOSE

Conformal radiotherapy in the head and neck region requires precise and reproducible patient setup. The definition of safety margins around the clinical target volume has to take into account uncertainties of fixation and positioning. Data are presented to quantify the involved uncertainties for the system used.

METHODS AND MATERIALS

Interfractional reproducibility of fixation and positioning of a target point in the brain was evaluated by biplanar films. 118 film pairs obtained at 52 fractions in 4 patients were analyzed. The setup was verified at the actual treatment table position by diagnostic X-ray units aligned to the isocenter and by a stereotactic X-ray localization technique. The stereotactic coordinates of the treated isocenter, of fiducials on the mask, and of implanted internal markers within the patient were measured to determine systematic and random errors. The data are corrected for uncertainty of the localization method.

RESULTS

Displacements in target point positioning were 0.35 +/- 0.41 mm, 1.22 +/- 0.25 mm, and -0.74 +/- 0.32 mm in the x, y, and z direction, respectively. The reproducibility of the fixation of the patient's head within the mask was 0.48 mm (x), 0.67 mm (y), and 0.72 mm (z). Rotational uncertainties around an axis parallel to the x, y, and z axis were 0.72 degrees, 0.43 degrees, and 0.70 degrees, respectively. A simulation, based on the acquired data, yields a typical radial overall uncertainty for positioning and fixation of 1.80 +/- 0.60 mm.

CONCLUSIONS

The applied setup technique showed to be highly reproducible. The data suggest that for the applied technique, a safety margin between clinical and planning target volume of 1-2 mm along one axis is sufficient for a target at the base of skull.

Relation between carbon ion ranges and x-ray CT numbers

Measurements of carbon ion ranges in various phantom materials and real bones are presented. Together with measured Hounsfield values, an empirical relation between ranges and Hounsfield units is derived, which is an important prerequisite for treatment planning in carbon ion therapy.

Treatment planning for heavy ion radiotherapy: clinical implementation and application

The clinical implementation and application of a novel treatment planning system (TPS) for scanned ion beams is described, which is in clinical use for carbon ion treatments at the German heavy ion facility (GSI). All treatment plans are evaluated on the basis of biologically effective dose distributions. For therapy control, in-beam positron emission tomography (PET) and an online monitoring system for the beam intensity and position are used. The absence of a gantry restricts the treatment plans to horizontal beams. Most of the treatment plans consist of two nearly opposing lateral fields or sometimes orthogonal fields. In only a very few cases a single beam was used. For patients with very complex target volumes lateral and even distal field patching techniques were applied. Additional improvements can be achieved when the patient's head is fixed in a tilted position, in order to achieve sparing of the organs at risk. In order to test the stability of dose distributions in the case of patient misalignments we routinely simulate the effects of misalignments for patients with critical structures next to the target volume. The uncertainties in the range calculation are taken into account by a margin around the target volume of typically 2-3 mm, which can, however, be extended if the simulation demonstrates larger deviations. The novel TPS developed for scanned ion beams was introduced into clinical routine in December 1997 and was used for the treatment planning of 63 patients with head and neck tumours until July 2000. Planning strategies and methods were developed for this tumour location that facilitate the treatment of a larger number of patients with the scanned heavy ion beam in a clinical setting. Further developments aim towards a simultaneous optimization of the treatment field intensities and more effective procedures for the patient set-up. The results demonstrate that ion beams can be integrated into a clinical environment for treatment planning and delivery.

A method for determining the alignment accuracy of the treatment table axis at an isocentric irradiation facility

At an isocentric irradiation facility, the rotation axis of the treatment table has to be accurately aligned in vertical orientation to the isocentre, which is usually marked by three perpendicular laser planes. In particular, high precision radiotherapy techniques, such as radiosurgery or intensity modulated radiotherapy, require a higher alignment accuracy of the table axis than routinely specified by the manufacturers. A simple and efficient method is presented to measure the direction and the size of the displacement of the table axis from the isocentre as marked by the lasers. In addition, the inclination of the table axis against the vertical direction can be determined. The measured displacement and inclination provide the required data to correct for possible misalignments of the treatment table axis and to maintain its alignment. Measurements were performed over a period of two years for a treatment table located at the German heavy ion therapy facility. The mean radial distance between the table axis and the isocentre was found to be 0.25 +/- 0.25 mm. The mean inclination of the table axis in the XZ- and YZ-planes was measured to be -0.03 +/- 0.02 degrees and -0.04 +/- 0.01 degrees, respectively. The measurements demonstrate the good alignment of the treatment table over the analysed time period. The described method can be applied to any isocentric irradiation facility, especially including isocentric linear accelerators used for radiosurgery or other high precision irradiation techniques.

Determination of tolerance dose uncertainties and optimal design of dose response experiments with small animal numbers

BACKGROUND

Dose response experiments aim to determine the complication probability as a function of dose. Adjusting the parameters of the frequently used dose response model P(D) = 1/[1 + (D50/D)k] to the experimental data, 2 intuitive quantities are obtained: the tolerance dose D50 and the slope parameter k. For mathematical reasons, however, standard statistic software uses a different set of parameters. Therefore, the resulting fit parameters of the statistic software as well as their standard errors have to be transformed to obtain D50 and k as well as their standard errors.

MATERIAL AND METHODS

The influence of the number of dose levels on the uncertainty of the fit parameters is studied by a simulation for a fixed number of animals. For experiments with small animal numbers, statistical artifacts may prevent the determination of the standard errors of the fit parameters. Consequences on the design of dose response experiments are investigated.

RESULTS

Explicit formulas are presented, which allow to calculate the parameters D50 and k as well as their standard errors from the output of standard statistic software. The simulation shows, that the standard errors of the resulting parameters are independent of the number of dose levels, as long as the total number of animals involved in the experiment, remains constant.

CONCLUSION

Statistical artifacts in experiments containing small animal numbers may be prevented by an adequate design of the experiment. For this, it is suggested to select a higher number of dose levels, rather than using a higher number of animals per dose level.

Quality assurance for a treatment planning system in scanned ion beam therapy

Conformal radiation therapy using dynamic beam delivery systems like scanned ion beams requires concise quality assurance procedures for the complete treatment planning process. For the heavy ion therapy facility at GSI, Darmstadt, a quality assurance program for the treatment planning system (TPS) has been developed. It covers the development and updating of software, data protection and safety, and the application of soft- and hardware. The tests also apply to the geometrical precision of imaging devices and the geometrical and dosimetrical verification of dose distributions in different phantoms. The quality assurance program addresses acceptance and constancy tests of the treatment planning program. Results of the acceptance tests served as a basis for its governmental approval. Two main results of the acceptance tests are representative for the overall performance of the system. (1) The geometrical uncertainty that could be achieved for the target point definition, setup accuracy, field contouring, and field alignment is typically 1.5 mm. The uncertainty for the setup verification using digitally reconstructed radiographs (DRR's) is limited to 2 mm. (2) The mean deviations between measured and planned dose values is 3% for standardized cases in a water phantom and up to 6% for more complicated treatment configurations.

Quality management of medical physics issues at the German heavy ion therapy project

For the commissioning and operation of the German Heavy Ion Therapy Project a quality assurance program was developed and successfully applied. The complete radiotherapy process using heavy ions was carefully analyzed and divided into three areas related to beam delivery and control, safety-interlock system and medical physics issues. In this paper, the medical physics issues are addressed. Since the irradiation with heavy ions is a nonstandard modality, new concepts and ion specific tests were developed. As far as possible, national and international standard specifications for radiotherapy were adopted. For each aspect, a performance characteristic and a corresponding acceptance test were introduced. In addition, test characteristics for the constancy tests were established. For all tests, intervention thresholds and test frequencies were specified. Using the described protocol of acceptance tests, the commissioning was passed successfully. The heavy ion irradiation facility was approved by the governmental authorities on the basis of these test results. During clinical operation, constancy tests are performed at the beginning of each treatment period, in order to maintain the quality found during the acceptance tests. Up to now, 48 patients have been treated within 6 treatment periods of 4 weeks each. The concepts used and the tests developed for the quality assurance program may serve as an example of how to introduce systematically a quality assurance program for a new treatment modality.

Quantification of microglial late reaction to stereotactic irradiation of the rat brain using computer-aided image analysis

The role of microglial cells in the late delayed reaction following radiotherapy of brain tumors has not been elucidated. To investigate the late delayed response of microglial cells to radiation, we stereotactically irradiated spherical treatment volumes in the right frontal lobe of rat brains. Doses of 20, 30, 40, and 50 Gy were used in combination with two different collimators. The response of microglial cells at 10 and 19 months after irradiation was determined with Anti-CD 11 b/c (Ox 42) as an immunohistochemical marker. For evaluation of immunostaining, we developed a method using computer-aided image analysis in which the ratio of the area of stained cells to that of nonstained brain tissue is calculated. In addition, quantification of Ox-42+ cells per microscopic field was performed. Animals treated with 30 Gy or more had significantly increased total areas of staining at both time points studied. In contrast, the number of stained cells at 10 months increased significantly only in animals treated with 30 or 40 Gy. Likewise, at 19 months, this number increased significantly only in animals treated with 40 Gy or more. These results indicate that computer-aided determination of the area of stained cells is more sensitive than the counting of stained cells. We have demonstrated that microglial cells respond to stereotactic irradiation in a dose-dependent fashion. The image analysis we employed for this purpose is a systematic method to evaluate immunohistochemical staining.

A system for three-dimensional dosimetric verification of treatment plans in intensity-modulated radiotherapy with heavy ions

The introduction of dynamic intensity modulation into radiotherapy using conventional photon beams or scanning particle beams requires additional and efficient methods of dose verification. Dose measurements in dynamically generated dose distributions with a single ionization chamber require a complete application of the treatment field for each single measurement. Therefore measurements are performed by simultaneous use of multiple ionization chambers. The measurement is performed by a computer controlled system and is comprised of the following steps: (a) automated positioning of the ionization chambers, (b) measurement at these points, (c) a comparison with the calculated dose from the treatment planning system, and (d) documentation of the measurement. The ionization chambers are read out by a multichannel electrometer and are densely packed into a mounting of polymethylmetacrylate, which is attached to the arm of a three-dimensional motor-driven water phantom. The measured and planned dose values are displayed numerically as well as graphically. The mean deviation between measured and planned doses as well as their standard deviation are calculated and displayed. Through printouts complete documentation of the measurement is obtained and a quick decision can be made whether the dose distribution is acceptable for the patient. The system is now routinely used for dose verification at the heavy ion therapy project at the Gesellschaft für Schwerionenforschung in Darmstadt. Up to now 242 measurements have been performed for heavy ion treatment of 30 patients. The system allows efficient verification and documentation of carbon ion fields and is in principle also applicable to intensity-modulated photon beams.

Delayed vascular injury after single high-dose irradiation in the rat brain: histologic immunohistochemical, and angiographic studies

PURPOSE

To investigate structural and functional changes in rats after focal brain irradiation by using histologic, immunohistochemical, and angiographic methods.

MATERIALS AND METHODS

Sixty rats were irradiated stereotactically with photons from a 15-MeV linear accelerator. Two collimators and single doses ranging from 20 to 100 Gy were used to treat stereotactically defined areas of 3.7- and 4.7-mm cross section (80% isodose) in the right frontal lobe. The dose-response relationship for the end-point necrosis at 19 months revealed a mean tolerance dose (D50) of 34.2 Gy (standard errors: +4.1, -3.7 Gy). Histologic, immunohistochemical, and angiographic examinations were performed to evaluate delayed radiation effects.

RESULTS

All animals irradiated with 100 Gy developed radiation necrosis after 9 months. Microangiography and immunohistochemical fluorescence staining of the endothelial cells revealed dose-dependent vascular dilatation and rarefaction. Animals irradiated with 20-50 Gy showed no morphologic changes after 9 months. With irradiation of 30-50 Gy, histologic vascular changes that increased with dose were found after 19 months. At that time, no changes were detected after irradiation with 20 Gy with both field sizes and after irradiation with 30 Gy and the 2-mm collimator. Radiation-induced functional disturbances of the brain vasculature could be demonstrated by extravasation of contrast medium by using a microangiographic technique.

CONCLUSION

The observed effect had a definite dependence on dose, volume, and time after treatment.

Treatment planning for the heavy-ion facility at GSI

A new treatment planning program was developed for the heavy-ion therapy facility at GSI. In addition, a concise quality standard for treatment planning has been set up. It covers acceptance and constancy checks of all critical aspects in treatment planning. Dose verification measurements done during the commissioning phase show an overall good agreement with the treatment planning calculations.

Determination of water absorbed dose in a carbon ion beam using thimble ionization chambers

The method to measure absorbed dose to water in a field of carbon ions as applied for the heavy ion therapy project at the Heavy Ion Research Laboratory in Darmstadt (GSI), Germany, is described in detail. Thimble ionization chambers with a water absorbed calibration factor are applied. The dose obtained with this method was compared with that obtained at the heavy ion therapy facility HIMAC at the National Institute of Radiological Sciences in Chiba, Japan, using the Japanese code of practice. The agreement found was better than 1%. The combined uncertainty of the determination of absorbed dose to water was estimated to amount to 5%.

Dose-response relationship for late functional changes in the rat brain after radiosurgery evaluated by magnetic resonance imaging

PURPOSE

Only few quantitative data are available on late effects in the healthy brain after radiosurgery. An animal model can contribute to systematically investigate such late effects. Therefore, a model applying radiosurgery at the rat brain was established. A long-term (19 months) follow up study with 66 animals after radiosurgery was carried out.

METHODS AND MATERIALS

In 60 animals, an area in the frontal lobe of the brain was irradiated stereotactically with a 15 MV linac. Different doses of 20, 30, 40, 50, and 100 Gy with two field sizes (3.9 and 5.9 mm collimator) were selected, using the integrated logistic formula with input parameters from human brain. The induced alteration of the blood-brain barrier permeability was investigated by means of contrast enhanced magnetic resonance imaging.

RESULTS

A first intracranial signal enhancement was observed in one animal 160 days after irradiation with 100 Gy. Beginning at 5 months all animals in the two 100 Gy groups homogeneously showed contrast enhancement, but none of the other groups. This remained until 13 months after irradiation. The volume of contrast enhancement as well as the increase of signal intensity were different between the two 100 Gy groups. After 19 months, the animals irradiated with lower doses also showed contrast enhancements, although not uniformly throughout one group. A maximum likelihood fit of the logistic formula P(D) = 1/[1 + (D50/D)k] to the incidence of late effects for the 5.9 mm collimator at 19 months after irradiation results in the parameters D50 = 37.4(-5.2,+6.1) Gy and k = 4.7 +/- 2.4.

CONCLUSIONS

An animal model was established to study late normal brain tissue response. The observed late effects appeared very similar to the estimation of the integrated logistic formula for human brain. Based on these radiosurgery techniques, future experiments will focus on modifications in the irradiation modalities, i.e., irregular volumes, radiation quality or fractionation.

A system for stereotactic irradiation and magnetic resonance evaluations in the rat brain

PURPOSE

A stereotactic fixation and localization device developed for small animal stereotactic radiosurgery is described.

METHODS AND MATERIALS

Irradiated volumes of spherical shape down to 1.7 mm in diameter at the 80% isodose level are attainable. The fixation device can also be used for magnetic resonance imaging (MRI) and allows target localization during magnetic resonance (MR) image content measurement. The capabilities of the entire system were investigated using a phantom that permitted measurement and localization of the three-dimensional dose distribution. Localization of target isocenter coordinates in MR images was also checked with the phantom.

RESULTS

An overall spatial error of about 1 mm for subsequent stereotactic irradiation was obtained.

CONCLUSIONS

The accuracy of the fixation and localization techniques is adequate to investigate radiation-induced changes in the rat brain.