A joint physics and radiobiology DREAM team vision - Towards better response prediction models to advance radiotherapy.

Radiotherapy developed empirically through experience balancing tumour control and normal tissue toxicities. Early simple mathematical models formalized this practical knowledge and enabled effective cancer treatment to date. Remarkable advances in technology, computing, and experimental biology now create opportunities to incorporate this knowledge into enhanced computational models. The ESTRO DREAM (Dose Response, Experiment, Analysis, Modelling) workshop brought together experts across disciplines to pursue the vision of personalized radiotherapy for optimal outcomes through advanced modelling. The ultimate vision is leveraging quantitative models dynamically during therapy to ultimately achieve truly adaptive and biologically guided radiotherapy at the population as well as individual patient-based levels. This requires the generation of models that inform response-based adaptations, individually optimized delivery and enable biological monitoring to provide decision support to clinicians. The goal is expanding to models that can drive the realization of personalized therapy for optimal outcomes. This position paper provides their propositions that describe how innovations in biology, physics, mathematics, and data science including AI could inform models and improve predictions. It consolidates the DREAM team's consensus on scientific priorities and organizational requirements. Scientifically, it stresses the need for rigorous, multifaceted model development, comprehensive validation and clinical applicability and significance. Organizationally, it reinforces the prerequisites of interdisciplinary research and collaboration between physicians, medical physicists, radiobiologists, and computational scientists throughout model development. Solely by a shared understanding of clinical needs, biological mechanisms, and computational methods, more informed models can be created. Future research environment and support must facilitate this integrative method of operation across multiple disciplines.

Optimal beam angle selection and knowledge-based planning significantly reduces radiotherapy dose to organs at risk for lung cancer patients.

BACKGROUND: Lung cancer patients struggle with high toxicity rates. This study investigates if IMRT plans with individually set beam angles or uni-lateral VMAT plans results in dose reduction to OARs. We investigate if introduction of a RapidPlan model leads to reduced dose to OARs. Finally, the model is validated prospectively. MATERIAL AND METHODS: Seventy-four consecutive lung cancer patients treated with IMRT were included. For all patients, new IMRT plans were made by an experienced dose planner re-tuning beam angles aiming for minimized dose to the lungs and heart. Additionally, VMAT plans were made. The IMRT plans were selected as input for a RapidPlan model, which was used to generate 74 new IMRT plans. The new IMRT plans were used as input for a second RapidPlan model. This model was clinically implemented and used for generation of clinical treatment plans. Dosimetric parameters were compared using a Wilcoxon signed rank test or a 1-sided student's t-test. p < .05 was considered significant. RESULTS: IMRT plans significantly reduced mean doses to lungs (MLD) and heart (MHD) by 1.6 Gy and 1.7 Gy in mean compared to VMAT plans. MLD was significantly (p < .001) reduced from 10.8 Gy to 9.4 Gy by using the second RapidPlan model. MHD was significantly (p < .001) reduced from 4.9 Gy to 3.9 Gy. The model was validated in prospectively collected treatment plans showing significantly lower MLD after the implementation of the second RapidPlan model. CONCLUSION: Introduction of RapidPlan and beam angles selected based on the target and OARs position reduces dose to OARs.

Energy layer optimization strategies for intensity-modulated proton therapy of lung cancer patients.

