Quantification of Artifacts and Image Distortions in 1.5 Tesla Magnetic Resonance Images of a Commercial Multi-Channel Vaginal Cylinder Brachytherapy Applicator Set

Background

The BEBIG Portio multi-channel applicator provides better target dose coverage and sparing organs-at-risk compared to a single-channel cylinder. However, artifacts and distortions of Portio in magnetic resonance images (MRI) have not yet been reported.

Objective

We aimed to quantify the artifacts and distortions in its 1.5-Tesla MR images before clinical use.

Material and Methods

In this experimental study, we employed a gelatin-filled phantom to conduct our measurements. T2-weighted (T2W) images were examined for artifacts and distortions. Computed tomography (CT) images were used as a reference to assess image distortions. Artifact severity was measured by recording the full-width-at-half-maximum (FWHM) image pixel values at various positions along the length of the applicator/channels. CT and MRI-based applicator reconstruction accuracy were then compared, and signal-to-noise ratio (SNR) and contrast were also determined for the applicator images.

Results

The applicator distortion level for the Portio applicator was less than the image spatial resolution (0.5±0.5 pixels). The average FWHM for the tandem applicator images was 5.23±0.39 mm, while it was 3.21±0.37 mm for all channels (compared to their actual diameters of 5.0 mm and 3.0 mm, respectively). The average applicator reconstruction difference between CT and MR images was 0.75±0.30 mm overall source dwell positions. The image SNR and contrast were both acceptable.

Conclusion

These findings indicate that the Portio applicator has a satisfactory low level of artifacts and image distortions in 1.5-Tesla, T2W images. It may, therefore, be a promising option for MRI-guided multi-channel vaginal brachytherapy.

Metal artifact reduction in cervix brachytherapy with titanium applicators using dual-energy CT through virtual monoenergetic images and an iterative algorithm: A phantom study

PURPOSE

To evaluate an iterative metal-artifact reduction (iMAR) algorithm, dual-energy CT (DECT) through virtual monoenergetic images (VMI), and a combination of iMAR and DECT for reducing metal artifact severity (AS) induced by Fletcher titanium applicators used in cervix brachytherapy, the efficacy of which are hitherto unreported.

METHODS AND MATERIALS

120 kV_p single-energy CT (SECT) (Siemens) of BEBIG tandem applicators, varying in shape (straight or curved) and diameter (3.5 mm or 5 mm) in a custom-made water-filled phantom, and their DECT images obtained from extrapolation of 80 kVp and 140 kVp, were reconstructed using four methods: DECT through VMI±iMAR, and SECT±iMAR. The DECT images were reconstructed monoenergetically at 70, 150, and 190 keV. AS was evaluated using measured values and statistical analysis.

RESULTS

iMAR, DECT, and combined DECT and iMAR reduced AS (p < 0.05). DECT had a lower AS than SECT, even without iMAR (p < 0.025). SECT+iMAR was more effective than DECT-iMAR with VMI at 70 and 190 keV (p < 0.05), whereas showing no statistically significant difference at 150 keV. With DECT and iMAR combined, AS was reduced more effectively compared to the SECT+iMAR or DECT alone. It also reduced the mean interobserver uncertainty by 0.2 mm.

CONCLUSIONS

These findings indicate that iMAR reduces the AS caused by Fletcher titanium applicators for both SECT and DECT, a combination of iMAR and DECT is superior to either strategy alone, and at low energies, DECT+iMAR also produces similar artifact reduction. These practical strategies promise more accurate source-position and structure definitions in CT-based gynecological brachytherapy treatment planning.

Synthesizing CT images from MR images with deep learning: model generalization for different datasets through transfer learning

PURPOSE

Synthetic CT generation is the focus of many studies, however, only models on data with the same dataset were tested. Therefore, how well the trained model will work for data from different hospitals and MR protocols is still unknown. In this study, we addressed the model generalization problem for the MR-to-CT conversion task.

METHODS

