Search

Objectives. To derive a collection efficiency formula,fGauss, for cylindrical ionization chambers in pulsed radiation beams from a volume recombination model of Boag et al (1996 Phys. Med. Biol. 41 885-97) including free electrons. To validatefGaussand a parallel plate chamber formulafexpusing an ion transport code and calculate changes in collection efficiencies caused by electric field charge screening at 0.1-100 mGy doses-per-pulse. And to determine collection efficienciesCE∞predicted at infinite voltage in the absence of avalanche effects by fitting scaled formulae to efficiencies computed for 100-400 V chamber voltages and 10 and 100 mGy doses-per-pulse. Approach. Calculations were performed for an idealized parallel plate chamber with 2 mm electrode separationd, and for an idealized cylindrical chamber with 0.5 and 2.333 mm inner and electrode radiirinandrout. Main results. fGaussandfexppredict the same collection efficiencies for cylindrical and parallel plate chambers satisfyingd2=(rout2-rin2)ln(rout/rin)/2, an equivalence condition met by the chambers studied. Without charge screening, efficiencies computed using the code equalledfGaussandfexp. With screening, efficiencies changed by ⩽0.03%, ⩽1.1% and ⩽21.3% at 1, 10 and 100 mGy doses-per-pulse, and differed between the chambers by ⩽0.9% and ⩽19.6% at ⩽10 and 100 mGy dose-per-pulse. For fits offexpandfGauss,CE∞values were ⩽1.2% and ⩽17.6% from unity at 10 and 100 mGy per pulse respectively, closer than for other formulae tested. Significance. Allowing for screening,fGaussandfexpdescribed computed collection efficiencies to within 0.03%, 1.1% and 21.3% at doses-per-pulse ⩽1, 10 and 100 mGy. Equivalence of the two chambers broke down at 100 mGy per pulse. Departures ofCE∞values from unity suggest that collection efficiencies determined experimentally by fittingfGaussorfexpto readings made at multiple voltages will be accurate to within 1.2% and 17.6% at 10 and 100 mGy per pulse respectively., (Creative Commons Attribution license.)

This Roadmap paper covers the field of precision preclinical x-ray radiation studies in animal models. It is mostly focused on models for cancer and normal tissue response to radiation, but also discusses other disease models. The recent technological evolution in imaging, irradiation, dosimetry and monitoring that have empowered these kinds of studies is discussed, and many developments in the near future are outlined. Finally, clinical translation and reverse translation are discussed., (Creative Commons Attribution license.)

Purpose: We have experimentally and computationally characterized the PTW microSilicon 60023-type diode's performance in 6 and 15 MV photon fields ≥5 × 5 mm 2 projected to isocenter. We tested the detector on- and off-axis at 5 and 15 cm depths in water, and investigated whether its response could be improved by including within it a thin airgap., Methods: Experimentally, detector readings were taken in fields generated by a Varian TrueBeam linac and compared with doses-to-water measured using Gafchromic film and ionization chambers. An unmodified 60023-type diode was tested along with detectors modified to include 0.6, 0.8, and 1.0 mm thick airgaps. Computationally, doses absorbed by water and detectors' sensitive volumes were calculated using the EGSnrc/BEAMnrc Monte Carlo radiation transport code. Detector response was characterized using k Q c l i n , 4 cm f c l i n , 4 cm , a factor that corrects for differences in the ratio of dose-to-water to detector reading between small fields and the reference condition, in this study 5 cm deep on-axis in a 4 × 4 cm 2 field., Results: The greatest errors in measurements of small field doses made using uncorrected readings from the unmodified 60023-type detector were over-responses of 2.6% ± 0.5% and 5.3% ± 2.0% determined computationally and experimentally, relative to the reading-per-dose in the reference field. Corresponding largest errors for the earlier 60017-type detector were 11.9% ± 0.6% and 11.7% ± 1.4% over-responses. Adding even the thinnest, 0.6 mm, airgap to the 60023-type detector over-corrected it, leading to under-responses of up to 4.8% ± 0.6% and 5.0% ± 1.8% determined computationally and experimentally. Further, Monte Carlo calculations indicate that a detector with a 0.3 mm airgap would read correctly to within 1.3% on-axis. The ratio of doses at 15 and 5 cm depths in water in a 6 MV 4 × 4 cm 2 field was measured more accurately using the unmodified 60023-type detector than using the 60017-type detector, and was within 0.3% of the ratio measured using an ion chamber. The 60023-type diode's sensitivity also varied negligibly as dose-rate was reduced from 13 to 4 Gy min -1 by decreasing the linac pulse repetition frequency, whereas the sensitivity of the 60017-type detector fell by 1.5%., Conclusions: The 60023-type detector performed well in small fields across a wide range of beam energies, field sizes, depths, and off-axis positions. Its response can potentially be further improved by adding a thin, 0.3 mm, airgap., (© 2021 American Association of Physicists in Medicine.)