PURPOSE: When treating lung cancer patients with intensity-modulated proton therapy (IMPT), target coverage can only be guaranteed when utilizing motion mitigation. The three motion mitigation techniques, gating, breath-hold, and dose repainting, all benefit from a more rapid application of the treatment plan. A lower limit for the ungated treatment time is defined by the number of energy layers in the IMPT plan. By limiting this number during treatment planning, IMPT could become more viable for lung cancer patients. We investigate to what extend the number of layers can be reduced in single-field optimization (SFO) and multifield optimization (MFO) plans and which implications it has on the plan quality and robustness. METHODS: We have implemented three distinct layer-reducing strategies in the treatment planning system Hyperion; constant energy steps, exponential energy steps, and an adaptive strategy, where the spot weights are exposed to a group sparsity penalty in combination with layer exclusion during optimization. Four levels of increasing layer removal are planned for each strategy. SFO and MFO plans with three treatment fields are created for eleven locally advanced NSCLC patients on the midventilation 4DCT phase to simulate a breath-hold. A minimum dose to the target is ensured for each degree of layer reduction, reflecting the plan quality in the homogeneity index (HI). Plan quality was also assessed by a robustness evaluation, where the patient setup was shifted 2 mm or 4 mm in six directions. RESULTS: The three strategies result in very similar target coverages and robustness levels as a function of removed layers. The HI increases unacceptably for all the SFO plans after 50% layer removal as compared to the reference plan, while all the MFO plans are clinically acceptable with up to a highest removed percentage of 75%. The robustness level is constant as a function of removed layers. The SFO plans are significantly more robust than the MFO plans with all P-values below 0.001 (Wilcoxon signed-rank). The overall mean D98% CTV dose difference is at 2-mm setup error amplitude: 0.7 Gy (SFO) and 1.9 Gy (MFO), and at 4 mm: 3.2 Gy (SFO) and 5.4 Gy (MFO), respectively. CONCLUSIONS: The number of layers in MFO plans can be reduced substantially more than in SFO plans without compromising plan quality. Furthermore, as the robustness is independent of the number of layers, it follows that if the level of robustness is acceptable or enforced via robust optimization, MFO plans could be candidates for treatment time reductions via energy layer reductions.

Cadherin composition and multicellular aggregate invasion in organotypic models of epithelial ovarian cancer intraperitoneal metastasis.

During epithelial ovarian cancer (EOC) progression, intraperitoneally disseminating tumor cells and multicellular aggregates (MCAs) present in ascites fluid adhere to the peritoneum and induce retraction of the peritoneal mesothelial monolayer prior to invasion of the collagen-rich submesothelial matrix and proliferation into macro-metastases. Clinical studies have shown heterogeneity among EOC metastatic units with respect to cadherin expression profiles and invasive behavior; however, the impact of distinct cadherin profiles on peritoneal anchoring of metastatic lesions remains poorly understood. In the current study, we demonstrate that metastasis-associated behaviors of ovarian cancer cells and MCAs are influenced by cellular cadherin composition. Our results show that mesenchymal N-cadherin-expressing (Ncad+) cells and MCAs invade much more efficiently than E-cadherin-expressing (Ecad+) cells. Ncad+ MCAs exhibit rapid lateral dispersal prior to penetration of three-dimensional collagen matrices. When seeded as individual cells, lateral migration and cell-cell junction formation precede matrix invasion. Neutralizing the Ncad extracellular domain with the monoclonal antibody GC-4 suppresses lateral dispersal and cell penetration of collagen gels. In contrast, use of a broad-spectrum matrix metalloproteinase (MMP) inhibitor (GM6001) to block endogenous membrane type 1 matrix metalloproteinase (MT1-MMP) activity does not fully inhibit cell invasion. Using intact tissue explants, Ncad+ MCAs were also shown to efficiently rupture peritoneal mesothelial cells, exposing the submesothelial collagen matrix. Acquisition of Ncad by Ecad+ cells increased mesothelial clearance activity but was not sufficient to induce matrix invasion. Furthermore, co-culture of Ncad+ with Ecad+ cells did not promote a 'leader-follower' mode of collective cell invasion, demonstrating that matrix remodeling and creation of invasive micro-tracks are not sufficient for cell penetration of collagen matrices in the absence of Ncad. Collectively, our data emphasize the role of Ncad in intraperitoneal seeding of EOC and provide the rationale for future studies targeting Ncad in preclinical models of EOC metastasis.

Reliability of dose volume constraint inference from clinical data.