Brain T2 MR and corresponding CT images were collected from one hospital and brain T1-FLAIR, T1-POST MR, and corresponding CT images were collected from another hospital. To investigate the model's generalizability ability, four potential solutions were proposed: source model, target model, combined model, and adapted model. All models were trained using the CycleGAN network. The source model was trained with a source domain dataset from scratch and tested with a target domain dataset. The target model was trained with a target domain dataset and tested with a target domain dataset. The combined model was trained with both source domain and target domain datasets, and tested with the target domain dataset. The adapted model used a transfer learning strategy to train a CycleGAN model with a source domain dataset and retrain the pre-trained model with a target domain dataset. MAE, RMSE, PSNR, and SSIM were used to quantitatively evaluate model performance on a target domain dataset.

RESULTS

The adapted model achieved best quantitative results of 74.56±8.61, 193.18±17.98, 28.30±0.83, and 0.84±0.01 for MAE, RMSE, PSNR, and SSIM using the T1-FLAIR dataset and 74.89±15.64, 195.73±31.29, 27.72±1.43, and 0.83±0.04 for MAE, RMSE, PSNR, and SSIM using the T1-POST dataset. The source model had the poorest performance.

CONCLUSIONS

This work indicates high generalization ability to generate synthetic CT images from small training datasets of MR images using pre-trained CycleGAN. The quantitative results of the test data, including different scanning protocols and different acquisition centers, indicated the proof of this concept.

Comprehensive methodology for commissioning modern 3D-image-based treatment planning systems for high dose rate gynaecological brachytherapy: A review

PURPOSE

Correct commissioning of treatment planning systems (TPSs) is important for reducing treatment failure events. There is currently no comprehensive and robust methodology available for TPS commissioning in modern brachytherapy. This review aimed to develop a comprehensive template for commissioning modern 3D-image-based brachytherapy TPSs for high dose rate (HDR) gynaecological applications.

METHODS

The literature relevant to TPS commissioning, including both external beam radiation therapy (EBRT) and brachytherapy, as well as guidelines by the International Atomic Energy Agency (IAEA), the American Association of Physicists in Medicine (AAPM), and the European Society for Radiotherapy and Oncology (ESTRO) were searched, studied and appraised. The applied relevant EBRT TPS commissioning tests were applied to brachytherapy. The developed template aimed to cover all dosimetric and non-dosimetric issues.

RESULTS

The essential commissioning items could be categorized into six parts: geometry, dose calculation, plan evaluation tools, plan optimization, TPS output, and end-to-end verification. The final template consists of 43 items. This paper presents the purpose and role of each test, as well as tolerance limits, to facilitate the use of the template.

CONCLUSION

The information and recommendations available in a collection of publications over many years have been reviewed in order to develop a comprehensive template for commissioning complex modern 3D-image-based brachytherapy TPSs for HDR gynaecological applications. The up-to-date and concise information contained in the template can aid brachytherapy physicists during TPS commissioning as well as devising a regular quality assurance program and allocation of time and resources.

MRI-only treatment planning: benefits and challenges

Over the past decade, the application of magnetic resonance imaging (MRI) has increased, and there is growing evidence to suggest that improvements in the accuracy of target delineation in MRI-guided radiation therapy may improve clinical outcomes in a variety of cancer types. However, some considerations should be recognized including patient motion during image acquisition and geometric accuracy of images. Moreover, MR-compatible immobilization devices need to be used when acquiring images in the treatment position while minimizing patient motion during the scan time. Finally, synthetic CT images (i.e. electron density maps) and digitally reconstructed radiograph images should be generated from MRI images for dose calculation and image guidance prior to treatment. A short review of the concepts and techniques that have been developed for implementation of MRI-only workflows in radiation therapy is provided in this document.

Technical Report: Evaluation of peripheral dose for flattening filter free photon beams

PURPOSE

To develop a comprehensive peripheral dose (PD) dataset for the two unflattened beams of nominal energy 6 and 10 MV for use in clinical care.

METHODS