Small-field dosimetry used in advance treatment technologies poses challenges due to loss of lateral charged particle equilibrium (LCPE), occlusion of the primary photon source, and the limited choice of suitable radiation detectors. These challenges greatly influence dosimetric accuracy. Many high-profile radiation incidents have demonstrated a poor understanding of appropriate methodology for small-field dosimetry. These incidents are a cause for concern because the use of small fields in various specialized radiation treatment techniques continues to grow rapidly. Reference and relative dosimetry in small and composite fields are the subject of the International Atomic Energy Agency (IAEA) dosimetry code of practice that has been published as TRS-483 and an AAPM summary publication (IAEA TRS 483; Dosimetry of small static fields used in external beam radiotherapy: An IAEA/AAPM International Code of Practice for reference and relative dose determination, Technical Report Series No. 483; Palmans et al., Med Phys 45(11):e1123, 2018). The charge of AAPM task group 155 (TG-155) is to summarize current knowledge on small-field dosimetry and to provide recommendations of best practices for relative dose determination in small megavoltage photon beams. An overview of the issue of LCPE and the changes in photon beam perturbations with decreasing field size is provided. Recommendations are included on appropriate detector systems and measurement methodologies. Existing published data on dosimetric parameters in small photon fields (e.g., percentage depth dose, tissue phantom ratio/tissue maximum ratio, off-axis ratios, and field output factors) together with the necessary perturbation corrections for various detectors are reviewed. A discussion on errors and an uncertainty analysis in measurements is provided. The design of beam models in treatment planning systems to simulate small fields necessitates special attention on the influence of the primary beam source and collimating devices in the computation of energy fluence and dose. The general requirements for fluence and dose calculation engines suitable for modeling dose in small fields are reviewed. Implementations in commercial treatment planning systems vary widely, and the aims of this report are to provide insight for the medical physicist and guidance to developers of beams models for radiotherapy treatment planning systems., (© 2021 American Association of Physicists in Medicine.)

Purpose: In small megavoltage photon fields, the accuracies of an unmodified PTW 60017-type diode dosimeter and six diodes modified by adding airgaps of thickness 0.6-1.6 mm and diameter 3.6 mm have been comprehensively characterized experimentally and computationally. The optimally thick airgap for density compensation was determined, and detectors were micro-CT imaged to investigate differences between experimentally measured radiation responses and those predicted computationally., Methods: Detectors were tested on- and off-axis, at 5 and 15 cm depths in 6 and 15 MV fields ≥ 0.5 × 0.5 cm 2 . Computational studies were carried out using the EGSnrc/BEAMnrc Monte Carlo radiation transport code. Experimentally, radiation was delivered using a Varian TrueBeam linac and doses absorbed by water were measured using Gafchromic EBT3 film and ionization chambers, and compared with diode readings. Detector response was characterized via the [Formula: see text] formalism, choosing a 4 × 4 cm 2 reference field., Results: For the unmodified 60017 diode, the maximum error in small field doses obtained from diode readings uncorrected by [Formula: see text] factors was determined as 11.9% computationally at +0.25 mm off-axis and 5 cm depth in a 15 MV 0.5 × 0.5 cm 2 field, and 11.7% experimentally at -0.30 mm off-axis and 5 cm depth in the same field. A detector modified to include a 1.6 mm thick airgap performed best, with maximum computationally and experimentally determined errors of 2.2% and 4.1%. The 1.6 mm airgap deepened the modified dosimeter's effective point of measurement by 0.5 mm. For some detectors significant differences existed between responses in small fields determined computationally and experimentally, micro-CT imaging indicating that these differences were due to within-tolerance variations in the thickness of an epoxy resin layer., Conclusions: The dosimetric performance of a 60017 diode detector was comprehensively improved throughout 6 and 15 MV small photon fields via density compensation. For this approach to work well with good detector-to-detector reproducibility, tolerances on dense component dimensions should be reduced to limit associated variations of response in small fields, or these components should be modified to have more water-like densities.