Dose volume histogram points (DVHPs) frequently serve as dose constraints in radiotherapy treatment planning. An experiment was designed to investigate the reliability of DVHP inference from clinical data for multiple cohort sizes and complication incidence rates. The experimental background was radiation pneumonitis in non-small cell lung cancer and the DVHP inference method was based on logistic regression. From 102 NSCLC real-life dose distributions and a postulated DVHP model, an 'ideal' cohort was generated where the most predictive model was equal to the postulated model. A bootstrap and a Cohort Replication Monte Carlo (CoRepMC) approach were applied to create 1000 equally sized populations each. The cohorts were then analyzed to establish inference frequency distributions. This was applied to nine scenarios for cohort sizes of 102 (1), 500 (2) to 2000 (3) patients (by sampling with replacement) and three postulated DVHP models. The Bootstrap was repeated for a 'non-ideal' cohort, where the most predictive model did not coincide with the postulated model. The Bootstrap produced chaotic results for all models of cohort size 1 for both the ideal and non-ideal cohorts. For cohort size 2 and 3, the distributions for all populations were more concentrated around the postulated DVHP. For the CoRepMC, the inference frequency increased with cohort size and incidence rate. Correct inference rates >[Formula: see text] were only achieved by cohorts with more than 500 patients. Both Bootstrap and CoRepMC indicate that inference of the correct or approximate DVHP for typical cohort sizes is highly uncertain. CoRepMC results were less spurious than Bootstrap results, demonstrating the large influence that randomness in dose-response has on the statistical analysis.

Compression-induced structural and mechanical changes of fibrin-collagen composites.

Fibrin and collagen as well as their combinations play an important biological role in tissue regeneration and are widely employed in surgery as fleeces or sealants and in bioengineering as tissue scaffolds. Earlier studies demonstrated that fibrin-collagen composite networks displayed improved tensile mechanical properties compared to the isolated protein matrices. Unlike previous studies, here unconfined compression was applied to a fibrin-collagen filamentous polymer composite matrix to study its structural and mechanical responses to compressive deformation. Combining collagen with fibrin resulted in formation of a composite hydrogel exhibiting synergistic mechanical properties compared to the isolated fibrin and collagen matrices. Specifically, the composite matrix revealed a one order of magnitude increase in the shear storage modulus at compressive strains>0.8 in response to compression compared to the mechanical features of individual components. These material enhancements were attributed to the observed structural alterations, such as network density changes, an increase in connectivity along with criss-crossing, and bundling of fibers. In addition, the compressed composite collagen/fibrin networks revealed a non-linear transformation of their viscoelastic properties with softening and stiffening regimes. These transitions were shown to depend on protein concentrations. Namely, a decrease in protein content drastically affected the mechanical response of the networks to compression by shifting the onset of stiffening to higher degrees of compression. Since both natural and artificially composed extracellular matrices experience compression in various (patho)physiological conditions, our results provide new insights into the structural biomechanics of the polymeric composite matrix that can help to create fibrin-collagen sealants, sponges, and tissue scaffolds with tunable and predictable mechanical properties.

Vascular perfusion of the dorsal and palmar condyles of the equine third metacarpal bone.

REASONS FOR PERFORMING STUDY: Palmar osteochondral disease (POD) is an overload arthrosis that commonly affects fetlock joints of racing Thoroughbreds (TB) but the aetiopathogenesis of the disease has not been well defined. OBJECTIVES: The aim of this study was to compare India ink perfusion in the dorsal and palmar condyles of the third metacarpal bone (McIII) in both passively flexed and maximally extended fetlock joints from paired equine cadaver limbs. STUDY DESIGN: Descriptive cadaver study comparing perfusion of condyles of McIII in paired cadaver limbs in flexion (control group) and maximal extension (intervention group). METHODS: Pairs of forelimbs were acquired from 5 TB horses subjected to euthanasia for reasons unrelated to lameness. Limb pairs were perfused intra-arterially with India ink and then randomly assigned to passive flexion or maximal extension of the fetlock joint. Limbs were sectioned sagittally in 3 mm sections through the fetlock and 12 sections per limb processed using a modified tissue-clearing technique. Sections were subsequently digitally imaged and bone perfusion evaluated with image analysis software. RESULTS: Greater perfusion of the dorsal condyle than of palmar condyle was observed in 78% of sections from limbs in passive flexion and 92% of maximally extended sections. Perfusion to the palmar aspect of the condyle was significantly decreased (P < 0.0001) when the limbs were placed in maximal extension compared to passive flexion. CONCLUSIONS: The palmar condyle of McIII had less perfusion than the dorsal condyle when the fetlock joint was in passive flexion and this difference was exacerbated by maximal extension. Based on the anatomical location of POD lesions, perfusion differences between the dorsal and palmar condyles of McIII may be associated with development of these lesions.