Measurements were made in a 40 × 120 × 20 cm(3) (width × length × depth) stack of solid water using an ionization chamber at varying depths (dmax, 5, and 10 cm), field sizes (3 × 3 to 30 × 30 cm(2)), and distances from the field edge (5-40 cm). The effects of the multileaf collimator (MLC) and collimator rotation were also evaluated for a 10 × 10 cm(2) field. Using the same phantom geometry, the accuracy of the analytic anisotropic algorithm (AAA) and Acuros dose calculation algorithm was assessed and compared to the measured values.

RESULTS

The PDs for both the 6 flattening filter free (FFF) and 10 FFF photon beams were found to decrease with increasing distance from the radiation field edge and the decreasing field size. The measured PD was observed to be higher for the 6 FFF than for the 10 FFF for all field sizes and depths. The impact of collimator rotation was not found to be clinically significant when used in conjunction with MLCs. AAA and Acuros algorithms both underestimated the PD with average errors of -13.6% and -7.8%, respectively, for all field sizes and depths at distances of 5 and 10 cm from the field edge, but the average error was found to increase to nearly -69% at greater distances.

CONCLUSIONS

Given the known inaccuracies of peripheral dose calculations, this comprehensive dataset can be used to estimate the out-of-field dose to regions of interest such as organs at risk, electronic implantable devices, and a fetus. While the impact of collimator rotation was not found to significantly decrease PD when used in conjunction with MLCs, results are expected to be machine model and beam energy dependent. It is not recommended to use a treatment planning system to estimate PD due to the underestimation of the out-of-field dose and the inability to calculate dose at extended distances due to the limits of the dose calculation matrix.

Revisiting fetal dose during radiation therapy: evaluating treatment techniques and a custom shield

To create a comprehensive dataset of peripheral dose (PD) measurements from a new generation of linear accelerators with and without the presence of a newly designed fetal shield, PD measurements were performed to evaluate the effects of depth, field size, distance from the field edge, collimator angle, and beam modi-fiers for common treatment protocols and modalities. A custom fetal lead shield was designed and made for our department that allows external beam treatments from multiple angles while minimizing the need to adjust the shield during patient treatments. PD measurements were acquired for a comprehensive series of static fields on a stack of Solid Water. Additionally, PDs from various clinically relevant treatment scenarios for pregnant patients were measured using an anthropomorphic phantom that was abutted to a stack of Solid Water. As expected, the PD decreased as the distance from the field edge increased and the field size decreased. On aver-age, a PD reduction was observed when a 90° collimator rotation was applied and/or when the tertiary MLCs and jaws defined the field aperture. However, the effect of the collimator rotation (90° versus 0°) in PD reduction was not found to be clini-cally significant when the tertiary MLCs were used to define the field aperture. In the presence of both the MLCs and the fetal shield, the PD was reduced by 58% at a distance of 10 cm from the field edge. The newly designed fetal shield may effectively reduce fetal dose and is relatively easy to setup. Due to its design, we are able to use a broad range of treatment techniques and beam angles. We believe the acquired comprehensive PD dataset collected with and without the fetal shield will be useful for treatment teams to estimate fetal dose and help guide decisions on treat-ment techniques without the need to perform pretreatment phantom measurements.

A surface energy spectral study on the bone heterogeneity and beam obliquity using the flattened and unflattened photon beams

AIM

Using flattened and unflattened photon beams, this study investigated the spectral variations of surface photon energy and energy fluence in the bone heterogeneity and beam obliquity.

BACKGROUND

Surface dose enhancement is a dosimetric concern when using unflattened photon beam in radiotherapy. It is because the unflattened photon beam contains more low-energy photons which are removed by the flattening filter of the flattened photon beam.

MATERIALS AND METHODS

We used a water and bone heterogeneity phantom to study the distributions of energy, energy fluence and mean energy of the 6 MV flattened and unflattened photon beams (field size = 10 cm × 10 cm) produced by a Varian TrueBEAM linear accelerator. These elements were calculated at the phantom surfaces using Monte Carlo simulations. The photon energy and energy fluence calculations were repeated with the beam angle turned from 0° to 15°, 30° and 45° in the water and bone phantom.

RESULTS