Andreo and Benmakhlouf (2017 Phys. Med. Biol. 62 1518-32) have disputed a finding of Scott et al (2012 Phys. Med. Biol. 57 4461-76) that the variation with field-size of the response of small ion chambers and solid-state dosimeters in small megavoltage photon radiation fields can largely be attributed to density. Further evidence for this finding was provided by Fenwick et al (2018 Phys. Med. Biol. 63 125003), but Andreo and Benmakhlouf (2018 Phys. Med. Biol. 63 125003) have now challenged the methodology used in that study. Specifically, Andreo and Benmakhlouf suggest that mass stopping-powers of fictitious materials used in Monte Carlo radiation transport calculations should be adjusted with material density according to the polarization effect, as if the materials were real and created by compressing other real materials. In this reply, we observe that fictitious materials are not real, and therefore their densities, mass stopping-powers and microscopic radiation interaction cross-sections can be freely and independently chosen to provide the clearest answers to the questions being studied. And we note that the key role played by density in small field detector response was further confirmed by our group back in 2013, using fictitious materials in which mass stopping-powers were adjusted with density, as preferred by Andreo and Benmakhlouf, as well as being held fixed, with very similar results being obtained in both circumstances (Underwood et al 2013a Med. Phys. 40 082102).

In small photon fields ionisation chambers can exhibit large deviations from Bragg-Gray behaviour; the EGSnrc Monte Carlo (MC) code system has been employed to investigate this 'Bragg-Gray breakdown'. The total electron (+positron) fluence in small water and air cavities in a water phantom has been computed for a full linac beam model as well as for a point source spectrum for 6 MV and 15 MV qualities for field sizes from 0.25  ×  0.25 cm(2) to 10  ×  10 cm(2). A water-to-air perturbation factor has been derived as the ratio of total electron (+positron) fluence, integrated over all energies, in a tiny water volume to that in a 'PinPoint 3D-chamber-like' air cavity; for the 0.25  ×  0.25 cm(2) field size the perturbation factors are 1.323 and 2.139 for 6 MV and 15 MV full linac geometries respectively. For the 15 MV full linac geometry for field sizes of 1  ×  1 cm(2) and smaller not only the absolute magnitude but also the 'shape' of the total electron fluence spectrum in the air cavity is significantly different to that in the water 'cavity'. The physics of this 'Bragg-Gray breakdown' is fully explained, making reference to the Fano theorem. For the 15 MV full linac geometry in the 0.25  ×  0.25 cm(2) field the directly computed MC dose ratio, water-to-air, differs by 5% from the product of the Spencer-Attix stopping-power ratio (SPR) and the perturbation factor; this 'difference' is explained by the difference in the shapes of the fluence spectra and is also formulated theoretically. We show that the dimensions of an air-cavity with a perturbation factor within 5% of unity would have to be impractically small in these highly non-equilibrium photon fields. In contrast the dose to water in a 0.25  ×  0.25 cm(2) field derived by multiplying the dose in the single-crystal diamond dosimeter (SCDDo) by the Spencer-Attix ratio is within 2.9% of the dose computed directly in the water voxel for full linac geometry at both 6 and 15 MV, thereby demonstrating that this detector exhibits quasi Bragg-Gray behaviour over a wide range of field sizes and beam qualities.