Accelerated evaluation of the robustness of treatment plans against geometric uncertainties by Gaussian processes.

In order to provide a consistently high quality treatment, it is of great interest to assess the robustness of a treatment plan under the influence of geometric uncertainties. One possible method to implement this is to run treatment simulations for all scenarios that may arise from these uncertainties. These simulations may be evaluated in terms of the statistical distribution of the outcomes (as given by various dosimetric quality metrics) or statistical moments thereof, e.g. mean and/or variance. This paper introduces a method to compute the outcome distribution and all associated values of interest in a very efficient manner. This is accomplished by substituting the original patient model with a surrogate provided by a machine learning algorithm. This Gaussian process (GP) is trained to mimic the behavior of the patient model based on only very few samples. Once trained, the GP surrogate takes the place of the patient model in all subsequent calculations.The approach is demonstrated on two examples. The achieved computational speedup is more than one order of magnitude.

Toxicity and outcome of pelvic IMRT for node-positive prostate cancer.

BACKGROUND AND PURPOSE: This study reports on the treatment techniques, toxicity, and outcome of pelvic intensity-modulated radiotherapy (IMRT) for lymph node-positive prostate cancer (LNPPC, T1-4, c/pN1 cM0). PATIENTS AND METHODS: Pelvic IMRT to 45-50.4 Gy was applied in 39 cases either after previous surgery of involved lymph nodes (n = 18) or with a radiation boost to suspicious nodes (n = 21) with doses of 60-70 Gy, usually combined with androgen deprivation (n = 37). The prostate and seminal vesicles received 70-74 Gy. In cases of previous prostatectomy, prostatic fossa and remnants of seminal vesicles were given 66-70 Gy. Treatment-related acute and late toxicity was graded according to the RTOG criteria. RESULTS: Acute radiation-related toxicity higher than grade 2 occurred in 2 patients (with the need for urinary catheter/subileus related to adhesions after surgery). Late toxicity was mild (grade 1-2) after a median follow-up of 70 months. Over 50% of the patients reported no late morbidity (grade 0). PSA control and cancer-specific survival reached 67% and 97% at over 5 years. CONCLUSION: Pelvic IMRT after the removal of affected nodes or with a radiation boost to clinically positive nodes led to an acceptable late toxicity (no grade 3/4 events), thus justifying further evaluation of this approach in a larger cohort.

Dosimetric treatment course simulation based on a statistical model of deformable organ motion.

We present a method of modeling dosimetric consequences of organ deformation and correlated motion of adjacent organ structures in radiotherapy. Based on a few organ geometry samples and the respective deformation fields as determined by deformable registration, principal component analysis (PCA) is used to create a low-dimensional parametric statistical organ deformation model (Söhn et al 2005 Phys. Med. Biol. 50 5893-908). PCA determines the most important geometric variability in terms of eigenmodes, which represent 3D vector fields of correlated organ deformations around the mean geometry. Weighted sums of a few dominating eigenmodes can be used to simulate synthetic geometries, which are statistically meaningful inter- and extrapolations of the input geometries, and predict their probability of occurrence. We present the use of PCA as a versatile treatment simulation tool, which allows comprehensive dosimetric assessment of the detrimental effects that deformable geometric uncertainties can have on a planned dose distribution. For this, a set of random synthetic geometries is generated by a PCA model for each simulated treatment course, and the dose of a given treatment plan is accumulated in the moving tissue elements via dose warping. This enables the calculation of average voxel doses, local dose variability, dose-volume histogram uncertainties, marginal as well as joint probability distributions of organ equivalent uniform doses and thus of TCP and NTCP, and other dosimetric and biologic endpoints. The method is applied to the example of deformable motion of prostate/bladder/rectum in prostate IMRT. Applications include dosimetric assessment of the adequacy of margin recipes, adaptation schemes, etc, as well as prospective 'virtual' evaluation of the possible benefits of new radiotherapy schemes.

SU-E-T-490: Comparison of XVMC Monte Carlo Dose Calculations with Eclipse AAA Calculations for RapidArc Plans.