Spectral results at the phantom surfaces showed that the unflattened photon beams contained more photons concentrated mainly in the low-energy range (0-2 MeV) than the flattened beams associated with a flattening filter. With a bone layer of 1 cm under the phantom surface and within the build-up region of the 6 MV photon beam, it is found that both the flattened and unflattened beams had slightly less photons in the energy range <0.4 MeV compared to the water phantom. This shows that the presence of the bone decreased the low-energy photon backscatters to the phantom surface. When both the flattened and unflattened photon beams were rotated from 0° to 45°, the number of photon and mean photon energy increased. This indicates that both photon beams became more hardened or penetrate when the beam angle increased. In the presence of bone, the mean energies of both photon beams increased. This is due to the absorption of low-energy photons by the bone, resulting in more beam hardening.

CONCLUSIONS

This study explores the spectral relationships of surface photon energy and energy fluence with bone heterogeneity and beam obliquity for the flattened and unflattened photon beams. The photon spectral information is important in studies on the patient's surface dose enhancement using unflattened photon beams in radiotherapy.

Clinical implementation of MR-guided vaginal cylinder brachytherapy

We present an institutional experience on the clinical implementation of magnetic resonance (MR)‐guided vaginal brachytherapy using commercially available solid applicator models. To test the fidelity of solid applicator models to digitize vaginal cylinder applicators, three datasets were evaluated. The first included 15 patients who were simulated with CT alone. Next, a water phantom was used to evaluate vaginal cylinders ranging from 20 to 35 mm in diameter with CT and MR. Finally, three patients undergoing HDR brachytherapy with vaginal cylinders that were simulated with both CT and MR were evaluated. In these assessments, the solid applicator models were aligned based on the outline of the applicators on the corresponding volumetric image, and deviations between the central source positions defined based on X‐ray markers (on CT) and solid applicator models (on CT and MR), and the percent dose difference between select reference points were calculated. The mean central source deviation defined based on X‐ray markers (on CT) and solid applicator models (on CT and MR) for the 15‐patient cohort, the phantom, and the 3‐patient cohort is 0.6 mm, 0.6 mm, and 1.2 mm, respectively. The average absolute percent dose difference for the bladder, rectum, prescription, and inferior reference points were 2.2%, 2.3%, 2.2%, and 2.4%, respectively, for the 15 patient cohort. For the phantom study, the average, absolute percent dose difference for the prescription and inferior reference points are 2.0% and 2.1% for the CT, 2.3% and 2.2% for the T1W, and 2.8% and 3.0% for the T2W images. For the three patient cohort, the average absolute percent dose difference for the bladder, rectum, prescription, and inferior reference points are 2.9%, 2.6%, 3.0%, and 4.2% for the CT, 6.5%, 1.6%, 2.5%, and 4.7% for the T1W, and 6.0%, 7.4%, 2.6, and 2.0% for the T2W images. Based on the current study, aligning the applicator model to MR images provides a practical, efficient approach to perform MR‐based brachytherapy planning.

PACS numbers: 87.53.Jw, 87.61.Tg

Semi-automated scoring of pulmonary emphysema from X-ray CT: trainee reproducibility and accuracy

OBJECTIVE

We developed a semi-automated tool to quantify emphysema from thoracic X-ray multi-detector (64-slice) computed tomography (CT) for training purposes and multi-reader studies.

MATERIALS AND METHODS

Thoracic X-ray CT was acquired in 93 ex-smokers, who were evaluated by six trainees with little or no expertise (trainees) and a single experienced thoracic radiologist (expert). A graphic user interface (GUI) was developed for emphysema quantification based on the percentile of lung where a score of 0=no abnormalities, 1=1-25%, 2=26-50%, 3=51-75% and 4=76-100% for each lung side/slice. Trainees blinded to subject characteristics scored randomized images twice; accuracy was determined by comparison to expert scores, density histogram 15th percentile (HU 15), relative area at -950 HU (RA(950)), low attenuation clusters at -950 HU (LAC(950)), -856 HU (LAC(856)) and the diffusing capacity for carbon monoxide (DL(CO%pred)). Intra- and inter-observer reproducibility was evaluated using coefficients-of-variation (COV), intra-class (ICC) and Pearson correlations.