Standard commercial diode detectors over-respond within small radiation fields, an effect largely attributable to the relatively high mass-density of silicon. However, Monte Carlo studies can be used to optimise dosimeter designs and have demonstrated that 'mass-density compensation'-for example, introducing a low-density air-gap upstream of a diode's high-density silicon volume-can substantially improve instrument response. In this work we used egs_chamber Monte Carlo simulations to predict the ideal air-gap thickness for a PTW 60017 unshielded diode detector. We then developed a prototype instrument incorporating that air-gap and, for a 6 MV linac, tested it experimentally against EBT3 film. We also tested a further three prototypes with different air-gap thicknesses. Our results demonstrate that for a 10 × 10 cm(2) reference field the DiodeAir, a PTW 60017 diode with a built-in air-gap of 1 mm, has on-axis correction factors near unity. Laterally the DiodeAir performs very well off-axis and reports FWHM and penumbra values consistent with those measured using EBT3. For PDD measurement, the performance of the DiodeAir matches that of the original PTW 60017. The experimental focus of this work was 6 MV but we also simulated the on-axis response of the DiodeAir within 15 MV beams and found that our modification proved robust to this substantial increase in beam energy. However, the original diode 60017 does exhibit low energy scatter dependencies and may over-respond to high linac dose-rates such that applying the mass-density compensation method to an alternative instrument (particularly a diamond detector) could ultimately take us even closer to the small-field ideal.

Dosimeters often consist of several components whose mass densities differ substantially from water. These components cause small-field correction factors to vary significantly as lateral electronic equilibrium breaks down. Even amongst instruments designed for small-field dosimetry, inter-detector variation in the correction factors associated with very small (∼0.5 cm) fields can amount to tens of per cent. For a given dosimeter, small-field correction factors vary not only with field size but also with detector azimuthal angle and position within the field. Furthermore the accurate determination of these factors typically requires time-intensive Monte Carlo simulations. Thus, if achievable, 'correction factor free' small-field dosimetry would be highly desirable. This study demonstrates that a new generation of mass-density compensated detectors could take us towards this goal. Using a 6 MV beam model, it shows that 'mass-density compensation' can be utilized to improve the performance of a range of different detectors under small-field conditions. Non-sensitive material of appropriate mass-density is incorporated into detector designs in order to make the instruments behave as if consisting only of water. The dosimeter perturbative effects are then reduced to those associated with volume averaging. An even better solution-which modifies detectors to obtain profiles that look like those measured by a point-like water structure-is also considered. Provided that adequate sensitivity can be achieved for a small measurement volume, this study shows that it may be possible to use mass-density compensation (and Monte Carlo-driven design) to produce a solid-state dosimeter/ionization chamber with a near-perfect non-equilibrium response.

Purpose: The Alfonso et al. [Med. Phys. 35, 5179-5186 (2008)] formalism for small field dosimetry proposes a set of correction factors (kQclin,Qmsrfclin,fmsr) which account for differences between the detector response in nonstandard (clinical) and machine-specific-reference fields. In this study, the Monte Carlo method was used to investigate the viability of such small field correction factors for four different detectors irradiated under a variety of conditions. Because kQclin,Qmsrfclin,fmsr values for single detector position measurements are influenced by several factors, a new theoretical formalism for integrated-detector-position [dose area product (DAP)] measurements is also presented and was tested using Monte Carlo simulations., Methods: A BEAMnrc linac model was built and validated for a Varian Clinac iX accelerator. Using the egs++ geometry package, detailed virtual models were built for four different detectors: a PTW 60012 unshielded diode, a PTW 60003 Diamond detector, a PTW 31006 PinPoint (ionization chamber), and a PTW 31018 MicroLion (liquid-filled ionization chamber). The egs_chamber code was used to investigate the variation of kQclin,Qmsrfclin,fmsr with detector type, detector construction, field size, off-axis position, and the azimuthal angle between the detector and beam axis. Simulations were also used to consider the DAP obtained by each detector: virtual detectors and water voxels were scanned through high resolution grids of positions extending far beyond the boundaries of the fields under consideration., Results: For each detector, the correction factor (kQclin,Qmsrfclin,fmsr) was shown to depend strongly on detector off-axis position and detector azimuthal angle in addition to field size. In line with previous studies, substantial interdetector variation was also observed. However, it was demonstrated that by considering DAPs rather than single-detector-position dose measurements the high level of interdetector variation could be eliminated. Under small field conditions, mass density was found to be the principal determinant of water equivalence. Additionally, the mass densities of components outside the sensitive volumes were found to influence the detector response., Conclusions: kQclin,Qmsrfclin,fmsr values for existing detector designs depend on a host of variables and their calculation typically relies on the use of time-intensive Monte Carlo methods. Future moves toward density-compensated detector designs or DAP based protocols may simplify the methodology of small field dosimetry.