PURPOSE: The consistency between the AAA and XVMC algorithm in the treatment planning for RapidArc is investigated. While the majority of the radiation field is blocked by the MLC system, multiple small dose islands with MLC opened only slightly can be observed in one control point. This raises questions on how accurate the clinically used AAA algorithm in Eclipse is able to calculate RapidArc dose distributions. The fast Monte Carlo Code XVMC was used as a benchmark to test the AAA algorithm. METHODS: RadpidArc plans of 25 patients were calculated with AAA and XVMC. The patient cohort consisted of 4 different cancer sites (H&N, upper abdominal, lung, prostate). Dose distributions, PTV and OAR coverage were compared looking at the PTV mean dose Dmean, the volume V95% of the PTV receiving 95% of the prescribed dose, the dose D95% delivered to 95% of the PTV Volume, the percentage PTV mean dose with respect to the prescribed dose Dmean/prescr and OAR mean dose. RESULTS: The recalculation of RapidArc plans yielded good agreement of both calculation algorithms for treatment plans of all four cancer sites. PTV mean dose differences of AAA and XVMC were found to be in between -0.11% and 4.89% of the prescribed dose. The mean dose difference found was 0.48±0.77 Gy. Local dose differences were found when comparing dose distributions in regions of big mass density differences and in high dose regions. One head and neck plan and one prostate plan revealed significant differences in PTV coverage (ΔDmean=3.25 Gy) and OAR mean dose (prostate mean dose -13.71 Gy) respectively. CONCLUSIONS: The vast majority of treatment plans calculated with the AAA algorithm were found to agree within the expected and acceptable tolerances compared to XVMC results. Nevertheless in some cases dose differences were observed that could be of clinical significance. This work was funded by a Varian grant. Wolfram Laub is working in the physics group of CMS.

[Neonatal intracerebral bleeding as initial symptom of gliomatosis cerebri WHO III--neurological outcome after partial hemispherectomy]. / Eine neonatale intrazerebrale Blutung als Initialsymptom einer Gliomatosis cerebri WHO III--neurologisches Outcome nach partieller Hemisphärektomie.

INTRODUCTION: Anaplastic astrocytomas presenting as gliomatosis cerebri in neonates are extremely rare. Tumours in newborns are mostly of neuroectodermal origin. CASE REPORT: We report about a female newborn at term [birth weight 3 600 g (P 90), head circumference 35 cm (P 95) APGAR 9/10/10] with an intracerebral partially clotted bleeding in the left parieto-occipital region. The bleeding was expansive leading to axial and lateral cerebral herniation. The intracerebral bleeding in the left occipital region was surgically removed. Macroscopically no solid tumour was seen, but small fragments of an anaplastic astrocytic tumour (WHO grade III) were diagnosed histologically. After surgery, no remaining tumour was visible in the MRI. 6 weeks later, a recurrent tumour (4×4 cm) was found in the area of the initial bleeding. Further treatment was initially refused by the parents. The child was readmitted to our hospital at the age of 11 months in good clinical condition and presented with left-sided hemiparesis, right-sided hemianopsia and intermittent strabismus convergens alternans. Because of the good clinical condition further therapeutic treatment was initiated. Due to the final extension of the tumour into the temporal, parietal and occipital lobes, a gliomatosis cerebri WHO III was diagnosed. An extended partial hemispherectomy was done. After neurosurgery, no further neurological failures occurred. In the follow-up examination, MRI showed no relapse of the tumour. Chemotherapy according to the HIT SKK protocol was initiated. A relapse did not occur over a follow-up of 2 years. CONCLUSION: This is a rare case report of a congenital gliomatosis cerebri WHO grade III, treated with partial hemispherectomy, leading to a good clinical and neurological long-term outcome.

On expedient properties of common biological score functions for multi-modality, adaptive and 4D dose optimization.