RESULTS

Trainee-expert correlations were significant (r=0.85-0.97, p<0.0001) and a significant trainee bias (0.15 ± 0.22) was observed. Emphysema score was correlated with RA(950) (r=0.88, p<0.0001), HU 15 (r=-0.77, p<0.0001), LAC(950) (r=0.76, p<0.0001), LAC(856) (r=0.74, p=0.0001) and DLCO%pred (r=-0.71, p<0.0001). Intra-observer reproducibility (COV=4-27%; ICC=0.75-0.94) was moderate to high for trainees; intra- and inter-observer COV were negatively and non-linearly correlated with emphysema score.

CONCLUSION

We developed a GUI for rapid and interactive emphysema scoring that allows for comparison of multiple readers with clinical and radiological standards.

Computed tomography density histogram analysis to evaluate pulmonary emphysema in ex-smokers

RATIONALE AND OBJECTIVES

High-resolution computed tomography (CT) measurements of emphysema typically use Hounsfield unit (HU) density histogram thresholds or observer scores based on regions of low x-ray attenuation. Our objective was to develop an automated measurement of emphysema using principal component analysis (PCA) of the CT density histogram.

MATERIALS AND METHODS

Ninety-seven ex-smokers, including 53 subjects with chronic obstructive pulmonary disease (COPD) and 44 asymptomatic subjects (AEs), provided written informed consent to imaging as well as plethysmography and spirometry. We applied PCA to the CT density histogram to generate whole lung and regional density histogram principal components including the first and second components and the sum of both principal components (density histogram principal component score [DHPCS]). Significant relationships for DHPCS with single HU thresholds, pulmonary function measurements, an expert's emphysema score, and hyperpolarized (3)He magnetic resonance imaging apparent diffusion coefficients (ADCs) were determined using linear regression and Pearson coefficients. Receiver operator characteristics analysis was performed using forced expiratory volume in 1 second (FEV1)/forced vital capacity (FVC) as the independent diagnostic.

RESULTS

There was a significant difference (P < .0001) between AE and COPD subjects for DHPCS; FEV1/FVC; diffusing capacity of lung for carbon monoxide%predicted; attenuation values below -950, -910, and -856 HU; and (3)He ADCs. There were significant correlations for DHPCS with FEV1/FVC (r = -0.85, P < .0001); diffusing capacity of lung for carbon monoxide%predicted (r = -0.67, P < .0001); attenuation values below -950/-910/-856 HU (r = 0.93/0.96/0.76, P < .0001); and (3)He ADCs (r = 0.85, P < .0001). Receiver operator characteristics analysis showed a 91% classification rate for DHPCS.

CONCLUSIONS

We generated an automated emphysema score using PCA of the CT density histogram with a 91% COPD classification rate that showed strong and significant correlations with pulmonary function tests, single HU thresholds, and (3)He magnetic resonance imaging ADCs.

Dosimetry of small bone joint calculated by the analytical anisotropic algorithm: a Monte Carlo evaluation using the EGSnrc

