Implementation of Silicon Based Dosimeters, the Dose Magnifying
University of Wollongong Research Online
University of Wollongong Thesis Collection University of Wollongong Thesis Collections
2011 Implementation of silicon based dosimeters, the dose magnifying glass and magic plate for the dosimetry of modulated radiation therapy Jeannie Hsiu Ding Wong University of Wollongong
Recommended Citation Wong, Jeannie Hsiu Ding, Implementation of silicon based dosimeters, the dose magnifying glass and magic plate for the dosimetry of modulated radiation therapy, Doctor of Philosophy thesis, Centre for Medical Radiation Physics, Engineering Physics, University of Wollongong, 2011. http://ro.uow.edu.au/theses/3348
Research Online is the open access institutional repository for the University of Wollongong. For further information contact Manager Repository Services: [email protected].
IMPLEMENTATION OF SILICON BASED DOSIMETERS, THE DOSE MAGNIFYING GLASS AND MAGIC PLATE FOR THE DOSIMETRY OF MODULATED RADIATION THERAPY
A Thesis Submitted in Fulfilment of the Requirements for the Award of the Degree of
Doctor of Philosophy
from
UNIVERSITY OF WOLLONGONG
by
Jeannie Hsiu Ding Wong BBiomed. Eng., MMed.Phys.
Centre for Medical Radiation Physics, Engineering Physics Faculty of Engineering
2011 © Copyright 2011 by Jeannie Hsiu Ding Wong ALL RIGHTS RESERVED CERTIFICATION
I, Jeannie Hsiu Ding Wong, declare that this thesis, submitted in fulfilment of the requirements for the award of Doctor of Philosophy, in the Centre for Medical Radia- tion Physics, Engineering Physics, Faculty of Engineering, University of Wollongong, is wholly my own work unless otherwise referenced or acknowledged. The document has not been submitted for qualifications at any other academic institution.
(Signature Required) Jeannie Hsiu Ding Wong 30 June 2011 Table of Contents
ListofTables...... vi List of Figures/Illustrations ...... xii
ABSTRACT...... xiii Acknowledgements...... xv Contribution of Collaborators ...... xviii Publication List ...... xx Conferences...... xxi
1 Introduction 1 1.1 Project aim ...... 2
2 Literature review 4 2.1Cancerstatistics...... 4 2.2Radiationtreatmenttrend...... 4 2.3 Intensity Modulated Radiation Therapy ...... 6 2.3.1 StepandshootIMRTdelivery...... 6 2.3.2 Dynamic sliding window IMRT delivery ...... 7 2.4 Volumetric Modulated Arc Therapy ...... 7 2.5StereotacticRadiosurgery/Radiotherapy...... 8 2.6HelicalTomotherapy...... 8 2.7Currentqualityassuranceanddosimetricapproaches...... 9 2.7.1 Ionisationchamber...... 10 2.7.2 Filmdosimetry...... 10 2.7.3 Silicon diode ...... 11 2.7.4 Diamonddetectors...... 14 2.7.5 Geldosimetry...... 14 2.7.6 Electronicportalimagingdevice...... 15 2.8Highspatialresolutiondosimetry...... 16 2.8.1 Concept of silicon strip detector ...... 16 2.8.2 Application of high spatial resolution dosimeters in medical ra- diationtherapy...... 19 2.9 Concept of a two dimensional array detector ...... 20 2.9.1 Pixelateddetector...... 24 2.9.2 MAESTRO project ...... 25
i TABLE OF CONTENTS ii
3 Methodology 26 3.1DoseMagnifyingGlass...... 26 3.1.1 Designandfabrication...... 26 3.1.2 Detectorpackaging...... 28 3.2MagicPlate...... 28 3.2.1 Epitaxialdiodes...... 29 3.3TERAreadoutsystem...... 31 3.3.1 HowTERAworks...... 33 3.3.2 Amplitude and timing ...... 37 3.3.3 Charge collection in silicon strip detector ...... 38 3.4Filmdosimetry...... 39 3.4.1 Radiographic film ...... 39 3.4.2 Radiochromicfilm...... 40
4 Radiation response and basic characterisation of the Dose Magnifying Glass 43 4.1Introduction...... 43 4.2Materialsandmethods...... 44 4.2.1 Percentdepthdosemeasurement...... 44 4.2.2 Dose per pulse response measurement ...... 45 4.2.3 Stem effect measurement ...... 48 4.2.4 Dose linearity measurement ...... 49 4.2.5 Energyresponsemeasurement...... 49 4.2.6 Angularresponsemeasurement...... 50 4.3Results...... 52 4.3.1 Percentdepthdose...... 52 4.3.2 Dose per pulse response ...... 54 4.3.3 Stem effect ...... 58 4.3.4 Dose linearity ...... 58 4.3.5 Energyresponse...... 59 4.3.6 Angularresponse...... 61 4.4Conclusion...... 63
5 Application of the Dose Magnifying Glass in the dosimetric verifica- tion of an intensity modulated radiation therapy treatment delivery 65 5.1Introduction...... 65 5.2Materialsandmethods...... 66 5.2.1 DoseMagnifyingGlass...... 66 5.2.2 Uniformitymeasurement...... 67 5.2.3 Penumbraresponsemeasurement...... 67 5.2.4 Clinical application in IMRT fields ...... 68 5.3Resultsanddiscussions...... 69 5.3.1 Uniformity...... 69 5.3.2 Penumbrameasurement...... 69 TABLE OF CONTENTS iii
5.3.3 Clinical application in IMRT fields ...... 70 5.4Conclusion...... 76
6 Application of the Dose Magnifying Glass in the dosimetric verifica- tion of a stereotactic radiosurgery treatment delivery 78 6.1Introduction...... 78 6.2Materialsandmethods...... 82 6.2.1 DoseMagnifyingGlass...... 82 6.2.2 SRSphantom...... 82 6.2.3 Detectorrelativesensitivityfactormeasurement...... 83 6.2.4 Angulardependencecorrection...... 85 6.2.5 Determination of the center of rotation measurement ...... 86 6.2.6 SRSconeprofilesandtotalscatterfactormeasurement..... 87 6.2.7 Clinical stereotactic arc measurement ...... 88 6.3Resultsanddiscussions...... 89 6.3.1 Uniformity...... 89 6.3.2 Determination of center of rotation ...... 89 6.3.3 SRSconeprofilesandpenumbra...... 91 6.3.4 SRSconetotalscatterfactor...... 91 6.3.5 Clinical SRS application ...... 93 6.3.6 On the volume averaging effect of small field dosimetry ..... 94 6.4Conclusion...... 102
7 Application of the Dose Magnifying Glass in the quality assurance of Helical Tomotherapy 105 7.1Introduction...... 105 7.2Methodsandmaterials...... 107 7.2.1 DoseMagnifyingGlass...... 107 7.2.2 Multileaf collimator (MLC) alignment measurement ...... 107 7.2.3 Leafopentimemeasurement...... 109 7.2.4 Leaffluenceoutputfactormeasurement...... 111 7.3Resultsanddiscussions...... 113 7.3.1 Multileaf collimator alignment ...... 113 7.3.2 Leafopentimethreshold...... 113 7.3.3 Leaffluenceoutputfactor...... 118 7.4Conclusion...... 119
8 Radiation response and basic characterisation of the Magic Plate 121 8.1Introduction...... 121 8.2Materialsandmethods...... 122 8.2.1 PackagingoftheMagicPlate...... 122 8.2.2 Radiation damage studies ...... 123 8.2.3 Dose per pulse response measurement ...... 126 8.2.4 Percentdepthdosemeasurement...... 129 TABLE OF CONTENTS iv
8.2.5 Dose linearity measurement ...... 129 8.2.6 Energydependencemeasurement...... 130 8.2.7 Temperaturedependencemeasurement...... 130 8.2.8 Fieldsizedependencemeasurement...... 132 8.2.9 Angularresponsemeasurement...... 133 8.2.10Beamperturbationmeasurements...... 135 8.3Results...... 137 8.3.1 Radiationdamage...... 137 8.3.2 Dose per pulse response ...... 139 8.3.3 Percentdepthdose...... 146 8.3.4 Dose linearity ...... 146 8.3.5 Energydependence...... 148 8.3.6 Temperaturedependence...... 150 8.3.7 Fieldsizedependence...... 152 8.3.8 Angulardependence...... 153 8.3.9 Beamperturbationstudy...... 154 8.4Conclusion...... 156
9 Investigation of the Magic Plate in clinical application 158 9.1Introduction...... 158 9.2Materialandmethods...... 161 9.2.1 Uniformity and absolute dose calibration ...... 161 9.2.2 Dosemeasurementinsolidwaterphantom...... 164 9.2.3 Fluencemeasurementsintransmissionmode...... 165 9.2.4 Performance index used for the dose distribution comparison . . 166 9.3Resultsanddiscussion...... 166 9.3.1 Uniformitycorrection...... 166 9.3.2 Dose distribution comparison ...... 167 9.4Conclusion...... 174
10 Conclusion 176 10.1DoseMagnifyingGlass...... 176 10.2MagicPlate...... 179
Appendices 183
A SolidWorks drawing 183 A.1 The 2nd generationDoseMagnifyingGlassdetectorholder...... 183 A.2SRSphantom...... 183 A.3MagicPlateY-shapeperspexframe...... 186 A.4MagicPlateholdertobeusedwiththeI’mRTphantom...... 186 TABLE OF CONTENTS v
B Matlab Scripts 189 B.1Uniformitycorrectionscript...... 189 B.2Gammaanalysisscript...... 192
References 219 List of Tables
2.1Specificationofcommercial2Darrays...... 22
3.1 Set reference values for the CMRP TERA 03 board...... 35
5.1 Comparison of the 80-20% penumbra width measurements for a 6 MV beam at 1.5 cm and 10.0 cm depth between the DMG, Gafchromic EBT films, and other published literature. The uncertainty reported represent the 95% CI of the mean for three sets of measurements. . . . 71
6.1 The penumbra and FWHM measurements of the dose profiles of 5 mm to 20 mm cone diameter at SDD 100 cm, comparing the DMG and EBT2measurements...... 92 6.2 Dose calculated based on the analytical model showing the effect of dose averaging effect with different detector sizes...... 98
7.1 Programmed leaf opening configuration for LFOF measurement for LOI = 33. The same delivery sinogram was also used with leaf 47 and 62 as the LOI. The symbol “|”denotesanopenleaf...... 112 7.2 LFOF measurements by DMG and tomo detectors. The uncertainties represent 1 s.d. of the mean...... 119
8.1 Measurement setup for the dose per pulse measurement...... 126 8.2Specificationsofthedosimetersused...... 128 8.3Energydependencemeasurementconfigurations...... 131 8.4 Dose per pulse dependence sections...... 141 8.5 Definition of symbols...... 142 8.6 Difference between the surface doses measured with open fields and a MP fields in percent. The values in the parentheses are the measured surfacedosenormalisedtothemaximumdose...... 155
9.1 MP positioning for the uniformity calibration of the diodes...... 163 9.2 Gamma analysis (3% DD,3 mm DTA) of IMRT dose distributions com- paringTPScalculations,MPandEBT2measurements...... 172
vi List of Figures
2.1 Silicon/water stopping power ratio...... 13 2.2 Schematic of a silicon p-n junction diode (Shi et al., 2003)...... 13 2.3 Operation principle of a silicon strip detector [modified from (Damerell, 1995)]...... 18
3.1SchematicdiagramoftheDoseMagnifyingGlass(DMG)...... 27 3.2 Prototype version of DMG (left) and the 2nd generation DMG (right). . 28 3.3TheMagicPlate(MP)2Darraydetector...... 29 3.4 The Magic Plate mounted on the linac gantry. The mounting plate was a modified TBI applicator...... 30 3.5Schematicdiagramofanepitaxialdiode(withguardring)...... 32 3.6Schematicdiagramofanepitaxialdiode(withoutguardring)...... 32 3.7 Schematic diagram of the TERA readout [modified from (Mazza et al., 2005)]...... 33 3.8 Charge subtraction waveform (Mazza et al., 2005)...... 35 3.9 Physical structure of Gafchromic (a) EBT and (b) EBT2 film (ISP, 2006, 2009)...... 42
4.1SketchofthesolidwaterencapsulationfortheprototypeDMG..... 45 4.2 The 2nd generation DMG in the solid water holder and the custom made phantom which has a cylindrical body and hemispherical head...... 53 4.3 Definition of polar and azimuthal angle...... 53 4.4 Depth dose curve for a 10 × 10 cm2 field size of a 6 MV photon beam comparing DMG with Farmer ion chamber (NE-2571) measurements in solidwater...... 54 4.5 Dose per pulse response of selected DMGs with different resistivity and preirradiation conditions for the dose per pulse range of 9.45 × 10−5 to 2.72 × 10−4 Gy/pulse. The error bars represent the 1 standard deviationofthemeanofallthechannelsinthestripdetector...... 55 4.6 Dose per pulse response for a low resistivity, preirradiated device for the dose per pulse range of 6.42 × 10−7 to 2.92 × 10−4 Gy/pulse. The dose per pulse responses were normalised to the dose per pulse of 2.2 × 10−4 Gy/pulse...... 57 4.7 Dose per pulse corrected response...... 58
vii LIST OF FIGURES viii
4.8 Dose linearity of the DMG...... 59 4.9 Energy response of the strip detector, normalized to 1 at a 6 MV photon energy. Error bars represent the combined uncertainties of the 95% CI of the mean for three sets of measurements and the 0.5% calibration uncertainty...... 60 4.10 Angular response correction factor generated as the measured dose by DMG/CC13 ion chamber measurements. As the angle of the incident beam increased from perpendicular beam (θ =0◦)toparallelbeam(θ = 90◦), the strip detector under-responded compared to the dose measured by CC13 ion chamber. The maximum angular response 28.1% ± 0.1% was found to be at the parallel incident beam at the gantry angle of 90◦. The difference between the response of the detector channels at the position of ± 1 cm away from the central axis (mid channels) were <1%...... 62 4.11 Angular response for the 2nd generationDMG...... 63
5.1 Detector response before and after uniformity correction. The error bars represent the 95% CI of the mean for three sets of measurements. 70 5.2 Dose profiles showing the penumbra of a 6 MV beam at 1.5 cm depth for the secondary X-jaw and the rounded leaf ends of a multileaf collimator (MLC) using the strip detector and Gafchromic EBT films. The error bars for the EBT profiles represent the 95% CI of the mean of three sets of measurements while the error bars for the DMG measurements represents the 2% reproducibility uncertainty...... 71 5.3 IMRT dose profiles for each gantry angle and the total summation of all the beams. The error bars represent the reproducibility uncertainty ofthemeasurementsatthelevelof95%CIofthemean...... 72 5.4 Comparison between Pinnacle predicted dose profiles, measurements with Gafchromic EBT film, and DMG. The error bars represent the reproducibility of the measurements which are in the order of ± 2% for DMG and 3% for EBT film, respectively. On average, the difference between the dose measured using the DMG with the dose points from the Pinnacle dose matrix and the dose measured with the EBT film were 1.1%± 1.8% (1 s.d.) and 1.0%± 1.6% (1 s.d.)...... 73 5.5 (a) 3D dose profiles of a single step and shoot IMRT field, depicting the modulation of the dose during each segment, within the width of 2.56 cm. (b) Dose profiles of a full IMRT delivery, each gantry angle and modulation in each IMRT segments can be shown temporally. The measurements were performed with the acquisition pulse width of 0.1 s. (c) The temporal cumulative 3D dose profiles of a full IMRT delivery. (d) Temporal dose rate pattern of an IMRT delivery for channel no. 8 and channel no. 127...... 74 LIST OF FIGURES ix
6.1 SRS phantom mounted on the Radionics SRS couch mount. A plate with flat top and matching concave bottom was used to achieve a flat calibration phantom. Solid water spacers were used to allow the lateral shiftingoftheDMGwithinthephantom...... 83 6.2 Detector response of the DMG before (open circle) and after uniformity correction(closedcircle)...... 89 6.3 Two profiles of a 5 mm diameter cone measured with the DMG showing discrepancy in the center of rotation due to couch rotation from 90◦ - 270◦...... 91 6.4 SRS cones total scatter factor comparing DMG with EBT2, CC04 mea- surements and Monte Carlo calculation. The measurements were taken within the spatial distance of ± 0.2 mm from the center axis. The CC04 ion chamber was used to measured the Scp for the cones with 15 mm - 40 mm. The Monte Carlo calculation was made with voxel size of 1 × 1 × 2mm3 for cone diameter ≤ 10 mm while voxel size of 2 × 2 × 2 mm3 was used to calculate for the cone diameter >10cm...... 93 6.5 SRS arc (0◦ - 180◦)relativeintensityprofile...... 94 6.6 (a) Pseudo colored image from EBT2 film measurement of a 5 mm diameter beam, normalized to the center of the beam, (b) isodose lines for the 5 mm diameter beam, (c) signal profile alone a horizontal line drawn across the middle of the beam and (d) zoom in to the top beam profile, data points >0.3 were extracted and curve fitted with a 4th degree polynomial function (r2 = 0.998)...... 95 6.7 Schematic drawing of detector area superimposed on the radial beam contour...... 96 6.8 Schematic diagram representing the top view of the beam and the deno- tation of the r- and the τ-directions. The isodose lines are represented astheconcentriccircles...... 99 1 6.9 Dose gradient in (a) r- and (b) τ-direction for a detector with 2 longitu- dinal length of b = 0.5 mm and 1 mm. The x-axis is the radial direction extending from the center of the beam (r = 0 mm) to the edge of the beam (r = - 6 mm). The negative gradient at the region r = -0.5 mm to 0 mm is due to the noise in the EBT2 film picked up by the fitted curve...... 101 6.10 Dose map of (a) 2 × 10 mm2 and (b) 1 × 10 mm2...... 102 6.11 Numerical integration of dose for detector of 0.2 × 2mm2 and 1 × 1 mm2 area...... 103 LIST OF FIGURES x
7.1 (a) EDR2 film measured leaf alignment test. The distances from the two outside peaks to the central peak were X1 = 24.94 ± 0.09 mm and X2 = 23.53 ± 0.09 mm. The MLC alignment error was taken as half the difference between these two distances (dX = 0.71 ± 0.09 mm) and (b) DMG measured leaf alignment test. Half the distance between the centre of each profile gives the MLC alignment error (dX = 0.55 ± 0.10 mm)...... 114 7.2 Actual measured leaf open times plotted against programmed leaf open times for (a) 50 ms, (b) 100 ms, (c) 200 ms and (d) 303 ms projec- tion time. A linear curve has been fitted for guidance. The error bars represent 1 s.d. of the mean of the five projections within the same modulation bracket...... 115 7.3 Comparison of the DMG and MVCT measured leaf open time for leaf 32 (a) 50 ms, (b) 100 ms, (c) 200 ms and (d) 303 ms projection time. A linear curve has been fitted for guidance...... 116 7.4 DMG measurement of the leaf opening and closing for a 200 ms projec- tion time. Each row of the ribbon represents a 3 ms acquisition. Red arrowspointtothedirectionofleaftravel...... 118
8.1 MP packaging: (a) MP mounted on the Y-shaped Perspex frame. The readout electronics box is connected to the MP on the left. (b) MP sandwiched between two 5 mm solid water slabs for measurements in theI’mRTphantom...... 123 8.2 Schematic diagram of the MP packaging for the use in (a) transmission modeand(b)inanI’mRTphantom...... 124 8.3 Modified accessory tray to position the Pb attenuator closed to linac gantry...... 127 8.4 Schematic drawing of the ‘face-up’ and ‘face-down’ configurations for theMPdiode...... 130 8.5Experimentsetupforthetemperaturedependencemeasurement.....132 8.6 Sensitivity of the detector as a function of the accumulated radiation dose. The error bars represent the ± 1 standard deviation of the chan- nels located within the central 20-80% region of the 11 × 11 MP detec- tor. The left axis shows the detector sensitivity expressed as counts/1Gy while the right axis shows the detector sensitivity normalised to the first irradiationof1Gy...... 138 8.7 Reproducibility of the MP. The error bars represent ± 1 standard de- viationofthemeanofthreemeasurements...... 138 8.8 Dose per pulse response for the MP and the commercial diode, nor- malised to the dose per pulse of 2.78 × 10−4 Gy/pulse...... 140 8.9 Division of the dose per pulse measurements into four (A to D) sections. 140 8.10 Minority carrier lifetime as a function of ionising radiation dose rate 16 −3 in a highly doped n-type silicon (NA =10 cm ) [reproduced from (Alexander, 2003)]...... 143 LIST OF FIGURES xi
8.11 Depth dose curve for a 10 × 10 cm2 field size of a 6 MV photon energy measuredwithaCC13ionchamberandtheMP...... 147 8.12 Build up region of the depth dose curve, comparing the MP, CC13 and Attix chamber measurements (courtesy of Dr. Bradley Oborn, ICCC andCMRP)...... 147 8.13 The dose linearity of the MP...... 148 8.14 The energy dependence of the MP, measured with four different se- tups; comparing ‘face-up’/‘face-down’ configurations. The measure- ments were made in free air geometry and with a solid water phantom. 150 8.15 The energy response of the MP in ‘face-up’ and ‘face-down’ configu- rations for 6 and 10 MV photon energies. The detector response was normalised to 1 at the reference setup of SSD 100 cm and dmax of 1.5 cm and 2.0 cm for the 6 and 10 MV photon energies, respectively. . . . 151 8.16 Sensitivity variation with temperature (svwt) of the MP. The detector response was normalized to 1 at 23.32◦C. Error bars represent 1 s.d. of themeanofthreemeasurements...... 151 8.17 The field size dependence of MP at SDD of 101.5 and 58 cm compared with the standard data measured using the Farmer chamber measured at SDD of 101.5 cm. The data was normalised to 1 at the field size of 10 × 10 cm2...... 153 8.18 Mean angular response of the MP for the detectors located at the central column (0 cm) and at ± 5cmoffaxisdistance...... 155 8.19 Dose at the build up region for a 20 × 20 cm2 and 30 × 302 cm at (a) SSD of 90 cm and (b) SSD of 80 cm. The measurement uncertainty was ± 0.5%, representing 1 standard deviation of the mean of three measurements...... 156
9.1 Sketch of the MP coordinate system. The diode at coordinate [0,0] is thecentrediode...... 162 9.2 Mean cross plane and in plane profiles of the MP diodes before and after the uniformity correction. The error bars represent the ± 1 standard deviation of the mean of 11 detectors located within the same row or column...... 167 9.3 Dose measurement for a single IMRT field delivered with gantry angle set to 0◦. (I.a) shows the normalised cross plane dose profile at an off axis distance of 0 mm. (I.b to d) shows the three planar dose distribu- tions of the MP, EBT2 and the Pinnacle TPS. Gamma analysis of the (I.e) TPS versus MP (pass rate = 97.5%), and (I.f) TPS versus EBT2 (pass rate = 92.4%) dose distributions. (II) shows the horizontal and (III)verticaldoseprofilesforthethreedatasets...... 169 LIST OF FIGURES xii
9.4 Dose measurement for a composite IMRT plan delivered with the actual treatment gantry angles. (I.a) shows the normalised cross plane dose profile at an off axis distance of 0 mm. (I.b to d) shows the three planar dose distributions of the MP, EBT2 and the Pinnacle TPS. The MP measurements were corrected for angular dependence. Gamma analysis of the (I.e) TPS versus MP (pass rate = 80.2%), and (I.f) TPS versus EBT2 (pass rate = 82%) dose distributions. (II) shows the horizontal and(III)verticaldoseprofilesforthethreedatasets...... 171
A.1 The 2nd generationDoseMagnifyingGlassdetectorholder...... 184 A.2SRSphantom...... 185 A.3MagicPlateperspexframe...... 187 A.4SolidwaterplatesusedwithanI’mRTphantom...... 188
B.1 Matlab GUI for processing and comparing the three data sets, (i) Magic Plate, (ii) EBT2 film and (iii) Pinnacle predicted and (iv) Monte Carlo simulated dose distributions...... 190 Implementation of silicon based dosimeters, the Dose
Magnifying Glass and Magic Plate for the dosimetry of
modulated radiation therapy
Jeannie Hsiu Ding Wong
A Thesis for Doctor of Philosophy
Centre for Medical Radiation Physics, Engineering Physics University of Wollongong
ABSTRACT
New cutting edge radiation therapy techniques such as Intensity Modulated Radiation Therapy (IMRT), Stereotactic Radiosurgery (SRS), Helical TomoTherapy and most recently Volumetric Modulated Arc Therapy (VMAT) produce radiation dose maps with high dose modulation and tight gradients between the high and low dose region. Difficulties in the dosimetric verification of these new complex treatment methods us- ing existing dosimeters have led to the need for a new generation of fast responding real time dosimeters with submillimetre precision. This thesis describes two detector systems based on silicon substrates that were developed to address this need. The first detector system was a silicon strip detector called the “Dose Magnifying Glass” (DMG). It consisted of 128 detectors spaced 0.2 mm apart. It was coupled to a TERA ASIC chip that enabled simultaneous readout of multiple channels at high temporal resolution. The first part of the thesis involved investigation of the basic character- istics of the DMG followed by its application in the dosimetric verification of IMRT, SRS and Helical TomoTherapy treatment deliveries. The high spatial resolution of this device was ideal for the measurement of high dose gradients in IMRT and small fields encountered in SRS. When compared with film dosimetry, DMG measurements showed agreements within 3% for a SRS treatment plan. The DMG was also success- fully employed as an independent quality assurance tool for the verification of helical Tomotherapy machine binary MLC leaf parameters. The second detector was a two dimensional array detector so named the “Magic Plate” (MP). The diode was based on epitaxial technology and has a very thin sensitive volume of 50 μm. The MP comprised of 11 × 11 epitaxial diodes mounted on a 0.6 mm thick Kapton substrate. This detector was designed to be used either as a transmission detector or to measure dose distributions in a solid water phantom. Preliminary testing of the MP in a clin- ical IMRT treatment delivery was carried out. The MP measurements demonstrated good agreement (>80%) with conventional EBT2 film dosimetry and with treatment planning system predicted dose distributions using 3%/3mm gamma criteria.
KEYWORDS: IMRT, small field dosimetry, Dose Magnifying Glass, Magic Plate, Tomotherapy, high spatial resolution silicon detector, transmission detector, quality assurance Acknowledgements
On the 31st December 2007, I boarded a flight to Australia, embarked on the journey of a PhD study. For the next 3.5 years, I have had the opportunity to learn so much and met so many wonderful people. It was a wonderful and enriching 3.5 years, indeed. First of all, I would like to thank my three supervisors; Professor Anatoly Rosen- feld, Professor Peter Metcalfe and Dr. Martin Carolan, for without their guidance and advice, this thesis would not be possible. Professor Anatoly Rosenfeld is a wonderful mentor and supervisor. Despite having many other students, he still find time to pull me out of my research portholes. His vast knowledge and experience in physics and solid state dosimetry is inspirational and am humbled by his patience to teach and explain to me the basic every time I ask a question. I also thank him for giving me the opportunities to extend my professional horizon in international research collab- orations and the financial support that goes with it. I would like to thank Professor Peter Metcalfe for his advice and help with my thesis writing. Professor Peter Met- calfe and Dr. Martin Carolan, with their extensive experience in the clinical medical physics provided sound advice to me. I would like to thank Dr. Martin Carolan for the many late nights that he stayed back at the Illawarra Cancer Care Centre (ICCC), supervising and assisting me with the measurements. His generosity and kindness in giving his time and energy towards nurturing young researcher is humbling. It was a priviledge to have worked and learned under these three great people. I would like to thank Dr. Michael Lerch and Dr. Marco Petesecca for their kind
xv assistance and advice in the electronic readout system and explanation of the basic physics of solid state detectors. I thank them, as well as Dean, Sami, Iolanda and Paul for giving me lifts with the equipment to the ICCC during those measurements days. Special thanks to Iolanda and Paul who had worked with me in some of my measurements. In my work, I often have to make special phantoms or holders for the detectors. I was so glad the UoW workshop was there to help me out with this. Thank you to Ron Marshall, Martin, Stuart, Keith, Andrew and Doug for their help in making all the big and small devices that I needed. I would also like to thank Peter Ihnat for the technical support. Special thanks to Terry Braddock for his friendship and ever ready help whenever I come and look for him. I would also like to thank Karen Ford for the administrative assistance and for her amazing organizing skills. I would like to thank Dr. Martin Butson and Dr. Matthew Williams for the useful advice in the use of the EBT and EBT2 films. Thank you to Jo and Ab for helping out and staying with us during some of our measurement evenings. I would like to thank Dr. Tony Knittel for the experience working with him in the SRS project. I thank him for his generosity with his time, spending the many weekends working with me at the Prince of Wales hospital, Randwick, even until the late hours. His vast knowledge, experience and hands on approach taught me a lot of things, not only in medical physics, but in conducting research. I will always remember the advice given to me “Change only ONE variable at a time!” To the physicist at POWH: thank you to Soo Ming for bunking me for a night, Dean Inwood for helping with one of the measurement and Carl Chan for the Monte Carlo calculation of the
SRS cone Scp values. During my short trip to the University of Wisconsin-Madison, WI, I had the op- portunity to work with Dr. Nick Hardcastle, Professor Wolfgang Tom´e, Ranjini To-
xvi lakanahalli and Dr. Adam Bayliss. I thank all of them for the help and advice given to me, leading to the success of this project. I also thank Dr. Peter Hoban for the tour of TomoTherapy Inc. and the TomoTherapy engineers for the feedback on the MLC leaf open time. To my fellow office mates: Thank you Amir for teaching me Geant4, although I do hope one day I can put it to use. Sianne, Norlaili, Elise, Amir, Ashley and Amy, thanks for the friendship, cakes and banter. To my family: Dad, mom, Ah Ma, Mama, Evelyn, James and Linda, thank you for your encouragement, support, prayers and for just being my family. To my fianc´e, Kong Sih Ying, I am ever so glad that I found you. Thank you for your support and for being there for me. I was never alone because I had you. To my friends outside the CMRP, my housemates, the Verbum Dei community, my church choir members, it was a wonderful 3.5 years in Wollongong and you guys have made my life so much richer and memorable. Special thanks to Chris for proof reading part of my thesis.
xvii Contribution of Collaborators
Professor Anatoly Rosenfeld is the inventor of the Dose Magnifying Glass (DMG) and the Magic Plate (MP). He also advised on the experiments design and results discussion. Professor Peter Metcalfe provided advice on the clinical aspects of the experiments and writing of the thesis. Dr. Martin Carolan provided supervision and assistance with the measurements conducted at the Illawarra Cancer Care Centre (ICCC), Wollongong, advised on the experiment designs and the writing of the thesis. Dr. Bradley Oborn provided the build up region depth dose data for a 10 × 10 cm2 field of a 6 MV photon beam. Dr. Michael Lerch and Dr. Marco Petasecca provided technical advice and assis- tance on the electronic readout system. Mr. Sam Khanna wrote the initial Labview software that was used for the data acquisition. Dr. V. Perevertaylo manufactured the DMG and MP detectors. Dr. Dean Cutajar provided advice and assistance on the use of DMG in brachytherapy. Miss Iolanda Fuduli assisted in the electronic readout and the measurements of the Magic Plate work. Mr. Paul Geenty assisted with the Magic Plate measurements. Professor Wolfgang Tom´e, Dr. Nicholas Hardcastle, Dr. Adam Bayliss and Ms Ranjini Tolakanahalli provided assistance with experimental design and data analysis in the Tomotherapy project. Dr. Tony Knittel provided advice and technical assistance with the design and machining of the SRS phantom, assisted in the measurement, and result discussion
xviii in the stereotactic radiosurgery project. Dr. Simon Downes and Dr. Michael Jack- son facilitated the experimental sessions at the Prince of Wales Hospital, Randwick, Sydney.
xix Publications
Wong,J.H.D., Carolan, M., Lerch, M. L. F., Petasecca, M., Khanna, S., Perever- taylo, V. L., Metcalfe, P. and Rosenfeld, A. B. (2010). A silicon strip detector dose magnifying glass for IMRT dosimetry. Medical Physics 37(2): 427-439.
Wong,J.H.D., Hardcastle, N., Tome, W. A., Bayliss, A., Tolakanahalli, R., Lerch, M. L. F., Petasecca, M., Carolan, M., Metcalfe, P. and Rosenfeld, A. B. (2011). Inde- pendent quality assurance of a helical tomotherapy machine using the dose magnifying glass. Medical Physics 38(4): 2256-2264.
Wong,J.H.D., Knittel, T., Downes, S., Carolan, M., Lerch, M. L. F., Petasecca, M., Perevertaylo, V. L., Metcalfe, P., Jackson, M. and Rosenfeld, A. B. (2011). The use of a silicon strip detector dose magnifying glass in stereotactic radiotherapy QA and dosimetry. Medical Physics 38(3): 1226-1238.
Wong,J.H.D., Cutajar, D., Lerch, M. L. F., Petasecca, M., Knittel, T., Carolan, M., Perevertaylo, V. L., Metcalfe, P. and Rosenfeld, A. B. (2011). From HEP to Med- ical Radiation Dosimetry - the silicon strip detector Dose Magnifying Glass. Radiation Measurements. DOI:10.1016/j.radmeas.2011.06.031(Accepted, 14 June 2011).
xx Conferences
Wong,J.H.D., Carolan, M., Lerch, M. L. F., Petasecca, M., Khanna, S., Perever- taylo, V. L., Metcalfe, P. and Rosenfeld, A. B. (2010). A Silicon Strip Detector Dose Magnifying Glass for IMRT Dosimetry, Abstracts in: EPSM-ABEC 2009: Engineering and Physical Sciences in Medicine & the Australian Biomedical Engineering College Conference, 8-12 November 2009, Hotel Realm, Canberra, ACT.” Australas Phys Eng Sci Med. 33(1): 68.
Wong,J.H.D., Carolan, M., Lerch, M. L. F., Petasecca, M., Khanna, S., Perever- taylo, V. L., Metcalfe, P. and Rosenfeld, A. B. (2009). A Silicon Strip Detector Dose Magnifying Glass For IMRT Dosimetry. Poster presented at the AIP Physics in In- dustry Day 2009. CSIRO Materials Science and Engineering, NSW, Australia.
Wong,J.H.D., Cutajar, D., Lerch, M. L. F., Petasecca, M., Knittel, T., Carolan, M., Perevertaylo, V. L., Metcalfe, P. and Rosenfeld, A. B. (2010)Silicon strip detector dose magnifying glass on stereotactic radiotherapy QA and dosimetry. Presented at the Solid State Dosimetry 16th International Conference, Sydney, 19-24 September 2010.
Wong,J.H.D., Hardcastle, N., Tome, W. A., Bayliss, A., Tolakanahalli and Rosen- feld, A. B. (2010). Independent quality assurance of a helical tomotherapy machine using the dose magnifying glass. Presented at the Engineering and Physical Sciences in Medicine and the Australian Biomedical Engineering Conference (EPSM-ABEC) 2010, Melbourne, Australia, 5 - 9 December 2010.
Wong,J.H.D., P Geenty, M Carolan, M L F Lerch, M Petasecca, V Perevertaylo, P Metcalfe, and A B Rosenfeld (2010), A novel 2D array silicon detector ‘Magic Plate” for IMRT dose verification. Poster presented at the Engineering and Physical Sciences in Medicine and the Australian Biomedical Engineering Conference (EPSM- ABEC) 2010, Melbourne, Australia, 5 - 9 December 2010.
Wong,J.H.D., M Carolan, I Fuduli, M L F Lerch, M Petasecca, P Metcalfe, and A B Rosenfeld (2010), Development of a 2D array silicon detector Magic Plate for the dosimetric verification of IMRT treatment delivery. Abstract accepted for the Engi- neering and Physical Sciences in Medicine and the Australian Biomedical Engineering
xxi Conference (EPSM-ABEC) 2011, Darwin, Australia, 14 - 18 August 2011.
xxii 1 Chapter 1
2 Introduction
3 The advancement of the treatment of cancer patients with radiation therapy is driven
4 by the intention to achieve higher tumour control and fewer side effects to critical
5 organs. Historical studies show a higher dose enhances tumour control and a lesser
6 dose reduces normal tissue complication (Pollack et al., 2000). In order to provide a
7 higher targeted dose to tumour and a lesser dose to critical organs, modern radiation
8 therapy delivery techniques are continually evolving and becoming more complex.
9 Delivery techniques such as intensity modulated radiation therapy (IMRT) challenge
10 the conventional mindset of boxed 3D conformal radiation therapy, while the issues
11 pertaining to the verification of small field dosimetry are still not completely resolved.
12 Conventional quality assurance (QA) approaches to dosimetry are used with vary-
13 ing degrees of success. Their adequacy in addressing the new dosimetric challenges
14 presented by modulated and small field radiation therapy delivery techniques in par-
15 ticular, lend themselves to development of new detector systems that complement the
16 existing detectors. In particular it is hoped these new detector systems will provide
17 higher spatial and temporal resolution than the existing dosimeters. This has led to the
18 development of new purpose-specific dosimeters to address the various requirements
19 of dosimetric verification.
1 1.1. Project aim 2
20 The technological advancement of the medical radiation physics tools has been
21 complemented somewhat by the rapid advancement in semiconductor and computer
22 technology. This enabled the development of device control systems such as the dy-
23 namically driven multi leaf collimators (MLCs). The same technology evolution has
24 also helped push the frontiers of semiconductor detectors.
25 1.1 Project aim
26 This thesis describes the development and characterisation of two new detector systems
27 based on silicon substrates. Chapter 3 describes the design and fabrication of these
28 two systems.
29 The thesis then describes the use of a silicon strip detector as a radiation dosime-
30 ter. The device is referred to in this thesis as the “Dose Magnifying Glass” (DMG).
31 This device consists of a linear array of silicon strips with high spatial resolution. It
32 was coupled to an electronic readout system capable of high temporal resolution ac-
33 quisition. This device has been successfully employed and characterised under simple
34 open medical linear accelerator x-ray fields (chapter 4). The device was then used to
35 characterise IMRT beams (chapter 5). It was subsequently used to characterise the
36 spatial and temporal properties of extremely small treatment fields as used in stereo-
37 tactic radiosurgery (SRS), see chapter 6. Then the DMG was used to verify temporal
38 and spatial properties of extremely fast moving binary MLC systems used in a helical
39 tomotherapy system, see (chapter 7).
40 A second detector system which is a two dimensional array detector so named
41 the “Magic Plate” (MP) is then described in subsequent chapters. This detector
42 prototype has been designed to be used as a transmission detector to verify dose from
43 IMRT and volumetric modulated arc therapy (VMAT) treatment delivery. The basic
44 characterisation of this system under simple open field linac x-ray beams is described 1.1. Project aim 3
45 in chapter 8. The clinical implementation and testing under more complex clinical
46 IMRT treatment beams is described in chapter 9.
47 Chapter 10 presents the summary of the main findings of this thesis. The ad-
48 vantages, limitations and the future direction of the development of these two novel
49 detector system prototypes are also discussed. 50 Chapter 2
51 Literature review
52 2.1 Cancer statistics
53 According to a report by the Australian Institute of Health and Welfare & Australasian
54 Association of Cancer Registries (AIHW & AACR, 2010), the highest number of re-
55 ported cancer incidence is prostate cancer (19403 cases per annum), followed by breast
56 cancer (12670 cases) for the Australian male and female population, respectively. This
57 is followed by melanoma of the skin (10340 cases) and lung cancer (9703 cases). One
58 in every ten deaths in Australia is caused by cancer, second only to cardiovascular
59 diseases. However, between 1982 and 2007, rate of deaths decreased for most cancer
60 sites. The fall in the mortality rate is attributed to improved treatment and early
61 detection of cancer.
62 2.2 Radiation treatment trend
63 Radiation therapy has long been recognised as an effective method in the treatment of
64 cancer. Radiation can be delivered to patients in three ways, externally using medical
65 linear accelerators, also called external beam radiation therapy (EBRT), internally via
4 2.2. Radiation treatment trend 5
66 brachytherapy, or systemically using radioisotope therapy.
67 Since the invention of the computed tomography (CT) scanner by Sir Godfrey
68 Hounsfield in 1971, this has been incorporated as an essential part of the treatment
69 planning work chain allowing three dimensional radiation therapy treatment tech-
70 niques (Battista et al., 1980). The use of blocking material to create irregularly shaped
71 fields and subsequently the use of multileaf collimators (MLC) enabled the pathway
72 to 3D conformal radiation therapy. The idea of 3D conformal radiation therapy is
73 essentially a forward planning method whereby uniform radiation fields are delivered
74 from multiple angles to the planning treatment volume (PTV), guided by the 3D CT
75 images, providing a homogeneous dose to the target while sparing surrounding nor-
76 mal tissues. This technique was largely successful in most cancer sites and reduced
77 radiation toxicity to the surrounding normal tissues.
78 In 1982, Brahme et al. (1982) published their paper which was later regarded as the
79 conceptual paper that introduced the concept of a treatment technique we now know
80 as intensity modulated radiation therapy (IMRT). They presented the problem of the
81 treatment planning as an inverse problem, as opposed to the forward planning method
82 that the experts of that time were used to. They also proposed that in order to provide
83 a uniform dose to a circular target structure while sparing the critical structure located
84 at the centre of the circular target, one needs to provide a radiation beam which is
85 non uniform. The advent of IMRT treatment technique increased the dose gradient
86 between tumour and close neighbour re-entrant normal structures. This allowed the
87 potential for dose escalation to PTVs as well as produced steep dose gradients between
88 PTV and some adjacent critical organs.
89 This chapter presents a brief introduction of some of the modern radiation therapy
90 delivery techniques used to deliver IMRT. It focuses on methods of delivery that are
91 subsequently used for dose delivery during dosimetric measurements. The two IMRT 2.3. Intensity Modulated Radiation Therapy 6
92 techniques include linac IMRT and Tomotherapy. In addition, the chapter discusses
93 stereotactic radiotherapy (SRS) as this is the method used for small field dosimetry.
94 Dosimeters that are currently used for the dosimetric verification for these deliveries
95 are also discussed as this provides a focus for comparison. This is followed by an
96 introduction to silicon based dosimeters for high spatial resolution dosimetry and two
97 dimensional array detectors.
98 2.3 Intensity Modulated Radiation Therapy
99 Intensity modulated radiation therapy is a complex radiation delivery. It utilizes an
100 MLC to modulate the radiation beam to achieve a high dose to the tumour while
101 sparing the critical organs close to the irradiated target. IMRT plans tend to produce
102 steep dose fall off regions at the interface of the target and organ-at-risk. In an
103 IMRT plan, a single field is made up of several segments comprised of different MLC
104 shapes. The sum of these segments will produce a non-uniformed radiation field. Two
105 techniques that are normally used on conventional linear accelerators are the step and
106 shoot delivery and sliding window delivery. Both techniques deliver IMRT treatment
107 with static gantry angles.
108 2.3.1 Step and shoot IMRT delivery
109 In the step and shoot technique, the radiation beam is only switched on when the
110 leaves are stationary. When the leaves are changing from one segment to the next, the
111 radiation beam is switched off (Boyer et al., 2001). 2.4. Volumetric Modulated Arc Therapy 7
112 2.3.2 Dynamic sliding window IMRT delivery
113 In the dynamic sliding window technique, the radiation beam is kept on all the time
114 as the MLC leaves sweep across the field. This technique is sometimes affected by the
115 mechanical limitation of the MLC whereby the leaves may not have reached the right
116 leaf positions, and therefore the radiation beam has to be halted temporarily to wait
117 for the leaf to come into the correct position (Boyer et al., 1992).
118 2.4 Volumetric Modulated Arc Therapy
119 The original concept is based on a paper by Yu (1995) called intensity modulated
120 arc therapy (IMAT). This concept did not include dose rate modulation. Volumetric
121 Modulated Arc Therapy (VMAT) is a special type of IMRT delivery technique where
122 treatment is delivered in a single dynamically modulated arc. In other words, the
123 radiation is delivered as the linac gantry rotates around the patient (Otto, 2008).
124 In VMAT, the dose rate and the gantry rotation speed is also varied to achieve the
125 required beam modulation. Because VMAT essentially uses all available angles in
126 the inverse plan domain of parameters, it has the potential to produce superior dose
127 distributions than IMRT while using the same inverse planning dose objective approach
128 as IMRT. Planning comparisons show some small gains over IMRT (Bzdusek et al.,
129 2009; Matuszak et al., 2010; Yan et al., 2010). The time advantages of not stopping and
130 starting the gantry at fixed gantry locations also suggest VMAT is faster to deliver
131 than IMRT, representing a gain of 30% in patient in room time when all in room
132 procedures are included (Hardcastle et al., n.d.). 2.5. Stereotactic Radiosurgery / Radiotherapy 8
133 2.5 Stereotactic Radiosurgery / Radiotherapy
134 Stereotaxy is a method by which a point is defined within the patient’s body by
135 an external three-dimensional coordinate system (Grosu et al., 2006). Stereotactic
136 radiosurgery (SRS) is the use of radiation ablation in place of conventional surgical
137 excision to remove or modify a benign lesion in the body (Metcalfe et al., 2007).
138 Traditionally, this method requires a delivery of a large dose in a single treatment,
139 resembling a surgical procedure. In malignancy cases, the treatment is carried out
140 using a series of equal dose fractions; this is then termed stereotactic radiotherapy
141 (SRT). The intention of the SRT/SRS treatment is to deliver a concentrated dose to
142 a small volume of tumour tissue, usually located in close proximity to critical organs.
143 Or, in other cases, it is intended as a boost dose to the target volume. The nature of
144 SRT/SRS treatment requires it to have a very high geometric precision, i.e. a tight
145 margin for the planning target volume (PTV) and a sharp dose fall off (Grosu et al.,
146 2006). Due to the high dose delivery and tight margins required in this technique, the
147 planning and delivery of the treatment requires great precision and accuracy.
148 2.6 Helical Tomotherapy
149 Tomotherapy, literary translated as “slice-therapy” was first introduced using the
150 PEACOCK (NOMOS) binary MLC (MIMiC) and retrofitted on conventional linacs
151 (Carol, 1995). On this device, the slices were treated with a single gantry rotation.
152 Then the patient was stepped to the next superior-inferior position.
TM 153 Helical Tomotherapy was born in the late 1980s at the University of Wisconsin,
154 Madison. It was designed with the aim to produce non-uniform intensity fields to
155 achieve a highly conformal dose distribution with a central avoidance structure (Mackie
156 et al., 1993; Mackie, 2006). The concept of helical Tomotherapy is to deliver intensity 2.7. Current quality assurance and dosimetric approaches 9
157 modulated radiation therapy in a continuous rotational manner with the linac mounted
158 on a slip ring gantry, while the patient is translated through the gantry bore. The
159 translation is analogous to a diagnostic computed tomography (CT) scanner, hence the
160 term Tomotherapy. The radiation source used is a fan beam 6 MV linear accelerator.
161 The beam modulation was achieved by a binary multileaf collimator (MLC) seated
162 at the distal end of the linear accelerator. The binary MLC consists of 64 tungsten
163 leaves extending across a field of 40 cm in the x-direction at the isocenter. In contrast
164 to the field shaping MLC in conventional linacs, the binary MLC has only two states,
165 either closed or opened. Hence, beam modulation is achieved by varying the leaf open
166 time, i.e. temporal modulation. The helical Tomotherapy device also has a linear
167 array of 738 xenon filled ion chambers that sits opposite the linac source on the gantry
168 ring. This detector is not only used for patient localisation but also allows the helical
169 Tomotherapy to perform MVCT image guided radiation therapy. The original design
170 (Mackie et al., 1993) had a kV x-ray source offset from the MV linac to provide kV
171 CT scans. This as yet has not been implemented and the MV CT scans are used for
172 IGRT.
173 2.7 Current quality assurance and dosimetric ap-
174 proaches
175 Currently, dosimetric verification of IMRT plans are patient based. Techniques used
176 are commonly based on a point dose measurement using ion chambers and planar
177 dose verification. Films, electronic portal imaging devices (EPIDs), and two dimen-
178 sional diode or ion chamber arrays are often used for planar dose verifications. Other
179 dosimeters such as diamond detectors, thermoluminescence detectors (TLDs) and gel
180 dosimeter are often used as complimentary dosimeters in IMRT delivery verification. 2.7. Current quality assurance and dosimetric approaches 10
181 2.7.1 Ionisation chamber
182 The ionisation chamber, the long established gold standard dosimeter used in radiation
183 therapy, is a highly reliable dosimeter. It measures the number of ion pairs produced in
184 a volume of air due to radiation. In the simplest arrangement, the ionisation chamber
185 exists as two electrode plates spaced apart in air. A large potential (100V - 400V)
186 is applied to the plates. When radiation traverses between the plates, it ionises the
187 air producing free negative electrons and positive ions. The positive and negative
188 charged particles are then swept towards the appropriate electrodes by the electric
189 field, producing a steady current flow in the external circuit, which can be measured
190 by an electrometer. The conversion from ionisation to dose is historically understood
191 (Bragg, 1910; Gray, 1929, 1936) and various protocols to correct for detector designs
192 are also well established (Almond et al., 1999; IAEA, 2006)
193 Due to the finite size of the air volume required for a sufficient signal, the ionisation
194 chambers are limited in physical dimension. When used in IMRT fields where high
195 dose gradients exist, the ionisation chamber tends to overestimate the penumbra width
196 due to its larger volume. Therefore it may not accurately represent areas with high
197 dose gradient fall off (Bucciolini et al., 2003). Exact output measurements in small
198 SRS fields are also a challenge (Das et al., 2000, 2008; Fan et al., 2009) due to the
199 volume averaging of the dose.
200 2.7.2 Film dosimetry
201 Radiographic film dosimetry has excellent spatial resolution, but it is non-tissue equiva-
202 lent, energy dependent and somewhat effected by processing conditions. Radiochromic
203 (Gafchromic) film is almost tissue-equivalent and self-developing, and its ease of han-
204 dling in room light makes this a convenient dosimeter for small field dosimetry. All
205 film methods provide a 2D pixel intensity map that can be converted to a 2D dose 2.7. Current quality assurance and dosimetric approaches 11
206 map given appropriate calibration. Film methods are non real time and are poten-
207 tially affected by polarization effects and non-uniformities in commercial film scanners
208 Butson et al. (2003b, 2006b); Paelinck et al. (2007).
209 2.7.3 Silicon diode
210 In medical radiation dosimetry, silicon is an attractive semiconductor material as a
211 radiation detector. Silicon diodes are able to achieve high signal response with a small
212 sensitive volume. Hence these detectors can measure to a high spatial resolution. This
213 is because, due to the higher density of a solid, it has 18000 times higher sensitivity
214 than an ionisation chamber with the same volume while the energy required to produce
215 an electron-hole pair is 3.6 eV, 10 times less than those required to ionise air. Another
216 feature that makes the silicon diode a good detector for radiation therapy application
217 is the constancy of the silicon to water electron stopping power ratio for the energy
218 range used in mega voltage radiation therapy (Figure 2.1). Other advantages of the
219 silicon diode are its robustness and mechanical stability, its excellent reproducibility,
220 and its ability to be used in passive mode (Saini & Zhu, 2002; Rosenfeld, 2006).
221 A silicon diode is primarily a p-n junction. A pure bulk silicon is first lightly doped
222 with phosphorus or boron to produce n-type or p-type silicon, respectively. Then, a
223 heavily doped impurity of the opposite type is implanted on the surface region to form
224 a p-n junction. The diode type is defined by the bulk material.
225 The principle of operation for a silicon diode as a radiation detector is illustrated
226 in Figure 2.2. At the p-n junction, the majority carriers of each side diffuses to the
227 opposite side, resulting in an electric field (or build-in potential, Ψo) which prevents
228 further diffusion of the majority carriers (Shi et al., 2003). This diffusion of the ma-
229 jority carriers also leave behind immobile negatively charged acceptors and positively
230 charge donor ions in the p-andn- side. This layer of the spatially charged region is 2.7. Current quality assurance and dosimetric approaches 12
231 called the ‘depletion layer’, W. The diode is said to be operating in passive mode when
232 no external bias voltage is applied. In this mode, no current will flow in the external
233 circuit. When ionising radiation traverses the diode, it generates electron-hole pairs
234 in the diode. The excess minority carriers (i.e. electrons in the p-side and holes in
235 the n-side) located within the diffusion lengths, Ln and Lp, will diffuse towards the
236 p-n junction and be swept by the built-in potential towards the electrodes. This flow
237 of charge results in a current flow in the external circuit, which is measurable by the
238 electrometer.
239 Silicon diodes are however not without limitations. The energy, angular, temper-
240 ature, and dose rate dependency require rigorous characterisation (Rice et al., 1987;
241 Beddar et al., 1994). Silicon diodes also over-respond at low photon energies (<150
242 keV) due to the increased photoelectric cross section (Rosenfeld, 2006; Wong et al.,
243 2010). Asymmetric packaging and the inherent anisotropy of the silicon diode will re-
244 sult in the detector’s directional dependence (Westermark et al., 2000; Higgins et al.,
245 2003; Marre & Marinello, 2004; Jursinic, 2009). The use of diodes for in-vivo dosimetry
246 needs to take into account possible temperature dependence (Van Dam et al., 1990;
247 Saini & Zhu, 2002). Dose rate dependence of the detector also needs to be understood,
248 particularly when it is used under pulsed dose rate such as in the linear accelerator
249 Rikner & Grusell (1983); Van Dam et al. (1990); Grusell & Rikner (1993); Shi et al.
250 (2003). Planar arrays of Si diodes are widely used in the quality assurance of IMRT.
4 251 Examples include the Sun Nuclear MapCHECK and the ScandiDos Delta bi-planar
252 system.
4 253 2.7.3.1 Delta system
4 254 The Delta system (Scandidos, Uppsala, Sweden) consists of two planar dosimeter
255 arrays mounted orthogonally. Dose reconstruction software can then generate 3D dose
256 maps. It has 1069 p-type silicon diodes on two orthogonal planes mounted within a 2.7. Current quality assurance and dosimetric approaches 13
Figure 2.1: Silicon/water stopping power ratio.
Figure 2.2: Schematic of a silicon p-n junction diode (Shi et al., 2003). 2.7. Current quality assurance and dosimetric approaches 14
257 polymethylmethacrylate (PMMA) phantom. The diodes are cylindrical in shape with
2 2 258 an active area of 0.78 mm . Diode pitched is 6 mm within the central 6 × 6cm
259 region, and 10 mm for the outer region.
260 Many 2D planar diode or ion chamber arrays were designed for normal incidence
261 beams while verification of composite dose distribution is achieved by summing mul-
262 tiple beams. This simplification of the QA method lacks the verification of gantry,
263 collimator and couch angles (Sadagopan et al., 2009).
4 264 The Delta system can be used in a mode that addresses an additional parameter
265 (i.e. gantry rotation) that exists in VMAT treatment delivery. It has an inclinometer
266 that is attached to the linac gantry providing independent measurements of the gantry
267 rotation for VMAT angle measurement.
268 2.7.4 Diamond detectors
269 Diamond detectors are lauded for their near tissue-equivalence and small sensitive
270 volume, which is particularly useful for small field dosimetry and measuring at high
271 dose gradient regions. The detectors were also found to be independent of photon beam
272 quality for the clinical range (Laub et al., 1997; De Angelis et al., 2002). However,
273 one needs to understand the dose-rate dependence and significant preirradiation dose
274 effect (Heydarian et al., 1996). Fidanzio et al. (2000) reported a sensitivity variation
275 of 1.8% for the dose per pulse range of 0.068 mGy to 0.472 mGy.
276 2.7.5 Gel dosimetry
277 Gel dosimeters are fabricated from radiosensitive chemicals dissolved in gelatine ma-
278 terial which upon irradiation, polymerised as a function of irradiated dose. Since the
279 first introduction of this dosimeter in the 1950s, the gel dosimeter undergoes several
280 modifications in its formulation to address various limitations of the dosimeter. The 2.7. Current quality assurance and dosimetric approaches 15
281 original ferrous sulphate chemical dosimeter was bog down with post irradiation ion
282 diffusion problem (Gore & Kang, 1984; Olsson et al., 1992) while the polymer gel
283 dosimeter is susceptible to atmospheric oxygen inhibiting the polymerisation process
284 (Maryanski et al., 1994). The recent nomoxic gel dosimeter (Fong et al., 2001) are
285 currently used by researchers in the field for clinical applications.
286 Various methods can be used to readout gel dosimeters such as the Magnetic Res-
287 onance Imaging (MRI), x-ray computed tomography (CT), optical CT or ultrasound
288 (Gore & Kang, 1984; Maryanski et al., 1994; Hilts et al., 2000; Mather et al., 2002).
289 Gel dosimetry has the advantage of having the same medium for dosimetry and dose
290 scattering. It is near water equivalent and have high spatial resolution in all three
291 dimensions. It is however affected by MRI- and gel-related non-uniformity and arte-
292 facts.(Vergote et al., 2004).
293 2.7.6 Electronic portal imaging device
294 The electronic portal imaging device (EPID) consists of a flat panel amorphous silicon
295 detector array mounted on a retractable arm opposite to the linac beam portal. It
296 measures the exit fluence of the radiation and is usually used for online portal imaging
297 for verification of patient localisation. In recent years, its use has been extended to the
298 fluence verification of IMRT deliveries. The EPIDs have excellent spatial resolution.
299 However, the use of the EPIDs as an IMRT dosimetric verification tool requires com-
300 plicated conversion of image to dose (Warkentin et al., 2003). Van Esch et al. (2004)
301 and Greer et al. (2007) have also achieve success in using this device as an IMRT
302 QA dosimeter device. Challenges in using the EPID as an online patient dosimeter
303 include its off axis energy response to spectral changes in the beam. Despite these
304 challenges, Partridge et al. (2002) have reported the use of EPID’s for online patient
305 dose verification. 2.8. High spatial resolution dosimetry 16
306 2.8 High spatial resolution dosimetry
307 One common feature shared by the various state of the art radiation therapy treatment
308 techniques described in the first part of this chapter is that the treatment plans gen-
309 erated with these new methods aim to achieve high target/tumour conformance while
310 avoiding the adjacent organ at risk. This generally produces a dramatic change in
311 dose intensity within a short spatial distance. Therefore, in order to measure the high
312 gradient dose regions generated by IMRT plans, the small field dosimetry SRS/SRT
313 plans and the non-conventional dosimetric considerations such that exist in helical To-
314 motherapy treatment plans, one needs to have a detector with high spatial resolution.
315 Hence, the sensitive volume of the detector has to be small in physical size (Kron
316 et al., 1993, 2002).
317 2.8.1 Concept of silicon strip detector
318 Research in high energy physics (HEP) pioneered the development of tracking vertex
319 detectors and data acquisition systems (DAQ). Experience gained from this research
320 was translated into the medical field. Strip detectors, which are position sensitive
321 detectors, were first used in high energy physics (HEP) experimental research for high
322 precision charge particle tracking in a vertex detector (Damerell, 1995; Rosenfeld et al.,
323 1993).
324 2.8.1.1 Theory of strip detector
325 In its generic form a strip detector is an array of p-n junction diodes. As an example,
326 for a p-type device, a low resistivity p-type wafer is used as the starting bulk material.
327 The p-n junctions are formed by ion implantation of the heavily doped impurities
+ 328 of the opposite side n strips onto the non-segmented wafer substrate. On the back
+ 329 end, a highly doped p implant is used to provide an ohmic contact connected to a 2.8. High spatial resolution dosimetry 17
330 negative voltage, VB (Figure 2.3). When a photon with sufficient energy traverses
331 through the silicon, photon interactions and electrons set in motion by the photon
332 creates electron-hole pairs along the tract. The free moving electrons and holes will
+ + 333 then move towards the p and n layers, respectively. A net current flow results from
334 this charge movement.
335 The current flow in the silicon is picked up by the metal contacts that are laid on
+ 336 top of the n strips. The metal contacts are isolated from the implanted strips by a
337 thin passivation layer silicon oxide (SiO2), so called field oxide. The strips are wire
338 bonded to a fan-out board connecting to an individual preamplifier for each detector.
339 The spatial resolutions of these devices are generally in the order of μm.
+ 340 The thin passivation layer of silicon oxide (SiO2) also exists between the n strips.
341 This passivation layer will inevitably collect positive charges which also increase with
342 accumulated radiation. The presence of this positive charge layer is compensated by a
343 thin layer of electrons in the bulk material. This may eventually create a low resistance
344 interstrip leakage path, effectively shorting the adjacent strips (Damerell 1995). One
+ 345 method to overcome this problem is to create floating p-type zones between the p
346 strips, also called p-stop layers. However, the addition of these p-stops resulted in
347 additional technical difficulties and costs. Another alternative to achieve interstrip
348 insulation is the use of a uniform blanket ion implant performed on the surface of the
349 silicon (p-spray) (Pellegrini et al., 2007).
350 Silicon strip detectors are usually fabricated with a thickness of approximately 300
351 μm. The thinness of a detector is limited by the loss of signal charge, made worse by
352 the poor signal-to-noise ratio due to the increased capacitance from strip to substrate
353 (Damerell, 1995). 2.8. High spatial resolution dosimetry 18
Figure 2.3: Operation principle of a silicon strip detector [modified from (Damerell, 1995)]. 2.8. High spatial resolution dosimetry 19
354 2.8.1.2 Readout electronics
355 The silicon strip detector typically comprises of 128 channels. The readout of all 128
356 channels simultaneously would be limited by the logistics of the readout electronics.
357 The evolution of the accompanying readout electronics for the silicon strip detectors
358 was driven by various factors:
359 (i) The rapid and diverse application of the detector by the HEP experimental re-
360 searches and the need to equip detectors to work in increasingly varied and
361 extreme hostile environments
362 (ii) The rapid development and commercialisation brought on by the integrated cir-
363 cuit industry
364 Current readout electronics mostly utilises Very Large Scale Integration Applica-
365 tion Specific Integration Circuit (VLSI ASIC) customised to perform the readout of
366 the silicon strip detector for the intended temporal and noise resolution.
367 2.8.2 Application of high spatial resolution dosimeters in med-
368 ical radiation therapy
369 Application of strip detectors in medical radiation dosimetry were reported by Pappas
2 370 et al. (2008) with their 128 channel strip detector (active area 0.25 × 0.25 mm and
371 pitch 0.3 mm) in stereotactic radiotherapy QA. The use of silicon strip detectors
372 with a pitch of 121 μm for dosimetric characterization and position imaging of Sr-
373 90 for cardiovascular brachytherapy was reported by Caccia et al. (2004). At the
374 Centre for Medical Radiation Physics (CMRP), a 0.2 mm pitch and a sensitive area
2 375 of 0.02 × 2mm high spatial resolution silicon strip detector, referred to as the Dose
376 Magnifying Glass (DMG) was designed, prototyped and used in IMRT, SRS, and
377 helical Tomotherapy QA verification (Wong et al., 2010, 2011a,b). 2.9. Concept of a two dimensional array detector 20
378 2.9 Concept of a two dimensional array detector
379 A strip detector provides excellent positional information in a single dimension, but
380 it lacks the capability to produce picture-like two-dimensional information. As IMRT
381 treatment becomes a common treatment delivery method in many radiation therapy
382 delivery centres, quality assurance (QA) of the safe and accurate delivery of IMRT
383 treatments has become ever more important. Routine patient specific QA is recom-
384 mended for IMRT treatments (Ezzell et al., 2003; Nelms et al., 2011).
385 In its very basic form, the QA of IMRT treatment involves
386 (i) point dose measurements, and
387 (ii) planar dose (or fluence) comparison.
388 For point dose measurement, the common practice is to measure a point dose at
389 a high dose-low gradient region. Planar dose comparison are made using films or two
390 dimensional detector arrays.
391 The use of electronic 2D detector arrays are gradually gaining widespread use.
392 They provide efficient means of measuring doses at multiple locations in the field,
393 providing real time feed back and performing planar dose comparison simultaneously
394 (Venkataraman et al., 2009). Examples of some of the commercially available 2D de-
395 tector arrays are the ionisation chamber based I’mRT MatriXX, COMPASS 2D trans-
396 mission detector (IBA, Germany), 2D-ARRAY Type 10024 (PTW Freiburg, Germany)
397 (Amerio et al., 2004; Spezi et al., 2005; Stasi et al., 2005; Poppe et al., 2006b) and
398 semiconductor based MapCHECK (Sun Nuclear, Melbourne, Fl) (Jursinic & Nelms,
399 2003; Letourneau et al., 2009). Efforts were also made to utilise the existing electronic
400 portal imaging device (EPID) on the linac as a 2D array detector (Greer et al., 2007).
401 The various existing 2D detector arrays varied in their respective detector/chamber
402 size and detector-to-detector spacing (Table 2.1). Poppe et al. (2007) studied the 2.9. Concept of a two dimensional array detector 21
403 volume averaging effect in large size detectors and detector array gaps. They proposed
404 the use of the Nyquist sampling theorem to estimate the optimal detector pitch for a
405 sample of IMRT plans. In their work, they found that for a moderately complicated
−1 406 head and neck IMRT plan, most of the spatial frequencies were <0.1 mm . Therefore
−1 407 the required sampling frequency should be 0.2 mm , corresponding to 5 mm detector
408 pitch.
409 In general, IMRT QA verification using 2D detector arrays for planar dose mea-
410 surement can be divided into three approaches:
411 (a) With the 2D detector array positioned on the linac couch or mounted onto the
412 gantry head with a special jig, in the absence of a patient
413 (b) With the 2D detector array positioned between the MLC collimator and the pa-
414 tient, acting as a transmission detector for on-line and in-vivo measurement, i.e.
415 during patient treatment
416 (c) With the 2D detector array positioned downstream of the patient, such as using
417 an EPID or mounting a detector array on the EPID.
418 The first approach, the use of 2D detector arrays on the linac couch for planar
419 dose measurement is one of the most common approaches. This type of QA is patient
420 specific and usually performed prior to the delivery of the first IMRT treatment. It
421 is used to check the ability of the MLC to achieve the intended dose and modulation
422 as planned by the treatment planning system (TPS). The measurement can only be
423 done on a phantom that was CT scanned and recalculated with the patient specific
424 IMRT plan. The underlying assumption is that IMRT delivery would be reproducible
425 throughout the whole patient treatment period.
426 The third approach, which utilised the EPID can be done with or without the
427 presence of the real patient. It requires accurate reconstruction and recalculation of 2.9. Concept of a two dimensional array detector 22 10 mm variable spacing (7.07 - 10 mm) jected to 10 mm at SDD 100 cm. 6.5 mm Detector pitch 7.5 mm 2 2 0.8 mm 5mm × × ter, 2 mm height tive size/volume mm diameter, 5 mm height) 40) 3.8 mm diame- 27) 5 32) cylindrical (4.5 × × × 371600 (40 4.4 mm pro- -type Si diode 445 0.8 Ion chamber 1027 (32 Ion chambern 729 (27 Multi wire ioni- sation chamber Planelel paral- chamber ionisation Table 2.1: Specification of commercial 2D arrays. Manufacturer Detector typetry,Germany No. of detectorsPTW Detector Freiburg, Germany sensi- Melbourne, Fl Germany IBA Dosimetry, Germany Commercial name I’mRT MatriXX IBA Dosime- 2D-Type 10024 ARRAY MapCHECKDAVID Sun system PTW Nuclear, Freiburg, COMPASS sys- tem 2.9. Concept of a two dimensional array detector 23
428 signals which include patient or phantom abrsoption and scattering. If it is being used
429 for in-vivo patient dose monitoring, daily patient variability could easily mask any
430 dose delivery error.
431 The ultimate aim of IMRT QA verification, although not fully realized, is to ensure
432 the safe, accurate and reproducible delivery of IMRT treatment onto the patient. This
433 not only means that IMRT QA has to be patient specific, correct MLC functioning
434 in the pretreatment QA verification, but also throughout the whole treatment period.
435 One way to achieve this is via the second approach, where the 2D detector array is
436 positioned upstream of the patient, usually on the linac accessory slot. The detector
437 has to be of the transmission type with minimal beam perturbation. With the on-
438 line and in-vivo measurements, daily and almost real time feedback on the accuracy
439 delivery of IMRT treatment is achievable.
440 The idea of a 2D array transmission detector was first conceptualised by Paliwal
441 et al. (1996). They used a commercial off the shelf dose area product (DAP) meter,
442 commonly used in a fluoroscopy unit in the radiology department to monitor the
443 radiation beam. The concept of a dose area product was further developed by Poppe
444 et al. (2006a) in the design of the DAVID system and by Islam et al. (2009) with the
445 Integral Quality Monitor (IQM). Both systems allow online comparison of the real
446 time, in-vivo measurement of dose area product due to the opened MLC leaves with
447 pre-recorded data. Both groups reported high sensitivity to MLC leaf positioning
448 error down to 1 mm. The DAVID system was eventually used in a clinic for daily
449 IMRT verifications (Poppe et al., 2010). Although these two systems were highly
450 sensitive to leaf positioning errors, the systems do not have 2D spatial information
451 of the delivered dose making point wise verification impossible. Consequently, the
452 detected MLC functioning errors could not be correlated with under- or overdosing of
453 anatomical features in the treatment planning system. In other words, the errors in 2.9. Concept of a two dimensional array detector 24
454 IMRT delivery cannot be correlated to clinical relevant impact or end points.
455 Up until this point, IMRT QA verifications involved primarily ensuring the correct
456 and reproducible MLC functioning, the stability of the radiation output and the cor-
457 rect setting of device related parameters. The comparison of the measured 2D dose
458 map were usually compared to the TPS generated dose map using performance metrics
459 such as the dose difference, distance to agreement (DTA) comparison or the Gamma
460 analysis (Low et al., 1998). Nelms et al. (2011) in their recently published paper have
461 quantified the adequacy of this performance metric (Gamma) and suggested it is not
462 always an ideal tool in the prediction of clinical relevant dose errors. Towards that
463 end, several commercial 2D detector arrays with advanced software reconstruction al-
464 gorithms attempt to close the gap between physics performance metrics with clinically
465 relevant dose errors. These systems measured the 2D dose maps and reconstructed the
466 3D dose distribution using the CT data input from the TPS. Examples of such systems
467 are the COMPASS system (IBA Dosimetry, Germany), DOSIMETRYCHECK (Math
468 Resolutions LLC) and 3DVH (Sun Nuclear Corporation).
469 2.9.1 Pixelated detector
470 On the extreme end of the spectrum of detector size, is the pixelated detector. “Pixel
471 detector” is taken to mean a device equipped with a two dimensional array of detectors.
472 This technology has taken flight in the commercial market particularly in camcorders
473 and digital cameras. The utilisation of a pixelated detector in the medical field is not
474 uncommon particularly in medical imaging, such as the positron emission tomography
475 (PET) and gamma cameras. 2.9. Concept of a two dimensional array detector 25
476 2.9.2 MAESTRO project
477 Pixelated silicon detectors were recently proposed and a prototype data acquisition
478 (DAQ) system was developed in the framework of the European project MAESTRO
+ 479 (MAESTRO, 2008). Each pixel is based on an n -p junction surrounded by a guard-
480 ring structure implanted on an epitaxial 50 μm thick p-type silicon layer grown on
481 a Czochralski (Sze, 2001) substrate, a monolithic silicon segmented module. The
482 best compromise between granularity and electronic complexity has been achieved by
2 483 choosing pixels with 2 × 2mm active area and 3 mm pitch. The discrete readout
484 electronics for 441 pixels/channels (21 × 21 pixels) was developed. The system demon-
485 strated good results but is quite bulky due to discrete readout electronics (Menichelli
486 et al., 2007; Talamonti et al., 2007). The next version will utilize nine TERA 06 ASIC
487 chips for readout of segmented detectors. The system in development offered better
488 resolution than 2D discrete diodes which have a spatial resolution of 3 mm. The use
489 of a TERA ASIC chip for the readout of a silicon pixelated detector or strip detector
490 was first proposed and utilised by the CMRP. 491 Chapter 3
492 Methodology
493 This chapter describes the two detector systems that were used in this thesis, the
494 Dose Magnifying Glass (DMG) and the Magic Plate (MP). The working principles
495 of the electronic readout system will also be described. In many parts of this thesis,
496 the film dosimetry was often used as the comparison dosimeter. The properties and
497 composition of the two types of films used in this thesis will also be briefly discussed.
498 3.1 Dose Magnifying Glass
499 3.1.1 Design and fabrication
+ 500 The Dose Magnifying Glass (DMG) is an array of 128 phosphor implanted n strips
+ 501 on a p-type silicon wafer. The sensitive area defined by a single n strip is 20 × 5000
2 2 nd 502 μm for the prototype version and 20 × 2000 μm for the 2 generation DMGs. The
503 thickness of the silicon wafer is 375 μm and the strip pitch is 200 μm.
504 Two types of silicon strip detectors based on p-Si were produced, one of high
0Part of this chapter has been published in Medical Physics: Wong, J. H. D., Carolan, M., Lerch, M. L. F., Petasecca, M., Khanna, S., Perevertaylo, V. L., Metcalfe, P. and Rosenfeld, A. B. (2010). A silicon strip detector dose magnifying glass for IMRT dosimetry. Medical Physics 37(2): 427-439.
26 3.1. Dose Magnifying Glass 27
Figure 3.1: Schematic diagram of the Dose Magnifying Glass (DMG).
505 resistivity (5 k Ω·cm) float zone silicon (Sze, 2001), and the other of low resistivity (10
506 Ω·cm) Chochralski silicon (Sze, 2001).
+ 507 The area between strips was implanted with Boron producing a p stop layer to
+ + 508 avoid shorting between adjacent n strips and between n strips with the common
509 electrode. The schematic diagram of the strip detector is presented in Figure 3.1.
+ 510 Aluminium was evaporated on top of the n areas. The detector was used in passive
511 mode and the readout was carried out with the detector configured as a planar detector
+ + 512 with a common electrode p from the same side as the n strips. 3.2. Magic Plate 28
Figure 3.2: Prototype version of DMG (left) and the 2nd generation DMG (right).
513 3.1.2 Detector packaging
514 The first version of DMG is mounted on a ceramic substrate. The ceramic substrate
3 515 was of a higher density (∼ 2.5 - 3.5 g/cm ) material than silicon. The presence of this
516 ceramic substrate in the mounting of the detector caused additional attenuation of
nd 517 the radiation and affected the angular dependence of the detector. The 2 generation
518 DMGs are mounted on a kapton substrate (Figure 3.2). The thin Kapton substrate
3 519 has a density of 1.42 g/cm , and has an effective atomic number, Zeff =6.6(Berger
520 et al., n.d.). This was more tissue-equivalent than ceramic (Al2O3,Zeff = 11.2) and
521 was expected to create less perturbation to the radiation beam. It was therefore
522 also expected to reduce the angular dependency of the detector due to the detector
523 packaging. The detector array was mounted at the end of a 170 mm long Kapton
524 pigtail. This enabled the positioning of the silicon detector in phantoms to provide for
525 sufficient scattering.
526 3.2 Magic Plate
527 The Magic Plate (MP) is a 2D array of 11 × 11 silicon diodes covering an area of 10 ×
2 528 10 cm (Figure 3.3). The diodes were mounted on a 0.64 mm thick Kapton substrate 3.2. Magic Plate 29
Figure 3.3: The Magic Plate (MP) 2D array detector.
529 using the ‘drop-in’ technology. This technology was proposed and developed at the
3 530 CMRP. The physical size of a single diode is 1.5 × 1.5 × 0.425 mm .TheMPwas
531 designed to be mounted on the linear accelerator head (accessory mount slot) in line
532 with the radiation beam (Figure 3.4). It was designed to operate as a transmission
533 detector measuring the 2D fluence map of a modulated radiation beam, hence the thin
534 Kapton substrate mounting. However, the MP can also be used as a planar 2D array
535 detector for phantom measurements.
536 3.2.1 Epitaxial diodes
537 The word epitaxy comes from the Greek word which means to arrange upon in an
538 orderly manner (Singh, 2001; Orton, 2004). This process involves growing a high
539 quality thin epitaxial film on top of a heavily doped bulk silicon wafer which acts as
540 a crystal seed for the epitaxial growth and later serves as a supporting structure.
541 An epitaxial diode is an attractive alternative to a conventional silicon detector 3.2. Magic Plate 30
Figure 3.4: The Magic Plate mounted on the linac gantry. The mounting plate was a modified TBI applicator. 3.3. TERA readout system 31
542 due to the thin epitaxial layer which is deemed to be more radiation hard. The ability
543 to reduce the detector sensitive area allows the detector to be made thinner compared
544 to the 300 μm thick silicon detector (Kramberger et al., 2003). The diodes that were
545 used in the MP were grown using the epitaxial-growth technique (Sze, 2001). For the
546 MP diodes, the p-epitaxial layer is a 50 μm thick p-Si layer grown on top of a 375
+ 547 μm thick high resistivity p substrate. The sensitive volume of the individual element
+ 3 548 defined by the n region is 0.5 × 0.5 × 0.05 mm , whilst the detector pitch is 1 cm.
549 The epitaxial diodes are of low resistivity type (100 Ω·cm). In the prototype design of
550 these epitaxial diodes, a 0.7 μm thick layer of SiO2 was grown on top of the epitaxial
551 layer. The silicon oxide (SiO2) layer is generally used as a mask to allow selective
+ 552 implantation of the n region as well as to serve as a protective layer for the silicon
+ 553 diode. The area surrounding the n region and the guard ring was implanted with
+ 11 3 554 Boron producing a p stop layer (0.1 μm and concentration = 10 Ω·cm )toprevent
+ 555 a positive charge building up at the p stop layer.
556 Two types of epitaxial diodes were manufactured, one with a guard ring (Figure
557 3.5) and one without a guard ring (Figure 3.6). In this thesis, the MP mounted with
558 epitaxial diodes that did not have a guard ring was used. The MP was operated in
559 passive mode.
560 3.3 TERA readout system
561 The TERA readout system was employed to readout the current signal from the DMG
562 and the MP. The TERA chip is a Very Large Scale Integration Application Specific
563 Integration Circuit (VLSI ASIC) that was designed by Istituto Nazionale di Fisica
564 Nucleare (INFN) - Torino Division and University of Torino microelectronics group
565 working on the readout of pixelated ionization or strip chambers for hadron therapy
566 (Bonazzola et al., 1998). It underwent several modifications during the last decade 3.3. TERA readout system 32
Figure 3.5: Schematic diagram of an epitaxial diode (with guard ring).
Figure 3.6: Schematic diagram of an epitaxial diode (without guard ring). 3.3. TERA readout system 33
Figure 3.7: Schematic diagram of the TERA readout [modified from (Mazza et al., 2005)].
567 (Mazza et al., 2005). Various versions of the TERA family of ASICSs have been
568 successfully implemented in different 2D and 3D “Magic Cube” ionization chambers
569 in hadron therapy (Brusasco et al., 1997; La Rosa et al., 2006, 2008). The latest
570 version (TERA 6.0) of the chip was used by the Scanditronix-Wellh¨ofer (IBA Group)
571 in commercial dosimeters used in Radiotherapy (I’mRT MatriXX and StarTrack) based
572 on gas ionization chambers.
573 3.3.1 How TERA works
574 The structure of the TERA ASIC is based on a current-to-frequency converter followed
575 by a digital counter. It operates based on the charge-balancing or recycling integrator
576 technique (Figure 3.7) (Gottschalk, 1983; Horowitz & Hill, 1989) whereby the circuit
577 counts the number of times a capacitor is charged by an input current and discharged
578 by the circuitry (Bonazzola et al., 1998).
579 When radiation impinges on the silicon detector, it forms electron-hole pairs in
580 the silicon. The free moving charge carriers generate a current which is detected by
581 the readout system. This is the input current, Iin. The input current charged up
582 the 600 fF integrating capacitor (Cint) in the operational transconductance amplifier
583 (OTA). The output waveform from the OTA appears as a voltage ramp, VA.Inthe 3.3. TERA readout system 34
584 comparator, the integrated charge is compared with a threshold voltage, Vth.When
585 the threshold voltage is reached (VA >Vth), the comparator fires a calibrated pulse
586 (VB) to the pulse generator (PG), which triggers the pulse generator to output two
587 pulses; one, to be sent to the digital counter and the second to the subtraction circuit.
588 The pulse issued by the PG to the subtraction circuit will charge up the 200fF
589 subtraction capacitor, Csub. The output responses of the Csub to Vsub is two δ-like
590 current pulses with equal charge but opposite in polarity (Figure 3.8). The amplitude
591 of the Vsub is defined by the difference of two externally set reference voltages, Vpulse+
592 and Vpulse− . The timing of the two current pulses coincides with the leading (δ+) and
593 trailing (δ-) edge of the voltage pulse. The current pulse with the same polarity to
594 the voltage output of the comparator is shorted to the OTA reference (Vref ), and the
595 other pulse is added algebraically to the input current to subtract a charge quantum
596 from the input current. This results in a sharp decrease of the charge across Cint,
597 proportional to the charge issued by the subtraction circuit (Bonazzola et al., 1998;
598 Mazza et al., 2005).
599 If at the end of the above process, the input voltage to the comparator is still above
600 the threshold voltage, the PG will continue to issue pulses to the counter and the
601 subtraction circuit. Once the integrated charge over Cint is below Vth, the comparator
602 trigger is reset and the PG will stop issuing pulses.
603 Using the charge-balancing technique, the integrated charge over Cint is “sub-
604 tracted” by the algebraic sum of the input current with a negative charge quantum.
605 This allows the Cint to continue to integrate while a fixed amount of charge is being
606 subtracted. Hence, no charge is lost. This is far superior to the conventional resetting
607 of the Cint, whereby dead time is introduced when the capacitor is reset.
608 The charge quantum, Qc is the unit charge needed for 1 count. It is defined by the 3.3. TERA readout system 35
Figure 3.8: Charge subtraction waveform (Mazza et al., 2005).
Table 3.1: Set reference values for the CMRP TERA 03 board. Value Vsupply 5V Vref 2V Vp+ 4.19 V Vp− 1.24 V Vth 2.50 V Cint 600 fF Csub 200 fF
609 following relationship (Eq.3.1),
Qc = Csub · Vsub (3.1)
610 where, Vsub =Vp+ -Vp−. The parameters Vp+ and Vp− are two reference voltages
611 set externally that determine the charge quantum. They can be varied to adjust the
612 Qc depending on the application. However, Qc should be kept within the range of 100
613 fC to 800 fC. The minimum limit of the charge quantum is set by the resolution of
614 the comparator, which is 100 fC. In this thesis, Qc = 590 fC (see Table 3.1 for detail
615 values).
616 The relationship between the input current and pulse frequency is described by 3.3. TERA readout system 36
617 Eq. 3.2:
I Pulse frequency, f = in (3.2) Qc
618 The maximum input current is limited by the pulse generator state machine (Mazza
619 et al., 2005). The current to frequency converter has a frequency limit of 5 MHz.
620 By which, the shortest time between two pulses issued sequentially is 0.2 μs. This
621 puts the limit of the maximum input current at 3 μA(forQc = 600 fC). However,
622 in the event of current overload, the subtraction circuit would not be able to keep
623 up with the integrated charge over Cint. The voltage input to the comparator, VA
624 would continue to remain at a level that is higher than Vth. This prompts the pulse
625 generator to continue to issue pulses at the maximum speed (both to the counter and
626 to the subtraction circuit) until the overload is removed and VA 627 In the case of pulsed linac radiation, a very high current (∼ 7 μA) is delivered 628 within a very short pulse width of 3.5 μs. The pulse period may range from 16 ms 629 to 2.7 ms (equivalent to 100 MU/min to 600 MU/min repetition rate). In this case, 630 within the initial pulse of 3.5 μs, the high input current charges the Cint way beyond 631 the Vth. A single Qc pulse will not be sufficient to bring the charge below the threshold 632 value. Hence, the pulse generator will continue to run at a maximum speed, issuing 633 pulses to the subtraction circuit until the overload is removed. Because the linac pulses 634 are brief, with a long pulse period, the subtraction circuit will have sufficient time to 635 remove the overload before the next pulse comes on. Hence, no error will result because 636 no charge will be lost (Gottschalk, 1983). However, one should be careful to ensure 637 that the input current producing voltage on the capacitor does not exceed the positive 638 rail voltage. In this case, it would be the power supply, Vsupply to the TERA board. 639 If the input current is driven up the positive rail, the output counts will approach 640 saturation where any increase in the input current cannot produce more counts. 3.3. TERA readout system 37 641 Intuitively, to resolve this problem, one could adjust the Qc to a larger value, 642 therefore subtracting the input current at a faster rate. This will be achieved at the 643 expense of a lower sensitivity of the counts produced. 644 Each TERA ASIC has 64 independent channels. In this thesis, two TERA chips 645 producing 128 channels were used. The channels are coupled to an individual digital 646 counter followed by a 16-bit register with a common load command. This allows the 647 counters to store the counts and readout at a specific time intervals. 648 3.3.2 Amplitude and timing 649 The accuracy of the comparator circuit defines the minimum equivalent Qc (100 fC 650 in this case) and hence the lower limit of the dynamic range. The upper limit of 651 the dynamic range is limited by the maximum voltage of the comparator circuit (the 652 supply voltage, 5 V), equating to a maximum charge Qc of 1 pC). The maximum 5 653 deviation from linearity for Qc = 600 fC is 0.5% for a dynamic range of 10 . 654 The data acquisition (DAQ) software was written using LabVIEW 8.6 (National 655 Instruments, USA) and allows on line and off line data analysis. Two twisted pair 656 ribbon cables, each 15 m long, connect the DAQ computer with the silicon strip 657 detector board. The FE4C front end board [produced by Physalus (Physalus S.r.l. 658 2001)] with two TERA 3.0 chips (PGA 144 package) mounted was controlled by a 659 Nuclear Instrument Card for logging data to the DAQ PC. The DAQ allows real time 660 measurements with a user defined readout frequency of between 250 Hz and 0.5 MHz 661 (effective input current integration time from 4 ms to 2 μs) 662 As mentioned above, a digital counter follows each channel output. The final out- 663 put of the chip is a 16-bit number that represents the number of times the integrating 664 capacitor has been discharged. In this case, one would expect to measure between 0 665 and 65536 Qc charge units depending on the level of dark current in the detector (zero 3.3. TERA readout system 38 666 in the case of a passive detector) that would give rise to a minimum >0 counts. After 667 each counter readout clock pulse, the counter is re-zeroed. If more than 65536 pulses 668 are received by the counter with one readout clock cycle, the counter rolls over and 669 starts counting from zero again. The user therefore has to have a systematic method 670 to determine a suitable readout clock cycle (referred to below as the pulse width) to 671 gather good counting statistics while ensuring that the counters do not roll over. The 672 check to determine that the counter has not rolled over can be easily performed by 673 acquiring data at different pulse widths and comparing the number of counts acquired. 674 For example, if the number of counts acquired with a 1-second pulse width doubled 675 the counts acquired with a 0.5-second pulse width, the counter has not rolled over at 676 1-second pulse width. In contrary, if the counter had rolled over, the counts acquired 677 with a 1-second pulse width would appear to be lower than those acquired with a 678 0.5-second pulse width. The user therefore should use a shorter pulse width for the 679 measurement. 680 3.3.3 Charge collection in silicon strip detector 681 Estimation of the charge collection in a single strip and the corresponding output 682 frequency of pulses from the TERA chip are important. Assume the linac provides 683 an average dose of 400 cGy/min corresponding to 6.67 cGy/s. This is made up of 684 200 separate 3 - 5 μs pulses every second (Rosenfeld, 2006). This corresponds to 685 0.33 mGy/pulse. For the estimation purpose, assume that the sensitive volume (100% 3 686 charge collection) of the single strip is at least 20 × 5000 × 300 μm . Taking the 3 3 −8 687 density of silicon to be 2.33 × 10 kg/m , the mass of silicon, m is 6.99 × 10 688 kg. Taking into account that the energy required to produce a electron-hole pair in 689 silicon, W = 3.6 eV and the dose in silicon and water are approximately similar for 690 linac photon energies,(Rosenfeld, 2006; Metcalfe et al., 2007) the expected charge is 3.4. Film dosimetry 39 691 (Eq.3.3), 692 Charge collected per linac beam pulse, D · m (3.3 × 10−4J) × (6.99 × 10−8kg) Q = = =6.32 × 10−12C (3.3) W/e 3.6J/C 693 where, D is the linac dose per pulse, m is the silicon mass, and 3.6eV/ (e-h pair) × 1.602 × 10−19J/eV W/e = =3.6J/C (3.4) 1.602 × 10−19C/e 694 For Qc = 600 fC/count on the TERA, the number of counts on the output of TERA 695 will be about 10 counts/linac pulse. That provides approximately 2000 counts/s for 696 the quoted dose rate. In the case where charge collection is <100% the number of 697 measured counts per second may be less. This justifies the use of the TERA chip for 698 the silicon strip detector in linac radiation fields. 699 3.4 Film dosimetry 700 3.4.1 Radiographic film 701 Radiographic film is made up of a clear polyester base coated with radiosensitive emul- 702 sion. This emulsion consists of silver bromide (AgBr) crystals embedded in gelatine 703 material. Absorption of photons or ionising radiation causes the silver bromide to 704 reduce to silver which appears as darkened areas on the film. The readout of the film 705 is done using an optical densitometer or commercial film digitizer such as the Vidar 706 scanner (Vidar Systems Corporation, Hendon, VA). The advantage of film is that it 707 has very high spatial resolution. Theoretically, it has resolution the size of a grain (0.2 708 -2μm) (Metcalfe et al., 2007). However, in reality, it is limited by the resolution of 709 the digitiser/readout mechanism. The disadvantages of radiographic film are that it is 3.4. Film dosimetry 40 710 a non real time dosimeter and it is energy dependent due to the high atomic number of 711 the silver. Other than radiation photons, films are also sensitive to light, which require 712 careful handling of the film. The reproducibility or reliability of the radiographic film 713 is also heavily dependent on the processor condition. 714 The radiographic film that was generally used in this thesis is the enhanced dynamic 715 range (EDR2) film (Eastman Kodak Company, Rochester, NY) which comes in a ready 716 packed form. It has a wide dynamic range of up to 600 cGy compared to its predecessor. 717 The EDR2 film is widely used in radiation therapy particularly in the QA verification 718 of IMRT deliveries (Zhu et al., 2002; Olch, 2002; Childress et al., 2005). 719 3.4.2 Radiochromic film 720 Radiochromic films are a relatively new dosimeter in medical radiation therapy. These 721 films are comprised of a polyester based coated with radiosensitive material which 722 polymerise upon irradiation, producing blue coloration of the film. The radiochromic 723 films used in this thesis were the Gafchromic EBT and EBT2 films (ISP, 2006, 2009). 724 The EBT2 film was rolled out in 2009 replacing its predecessor the EBT film. The 725 EBT2 film was distinguished from its predecessor (EBT) by its yellow color as well as 726 the film construction. The physical structure of the EBT and EBT2 film is shown in 727 Figure 3.9. The EBT2 film has a nominal thickness of 285 μm and is slightly thicker 728 than the EBT film. The advantages of using radiochromic films are, 729 (i) It is self developing and does not require the use of a conventional processor. 730 (ii) It has high spatial resolution which is only limited by the scanning resolution. 731 (iii) It is not sensitive to ambient light, hence it is easy to handle in room light. 732 (iv) Film composition has a lower Zeff = 6.84 (ISP, 2009) hence making it more 733 tissue equivalent compared to radiographic film. 3.4. Film dosimetry 41 734 (v) It is less energy dependent compared to the EBT film (Butson et al., 2010). 735 (vi) Readout can be performed using a professional film digitizer or an office flatbed 736 scanner. 737 (vii) It can be cut into any shapes or sizes. 738 (viii) It can be used for measurements in water. 739 The EBT2 film however also has its disadvantages: 740 (i) The film readout is non real time due to the effect of post-irradiation coloration. 741 Films are generally readout 24 hrs after exposure (Cheung et al., 2005). 742 (ii) It is susceptive to the scanner uniformity and artefacts (Butson et al., 2003b). 743 (iii) Due to the elongated grains, the scanning orientation is critical to ensure correct 744 dose readout (Butson et al., 2006a). 745 (iv) The film is susceptive to finger prints, dust, lint, scratches and moisture. This 746 will result in artefacts during readout. Therefore, careful handling using gloves 747 and lint free cloths are recommended. 748 (v) The EBT2 film suffers from film inhomogeneity up to 3.7% within a piece of 749 film, resulting in dose measurement uncertainty of 8.5% at 1 Gy (Hartmann 750 et al., 2010). 751 In short, radiochromic film is a useful comparative dosimeter and is used extensively 752 in this thesis. However, careful and proper handling of the film is important to ensure 753 correct and consistent readout. 3.4. Film dosimetry 42 Figure 3.9: Physical structure of Gafchromic (a) EBT and (b) EBT2 film (ISP, 2006, 2009). 754 Chapter 4 755 Radiation response and basic 756 characterisation of the Dose 757 Magnifying Glass 758 . 759 4.1 Introduction 760 Silicon diodes are widely used for dosimetry in radiation therapy. As with ion chamber 761 dosimetry, silicon diode dosimeters produce current or charge, which is proportional 762 to the dose rate or accumulated dose respectively. One advantage of silicon detectors 763 for radiotherapy is the constancy of the silicon-water electron stopping power ratio 764 over a wide energy range. However, silicon detectors have some disadvantages due 0Part of this chapter has been published in Medical Physics and Radiation Measurements: Wong, J. H. D., Carolan, M., Lerch, M. L. F., Petasecca, M., Khanna, S., Perevertaylo, V. L., Metcalfe, P. and Rosenfeld, A. B. (2010). A silicon strip detector dose magnifying glass for IMRT dosimetry. Medical Physics 37(2): 427-439. Wong, J. H. D., Cutajar, D., Lerch, M. L. F., Petasecca, M., Knittel, T., Carolan, M., Pere- vertaylo, V. L., Metcalfe, P. and Rosenfeld, A. B. (2011). From HEP to Medical Radia- tion Dosimetry - the silicon strip detector Dose Magnifying Glass. Radiation Measurements. DOI:10.1016/j.radmeas.2011.06.031.(Accepted, 14 June 2011) 43 4.2. Materials and methods 44 765 to their dose rate dependence, angular dependence, energy dependence, and radiation 766 damage. The dependence of the silicon diode on these parameters needs to be properly 767 characterised and understood prior to the clinical use of a silicon based detector. 768 In this chapter, the detector response of the Dose Magnifying Glass (DMG) is 769 characterised using similar tests to those carried out when testing the suitability of 770 detectors as dosimeters in the radiation therapy environment. The parameters tested 771 were: 772 (a) percent depth dose 773 (b) dose per pulse dependency 774 (c) stem effect 775 (d) dose linearity 776 (e) energy response dependency, and 777 (f) angular response dependency 778 4.2 Materials and methods 779 4.2.1 Percent depth dose measurement 780 The percent depth dose measures the dose drop off as a function of distance. This mea- 781 surement represents one of the principle measurements in determining the suitability 782 of a device as a radiation detector. Various literatures have shown that the depth dose 783 measurement made using silicon diodes agrees well with the ion chamber for most MV 784 photon energy beams. However, due to the non tissue equivalence of silicon, a slight 785 over estimation (<2.5%) is often observed at an increased depth in water/solid water 786 (Wilkins et al., 1997; Westermark et al., 2000; Bucciolini et al., 2003). 4.2. Materials and methods 45 Figure 4.1: Sketch of the solid water encapsulation for the prototype DMG. 787 The prototype DMG was sandwiched between two small slabs of solid water (2.5 3 788 × 5.0 × 0.5 cm ) machined to fit the sensitive area (Figure 4.1). The solid water was 789 then sandwiched inside another larger piece of polymethlmetacrylate (PMMA) (7.0 × 3 790 15.0 × 3.0 cm ), which was machined to accept the strip detector inside its solid water 791 encapsulation. The strip detector with the PMMA housing was placed on top of a 30 3 792 × 30 × 15 cm solid water block. Sections of solid water and PMMA were used to 793 provide scatter around the detector. The percent depth dose profiles were obtained 794 using solid water at the depth of 1.0 cm to 1.5 cm (at 0.1 cm intervals), 5 cm, 10 cm, 2 795 15 cm, and 20 cm for a 10 × 10 cm field size and at source to surface distance of 100 796 cm. Photon energy of 6 MV was used. The result was compared with a Farmer ion 797 chamber (NE-2571) in a solid water phantom. The active volume of the Farmer ion 798 chamber has an inner diameter of 6.3 mm and an internal length of 24.0 mm. 799 4.2.2 Dose per pulse response measurement 800 Dose per pulse dependency refers to the change of a detector’s sensitivity due to the 801 change of dose rate in a pulsed radiation beam such as those delivered by a medical 802 linear accelerator (Hoban et al., 1994). The dose during each pulse is sufficiently 803 high to cause a variation in the silicon detector response. It is therefore important to 804 characterize new detectors for dose per pulse sensitivity. This effect was first reported 4.2. Materials and methods 46 805 by Rikner & Grusell (1983). They found that n-type silicon diodes are more sensitive to 806 dose per pulse variation compared to p-type diodes, showing an increase in sensitivity 807 with increased dose per pulse. They also found that preirradiation of the diode also 808 reduces the dose per pulse dependence. The dose within a pulse can be changed by 809 changing the source to detector distance (SDD) and introducing attenuators such as 810 wedges or multileaf collimator (MLC) in the way of the radiation beam. However, by 811 the introduction of attenuators in the beam path, the dose per pulse response of the 812 silicon detector would be influenced by the alteration of the radiation beam energy 813 spectrum. The study of the dose per pulse response of the DMG was divided into two 814 sections; 815 (i) the effect of preirradiation condition and device resistivity at high dose rates, 816 and 817 (ii) the dose per pulse response of a low resistivity, preirradiated device for a larger 818 dose per pulse range. 819 In the later section, the influence of the alteration of the radiation beam energy 820 spectrum was also studied. The DMGs investigated in this thesis were produced with 821 two types of resistivity, 10 Ω·cm and 5 k Ω·cm. Of these, some of the DMGs were not 822 preirradiated while the others were preirradiated with 1 MeV electrons up to 15 kGy. 823 4.2.2.1 Effect of device resistivity and preirradiation condition on the dose 824 per pulse response at high dose rates 825 Four DMGs were used to study the effect of device resistivity and preirradiation con- 826 dition on the dose per pulse response. 827 (i) high resistivity, preirradiated, 828 (ii) low resistivity, preirradiated, 4.2. Materials and methods 47 829 (iii) high resistivity, unirradiated, and 830 (iv) low resistivity, unirradiated 831 The DMG was positioned at 9 cm depth in an I’mRT phantom and irradiated with 2 832 a10× 10 cm beam. Dose per pulse variation was obtained by varying the SDD from 833 91 cm to 158 cm. At each position, the detector was irradiated with a 6 MV photon 834 beam for a fixed amount of monitor units (MUs). The dose per pulse range measured −5 −4 835 was 9.45 × 10 to 2.72 × 10 Gy/pulse. 836 The actual dose at each measurement position was obtained with a CC13 ion 837 chamber (Scanditronix Wellh¨ofer) measurement with a similar set up. The ion chamber 838 dose rate dependence was measured to be <0.4% over the range of dose rates used. The 839 dose rate response was taken as the ratio of the DMG readings over the ion chamber 840 readings. 841 A Varian Clinac 21EX (Varian, Palo Alto, USA) linear accelerator was used in 842 these measurements. At the repetition rate of 600 MU/min, the linear accelerator 843 delivered radiation in short beam pulses of 3 μs at a frequency of 360 Hz. The dose 844 rate was defined by the amount of dose delivered by a radiation beam pulse. This was 845 estimated by dividing the dose measured by the ion chamber with the estimated total 846 number of beam pulses. At the source to detector distance of 100 cm the estimated −4 847 dose per pulse was 2.28 × 10 Gy/pulse. 848 4.2.2.2 Dose per pulse response of low resistivity, preirradiated device for 849 a large dose per pulse range 850 The low resistivity, preirradiated DMG was used to study the dose per pulse response −7 −4 851 of the device for a larger dose per pulse range of 6.42 × 10 to 2.92 × 10 Gy/pulse. 852 The DMG was positioned at 1.5 cm depth in solid water phantom and irradiated with 2 853 10 × 10 cm beam. The measurements were made in four conditions; 4.2. Materials and methods 48 2 854 (i) in a 10 × 10 cm open field, 855 (ii) with a 35 mm thick low melting point alloy (LMA) beam attenuator attached to 856 the linac accessory slot, 857 (iii) with a LMA and 50 mm lead (Pb) block beam attenuator, and 858 (iv) Under a closed MLC. 859 For each of these setups, measurements were made at several SDDs ranging from 860 88 cm to 150 cm. The radiation beam was delivered using a Siemens Oncor linear 861 accelerator with a 300 MU/min repetition rate. The linac beam pulse period is 4.4 862 ms and the estimated dose per pulse measured at SDD 101.5 cm and 1.5 cm depth in −4 863 solid water phantom is 2.2 × 10 Gy/pulse. 864 4.2.3 Stem effect measurement 865 Stem effect is the radiation induced conductivity in the cable or stem of the radiation 866 detector when exposed to a radiation beam (Beentjes & Garrett, 1966; Ibbott et al., 867 1975) 868 The stem effect of this device was investigated by irradiating the DMG with a 2 869 rectangular field size of 20 × 5cm at a 1.5 cm depth in solid water phantom. The 870 rectangular field size was chosen so that the narrow width of the field size was sufficient 871 to cover the detector area while avoiding the cable connectors when the long axis 872 was orthogonal to the cables. Measurements were made with the long axis of the 873 field parallel and then perpendicular to the cables. The difference between the two 874 measurements was deemed to be due to the stem effect. 4.2. Materials and methods 49 875 4.2.4 Dose linearity measurement 876 Dose linearity measures the linearity of the detector signal as a function of dose. A 877 good detector would require the detector’s response to be linearly proportional to a 878 given dose. A linear dose response also reduces the likelihood of errors and enables 879 faster processing (Metcalfe et al., 2007). 880 The DMG was inserted into the I’mRT phantom (IBA Dosimetry) and positioned 881 at 10 cm depth. PMMA sections were machined specifically to pack the detector 882 snugly in the I’mRT phantom. The set up was irradiated with 6 MV photons from 883 5 to 400 monitor units (MUs) delivered with a 400 MU/min repetition rate. The 884 corresponding doses delivered to the detector ranged from 3.89 cGy to 311.05 cGy. 885 4.2.5 Energy response measurement 886 Due to the higher atomic number of silicon, silicon diodes tend to display an enhanced 887 response to photoelectric effect at low photon energy <150 keV. This is usually mani- 888 fested in the over-estimation of dose at greater depths in water or solid water phantom 889 and with larger field sizes (Westermark et al., 2000; Yin et al., 2004). The energy de- 890 pendence of a silicon diode may also be due to the material around the diode such 891 as the electrode attachment, protective housing and build up material which contain 892 high Z materials (Saini & Zhu, 2007; Jursinic, 2009). 893 To study the energy response of the DMG, the strip detector was irradiated under 894 an orthovoltage machine (Guymay DX 3300) which operated at the range of 50 kV 895 - 250 kV. The orthovoltage machine was calibrated following IPEMB 1996 protocol 896 (IPEMB et al., 1996; IPEM et al., 2005). The absolute calibration accuracy is 1.9% 897 whereas the reproducibility of the calibration is better than 0.5%. The detector was 898 irradiated without any build up. For 50 - 150 kV tube voltage, a 100 mm diameter 899 circular applicator was used and the detector was placed at a focus to surface distance 4.2. Materials and methods 50 900 of 312 mm. For the 200 and 250 kV tube voltage, a 100 mm diameter circular applicator 901 was used with the detector positioned at focus to surface distance of 512 mm. For 902 high energy photons, the detector was irradiated using 6 MV and 10 MV beams from 903 a Varian Clinac 21EX linear accelerator (Varian, Palo Alto, USA) (dmax of 1.5 cm for 904 6 MV and 2.0 cm for 10 MV). The linac was calibrated according to IAEA TRS-398 905 protocol (IAEA, 2006).The relative standard uncertainty of the absolute calibration 906 is approximately 1.7%. 1 Gy water equivalent dose was delivered at all energies. The 907 readings are normalized to 1 at the energy of 6 MV. 908 4.2.6 Angular response measurement 909 Most silicon diodes display angular dependence (Westermark et al., 2000; Higgins 910 et al., 2003) due to the asymmetry in their mounting using high Z elements (Cu), 911 anisothropy in the silicon substrate and the packaging of the detector. Methods have 912 been devised to improve the angular dependence of the silicon diodes (Jursinic, 2009). 913 The angular response of the DMGs was studied to determine the effect of detector nd 914 packaging for the prototype version and the 2 generation DMG. 915 4.2.6.1 Prototype DMG 916 The design of the prototype detector with its ceramic mounting is expected to cause 917 some angular response to the DMG. The angular response of the prototype DMG was 918 studied in the I’mRT phantom. The whole DMG was inserted into the middle of the 919 I’mRT phantom and positioned at the isocenter. Radiation beams with the field size 2 ◦ 920 of 10 × 10 cm were delivered with static gantry angles at the intervals of 15 from ◦ ◦ 921 0 to 90 . The readings were then compared with CC13 ion chamber measurements 922 at the isocenter and at the positions ± 1 cm lateral from the isocenter. The angular 923 response is defined as the ratio of the signal measured by the DMG to the ion chamber 4.2. Materials and methods 51 924 measurements at each point of measurement. For the DMG measurement, the average 925 reading was taken of the detector channels within the distance of ± 1.5 mm from the 926 measurement point. This width is comparable to half of the inner diameter of the 927 CC13 ion chamber. The angular response taken at the 1 cm spatial distance from 928 the center of the detector enabled the investigation of the response of the detector 929 channels as a function of the distance from the center channel. nd 930 4.2.6.2 2 generation DMG nd 931 The packaging of the 2 generation DMG was modified by mounting the silicon wafer 932 on a 0.12 mm thick Kapton substrate. The thin Kapton substrate has a density of 3 933 1.42 g/cm , and has the effective atomic number, Zeff = 6.6 (Hubbell & Seltzer, 934 2004). This is more tissue-equivalent (Zeff = 7.5 for soft tissue) than ceramic (Al2O3, 935 Zeff = 11.2) and is expected to create less perturbation to the radiation beam. It 936 is therefore expected to reduce the angular dependency of the detector due to the 937 detector packaging. The whole detector with the Kapton was encapsulated in a solid 938 water holder which was machined to fit snugly with the device, reducing the presence 939 of air gaps between the solid water holder and the detector (Figure 4.2). The angular nd 940 response of the 2 generation DMG was investigated in a custom made phantom (SRS 941 phantom) which has a cylindrical body and a hemispherical head. The custom made 942 phantom design will be described in detail in chapter 6. The detector array is mounted 943 on the end of a 170 mm long Kapton pigtail. This enables the positioning of the silicon 944 detector in phantoms to provide for sufficient scattering. The angular response of the 945 DMG mounted on Kapton substrate was investigated in two orthogonal directions, 946 i.e. the azimuth angles and polar angles. Refer to Figure 4.3 for pictorial definition of 947 polar and azimuth angles. 948 For the angular response measurements, the DMG was positioned at SDD 100 cm, 949 9 cm depth of the SRS phantom. For azimuth angle measurements, the couch angle 4.3. Results 52 ◦ 2 950 was set to 0 andasmallfieldsizeof2.5× 2.5 cm was used. The gantry was rotated ◦ ◦ ◦ 951 from 0 to 350 with stepping angles of 10 intervals and 100 MUs were delivered at 952 each angle. ◦ 953 For measurement of polar angles, the couch was rotated to 270 . Due to the 954 design of the detector mounting on one end of the Kapton pigtail and the logistics 955 of the readout connection/cables at the other end of the Kapton pigtail, polar angle ◦ ◦ 956 measurements were limited to gantry angles of 320 to 180 , clockwise. Measurements ◦ 2 957 were taken in 10 increments, with a 2.5 × 2.5 cm radiation field size. 958 The DMG detector response was taken as the average response of the silicon strips 959 located within ± 1 mm of the axis of rotation. The angular dependence of the DMG 960 was taken as a ratio of the DMG response at different angles to the DMG response 961 when the radiation beam was perpendicular to the detector plane, corresponding to ◦ 962 the gantry angle 0 (Eq.4.1). Rθ Angular dependence, Aθ = (4.1) R◦ ◦ 963 where, Rθ is the response at gantry angle θ and R◦ is the response at gantry angle ◦ 964 0 . 965 4.3 Results 966 4.3.1 Percent depth dose 967 Figure 4.4 shows the measured percent depth dose profile using the strip detector 968 compared with a Farmer ion chamber (NE-2571) measurement, using solid water for a 2 969 10 × 10 cm field size at a source to surface distance of 100 cm. The photon energy was 970 6 MV. The DMG was corrected for the dose per pulse effect using the method described 971 in section 4.3.2.3. The resulting depth dose curve agrees with the Farmer ion chamber 4.3. Results 53 Figure 4.2: The 2nd generation DMG in the solid water holder and the custom made phantom which has a cylindrical body and hemispherical head. Figure 4.3: Definition of polar and azimuthal angle. 4.3. Results 54 Figure 4.4: Depth dose curve for a 10 × 10 cm2 field size of a 6 MV photon beam comparing DMG with Farmer ion chamber (NE-2571) measurements in solid water. 972 measurements within 0.8% up to 20 cm depth in solid water. The over-response of 973 the DMG measurements at greater depths may be due to the energy response of the 974 silicon detector with lower energy electrons. The D20/10 for the strip detector and the 975 Farmer ion chamber was 0.573 and 0.569, respectively. 976 4.3.2 Dose per pulse response 977 4.3.2.1 Effect of device resistivity and preirradiation condition on the dose 978 per pulse response at high dose rates 979 Figure 4.5 shows the dose rate dependence (expressed as dose per pulse) of the detector −4 980 normalized to the dose per pulse of 2.28 × 10 Gy/pulse measured at 100 cm SDD. 981 By changing the SDD, the dose rate was varied by a factor of 2.9 (corresponding to −5 −4 982 dose per pulse of 9.45 × 10 to 2.72 × 10 Gy/pulse). The low resistivity devices 983 appear to have less dose rate dependence compared to the high resistivity devices 984 whilst preirradiating the detector also improved the dose rate dependence. This is in 985 agreement with the published literature (Rikner & Grusell, 1983; Jursinic, 2009). The 4.3. Results 55 Figure 4.5: Dose per pulse response of selected DMGs with different resistivity and preirradiation conditions for the dose per pulse range of 9.45 × 10−5 to 2.72 × 10−4 Gy/pulse. The error bars represent the 1 standard deviation of the mean of all the channels in the strip detector. 986 unirradiated, high resistivity device showed the largest dose rate dependence of 8.5% 987 while the preirradiated, high resistivity device showed a slight improvement of 2%. 988 For the low resistivity, unirradiated devices, the maximum dose rate change was found 989 to be 4.1% while its preirradiated counterpart showed a maximum variation of 1.1%. 990 4.3.2.2 Dose per pulse response of a low resistivity, preirradiated device 991 for a large dose per pulse range 992 Figure 4.6 shows the dose per pulse response of a low resistivity, preirradiated device for −7 −4 993 a dose per pulse range of 6.42 × 10 to 2.92 × 10 Gy/pulse, a dose per pulse change 994 of 454-fold. This large dose per pulse range was achieved by placing beam attenuators 995 such as LMA, Pb block and closing the MLC. For each setup, measurements were 996 made at multiple SDDs creating a dynamic range of dose per pulse that overlapped 997 at either ends. 4.3. Results 56 −4 998 For the dose per pulse >1.01 × 10 Gy/pulse, this device appears to be dose −5 −5 999 per pulse independent. For the dose per pulse range of 3.92 × 10 to 2.66 × 10 1000 Gy/pulse, the DMG’s sensitivity increases as the dose per pulse increases. At the lower 1001 end of this range (between the measurements setup using the LMA + Pb block and −5 1002 closed MLC), there is a discontinuity of the data points. At 6.0 × 10 Gy/pulse, the 1003 normalised ratio of DMG/ion chamber gives 0.758 with the LMA + Pb block whilst 1004 under closed MLC condition, the DMG/ion chamber ratio is 0.364. This may indicate 1005 that with the closed MLC setup, there is a significant change in the beam spectrum 1006 affecting the silicon response. By changing the SDDs, the dose per pulse was changed 1007 by inversed square distance with the source; the beam spectrum remains the same. 1008 However, when attenuators were introduced, not only was the dose per pulse reduced 1009 but the radiation beam becomes harder as most of the soft x-rays were absorbed by 1010 the MLC. Silicon diode are known to over response with low energy photons, hence 1011 the hardening of the beam resulted in a generally much lower signal to the silicon. 1012 The dose per pulse response under closed MLC was 58%± 6%. This only affects 1013 the measurements when the radiation is highly attenuated, particularly the out of 1014 field dose. However, it should be noted that the dose under a closed MLC should 1015 comprised of leaf transmission leakage (which is usually in the order of 2% of an open 1016 field (Metcalfe et al., 2007)) and scattered dose from the phantom. 1017 The error bars represent the 1 standard deviation of the mean of all the channels 1018 in the strip detector. The large error bars for the LMA and Pb blocks were partly 1019 attributed to the slight inhomogeneity of the LMA block across the silicon strips, 1020 while the large error bars in the measurements under closed MLCs may be due to 1021 the combined effect of the low signal to noise ratio and the MLC leakage through the 1022 tongue and groves. 4.3. Results 57 Figure 4.6: Dose per pulse response for a low resistivity, preirradiated device for the dose per pulse range of 6.42 × 10−7 to 2.92 × 10−4 Gy/pulse. The dose per pulse responses were normalised to the dose per pulse of 2.2 × 10−4 Gy/pulse. 1023 4.3.2.3 Dose per pulse correction 1024 From the dose per pulse response measured in the previous sections, a correction 1025 function can be derived to correct for the dose per pulse dependence. This is achieved 1026 by fitting a best fit curve to the data points. The dose per pulse effect (α) is a function 1027 of the measured signal (Ni) to the signal at a calibrated dose per pulse, e.g. 2.2 × −4 1028 10 Gy/pulse (N◦) ratio (Eq.4.2).