Frequently, radiotherapy treatments are comprised of several dose distributions computed or optimized in different patient geometries. Therefore, the need arises to compute the comprehensive biological effect or physical figure of merit of the combined dose of a number of distinct geometry instances. For that purpose the dose is typically accumulated in a reference geometry through deformation fields obtained from deformable image registration. However, it is difficult to establish precise voxel-by-voxel relationships between different anatomical images in many cases. In this work, the mathematical properties of commonly used score functions are exploited to derive an upper boundary for the maximum effect for normal tissue and a lower boundary for the minimum effect for the target of accumulated doses on multiple geometry instances.

Robust optimization based upon statistical theory.

PURPOSE: Organ movement is still the biggest challenge in cancer treatment despite advances in online imaging. Due to the resulting geometric uncertainties, the delivered dose cannot be predicted precisely at treatment planning time. Consequently, all associated dose metrics (e.g., EUD and maxDose) are random variables with a patient-specific probability distribution. The method that the authors propose makes these distributions the basis of the optimization and evaluation process. METHODS: The authors start from a model of motion derived from patient-specific imaging. On a multitude of geometry instances sampled from this model, a dose metric is evaluated. The resulting pdf of this dose metric is termed outcome distribution. The approach optimizes the shape of the outcome distribution based on its mean and variance. This is in contrast to the conventional optimization of a nominal value (e.g., PTV EUD) computed on a single geometry instance. The mean and variance allow for an estimate of the expected treatment outcome along with the residual uncertainty. Besides being applicable to the target, the proposed method also seamlessly includes the organs at risk (OARs). RESULTS: The likelihood that a given value of a metric is reached in the treatment is predicted quantitatively. This information reveals potential hazards that may occur during the course of the treatment, thus helping the expert to find the right balance between the risk of insufficient normal tissue sparing and the risk of insufficient tumor control. By feeding this information to the optimizer, outcome distributions can be obtained where the probability of exceeding a given OAR maximum and that of falling short of a given target goal can be minimized simultaneously. CONCLUSIONS: The method is applicable to any source of residual motion uncertainty in treatment delivery. Any model that quantifies organ movement and deformation in terms of probability distributions can be used as basis for the algorithm. Thus, it can generate dose distributions that are robust against interfraction and intrafraction motion alike, effectively removing the need for indiscriminate safety margins.

Estimation of protein and domain interactions in the switching motility system of Myxococcus xanthus.

The gram-negative myxobacterium Myxococcus xanthus is equipped with an interesting motility system that allows it to reverse direction on average every 8 minutes by switching the construction of two motility engines at the ends of this rod-shaped bacterium. While the mechanisms responsible for timing and engine construction/deconstruction are relatively well understood, there are several competing hypotheses as to how they are coupled together. In this paper we examine the evidence for protein interactions underlying these possible couplings using a novel framework consisting of a probabilistic model describing protein and domain interactions and a belief propagation inference algorithm. When provided with large amount of indirect pieces of information, such as high-throughput experiment results, and protein structures, we can reliably determine the relative likelihoods of these hypotheses, even though each individual piece of evidence by itself has very limited reliability. The same framework can be used to map large protein and domain interaction networks in myxobacteria and other organisms.

A virtual source model of electron contamination of a therapeutic photon beam.

The most efficient way of generating particles for Monte Carlo (MC) dose calculation is through a virtual source model (VSM) of the linear accelerator head. We have previously developed a VSM based on three sources: a primary photon source, a secondary photon source and an electron contamination source (Sikora et al 2007). In this work, we present an improvement of the electron contamination source. The VSM of contamination electrons (eVSM) is derived from a full MC simulation of the accelerator head with the BEAMnrc MC system. It comprises a Gaussian source located at the base of the flattening filter. The eVSM models two effects: an energy-dependent source diameter and an angular dependence of the particle fluence. The air scatter of the contamination electrons is approximated by energetic properties of the eVSM so that explicit in-air transport is not required during MC simulation of the dose distributions in the patient. The calculations of electron dose distributions were compared between the eVSM and the full MC simulation. Good agreement was achieved for various rectangular field sizes as well as for complex conformal segment shapes for the contamination electrons of 6 and 15 MV beams. The 3D dose evaluation of the surface dose in a CT-based patient geometry shows high accuracy (2%/2 mm) of the eVSM for both energies. The model has one tunable parameter, the mean energy of the spectrum at the patient surface. High accuracy and efficiency of particle generation make the eVSM a valuable virtual source of contamination electrons for MC treatment planning systems.