The dose imparted by a small non-equilibrium photon radiation field to the sensitive volume of a detector located within a water phantom depends on the density of the sensitive volume. Here this effect is explained using cavity theory, and analysed using Monte Carlo data calculated for schematically modelled diamond and Pinpoint-type detectors. The combined impact of the density and atomic composition of the sensitive volume on its response is represented as a ratio, Fw,det, of doses absorbed by equal volumes of unit density water and detector material co-located within a unit density water phantom. The impact of density alone is characterized through a similar ratio, Pρ -, of doses absorbed by equal volumes of unit and modified density water. The cavity theory is developed by splitting the dose absorbed by the sensitive volume into two components, imparted by electrons liberated in photon interactions occurring inside and outside the volume. Using this theory a simple model is obtained that links Pρ - to the degree of electronic equilibrium, see, at the centre of a field via a parameter Icav determined by the density and geometry of the sensitive volume. Following the scheme of Bouchard et al (2009 Med. Phys. 36 4654-63) Fw,det can be written as the product of Pρ -, the water-to-detector stopping power ratio [L[overline](Δ)/ρ](w)(det), and an additional factor Pfl -. In small fields [L[overline](Δ)/ρ](w)(det) changes little with field-size; and for the schematic diamond and Pinpoint detectors Pfl - takes values close to one. Consequently most of the field-size variation in Fw,det originates from the Pρ - factor. Relative changes in see and in the phantom scatter factor sp are similar in small fields. For the diamond detector, the variation of Pρ - with see (and thus field-size) is described well by the simple cavity model using an Icav parameter in line with independent Monte Carlo estimates. The model also captures the overall field-size dependence of Pρ - for the schematic Pinpoint detector, again using an Icav value consistent with independent estimates.

The impact of density and atomic composition on the dosimetric response of various detectors in small photon radiation fields is characterized using a 'density-correction' factor, F(detector), defined as the ratio of Monte Carlo calculated doses delivered to water and detector voxels located on-axis, 5 cm deep in a water phantom with a SSD of 100 cm. The variation of F(detector) with field size has been computed for detector voxels of various materials and densities. For ion chambers and solid-state detectors, the well-known variation of F(detector) at small field sizes is shown to be due to differences between the densities of detector active volumes and water, rather than differences in atomic number. However, associated changes in the measured shapes of small-field profiles offset these variations in F(detector), so that integral doses measured using the different detectors are quite similar, at least for slit fields. Since changes in F(detector) with field size arise primarily from differences between the densities of the detector materials and water, ideal small-field relative dosimeters should have small active volumes and water-like density.