This study compared a small bone joint dosimetry calculated by the anisotropic analytical algorithm (AAA) and Monte Carlo simulation using megavoltage (MV) photon beams. The performance of the AAA in the joint dose calculation was evaluated using Monte Carlo simulation, and dependences of joint dose on its width and beam angle were investigated. Small bone joint phantoms containing a vertical water layer (0.5‐2 mm) sandwiched by two bones (2×2×2cm3) were irradiated by the 6 and 15 MV photon beams with field size equal to 4×4 cm2. Depth doses along the central beam axis in a joint (cartilage) were calculated with and without a bolus (thickness=1.5cm) added on top of the phantoms. Different beam angles (0°‐15°) were used with the isocenter set to the center of the bone joint for dose calculations using the AAA (Eclipse treatment planning system) and Monte Carlo simulation (the EGSnrc code). For dosimetry comparison and normalization, dose calculations were repeated in homogeneous water phantoms with the bone substituted by water. Comparing the calculated dosimetry between the AAA and Monte Carlo simulation, the AAA underestimated joint doses varying with its widths by about 6%‐12% for 6 MV and 12%‐23% for 15 MV without bolus, and by 7% for 6 MV and 13%‐17% for 15 MV with bolus. Moreover, joint doses calculated by the AAA did not vary with the joint width and beam angle. From Monte Carlo results, there was a decrease in the calculated joint dose as the joint width increased, and a slight decrease as the beam angle increased. When bolus was added to the phantom, it was found that variations of joint dose with its width and beam angle became less significant for the 6 MV photon beams. In conclusion, dosimetry deviation in small bone joint calculated by the AAA and Monte Carlo simulation was studied using the 6 and 15 MV photon beam. The AAA could not predict variations of joint dose with its width and beam angle, which were predicted by the Monte Carlo simulations.

PACS numbers: 87.55.K‐; 87.53.Bn; 87.53.‐j

Surface dose reduction from bone interface in kilovoltage X-ray radiation therapy: a Monte Carlo study of photon spectra

This study evaluated the dosimetric impact of surface dose reduction due to the loss of backscatter from the bone interface in kilovoltage (kV) X‐ray radiation therapy. Monte Carlo simulation was carried out using the EGSnrc code. An inhomogeneous phantom containing a thin water layer (0.5–5 mm) on top of a bone (thickness=1 cm) was irradiated by a clinical 105 kVp photon beam produced by a Gulmay D3225 X‐ray machine. Field sizes of 2, 5, and 10 cm diameter and source‐to‐surface distance of 20 cm were used. Surface doses for different phantom configurations were calculated using the DOSXYZnrc code. Photon energy spectra at the phantom surface and bone were determined according to the phase‐space files at the particle scoring planes which included the multiple crossers. For comparison, all Monte Carlo simulations were repeated in a phantom with the bone replaced by water. Surface dose reduction was found when a bone was underneath the water layer. When the water thickness was equal to 1 mm for the circular field of 5 cm diameter, a surface dose reduction of 6.3% was found. The dose reduction decreased to 4.7% and 3.4% when the water thickness increased to 3 and 5 mm, respectively. This shows that the impact of the surface dose uncertainty decreased while the water thickness over the bone increased. This result was supported by the decrease in relative intensity of the lower energy photons in the energy spectrum when the water layer was with and over the bone, compared to without the bone. We concluded that surface dose reduction of 7.8%–1.1% was found when the water thickness increased from 0.5–5 mm for circular fields with diameters ranging from 2–10 cm. This decrease of surface dose results in an overestimation of prescribed dose at the patient's surface, and might be a concern when using kV photon beam to treat skin tumors in sites such as forehead, chest wall, and kneecap.

PACS number: 87.55.K‐; 87.55.ne; 87.57.uq

Pulmonary tumor measurements from x-ray computed tomography in one, two, and three dimensions

RATIONALE AND OBJECTIVES

We evaluated the accuracy and reproducibility of three-dimensional (3D) measurements of lung phantoms and patient tumors from x-ray computed tomography (CT) and compared these to one-dimensional (1D) and two-dimensional (2D) measurements.

MATERIALS AND METHODS

CT images of three spherical and three irregularly shaped tumor phantoms were evaluated by three observers who performed five repeated measurements. Additionally, three observers manually segmented 29 patient lung tumors five times each. Follow-up imaging was performed for 23 tumors and response criteria were compared. For a single subject, imaging was performed on nine occasions over 2 years to evaluate multidimensional tumor response. To evaluate measurement accuracy, we compared imaging measurements to ground truth using analysis of variance. For estimates of precision, intraobserver and interobserver coefficients of variation and intraclass correlations (ICC) were used. Linear regression and Pearson correlations were used to evaluate agreement and tumor response was descriptively compared.

RESULTS