Tools for the analysis of dose optimization: III. Pointwise sensitivity and perturbation analysis.

The major challenge in intensity-modulated radiotherapy planning is to find the right balance between tumor control and normal tissue sparing. The most desirable solution is never physically feasible, and a compromise has to be found. One possible way to approach this problem is constrained optimization. In this context, it is worthwhile to quantitatively predict the impact of adjustments of the constraints on the optimum dose distribution. This has been dealt with in regard to cost functions in a previous paper. The aim of the present paper is to introduce spatial resolution to this formalism. Our method reveals the active constraints in a target subvolume that was previously selected by the practitioner for its insufficient dose. This is useful if a multitude of constraints can be the cause of a cold spot. The response of the optimal dose distribution to an adjustment of constraints (perturbation) is predicted. We conclude with a clinical example.

A virtual photon source model of an Elekta linear accelerator with integrated mini MLC for Monte Carlo based IMRT dose calculation.

A dedicated, efficient Monte Carlo (MC) accelerator head model for intensity modulated stereotactic radiosurgery treatment planning is needed to afford a highly accurate simulation of tiny IMRT fields. A virtual source model (VSM) of a mini multi-leaf collimator (MLC) (the Elekta Beam Modulator (EBM)) is presented, allowing efficient generation of particles even for small fields. The VSM of the EBM is based on a previously published virtual photon energy fluence model (VEF) (Fippel et al 2003 Med. Phys. 30 301) commissioned with large field measurements in air and in water. The original commissioning procedure of the VEF, based on large field measurements only, leads to inaccuracies for small fields. In order to improve the VSM, it was necessary to change the VEF model by developing (1) a method to determine the primary photon source diameter, relevant for output factor calculations, (2) a model of the influence of the flattening filter on the secondary photon spectrum and (3) a more realistic primary photon spectrum. The VSM model is used to generate the source phase space data above the mini-MLC. Later the particles are transmitted through the mini-MLC by a passive filter function which significantly speeds up the time of generation of the phase space data after the mini-MLC, used for calculation of the dose distribution in the patient. The improved VSM model was commissioned for 6 and 15 MV beams. The results of MC simulation are in very good agreement with measurements. Less than 2% of local difference between the MC simulation and the diamond detector measurement of the output factors in water was achieved. The X, Y and Z profiles measured in water with an ion chamber (V = 0.125 cm(3)) and a diamond detector were used to validate the models. An overall agreement of 2%/2 mm for high dose regions and 3%/2 mm in low dose regions between measurement and MC simulation for field sizes from 0.8 x 0.8 cm(2) to 16 x 21 cm(2) was achieved. An IMRT plan film verification was performed for two cases: 6 MV head&neck and 15 MV prostate. The simulation is in agreement with film measurements within 2%/2 mm in the high dose regions (> or = 0.1 Gy = 5% D(max)) and 5%/2 mm in low dose regions (<0.1 Gy).

A finite size pencil beam algorithm for IMRT dose optimization: density corrections.

For beamlet-based IMRT optimization, fast and less accurate dose computation algorithms are frequently used, while more accurate algorithms are needed to recompute the final dose for verification. In order to speed up the optimization process and ensure close proximity between dose in optimization and verification, proper consideration of dose gradients and tissue inhomogeneity effects should be ensured at every stage of the optimization. Due to their speed, pencil beam algorithms are often used for precalculation of beamlet dose distributions in IMRT treatment planning systems. However, accounting for tissue heterogeneities with these models requires the use of approximate rescaling methods. Recently, a finite size pencil beam (fsPB) algorithm, based on a simple and small set of data, was proposed which was specifically designed for the purpose of dose pre-computation in beamlet-based IMRT. The present work describes the incorporation of 3D density corrections, based on Monte Carlo simulations in heterogeneous phantoms, into this method improving the algorithm accuracy in inhomogeneous geometries while keeping its original speed and simplicity of commissioning. The algorithm affords the full accuracy of 3D density corrections at every stage of the optimization, hence providing the means for density related fluence modulation like penumbra shaping at field edges.