Purpose: The purpose of this work is twofold: To further develop an approach to multiobjective optimization of rotational therapy treatments recently introduced by the authors [J. Pardo-Montero and J. D. Fenwick, "An approach to multiobjective optimization of rotational therapy," Med. Phys. 36, 3292-3303 (2009)], especially regarding its application to realistic geometries, and to study the quality (Pareto optimality) of plans obtained using such an approach by comparing them with Pareto optimal plans obtained through inverse planning., Methods: In the previous work of the authors, a methodology is proposed for constructing a large number of plans, with different compromises between the objectives involved, from a small number of geometrically based arcs, each arc prioritizing different objectives. Here, this method has been further developed and studied. Two different techniques for constructing these arcs are investigated, one based on image-reconstruction algorithms and the other based on more common gradient-descent algorithms. The difficulty of dealing with organs abutting the target, briefly reported in previous work of the authors, has been investigated using partial OAR unblocking. Optimality of the solutions has been investigated by comparison with a Pareto front obtained from inverse planning. A relative Euclidean distance has been used to measure the distance of these plans to the Pareto front, and dose volume histogram comparisons have been used to gauge the clinical impact of these distances. A prostate geometry has been used for the study., Results: For geometries where a blocked OAR abuts the target, moderate OAR unblocking can substantially improve target dose distribution and minimize hot spots while not overly compromising dose sparing of the organ. Image-reconstruction type and gradient-descent blocked-arc computations generate similar results. The Pareto front for the prostate geometry, reconstructed using a large number of inverse plans, presents a hockey-stick shape comprising two regions: One where the dose to the target is close to prescription and trade-offs can be made between doses to the organs at risk and (small) changes in target dose, and one where very substantial rectal sparing is achieved at the cost of large target underdosage. Plans computed following the approach using a conformal arc and four blocked arcs generally lie close to the Pareto front, although distances of some plans from high gradient regions of the Pareto front can be greater. Only around 12% of plans lie a relative Euclidean distance of 0.15 or greater from the Pareto front. Using the alternative distance measure of Craft ["Calculating and controlling the error of discrete representations of Pareto surfaces in convex multi-criteria optimization," Phys. Medica (to be published)], around 2/5 of plans lie more than 0.05 from the front. Computation of blocked arcs is quite fast, the algorithms requiring 35%-80% of the running time per iteration needed for conventional inverse plan computation., Conclusions: The geometry-based arc approach to multicriteria optimization of rotational therapy allows solutions to be obtained that lie close to the Pareto front. Both the image-reconstruction type and gradient-descent algorithms produce similar modulated arcs, the latter one perhaps being preferred because it is more easily implementable in standard treatment planning systems. Moderate unblocking provides a good way of dealing with OARs which abut the PTV. Optimization of geometry-based arcs is faster than usual inverse optimization of treatment plans, making this approach more rapid than an inverse-based Pareto front reconstruction.

Accurate characterization of small-field dosimetry requires measurements to be made with precisely aligned specialized detectors and is thus time consuming and error prone. This work explores measurement differences between detectors by using a Monte Carlo model matched to large-field data to predict properties of smaller fields. Measurements made with a variety of detectors have been compared with calculated results to assess their validity and explore reasons for differences. Unshielded diodes are expected to produce some of the most useful data, as their small sensitive cross sections give good resolution whilst their energy dependence is shown to vary little with depth in a 15 MV linac beam. Their response is shown to be constant with field size over the range 1-10 cm, with a correction of 3% needed for a field size of 0.5 cm. BEAMnrc has been used to create a 15 MV beam model, matched to dosimetric data for square fields larger than 3 cm, and producing small-field profiles and percentage depth doses (PDDs) that agree well with unshielded diode data for field sizes down to 0.5 cm. For fields sizes of 1.5 cm and above, little detector-to-detector variation exists in measured output factors, however for a 0.5 cm field a relative spread of 18% is seen between output factors measured with different detectors-values measured with the diamond and pinpoint detectors lying below that of the unshielded diode, with the shielded diode value being higher. Relative to the corrected unshielded diode measurement, the Monte Carlo modeled output factor is 4.5% low, a discrepancy that is probably due to the focal spot fluence profile and source occlusion modeling. The large-field Monte Carlo model can, therefore, currently be used to predict small-field profiles and PDDs measured with an unshielded diode. However, determination of output factors for the smallest fields requires a more detailed model of focal spot fluence and source occlusion.

In a previous paper, we described quality assurance procedures for Hi-Art helical tomotherapy machines. Here, we develop further some ideas discussed briefly in that paper. Simple helically generated dose distributions are modelled, and relationships between these dose distributions and underlying characteristics of Hi-Art treatment systems are elucidated. In particular, we describe the dependence of dose levels along the central axis of a cylinder aligned coaxially with a Hi-Art machine on fan beam width, couch velocity and helical delivery lengths. The impact on these dose levels of angular variations in gantry speed or output per linear accelerator pulse is also explored.