For spherical shaped phantoms, all measurements were highly accurate, but for irregularly shaped phantoms, only 3D measurements were in high agreement with ground truth measurements. All phantom and patient measurements showed high intra- and interobserver reproducibility (ICC >0.900). Over a 2-year period for a single patient, there was disagreement between tumor response classifications based on 3D measurements and those generated using 1D and 2D measurements.

CONCLUSION

Tumor volume measurements were highly reproducible and accurate for irregular, spherical phantoms and patient tumors with nonuniform dimensions. Response classifications obtained from multidimensional measurements suggest that 3D measurements provide higher sensitivity to tumor response.

Dependences of mucosal dose on photon beams in head-and-neck intensity-modulated radiation therapy: a Monte Carlo study

Dependences of mucosal dose in the oral or nasal cavity on the beam energy, beam angle, multibeam configuration, and mucosal thickness were studied for small photon fields using Monte Carlo simulations (EGSnrc-based code), which were validated by measurements. Cylindrical mucosa phantoms (mucosal thickness = 1, 2, and 3 mm) with and without the bone and air inhomogeneities were irradiated by the 6- and 18-MV photon beams (field size = 1 × 1 cm(2)) with gantry angles equal to 0°, 90°, and 180°, and multibeam configurations using 2, 4, and 8 photon beams in different orientations around the phantom. Doses along the central beam axis in the mucosal tissue were calculated. The mucosal surface doses were found to decrease slightly (1% for the 6-MV photon beam and 3% for the 18-MV beam) with an increase of mucosal thickness from 1-3 mm, when the beam angle is 0°. The variation of mucosal surface dose with its thickness became insignificant when the beam angle was changed to 180°, but the dose at the bone-mucosa interface was found to increase (28% for the 6-MV photon beam and 20% for the 18-MV beam) with the mucosal thickness. For different multibeam configurations, the dependence of mucosal dose on its thickness became insignificant when the number of photon beams around the mucosal tissue was increased. The mucosal dose with bone was varied with the beam energy, beam angle, multibeam configuration and mucosal thickness for a small segmental photon field. These dosimetric variations are important to consider improving the treatment strategy, so the mucosal complications in head-and-neck intensity-modulated radiation therapy can be minimized.

Solid water as phantom material for dosimetry of electron backscatter using low-energy electron beams: a Monte Carlo evaluation

This study evaluated the dosimetry of electron backscatter when Solid Water is used to substitute water as phantom in electron radiotherapy. Monte Carlo simulation (EGSnrc-based code) was employed to predict electron energy spectra and depth doses for the 0.5 and 1 cm of Solid Water and water slabs above 3 mm of lead (Pb) layers using electron beams with energies of 4 and 6 MeV. For comparison, Monte Carlo simulations were repeated with Pb layers taken out from the phantoms using the same experimental configuration. Analyses on electron energy spectra for the 4 and 6 MeV electron beams showed that deviations of electron energy distributions between the Solid Water and water phantom were more significant in the high-energy range (i.e., close to the maximal electron energy) than the lower range corresponding to the electron backscatter. These deviations of electron energy spectra varied with depth and were mainly due to the electron fluence or beam attenuation. Dosimetry results from Monte Carlo simulations showed that the Solid Water phantom had lower depth dose compared to water with the same experimental setup. For the 4 MeV electron beams with 0.5 cm of Solid Water, depth doses were 1.8%-3.9% and 2.3%-4.4% lower than those in water, with and without the Pb layer underneath, respectively. Thicker Solid Water of 1 cm resulted in different decreases in depth doses of 1.8%-4.6% (with Pb) and 2.3%-4.4% (without Pb) compared to water. For higher nominal electron beam energy of 6 MeV with 0.5 cm of Solid Water, depth doses decreased 1.7%-2.9% (with Pb) and 1.6%-2.1% (without Pb) compared to water. These decreases in depth doses changed to 1.7%-3.7% (with Pb) and 1.7%-3% (without Pb) when the thickness of Solid Water was increased to 1 cm. The dosimetry data in this study are useful in determining the correction factor when using Solid Water to substitute water for the electron backscatter measurement in electron radiotherapy.