A trial of nonescalated conformal versus conventional radiotherapy treatment of prostate cancer has been carried out at the Royal Marsden NHS Trust (RMH) and Institute of Cancer Research (ICR), demonstrating a significant reduction in the rate of rectal bleeding reported for patients treated using the conformal technique. The relationship between planned rectal dose-distributions and incidences of bleeding has been analyzed, showing that the rate of bleeding falls significantly as the extent of the rectal wall receiving a planned dose-level of more than 57 Gy is reduced. Dose-distributions delivered to the rectal wall over the course of radiotherapy treatment inevitably differ from planned distributions, due to sources of uncertainty such as patient setup error, rectal wall movement and variation in the absolute rectal wall surface area. In this paper estimates of the differences between planned and treated rectal dose-distribution parameters are obtained for the RMH/ICR nonescalated conformal technique, working from a distribution of setup errors observed during the RMH/ICR trial, movement data supplied by Lebesque and colleagues derived from repeat CT scans, and estimates of rectal circumference variations extracted from the literature. Setup errors and wall movement are found to cause only limited systematic differences between mean treated and planned rectal dose-distribution parameter values, but introduce considerable uncertainties into the treated values of some dose-distribution parameters: setup errors lead to 22% and 9% relative uncertainties in the highly dosed fraction of the rectal wall and the wall average dose, respectively, with wall movement leading to 21% and 9% relative uncertainties. Estimates obtained from the literature of the uncertainty in the absolute surface area of the distensible rectal wall are of the order of 13%-18%. In a subsequent paper the impact of these uncertainties on analyses of the relationship between incidences of bleeding and planned rectal dose-distributions is explored.

A trial of nonescalated conformal versus conventional radiotherapy treatment of prostate cancer has been carried out at the Royal Marsden NHS Trust (RMH) and Institute of Cancer Research (ICR), demonstrating a significant reduction in the rate of rectal bleeding reported for patients treated using the conformal technique. The rate of bleeding has been shown to fall significantly as the extent of rectal wall receiving a planned dose-level in excess of 57 Gy is reduced. Dose-distributions delivered to the rectal wall over the course of radiotherapy treatment inevitably differ from planned distributions. In a previous paper estimates were obtained of the uncertainties in some planned rectal dose-distribution parameters generated by patient setup error, rectal wall movement and the variable degree of rectal wall distension. Here these uncertainties are combined to obtain estimates of the total planning uncertainties in rectal dose-distribution parameters thought likely, on the basis of mechanistic biological modeling, to correlate strongly with the complication rate. Working from these totaled uncertainty values, together with values of patient-to-patient and technique-to-technique differences in planned dose-distribution parameters, it can be inferred that the rectal dose-distribution uncertainties: (i) Have only a marginal impact on fits of a normal tissue complication probability (ntcp) model to RMH/ICR dose-distribution and grade 1, 2, 3 bleeding data (slightly flattening observed volume-response curves); (ii) only slightly reduce the power of a 2 x 100 patient trial of conformal versus conventional prostate radiotherapy to detect a significantly lower rate of grade 1,2,3 rectal bleeding amongst conformally treated patients; (iii) do not diminish the information content of individual planned patient dose-distribution data to the point where the fitting of technique-averaged data would provide as sensitive a test of the existence of a volume effect as the fitting of individual patient data.

An analysis technique, based on simulated annealing, is described which is employed to derive megavoltage photon beam spectral information from narrow beam attenuation measurements. Megavoltage photon beam spectra have been determined using this technique for linear accelerators from different manufacturers, and different models from individual manufacturers at a range of energies from nominal 6 MV to nominal 25 MV. All of the photon beams included in the study are in routine clinical use. The subsequent effects on dosimetry of employing derived primary spectra to specify beam quality are examined. The results suggest that the quality index TPR(20)10 may be insensitive to beam quality changes for high-energy beams in the range of 15 MV to 25 MV. Although the quality index may be insensitive as a beam quality specifier at these higher qualities, the actual difference in the calculated dose delivered using derived spectra as the quality specifier rather than TPR(20)10 is likely to be small, the results obtained indicating a difference of between 0.2% and 0.7% in the calculated dose delivered.