PURPOSE: Ultrahypofractionated radiation therapy (UHRT) is an effective treatment for localized prostate cancer with an acceptable toxicity profile; boosting the visible intraprostatic tumor has been shown to improve biochemical disease-free survival with no significant effect on genitourinary (GU) and gastrointestinal (GI) toxicity. METHODS AND MATERIALS: HERMES is a single-center noncomparative randomized phase 2 trial in men with intermediate or lower high risk prostate cancer. Patients were allocated (1:1) to 36.25 Gy in 5 fractions over 2 weeks or 24 Gy in 2 fractions over 8 days with an integrated boost to the magnetic resonance imaging (MRI) visible tumor of 27 Gy in 2 fractions. A minimization algorithm with a random element with risk group as a balancing factor was used for participant randomization. Treatment was delivered on the Unity MR-Linac (Elekta AB) with daily online adaption. The primary endpoint was acute GU Common Terminology Criteria for Adverse Events version 5.0 toxicity with the aim of excluding a doubling of the rate of acute grade 2+ GU toxicity seen in PACE. Analysis was by treatment received and included all participants who received at least 1 fraction of study treatment. This interim analysis was prespecified (stage 1 of a 2-stage Simon design) for when 10 participants in each treatment group had completed the acute toxicity monitoring period (12 weeks after radiation therapy). RESULTS: Acute grade 2 GU toxicity was reported in 1 (10%) patient in the 5-fraction group and 2 (20%) patients in the 2-fraction group. No grade 3+ GU toxicities were reported. CONCLUSIONS: At this interim analysis, the rate of GU toxicity in the 2-fraction and 5-fraction treatment groups was found to be below the prespecified threshold (5/10 grade 2+) and continuation of the study to complete recruitment of 23 participants per group was recommended..

BACKGROUND: Real-time dose estimation is a key-prerequisite to enable online intra-fraction treatment adaptation in magnetic resonance (MR)-guided radiotherapy (MRgRT). It is an essential component for the assessment of the dosimetric benefits and risks of online adaptive treatments, such as multi-leaf collimator (MLC)-tracking. PURPOSE: We present a proof-of-concept for a software workflow for real-time dose estimation of MR-guided adaptive radiotherapy based on real-time data-streams of the linac delivery parameters and target positions. METHODS: A software workflow, combining our in-house motion management software DynaTrack, a real-time dose calculation engine that connects to a research version of the treatment planning software (TPS) Monaco (v.6.09.00, Elekta AB, Stockholm, Sweden) was developed and evaluated. MR-guided treatment delivery on the Elekta Unity MR-linac was simulated with and without MLC-tracking for three prostate patients, previously treated on the Elekta Unity MR-linac (36.25 Gy/five fractions). Three motion scenarios were used: no motion, regular motion, and erratic prostate motion. Accumulated monitor units (MUs), centre of mass target position and MLC-leaf positions, were forwarded from DynaTrack at a rate of 25 Hz to a Monte Carlo (MC) based dose calculation engine which utilises the research GPUMCD-library (Elekta AB, Stockholm, Sweden). A rigid isocentre shift derived from the selected motion scenarios was applied to a bulk density-assigned session MR-image. The respective electron density used for treatment planning was accessed through the research Monaco TPS. The software workflow including the online dose reconstruction was validated against offline dose reconstructions. Our investigation showed that MC-based real-time dose calculations that account for all linac states (including MUs, MLC positions and target position) were infeasible, hence states were randomly sampled and used for calculation as follows; Once a new linac state was received, a dose calculation with 106 photons was started. Linac states that arrived during the time of the ongoing calculation were put into a queue. After completion of the ongoing calculation, one new linac state was randomly picked from the queue and assigned the MU accumulated from the previous state until the last sample in the queue. The queue was emptied, and the process repeated throughout treatment simulation. RESULTS: On average 27% (23%-30%) of received samples were used in the real-time calculation, corresponding to a calculation time for one linac state of 148 ms. Median gamma pass rate (2%/3 mm local) was 100.0% (99.9%-100%) within the PTV volume and 99.1% (90.1%-99.4.0%) with a 15% dose cut off. Differences in PTVDmean , CTVDmean , RectumD2% , and BladderD2% (offline-online, % of prescribed dose) were below 0.64%. Beam-by-beam comparisons showed deviations below 0.07 Gy. Repeated simulations resulted in standard deviations below 0.31% and 0.12 Gy for the investigated volume and dose criteria respectively. CONCLUSIONS: Real-time dose estimation was successfully performed using the developed software workflow for different prostate motion traces with and without MLC-tracking. Negligible dosimetric differences were seen when comparing online and offline reconstructed dose, enabling online intra-fraction treatment decisions based on estimates of the delivered dose..

Objective.Respiratory motion of lung tumours and adjacent structures is challenging for radiotherapy. Online MR-imaging cannot currently provide real-time volumetric information of the moving patient anatomy, therefore limiting precise dose delivery, delivered dose reconstruction, and downstream adaptation methods.Approach.We tailor a respiratory motion modelling framework towards an MR-Linac workflow to estimate the time-resolved 4D motion from real-time data. We develop a multi-slice acquisition scheme which acquires thick, overlapping 2D motion-slices in different locations and orientations, interleaved with 2D surrogate-slices from a fixed location. The framework fits a motion model directly to the input data without the need for sorting or binning to account for inter- and intra-cycle variation of the breathing motion. The framework alternates between model fitting and motion-compensated super-resolution image reconstruction to recover a high-quality motion-free image and a motion model. The fitted model can then estimate the 4D motion from 2D surrogate-slices. The framework is applied to four simulated anthropomorphic datasets and evaluated against known ground truth anatomy and motion. Clinical applicability is demonstrated by applying our framework to eight datasets acquired on an MR-Linac from four lung cancer patients.Main results.The framework accurately reconstructs high-quality motion-compensated 3D images with 2 mm3isotropic voxels. For the simulated case with the largest target motion, the motion model achieved a mean deformation field error of 1.13 mm. For the patient cases residual error registrations estimate the model error to be 1.07 mm (1.64 mm), 0.91 mm (1.32 mm), and 0.88 mm (1.33 mm) in superior-inferior, anterior-posterior, and left-right directions respectively for the building (application) data.Significance.The motion modelling framework estimates the patient motion with high accuracy and accurately reconstructs the anatomy. The image acquisition scheme can be flexibly integrated into an MR-Linac workflow whilst maintaining the capability of online motion-management strategies based on cine imaging such as target tracking and/or gating..

Abstract Objective. Intensity modulated proton therapy (IMPT) is susceptible to uncertainties in patient setup and proton range. Robust optimization is employed in IMPT treatment planning to ensure sufficient coverage of the clinical target volume (CTV) in predefined scenarios, albeit at a price of increased planning times. We investigated a deep learning (DL) strategy for dose predictions in individual error scenarios in head and neck cancer IMPT treatment planning, enabling direct evaluation of plan robustness. The model is able to differentiate between scenarios by using embeddings of the scenario index. Approach. To accommodate resolution disparities in planning CT-scans and accommodate the setup error scenarios, we introduced scenario-specific isocentric distance maps as inputs to the DL models. For 392 H&N cancer patients, high-quality 9-scenario ground truth (GT) robust plans were generated with wish-list driven fully automated multi-criteria optimization. The scenario index is converted to one-hot-vector that is used to derive the scenarios embeddings through the training of the DL model, aiding the model to predict a scenario specific dose distribution. Main results. The model achieved within 1%-point of agreement with the GT the predicted V 95 % of the voxelwise minimum dose for CTV Low and CTV High for 96% and 75% respectively of the test patients. Considering all robustness scenarios, median differences were 0.035%-point for CTV High V 95 % , 0.11%-point for CTV Low V 95 % , 0.29 GyE for parotids D mean , 0.7 GyE for submandibular glands D mean and 0.9 GyE for oral cavity D mean . Prediction of full 3D dose distributions for all scenarios took around 14 s. Significance. Predicting individual scenarios for robust proton therapy using DL dose prediction is feasible, enabling direct robustness evaluation of the predicted scenario doses..

BACKGROUND: MR-Linac allows for daily online treatment adaptation to the observed geometry of tumor targets and organs at risk (OARs). Manual delineation for head and neck cancer (HNC) patients takes 45-75 minutes, making it unsuitable for online adaptive radiotherapy. This study aims to clinically and dosimetrically validate an in-house developed algorithm which automatically delineates the elective target volume and OARs for HNC patients in under a minute. METHODS: Auto-contours were generated by an in-house model with 2D U-Net architecture trained and tested on 52 MRI scans via leave-one-out cross-validation. A randomized selection of 684 automated and manual contours (split half-and-half) was presented to an oncologist to perform a blind test and determine the clinical acceptability. The dosimetric impact was investigated for 13 patients evaluating the differences in dosage for all structures. RESULTS: Automated contours were generated in 8 seconds per MRI scan. The blind test concluded that 114 (33%) of auto-contours required adjustments with 85 only minor and 15 (4.4%) of manual contours required adjustments with 12 only minor. Dosimetric analysis showed negligible dosimetric differences between clinically acceptable structures and structures requiring minor changes. The Dice Similarity coefficients for the auto-contours ranged from 0.66 ± 0.11 to 0.88 ± 0.06 across all structures. CONCLUSION: Majority of auto-contours were clinically acceptable and could be used without any adjustments. Majority of structures requiring minor adjustments did not lead to significant dosimetric differences, hence manual adjustments were needed only for structures requiring major changes, which takes no longer than 10 minutes per patient..

AIMS: Neoadjuvant chemoradiotherapy followed by surgery is the mainstay of treatment for patients with rectal cancer. Standard clinical target volume (CTV) to planning target volume (PTV) margins of 10 mm are used to accommodate inter- and intrafraction motion of target. Treating on magnetic resonance-integrated linear accelerators (MR-linacs) allows for online manual recontouring and adaptation (MRgART) enabling the reduction of PTV margins. The aim of this study was to investigate motion of the primary CTV (CTVA; gross tumour volume and macroscopic nodes with 10 mm expansion to cover microscopic disease) in order to develop a simultaneous integrated boost protocol for use on MR-linacs. MATERIALS AND METHODS: Patients suitable for neoadjuvant chemoradiotherapy were recruited for treatment on MR-linac using a two-phase technique; only the five phase 1 fractions on MR-linac were used for analysis. Intrafraction motion of CTVA was measured between pre-treatment and post-treatment MRI scans. In MRgART, isotropically expanded pre-treatment PTV margins from 1 to 10 mm were rigidly propagated to post-treatment MRI to determine overlap with 95% of CTVA. The PTV margin was considered acceptable if overlap was >95% in 90% of fractions. To understand the benefit of MRgART, the same methodology was repeated using a reference computed tomography planning scan for pre-treatment imaging. RESULTS: In total, nine patients were recruited between January 2018 and December 2020 with T3a-T4, N0-N2, M0 disease. Forty-five fractions were analysed in total. The median motion across all planes was 0 mm, demonstrating minimal intrafraction motion. A PTV margin of 3 and 5mm was found to be acceptable in 96 and 98% of fractions, respectively. When comparing to the computed tomography reference scan, the analysis found that PTV margins to 5 and 10 mm only acceptably covered 51 and 76% of fractions, respectively. CONCLUSION: PTV margins can be reduced to 3-5 mm in MRgART for rectal cancer treatment on MR-linac within an simultaneous integrated boost protocol..

PURPOSE: A scoping literature review was conducted to identify gastrointestinal (GI) factors most likely to influence prostate motion during radiotherapy. We proffer that patient specific measurement of these GI factors could predict motion uncertainty during radiotherapy, facilitating personalised care by optimising treatment technique e.g., daily adaption or via bespoke patient pre-habilitation and preparation. METHODS: The scoping review was undertaken as per JBI guidelines. Searches were conducted across four databases: Ovid Medline®, EMBASE, CINAHL and EBSCO discovery. Articles written in English from 2010-present were included. Those pertaining to paediatrics, biological women exclusively, infectious and post-treatment GI morbidity and diet were excluded.Common GI factors impacting men were identified and related symptoms, incidence and measurement tools examined. Prevalence among persons with prostate cancer was explored and suitable assessment tools discussed. RESULTS: A preliminary search identified four prominent GI-factors: mental health, co-morbidity and medication, physical activity, and pelvic floor disorder. The scoping search found 3644 articles; 1646 were removed as duplicates. A further 1249 were excluded after title and abstract screening, 162 remained subsequent to full text review: 42 mental health, 53 co-morbidity and medication, 39 physical activity and 28 pelvic floor disorder.Six GI factors prevalent in the prostate cancer population and estimated most likely to influence prostate motion were identified: depression, anxiety, diabetes, obesity, low physical activity, and pelvic floor disorder. Reliable, quick, and easy to use tools are available to quantify these factors. CONCLUSION: A comprehensive GI factor assessment package suitable to implement into the radiotherapy clinic has been created. Unveiling these GI factors upfront will guide improved personalisation of radiotherapy..

Hybrid systems that combine Magnetic Resonance Imaging (MRI) and linear accelerators are available clinically to guide and adapt radiotherapy. Vendor-approved MRI sequences are provided, however alternative sequences may offer advantages. The aim of this study was to develop a systematic approach for non-vendor sequence evaluation, to determine safety, accuracy and overall clinical application of two potential sequences for bladder cancer MRI guided radiotherapy. Non-vendor sequences underwent and passed clinical image qualitative review, phantom quality assurance, and radiotherapy planning assessments. Volunteer workflow tests showed the potential for one sequence to reduce workflow time by 27% compared to the standard vendor sequence..

BACKGROUND AND PURPOSE: Physiological motion impacts the dose delivered to tumours and vital organs in external beam radiotherapy and particularly in particle therapy. The excellent soft-tissue demarcation of 4D magnetic resonance imaging (4D-MRI) could inform on intra-fractional motion, but long image reconstruction times hinder its use in online treatment adaptation. Here we employ techniques from high-performance computing to reduce 4D-MRI reconstruction times below two minutes to facilitate their use in MR-guided radiotherapy. MATERIAL AND METHODS: Four patients with pancreatic adenocarcinoma were scanned with a radial stack-of-stars gradient echo sequence on a 1.5T MR-Linac. Fast parallelised open-source implementations of the extra-dimensional golden-angle radial sparse parallel algorithm were developed for central processing unit (CPU) and graphics processing unit (GPU) architectures. We assessed the impact of architecture, oversampling and respiratory binning strategy on 4D-MRI reconstruction time and compared images using the structural similarity (SSIM) index against a MATLAB reference implementation. Scaling and bottlenecks for the different architectures were studied using multi-GPU systems. RESULTS: All reconstructed 4D-MRI were identical to the reference implementation (SSIM > 0.99). Images reconstructed with overlapping respiratory bins were sharper at the cost of longer reconstruction times. The CPU  + GPU implementation was over 17 times faster than the reference implementation, reconstructing images in 60 ± 1 s and hyper-scaled using multiple GPUs. CONCLUSION: Respiratory-resolved 4D-MRI reconstruction times can be reduced using high-performance computing methods for online workflows in MR-guided radiotherapy with potential applications in particle therapy..

BACKGROUND: T2 * mapping can characterize tumor hypoxia, which may be associated with resistance to therapy. Acquiring T2 * maps during MR-guided radiotherapy could inform treatment adaptation by, for example, escalating the dose to resistant sub-volumes. PURPOSE: The purpose of this work is to demonstrate the feasibility of the accelerated T2 * mapping technique using model-based image reconstruction with integrated trajectory auto-correction (TrACR) for MR-guided radiotherapy on an MR-Linear accelerator (MR-Linac). MATERIALS AND METHODS: The proposed method was validated in a numerical phantom, where two T2 * mapping approaches (sequential and joint) were compared for different noise levels (0,0.1,0.5,1) and gradient delays ([1, -1] and [1, -2] in units of dwell time for x- and y-axis, respectively). Fully sampled k-space was retrospectively undersampled using two different undersampling patterns. Root mean square errors (RMSEs) were calculated between reconstructed T2 * maps and ground truth. In vivo data was acquired twice weekly in one prostate and one head and neck cancer patient undergoing treatment on a 1.5 T MR-Linac. Data were retrospectively undersampled and T2 * maps reconstructed, with and without trajectory corrections were compared. RESULTS: Numerical simulations demonstrated that, for all noise levels, T2 * maps reconstructed with a joint approach demonstrated less error compared to an uncorrected and sequential approach. For a noise level of 0.1, uniform undersampling and gradient delay [1, -1] (in units of dwell time for x- and y-axis, respectively), RMSEs for sequential and joint approaches were 13.01 and 9.32 ms, respectively, which reduced to 10.92 and 5.89 ms for a gradient delay of [1, 2]. Similarly, for alternate undersampling and gradient delay [1, -1], RMSEs for sequential and joint approaches were 9.80 and 8.90 ms, respectively, which reduced to 9.10 and 5.40 ms for gradient delay [1, 2]. For in vivo data, T2 * maps reconstructed with our proposed approach resulted in less artifacts and improved visual appearance compared to the uncorrected approach. For both prostate and head and neck cancer patients, T2 * maps reconstructed from different treatment fractions showed changes within the planning target volume (PTV). CONCLUSION: Using the proposed approach, a retrospective data-driven gradient delay correction can be performed, which is particularly relevant for hybrid devices, where full information on the machine configuration is not available for image reconstruction. T2 * maps were acquired in under 5 min and can be integrated into MR-guided radiotherapy treatment workflows, which minimizes patient burden and leaves time for additional imaging for online adaptive radiotherapy on an MR-Linac..

INTRODUCTION: Distended rectums on pre-radiotherapy scans are historically associated with poorer outcomes in patients treated with two-dimensional IGRT. Subsequently, strict rectal tolerances and preparation regimes were implemented. Contemporary IGRT, daily online registration to the prostate, corrects interfraction motion but intrafraction motion remains. We re-examine the need for rectal management strategies when using contemporary IGRT by quantifying rectal volume and its effect on intrafraction motion. MATERIALS AND METHODS: Pre and during radiotherapy rectal volumes and intrafraction motion were retrospectively calculated for 20 patients treated in 5-fractions and 20 treated in 20-fractions. Small (rectal volume at planning-CT ≤ median), and large (volume > median) subgroups were formed, and rectal volume between timepoints and subgroups compared. Rectal volume and intrafraction motion correlation was examined using Spearman's rho. Intrafraction motion difference between small and large subgroups and between fractions with rectal volume < or ≥ 90 cm3 were assessed. RESULTS: Median rectal volume was 74 cm3, 64 cm3 and 65 cm3 on diagnostic-MRI, planning-CT and treatment imaging respectively (ns). No significant correlation was found between patient's rectal volume at planning-CT and median intrafraction motion, nor treatment rectal volume and intrafraction motion for individual fractions. No significant difference in intrafraction motion between small and large subgroups presented and for fractions where rectal volume breached 90 cm3, motion during that fraction was not significantly greater. CONCLUSION: Larger rectal volumes before radiotherapy and during treatment did not cause greater intrafraction motion. Findings support the relaxation of strict rectal diameter tolerances and do not support the need for rectal preparation when delivering contemporary IGRT to the prostate..

INTRODUCTION: The Elekta Unity MR-Linac (MRL) has enabled adaptive radiotherapy (ART) for patients with head and neck cancers (HNC). Adapt-To-Shape-Lite (ATS-Lite) is a novel Adapt-to-Shape strategy that provides ART without requiring daily clinician presence to perform online target and organ at risk (OAR) delineation. In this study we compared the performance of our clinically-delivered ATS-Lite strategy against three Adapt-To-Position (ATP) variants: Adapt Segments (ATP-AS), Optimise Weights (ATP-OW), and Optimise Shapes (ATP-OS). METHODS: Two patients with HNC received radical-dose radiotherapy on the MRL. For each fraction, an ATS-Lite plan was generated online and delivered and additional plans were generated offline for each ATP variant. To assess the clinical acceptability of a plan for every fraction, twenty clinical goals for targets and OARs were assessed for all four plans. RESULTS: 53 fractions were analysed. ATS-Lite passed 99.9% of mandatory dose constraints. ATP-AS and ATP-OW each failed 7.6% of mandatory dose constraints. The Planning Target Volumes for 54 Gy (D95% and D98%) were the most frequently failing dose constraint targets for ATP. ATS-Lite median fraction times for Patient 1 and 2 were 40 mins 9 s (range 28 mins 16 s - 47 mins 20 s) and 32 mins 14 s (range 25 mins 33 s - 44 mins 27 s), respectively. CONCLUSIONS: Our early data show that the novel ATS-Lite strategy produced plans that fulfilled 99.9% of clinical dose constraints in a time frame that is tolerable for patients and comparable to ATP workflows. Therefore, ATS-Lite, which bridges the gap between ATP and full ATS, will be further utilised and developed within our institute and it is a workflow that should be considered for treating patients with HNC on the MRL..

AIMS: With interest in normal tissue sparing and dose-escalated radiotherapy in the treatment of inoperable locally advanced non-small cell lung cancer, this study investigated the impact of motion-managed moderate deep inspiration breath hold (mDIBH) on normal tissue sparing and dose-escalation potential and compared this to planning with a four-dimensional motion-encompassing internal target volume or motion-compensating mid-ventilation approach. MATERIALS AND METHODS: Twenty-one patients underwent four-dimensional and mDIBH planning computed tomography scans. Internal and mid-ventilation target volumes were generated on the four-dimensional scan, with mDIBH target volumes generated on the mDIBH scan. Isotoxic target dose-escalation guidelines were used to generate six plans per patient: three with a target dose cap and three without. Target dose-escalation potential, normal tissue complication probability and differences in pre-specified dose-volume metrics were evaluated for the three motion-management techniques. RESULTS: The mean total lung volume was significantly greater with mDIBH compared with four-dimensional scans. Lung dose (mean and V21 Gy) and mean heart dose were significantly reduced with mDIBH in comparison with four-dimensional-based approaches, and this translated to a significant reduction in heart and lung normal tissue complication probability with mDIBH. In 20/21 patients, the trial target prescription dose cap of 79.2 Gy was achievable with all motion-management techniques. CONCLUSION: mDIBH aids lung and heart dose sparing in isotoxic dose-escalated radiotherapy compared with four-dimensional planning techniques. Given concerns about lung and cardiac toxicity, particularly in an era of consolidation immunotherapy, reduced normal tissue doses may be advantageous for treatment tolerance and outcome..

Bladder tumour-focused magnetic resonance image-guided adaptive radiotherapy using a 1.5 Tesla MR-linac is feasible. A full online workflow adapting to anatomy at each fraction is achievable in approximately 30 min. Intra-fraction bladder filling did not compromise target coverage with the class solution employed..

AIMS: Prostate morphological changes during external beam radiotherapy are poorly understood. Excellent soft-tissue visualisation offered by magnetic resonance image-guided radiotherapy (MRIgRT) provides an opportunity to better understand such changes. The aim of this study was to quantify prostate volume and dimension changes occurring during extreme and moderately hypofractionated schedules. MATERIALS AND METHODS: Forty prostate cancer patients treated on the Unity 1.5 Tesla magnetic resonance linear accelerator (MRL) were retrospectively reviewed. The cohort comprised patients treated with 36.25 Gy in five fractions (n = 20) and 60 Gy in 20 fractions (n = 20). The volume of the delineated prostates on reference planning computed tomography (fused with MRI) and daily T2-weighted 2-min session images acquired on Unity were charted. Forty planning computed tomography and 500 MRL prostate volumes were evaluated. The mean absolute and relative change in prostate volume during radiotherapy was compared using a paired t-test (P value <0.01 considered significant to control for multiple comparisons). The maximum dimension of the delineated prostate was measured in three isocentric planes. RESULTS: Significant prostate volume changes, relative to MRL imaging fraction 1 (MRL#1), were seen at all time points for the five-fraction group. The peak mean relative volume increase was 21% (P < 0.001), occurring at MRL#3 and MRL#4 after 14.5 and 21.75 Gy, respectively. Prostate expansion was greatest in the superior-inferior direction; the peak mean maximal extension was 5.9 mm. The maximal extension in the left-right and anterior-posterior directions measured 1.1 and 2.2 mm, respectively. For the 20-fraction group, prostate volume increased relative to MRL#1, for all treatment time points. The mean relative volume increase was 11% (P < 0.001) at MRL#5 after 12 Gy, it then fluctuated between 8 and 13%. From MRL#5 to MRL#20, the volume increase was significant (P < 0.01) for 12 of 16 time points calculated. The peak mean maximal extension in the superior-inferior direction was 3.1 mm. The maximal extension in the left-right and anterior-posterior directions measured 1.7 and 3.7 mm, respectively. CONCLUSION: Significant prostate volume and dimension changes occur during extreme and moderately hypofractionated radiotherapy. The extent of change was greater during extreme hypofractionation. MRIgRT offers the opportunity to reveal, quantify and correct for this deformation..

PURPOSE: This study aims to develop and validate a simple geometric model of the accelerator head, from which a particle phase space can be calculated for application to fast Monte Carlo dose calculation in real-time adaptive photon radiotherapy. With this objective in view, the study investigates whether the phase space model can facilitate dose calculations which are compatible with those of a commercial treatment planning system, for convenient interoperability. MATERIALS AND METHODS: A dual-source model of the head of a Versa HD accelerator (Elekta AB, Stockholm, Sweden) was created. The model used parameters chosen to be compatible with those of 6-MV flattened and 6-MV flattening filter-free photon beams in the RayStation treatment planning system (RaySearch Laboratories, Stockholm, Sweden). The phase space model was used to calculate a photon phase space for several treatment plans, and the resulting phase space was applied to the Dose Planning Method (DPM) Monte Carlo dose calculation algorithm. Simple fields and intensity-modulated radiation therapy (IMRT) treatment plans for prostate and lung were calculated for benchmarking purposes and compared with the convolution-superposition dose calculation within RayStation. RESULTS: For simple square fields in a water phantom, the calculated dose distribution agrees to within ±2% with that from the commercial treatment planning system, except in the buildup region, where the DPM code does not model the electron contamination. For IMRT plans of prostate and lung, agreements of ±2% and ±6%, respectively, are found, with slightly larger differences in the high dose gradients. CONCLUSIONS: The phase space model presented allows convenient calculation of a phase space for application to Monte Carlo dose calculation, with straightforward translation of beam parameters from the RayStation beam model. This provides a basis on which to develop dose calculation in a real-time adaptive setting..

BACKGROUND: The prostate demonstrates inter- and intra- fractional changes and thus adaptive radiotherapy would be required to ensure optimal coverage. Daily adaptive radiotherapy for MRI-guided radiotherapy can be both time and resource intensive when structure delineation is completed manually. Contours can be auto-generated on the MR-Linac via a deformable image registration (DIR) based mapping process from the reference image. This study evaluates the performance of automatically generated target structure contours against manually delineated contours by radiation oncologists for prostate radiotherapy on the Elekta Unity MR-Linac. METHODS: Plans were generated from prostate contours propagated by DIR and rigid image registration (RIR) for forty fractions from ten patients. A two-dose level SIB (simultaneous integrated boost) IMRT plan is used to treat localised prostate cancer; 6000 cGy to the prostate and 4860 cGy to the seminal vesicles. The dose coverage of the PTV 6000 and PTV 4860 created from the manually drawn target structures was evaluated with each plan. If the dose objectives were met, the plan was considered successful in covering the gold standard (clinician-delineated) volume. RESULTS: The mandatory PTV 6000 dose objective (D98% > 5580 cGy) was met in 81 % of DIR plans and 45 % of RIR plans. The SV were mapped by DIR only and for all the plans, the PTV 4860 dose objective met the optimal target (D98% > 4617 cGy). The plans created by RIR led to under-coverage of the clinician-delineated prostate, predominantly at the apex or the bladder-prostate interface. CONCLUSION: Plans created from DIR propagation of prostate contours outperform those created from RIR propagation. In approximately 1 in 5 DIR plans, dosimetric coverage of the gold standard PTV was not clinically acceptable. Thus, at our institution, we use a combination of DIR propagation of contours alongside manual editing of contours where deemed necessary for online treatments..

A shift of the daily plan can mitigate target position changes that occur between daily MR acquisition and treatment for MR-linac radiotherapy, but increases the session time. We demonstrated that our workflow strategy and decision-making process, to determine whether a subsequent shift is necessary, is appropriate..

Radiation therapy is a major component of cancer treatment pathways worldwide. The main aim of this treatment is to achieve tumor control through the delivery of ionizing radiation while preserving healthy tissues for minimal radiation toxicity. Because radiation therapy relies on accurate localization of the target and surrounding tissues, imaging plays a crucial role throughout the treatment chain. In the treatment planning phase, radiological images are essential for defining target volumes and organs-at-risk, as well as providing elemental composition (e.g., electron density) information for radiation dose calculations. At treatment, onboard imaging informs patient setup and could be used to guide radiation dose placement for sites affected by motion. Imaging is also an important tool for treatment response assessment and treatment plan adaptation. MRI, with its excellent soft tissue contrast and capacity to probe functional tissue properties, holds great untapped potential for transforming treatment paradigms in radiation therapy. The MR in Radiation Therapy ISMRM Study Group was established to provide a forum within the MR community to discuss the unmet needs and fuel opportunities for further advancement of MRI for radiation therapy applications. During the summer of 2021, the study group organized its first virtual workshop, attended by a diverse international group of clinicians, scientists, and clinical physicists, to explore our predictions for the future of MRI in radiation therapy for the next 25 years. This article reviews the main findings from the event and considers the opportunities and challenges of reaching our vision for the future in this expanding field..

BACKGROUND: Size-based assessments are inaccurate indicators of tumor response in soft-tissue sarcoma (STS), motivating the requirement for new response imaging biomarkers for this rare and heterogeneous disease. In this study, we assess the test-retest repeatability of radiomic features from MR diffusion-weighted imaging (DWI) and derived maps of apparent diffusion coefficient (ADC) in retroperitoneal STS and compare baseline repeatability with changes in radiomic features following radiotherapy (RT). MATERIALS AND METHODS: Thirty patients with retroperitoneal STS received an MR examination prior to treatment, of whom 23/30 were investigated in our repeatability analysis having received repeat baseline examinations and 14/30 patients were investigated in our post-treatment analysis having received an MR examination after completing pre-operative RT. One hundred and seven radiomic features were extracted from the full manually delineated tumor region using PyRadiomics. Test-retest repeatability was assessed using an intraclass correlation coefficient (baseline ICC), and post-radiotherapy variance analysis (post-RT-IMS) was used to compare the change in radiomic feature value to baseline repeatability. RESULTS: For the ADC maps and DWI images, 101 and 102 features demonstrated good baseline repeatability (baseline ICC > 0.85), respectively. Forty-three and 2 features demonstrated both good baseline repeatability and a high post-RT-IMS (>0.85), respectively. Pearson correlation between the baseline ICC and post-RT-IMS was weak (0.432 and 0.133, respectively). CONCLUSIONS: The ADC-based radiomic analysis shows better test-retest repeatability compared with features derived from DWI images in STS, and some of these features are sensitive to post-treatment change. However, good repeatability at baseline does not imply sensitivity to post-treatment change..

The drive towards hypofractionated prostate radiotherapy is motivated by a low alpha/beta ratio for prostate cancer (1 to 3 Gy) compared to surrounding organs at risk, implying an improved therapeutic ratio with increasing dose per fraction. Early evidence from studies of ultrahypofractionated (UHF) prostate HDR brachytherapy has shown good tolerability in terms of normal tissue toxicities and clinical outcomes similar to conventional fractionation schedules. MR-guided stereotactic body radiotherapy (SBRT) with online plan adaptation and real-time tumour imaging may enable UHF doses to be delivered to the prostate safely, without the invasiveness of brachytherapy. The feasibility of UHF prostate treatment planning for the Unity MR-Linac (MRL, Elekta AB, Stockholm) was investigated for target prescriptions and planning constraints derived from the HDR brachytherapy and SBRT literature. Monaco 5.40 (Elekta) was used to generate MRL step-and-shoot IMRT plans for three dose fractionation protocols (5, 2 and 1 fractions), for ten randomly selected previously treated prostate cancer patients. Of the ten plans per UHF scheme, all clinical goals were met in all cases for 5 fractions, and in six cases for both 2 and 1 fraction schemes. PTV D95% was compromised by up to 6.4% and 3.9% of the associated target dose for 2 and 1 fraction plans respectively. There were two cases of PTV D95% compromise greater than a 5% dose decrease for the 2 fraction plans. The study suggests feasibility of the UHF treatment planning approaches if combined with real-time motion mitigation strategies..

PURPOSE: Online adaptive radiotherapy would greatly benefit from the development of reliable auto-segmentation algorithms for organs-at-risk and radiation targets. Current practice of manual segmentation is subjective and time-consuming. While deep learning-based algorithms offer ample opportunities to solve this problem, they typically require large datasets. However, medical imaging data are generally sparse, in particular annotated MR images for radiotherapy. In this study, we developed a method to exploit the wealth of publicly available, annotated CT images to generate synthetic MR images, which could then be used to train a convolutional neural network (CNN) to segment the parotid glands on MR images of head and neck cancer patients. METHODS: Imaging data comprised 202 annotated CT and 27 annotated MR images. The unpaired CT and MR images were fed into a 2D CycleGAN network to generate synthetic MR images from the CT images. Annotations of axial slices of the synthetic images were generated by propagating the CT contours. These were then used to train a 2D CNN. We assessed the segmentation accuracy using the real MR images as test dataset. The accuracy was quantified with the 3D Dice similarity coefficient (DSC), Hausdorff distance (HD), and mean surface distance (MSD) between manual and auto-generated contours. We benchmarked the approach by a comparison to the interobserver variation determined for the real MR images, as well as to the accuracy when training the 2D CNN to segment the CT images. RESULTS: The determined accuracy (DSC: 0.77±0.07, HD: 18.04±12.59mm, MSD: 2.51±1.47mm) was close to the interobserver variation (DSC: 0.84±0.06, HD: 10.85±5.74mm, MSD: 1.50±0.77mm), as well as to the accuracy when training the 2D CNN to segment the CT images (DSC: 0.81±0.07, HD: 13.00±7.61mm, MSD: 1.87±0.84mm). CONCLUSIONS: The introduced cross-modality learning technique can be of great value for segmentation problems with sparse training data. We anticipate using this method with any nonannotated MRI dataset to generate annotated synthetic MR images of the same type via image style transfer from annotated CT images. Furthermore, as this technique allows for fast adaptation of annotated datasets from one imaging modality to another, it could prove useful for translating between large varieties of MRI contrasts due to differences in imaging protocols within and between institutions..

AIMS: Target delineation uncertainty is arguably the largest source of geometric uncertainty in radiotherapy. Several factors can affect it, including the imaging modality used for delineation. It is accounted for by applying safety margins to the target to produce a planning target volume (PTV), to which treatments are designed. To determine the margin, the delineation uncertainty is measured as the delineation error, and then a margin recipe used. However, there is no published evidence of such analysis for recurrent gynaecological cancers (RGC). The aims of this study were first to quantify the delineation uncertainty for RGC gross tumour volumes (GTVs) and to calculate the associated PTV margins and then to quantify the difference in GTV, delineation uncertainty and PTV margin, between a computed tomography-magnetic resonance imaging (CT-MRI) and MRI workflow. MATERIALS AND METHODS: Seven clinicians delineated the GTV for 20 RGC tumours on co-registered CT and MRI datasets (CT-MRI) and on MRI alone. The delineation error, the standard deviation of distances from each clinician's outline to a reference, was measured and the required PTV margin determined. Differences between using CT-MRI and MRI alone were assessed. RESULTS: The overall delineation error and the resulting margin were 3.1 mm and 8.5 mm, respectively, for CT-MRI, reducing to 2.5 mm and 7.1 mm, respectively, for MRI alone. Delineation errors and therefore the theoretical margins, varied widely between patients. MRI tumour volumes were on average 15% smaller than CT-MRI tumour volumes. DISCUSSION: This study is the first to quantify delineation error for RGC tumours and to calculate the corresponding PTV margin. The determined margins were larger than those reported in the literature for similar patients, bringing into question both current margins and margin calculation methods. The wide variation in delineation error between these patients suggests that applying a single population-based margin may result in PTVs that are suboptimal for many. Finally, the reduced tumour volumes and safety margins suggest that patients with RGC may benefit from an MRI-only treatment workflow..

BACKGROUND AND PURPOSE: 4D and midposition MRI could inform plan adaptation in lung and abdominal MR-guided radiotherapy. We present deep learning-based solutions to overcome long 4D-MRI reconstruction times while maintaining high image quality and short scan times. METHODS: Two 3D U-net deep convolutional neural networks were trained to accelerate the 4D joint MoCo-HDTV reconstruction. For the first network, gridded and joint MoCo-HDTV-reconstructed 4D-MRI were used as input and target data, respectively, whereas the second network was trained to directly calculate the midposition image. For both networks, input and target data had dimensions of 256 × 256 voxels (2D) and 16 respiratory phases. Deep learning-based MRI were verified against joint MoCo-HDTV-reconstructed MRI using the structural similarity index (SSIM) and the naturalness image quality evaluator (NIQE). Moreover, two experienced observers contoured the gross tumour volume and scored the images in a blinded study. RESULTS: For 12 subjects, previously unseen by the networks, high-quality 4D and midposition MRI (1.25 × 1.25 × 3.3 mm3) were each reconstructed from gridded images in only 28 seconds per subject. Excellent agreement was found between deep-learning-based and joint MoCo-HDTV-reconstructed MRI (average SSIM ≥ 0.96, NIQE scores 7.94 and 5.66). Deep-learning-based 4D-MRI were clinically acceptable for target and organ-at-risk delineation. Tumour positions agreed within 0.7 mm on midposition images. CONCLUSION: Our results suggest that the joint MoCo-HDTV and midposition algorithms can each be approximated by a deep convolutional neural network. This rapid reconstruction of 4D and midposition MRI facilitates online treatment adaptation in thoracic or abdominal MR-guided radiotherapy..

Magnetic resonance imaging (MRI)-guided radiotherapy (RT) (MRIgRT) falls outside the scope of existing high energy photon therapy dosimetry protocols, because those protocols do not consider the effects of the magnetic field on detector response and on absorbed dose to water. The aim of this study is to evaluate and demonstrate the traceable measurement of absorbed dose in MRIgRT systems using alanine, made possible by the characterisation of alanine sensitivity to magnetic fields reported previously by Billaset al(2020Phys. Med. Biol.65115001), in a way which is compatible with existing standards and calibrations available for conventional RT. In this study, alanine is used to transfer absorbed dose to water to MRIgRT systems from a conventional linac. This offers an alternative route for the traceable measurement of absorbed dose to water, one which is independent of the transfer using ionisation chambers. The alanine dosimetry is analysed in combination with measurements with several Farmer-type chambers, PTW 30013 and IBA FC65-G, at six different centres and two different MRIgRT systems (Elekta Unity™ and ViewRay MRIdian™). The results are analysed in terms of the magnetic field correction factors, and in terms of the absorbed dose calibration coefficients for the chambers, determined at each centre. This approach to reference dosimetry in MRIgRT produces good consistency in the results, across the centres visited, at the level of 0.4% (standard deviation). Farmer-type ionisation chamber magnetic field correction factors were determined directly, by comparing calibrations in some MRIgRT systems with and without the magnetic field ramped up, and indirectly, by comparing calibrations in all the MRIgRT systems with calibrations in a conventional linac. Calibration coefficients in the MRIgRT systems were obtained with a standard uncertainty of 1.1% (Elekta Unity™) and 0.9% (ViewRay MRIdian™), for three different chamber orientations with respect to the magnetic field. The values obtained for the magnetic field correction factor in this investigation are consistent with those presented in the summary by de Pooteret al(2021Phys. Med. Biol.6605TR02), and would tend to support the adoption of a magnetic field correction factor which depends on the chamber type, PTW 30013 or IBA FC65-G..

PURPOSE: High-field magnetic resonance-linear accelerators (MR-Linacs), linear accelerators combined with a diagnostic magnetic resonance imaging (MRI) scanner and online adaptive workflow, potentially give rise to novel online anatomic and response adaptive radiation therapy paradigms. The first high-field (1.5T) MR-Linac received regulatory approval in late 2018, and little is known about clinical use, patient tolerability of daily high-field MRI, and toxicity of treatments. Herein we report the initial experience within the MOMENTUM Study (NCT04075305), a prospective international registry of the MR-Linac Consortium. METHODS AND MATERIALS: Patients were included between February 2019 and October 2020 at 7 institutions in 4 countries. We used descriptive statistics to describe the patterns of care, tolerability (the percentage of patients discontinuing their course early), and safety (grade 3-5 Common Terminology Criteria for Adverse Events v.5 acute toxicity within 3 months after the end of treatment). RESULTS: A total 943 patients participated in the MOMENTUM Study, 702 of whom had complete baseline data at the time of this analysis. Patients were primarily male (79%) with a median age of 68 years (range, 22-93) and were treated for 39 different indications. The most frequent indications were prostate (40%), oligometastatic lymph node (17%), brain (12%), and rectal (10%) cancers. The median number of fractions was 5 (range, 1-35). Six patients discontinued MR-Linac treatments, but none due to an inability to tolerate repeated high-field MRI. Of the 415 patients with complete data on acute toxicity at 3-month follow-up, 18 (4%) patients experienced grade 3 acute toxicity related to radiation. No grade 4 or 5 acute toxicity related to radiation was observed. CONCLUSIONS: In the first 21 months of our study, patterns of care were diverse with respect to clinical utilization, body sites, and radiation prescriptions. No patient discontinued treatment due to inability to tolerate daily high-field MRI scans, and the acute radiation toxicity experience was encouraging..

Target volume delineation uncertainty (DU) is arguably one of the largest geometric uncertainties in radiotherapy that are accounted for using planning target volume (PTV) margins. Geometrical uncertainties are typically derived from a limited sample of patients. Consequently, the resultant margins are not tailored to individual patients. Furthermore, standard PTVs cannot account for arbitrary anisotropic extensions of the target volume originating from DU. We address these limitations by developing a method to measure DU for each patient by a single clinician. This information is then used to produce PTVs that account for each patient's unique DU, including any required anisotropic component. We do so using a two-step uncertainty evaluation strategy that does not rely on multiple samples of data to capture the DU of a patient's gross tumour volume (GTV) or clinical target volume. For simplicity, we will just refer to the GTV in the following. First, the clinician delineates two contour sets; one which bounds all voxels believed to have a probability of belonging to the GTV of 1, while the second includes all voxels with a probability greater than 0. Next, one specifies a probability density function for the true GTV boundary position within the boundaries of the two contours. Finally, a patient-specific PTV, designed to account for all systematic errors, is created using this information along with measurements of the other systematic errors. Clinical examples indicate that our margin strategy can produce significantly smaller PTVs than the van Herk margin recipe. Our new radiotherapy target delineation concept allows DUs to be quantified by the clinician for each patient, leading to PTV margins that are tailored to each unique patient, thus paving the way to a greater personalisation of radiotherapy..

Microbeam radiotherapy (MRT) is a preclinical method of delivering spatially-fractionated radiotherapy aiming to improve the therapeutic window between normal tissue complication and tumour control. Previously, MRT was limited to ultra-high dose rate synchrotron facilities. The aim of this study was to investigate in vitro effects of MRT on tumour and normal cells at conventional dose rates produced by a bench-top X-ray source. Two normal and two tumour cell lines were exposed to homogeneous broad beam (BB) radiation, MRT, or were separately irradiated with peak or valley doses before being mixed. Clonogenic survival was assessed and compared to BB-estimated surviving fractions calculated by the linear-quadratic (LQ)-model. All cell lines showed similar BB sensitivity. BB LQ-model predictions exceeded the survival of cell lines following MRT or mixed beam irradiation. This effect was stronger in tumour compared to normal cell lines. Dose mixing experiments could reproduce MRT survival. We observed a differential response of tumour and normal cells to spatially fractionated irradiations in vitro, indicating increased tumour cell sensitivity. Importantly, this was observed at dose rates precluding the presence of FLASH effects. The LQ-model did not predict cell survival when the cell population received split irradiation doses, indicating that factors other than local dose influenced survival after irradiation..

PURPOSE: Flow-compensated (FC) diffusion-weighted MRI (DWI) for intravoxel-incoherent motion (IVIM) modeling allows for a more detailed description of tissue microvasculature than conventional IVIM. The long acquisition time of current FC-IVIM protocols, however, has prohibited clinical application. Therefore, we developed an optimized abdominal FC-IVIM acquisition with a clinically feasible scan time. METHODS: Precision and accuracy of the FC-IVIM parameters were assessed by fitting the FC-IVIM model to signal decay curves, simulated for different acquisition schemes. Diffusion-weighted acquisitions were added subsequently to the protocol, where we chose the combination of b-value, diffusion time and gradient profile (FC or bipolar) that resulted in the largest improvement to its accuracy and precision. The resulting two optimized FC-IVIM protocols with 25 and 50 acquisitions (FC-IVIMopt25 and FC-IVIMopt50 ), together with a complementary acquisition consisting of 50 diffusion-weighting (FC-IVIMcomp ), were acquired in repeated abdominal free-breathing FC-IVIM imaging of seven healthy volunteers. Intersession and intrasession within-subject coefficient of variation of the FC-IVIM parameters were compared for the liver, spleen, and kidneys. RESULTS: Simulations showed that the performance of FC-IVIM improved in tissue with larger perfusion fraction and signal-to-noise ratio. The scan time of the FC-IVIMopt25 and FC-IVIMopt50 protocols were 8 and 16 min. The best in vivo performance was seen in FC-IVIMopt50 . The intersession within-subject coefficients of variation of FC-IVIMopt50 were 11.6%, 16.3%, 65.5%, and 36.0% for FC-IVIM model parameters diffusivity, perfusion fraction, characteristic time and blood flow velocity, respectively. CONCLUSIONS: We have optimized the FC-IVIM protocol, allowing for clinically feasible scan times (8-16 min)..

For multimodality therapies such as the combination of hyperthermia and radiation, quantification of biological effects is key for dose prescription and response prediction. Tumour spheroids have a microenvironment that more closely resembles that of tumours in vivo and may thus be a superior in vitro cancer model than monolayer cultures. Here, the response of tumour spheroids formed from two established human cancer cell lines (HCT116 and CAL27) to single and combination treatments of radiation (0-20 Gy), and hyperthermia at 47 °C (0-780 CEM43) has been evaluated. Response was analysed in terms of spheroid growth, cell viability and the distribution of live/dead cells. Time-lapse imaging was used to evaluate mechanisms of cell death and cell detachment. It was found that sensitivity to heat in spheroids was significantly less than that seen in monolayer cultures. Spheroids showed different patterns of shrinkage and regrowth when exposed to heat or radiation: heated spheroids shed dead cells within four days of heating and displayed faster growth post-exposure than samples that received radiation or no treatment. Irradiated spheroids maintained a dense structure and exhibited a longer growth delay than spheroids receiving hyperthermia or combination treatment at (thermal) doses that yielded equivalent levels of clonogenic cell survival. We suggest that, unlike radiation, which kills dividing cells, hyperthermia-induced cell death affects cells independent of their proliferation status. This induces microenvironmental changes that promote spheroid growth. In conclusion, 3D tumour spheroid growth studies reveal differences in response to heat and/or radiation that were not apparent in 2D clonogenic assays but that may significantly influence treatment efficacy..

Reference dosimetry in the presence of a strong magnetic field is challenging. Ionisation chambers have shown to be strongly affected by magnetic fields. There is a need for robust and stable detectors in MRI-guided radiotherapy (MRIgRT). This study investigates the behaviour of the alanine dosimeter in magnetic fields and assesses its suitability to act as a reference detector in MRIgRT. Alanine pellets were loaded in a waterproof holder, placed in an electromagnet and irradiated by 60Co and 6 MV and 8 MV linac beams over a range of magnetic flux densities. Monte Carlo simulations were performed to calculate the absorbed dose, to water and to alanine, with and without magnetic fields. Combining measurements with simulations, the effect of magnetic fields on alanine response was quantified and a correction factor for the presence of magnetic fields on alanine was determined. This study finds that the response of alanine to ionising radiation is modified when the irradiation is in the presence of a magnetic field. The effect is energy independent and may increase the alanine/electron paramagnetic resonance (EPR) signal by 0.2% at 0.35 T and 0.7% at 1.5 T. In alanine dosimetry for MRIgRT, this effect, if left uncorrected, would lead to an overestimate of dose. Accordingly, a correction factor, [Formula: see text], is defined. Values are obtained for this correction as a function of magnetic flux density, with a standard uncertainty which depends on the magnetic field and is 0.6% or less. The strong magnetic field has a measurable effect on alanine dosimetry. For alanine which is used to measure absorbed dose to water in a strong magnetic field, but which has been calibrated in the absence of a magnetic field, a small correction to the reported dose is required. With the inclusion of this correction, alanine/EPR is a suitable reference dosimeter for measurements in MRIgRT..

BACKGROUND AND PURPOSE: Anatomical changes during external beam radiotherapy prevent the accurate delivery of the intended dose distribution. Resolving the delivered dose, which is currently unknown, is crucial to link radiotherapy doses to clinical outcomes and ultimately improve the standard of care. MATERIAL AND METHODS: In this study, we present a dose reconstruction workflow based on data routinely acquired during MR-guided radiotherapy. It employs 3D MR images, 2D cine MR images and treatment machine log files to calculate the delivered dose taking intrafractional motion into account. The developed pipeline was used to measure anatomical changes and assess their dosimetric impact in 89 prostate radiotherapy fractions delivered with a 1.5 T MR-linac at our institute. RESULTS: Over the course of radiation delivery, the CTV shifted 0.6 mm ± 2.1 mm posteriorly and 1.3 mm ± 1.5 mm inferiorly. When extrapolating the dose changes in each case to 20 fractions, the mean clinical target volume D98% and clinical target volume D50% dose-volume metrics decreased by 1.1 Gy ± 1.6 Gy and 0.1 Gy ± 0.2 Gy, respectively. Bladder D3% did not change (0.0 Gy ± 1.2 Gy), while rectum D3% decreased by 1.0 Gy ± 2.0 Gy. Although anatomical changes and their dosimetric impact were small in the majority of cases, large intrafractional motion caused the delivered dose to substantially deviate from the intended plan in some fractions. CONCLUSIONS: The presented end-to-end workflow is able to reliably, non-invasively and automatically reconstruct the delivered prostate radiotherapy dose by processing MR-linac treatment log files and online MR images. In the future, we envision this workflow to be adapted to other cancer sites and ultimately to enter widespread clinical use..

PURPOSE: Several studies have demonstrated potential improvements in treatment time through the use of dynamic arcs for delivery of stereotactic body radiation therapy (SBRT) on Cyberknife. However, the delivery system has a finite accuracy, so that potential exists for dosimetric uncertainties. This study estimates the expected dosimetric accuracy of dynamic delivery of SBRT, based on realistic estimates of the uncertainties in delivery parameters. METHODS: Five SBRT patient cases (prostate A - conventional, prostate B - brachytherapy-type, lung, liver, partial left breast) were retrospectively studied. Treatment plans were produced for a fixed arc trajectory using fluence optimization, segmentation, and direct aperture optimization. Dose rate uncertainty was modeled as a smoothly varying random fluctuation of ± 1.0%, ±2.0% or ± 5.0% over a time period of 10, 30 or 60 s. Multileaf collimator uncertainty was modeled as a lag in position of each leaf up to 0.25 or 0.5 mm. Robot pointing error was modeled as a shift of the target location, with the direction of the shift chosen as a random angle with respect to the multileaf collimator and with a random magnitude in the range 0.0-1.0 mm at the delivery nodes and with an additional random magnitude of 0.5-1.0 mm in between the delivery nodes. The impact of the errors was investigated using dose-volume histograms. RESULTS: Uncertainty in dose rate has the effect of varying the total monitor units delivered, which in turn produces a variation in mean dose to the planning target volume. The random sampling of dose rate error produces a distribution of mean doses with a standard deviation proportional to the magnitude of the dose rate uncertainty. A lag in multileaf collimator position of 0.25 or 0.5 mm produces a small impact on the delivered dose. In general, an increase in the PTV mean dose of around 1% is observed. An error in robot pointing of the order of 1 mm produces a small increase in dose inhomogeneity to the planning target volume, sometimes accompanied by an increase in mean dose by around 1%. CONCLUSIONS: Based upon the limited data available on the dose rate stability and geometric accuracy of the Cyberknife system, this study estimates that dynamic arc delivery can be accomplished with sufficient accuracy for clinical application. Dose rate variation produces a change in dose to the planning target volume according to the perturbation of total monitor units delivered, while multileaf collimator lag and robot pointing error typically increase the mean dose to the planning target volume by up to 1%..

Previous work has shown that PRESAGE® can be used successfully to perform 3D dosimetric measurements of complex radiotherapy treatments. However, measurements near the sample edges are known to be difficult to achieve. This is an issue when the doses at air-material interfaces are of interest, for example when investigating the electron return effect (ERE) present in treatments delivered by magnetic resonance (MR)-linac systems. To study this effect, a set of 3.5 cm-diameter cylindrical PRESAGE® samples was uniformly irradiated with multiple dose fractions, using either a conventional linac or an MR-linac. The samples were imaged between fractions using an optical-CT, to read out the corresponding accumulated doses. A calibration between TPS-predicted dose and optical-CT pixel value was determined for individual dosimeters as a function of radial distance from the axis of rotation. This data was used to develop a correction that was applied to four additional samples of PRESAGE® of the same formulation, irradiated with 3D-CRT and IMRT treatment plans, to recover significantly improved 3D measurements of dose. An alternative strategy was also tested, in which the outer surface of the sample was physically removed prior to irradiation. Results show that for the formulation studied here, PRESAGE® samples have a central region that responds uniformly and an edge region of 6-7 mm where there is gradual increase in dosimeter response, rising to an over-response of 24%-36% at the outer boundary. This non-uniform dose response increases in both extent and magnitude over time. Both mitigation strategies investigated were successful. In our four exemplar studies, we show how discrepancies at edges are reduced from 13%-37% of the maximum dose to between 2 and 8%. Quantitative analysis shows that the 3D gamma passing rates rise from 90.4, 69.3, 63.7 and 43.6% to 97.3, 99.9, 96.7 and 98.9% respectively..

INTRODUCTION: MR-guided adapted radiotherapy (MRgART) using a high field MR-linac has recently become available. We report the estimated delivered fractional dose of the first five prostate cancer patients treated at our centre using MRgART and compare this to C-Arm linac daily Image Guided Radiotherapy (IGRT). METHODS: Patients were treated using adapted treatment plans shaped to their daily anatomy. The treatments were recalculated on an MR image acquired immediately prior to treatment delivery in order to estimate the delivered fractional dose. C-arm linac non-adapted VMAT treatment plans were recalculated on the same MR images to estimate the fractional dose that would have been delivered using conventional radiotherapy techniques using a daily IGRT protocol. RESULTS: 95% and 93% of mandatory target coverage objectives and organ at risk dose constraints were achieved by MRgART and C-arm linac delivered dose estimates, respectively. Both delivery techniques were estimated to have achieved 98% of mandatory Organ At Risk (OAR) dose constraints whereas for the target clinical goals, 86% and 80% were achieved by MRgART and C-arm linac delivered dose estimates. CONCLUSIONS: Prostate MRgART can be delivered using the a high field MR-linac. Radiotherapy performed on a C-arm linac offers a good solution for prostate cancer patients who present with favourable anatomy at the time of reference imaging and demonstrate stable anatomy throughout the course of their treatment. For patients with critical OARs abutting target volumes on their reference image we have demonstrated the potential for a target dose coverage improvement for MRgART compared to C-arm linac treatment..

An MR-Linac can provide motion information of tumour and organs-at-risk before, during, and after beam delivery. However, MR imaging cannot provide real-time high-quality volumetric images which capture breath-to-breath variability of respiratory motion. Surrogate-driven motion models relate the motion of the internal anatomy to surrogate signals, thus can estimate the 3D internal motion from these signals. Internal surrogate signals based on patient anatomy can be extracted from 2D cine-MR images, which can be acquired on an MR-Linac during treatment, to build and drive motion models. In this paper we investigate different MRI-derived surrogate signals, including signals generated by applying principal component analysis to the image intensities, or control point displacements derived from deformable registration of the 2D cine-MR images. We assessed the suitability of the signals to build models that can estimate the motion of the internal anatomy, including sliding motion and breath-to-breath variability. We quantitatively evaluated the models by estimating the 2D motion in sagittal and coronal slices of 8 lung cancer patients, and comparing them to motion measurements obtained from image registration. For sagittal slices, using the first and second principal components on the control point displacements as surrogate signals resulted in the highest model accuracy, with a mean error over patients around 0.80 mm which was lower than the in-plane resolution. For coronal slices, all investigated signals except the skin signal produced mean errors over patients around 1 mm. These results demonstrate that surrogate signals derived from 2D cine-MR images, including those generated by applying principal component analysis to the image intensities or control point displacements, can accurately model the motion of the internal anatomy within a single sagittal or coronal slice. This implies the signals should also be suitable for modelling the 3D respiratory motion of the internal anatomy..

BACKGROUND AND PURPOSE: Retrieving quantitative parameters from magnetic resonance imaging (MRI), e.g. for early assessment of radiotherapy treatment response, necessitates contouring regions of interest, which is time-consuming and prone to errors. This becomes more pressing for daily imaging on MRI-guided radiotherapy systems. Therefore, we trained a deep convolutional neural network to automatically contour involved lymph nodes on diffusion-weighted (DW) MRI of head and neck cancer (HNC) patients receiving radiotherapy. MATERIALS AND METHODS: DW-images from 48 HNC patients (18 induction-chemotherapy + chemoradiotherapy; 30 definitive chemoradiotherapy) with 68 involved lymph nodes were obtained on a diagnostic 1.5 T MR-scanner prior to and 2-3 timepoints throughout treatment. A radiation oncologist delineated the lymph nodes on the b = 50 s/mm2 images. A 3D U-net was trained to contour involved lymph nodes. Its performance was evaluated in all 48 patients using 8-fold cross-validation and calculating the Dice similarity coefficient (DSC) and the absolute difference in median apparent diffusion coefficient (ΔADC) between the manual and generated contours. Additionally, the performance was evaluated in an independent dataset of three patients obtained on a 1.5 T MR-Linac. RESULTS: In the definitive chemoradiotherapy patients (n = 96 patients/lymphnodes/timepoints) the DSC was 0.87 (0.81-0.91) [median (1st-3rd quantiles)] and ΔADC was 1.9% (0.8-3.4%) and both remained stable throughout treatment. The network performed worse in the patients receiving induction-chemotherapy (n = 65), with DSC = 0.80 (0.71-0.87) and ΔADC = 3.3% (1.6-8.0%). The network performed well on the MR-Linac data (n = 8) with DSC = 0.80 (0.75-0.82) and ΔADC = 4.0% (0.6-9.1%). CONCLUSIONS: We established accurate automatic contouring of involved lymph nodes for HNC patients on diagnostic and MR-Linac DW-images..

.

Whole bladder magnetic resonance image-guided radiotherapy using the 1.5 Telsa MR-linac is feasible. Full online adaptive planning workflow based on the anatomy seen at each fraction was performed. This was delivered within 45 min. Intra-fraction bladder filling did not compromise target coverage. Patients reported acceptable tolerance of treatment..

PURPOSE: The POP-ART RT study aims to determine to what extent and how intrafractional real-time respiratory motion management (RRMM), and plan adaptation for interfractional anatomical changes (ART) are used in clinical practice and to understand barriers to implementation. Here we report on part II: ART using more than one plan per target per treatment course. MATERIALS AND METHODS: A questionnaire on the current practice of ART, wishes for expansion or implementation, and barriers to implementation was distributed worldwide. Four types of ART were discriminated: daily online replanning, online plan library, protocolled offline replanning (all three based on a protocol), and ad-hoc offline replanning. RESULTS: The questionnaire was completed by 177 centres from 40 countries. ART was used by 61% of respondents (31% with protocol) for a median (range) of 3 (1-8) tumour sites. CBCT/MVCT was the main imaging modality except for online daily replanning (11 users) where 10 users used MR. Two thirds of respondents wished to implement ART for a new tumour site; 40% of these had plans to do it in the next 2 years. Human/material resources and technical limitations were the main barriers to further use and implementation. CONCLUSIONS: ART was used for a broad range of tumour sites, mainly with ad-hoc offline replanning and for a median of 3 tumour sites. There was a large interest in implementing ART for more tumour sites, mainly limited by human/material resources and technical limitations. Daily online replanning was primarily performed on MR-linacs..

PURPOSE: The POP-ART RT study aims to determine to what extent and how intra-fractional real-time respiratory motion management (RRMM) and plan adaptation for inter-fractional anatomical changes (ART), are used in clinical practice and to understand barriers to implementation. Here we report on part I: RRMM. MATERIAL AND METHODS: A questionnaire was distributed worldwide to assess current clinical practice, wishes for expansion or new implementation and barriers to implementation. RRMM was defined as inspiration/expiration gating in free-breathing or breath-hold, or tracking where the target and the beam are continuously realigned. RESULTS: The questionnaire was completed by 200 centres from 41 countries. RRMM was used by 68% of respondents ('users') for a median (range) of 2 (1-6) tumour sites. Eighty-one percent of users applied inspiration breath-hold in at least one tumour site (breast: 96%). External marker was used to guide RRMM by 61% of users. KV/MV imaging was frequently used for liver and pancreas (with fiducials) and for lung (with or without fiducials). Tracking was mainly performed on robotic linacs with hybrid internal-external monitoring. For breast and lung, approximately 75% of respondents used or wished to implement RRMM, which was lower for liver (44%) and pancreas (27%). Seventy-one percent of respondents wished to implement RRMM for a new tumour site. Main barriers were human/financial resources and capacity on the machine. CONCLUSION: Sixty-eight percent of respondents used RRMM and 71% wished to implement RRMM for a new tumour site. The main barriers to implementation were human/financial resources and capacity on treatment machines..

Breathing motion is challenging for radiotherapy planning and delivery. This requires advanced four-dimensional (4D) imaging and motion mitigation strategies and associated validation tools with known deformations. Numerical phantoms such as the XCAT provide reproducible and realistic data for simulation-based validation. However, the XCAT generates partially inconsistent and non-invertible deformations where tumours remain rigid and structures can move through each other. We address these limitations by post-processing the XCAT deformation vector fields (DVF) to generate a breathing phantom with realistic motion and quantifiable deformation. An open-source post-processing framework was developed that corrects and inverts the XCAT-DVFs while preserving sliding motion between organs. Those post-processed DVFs are used to warp the first XCAT-generated image to consecutive time points providing a 4D phantom with a tumour that moves consistently with the anatomy, the ability to scale lung density as well as consistent and invertible DVFs. For a regularly breathing case, the inverse consistency of the DVFs was verified and the tumour motion was compared to the original XCAT. The generated phantom and DVFs were used to validate a motion-including dose reconstruction (MIDR) method using isocenter shifts to emulate rigid motion. Differences between the reconstructed doses with and without lung density scaling were evaluated. The post-processing framework produced DVFs with a maximum [Formula: see text]-percentile inverse-consistency error of 0.02 mm. The generated phantom preserved the dominant sliding motion between the chest wall and inner organs. The tumour of the original XCAT phantom preserved its trajectory while deforming consistently with the underlying tissue. The MIDR was compared to the ground truth dose reconstruction illustrating its limitations. MIDR with and without lung density scaling resulted in small dose differences up to 1 Gy (prescription 54 Gy). The proposed open-source post-processing framework overcomes important limitations of the original XCAT phantom and makes it applicable to a wider range of validation applications within radiotherapy..

Purpose: MR-guided Radiation Therapy (MRgRT) allows for high-precision radiotherapy under real-time MR visualization. This enables margin reduction and subsequent dose escalation which may lead to higher tumor control and less toxicity. The Unity MR-linac (Elekta AB, Stockholm, Sweden) integrates a linear accelerator with a 1.5T diagnostic quality MRI and an online adaptive workflow. A prospective international registry was established to facilitate the evidence-based implementation of the Unity MR-linac into clinical practice, to systemically evaluate long-term outcomes, and to aid further technical development of MR-linac-based MRgRT. Methods and Results: In February 2019, the Multi-OutcoMe EvaluatioN of radiation Therapy Using the MR-linac study (MOMENTUM) started within the MR-linac Consortium. The MOMENTUM study is an international academic-industrial partnership between several hospitals and industry partner Elekta. All patients treated on the MR-linac are eligible for inclusion in MOMENTUM. For participants, we collect clinical patient data (e.g., patient, tumor, and treatment characteristics) and technical patient data which is defined as information generated on the MR-linac during treatment. The data are captured, pseudonymized, and stored in an international registry at set time intervals up to two years after treatment. Patients can choose to provide patient-reported outcomes and consent to additional MRI scans acquired on the MR-linac. This registry will serve as a data platform that supports multicenter research investigating the MR-linac. Rules and regulations on data sharing, data access, and intellectual property rights are summarized in an academic-industrial collaboration agreement. Data access rules ensure secure data handling and research integrity for investigators and institutions. Separate data access rules exist for academic and industry partners. This study is registered at ClinicalTrials.gov with ID: NCT04075305 (https://clinicaltrials.gov/ct2/show/NCT04075305). Conclusion: The multi-institutional MOMENTUM study has been set up to collect clinical and technical patient data to advance technical development, and facilitate evidenced-based implementation of MR-linac technology with the ultimate purpose to improve tumor control, survival, and quality of life of patients with cancer..

Geometric uncertainties in radiotherapy are conventionally addressed by defining a safety margin around the radiotherapy target. Misappropriation of such margins could result in disease recurrence from geometric miss or unnecessary irradiation of normal tissue. Numerous quantitative organ motion studies in adults have been published, but the first paediatric-specific studies were only published in recent years. In the very near future, intensity-modulated proton beam therapy and magnetic resonance-guided radiotherapy will be clinically implemented in the UK. Such techniques offer the ability to deliver radiotherapy to the pinnacle of precision and accuracy, if geometric uncertainty relating to internal organ motion and deformation can be optimally managed. The optimal margin to account for internal organ motion in children remains largely undefined. Continuing efforts to characterise motion in children and young people is necessary to optimally define safety margins and to realise the full potential of intensity-modulated radiotherapy, magnetic resonance-guided radiotherapy and intensity-modulated proton beam therapy. This overview offers a timely review of published reports on paediatric organ motion, in anticipation of the increasing application of advanced radiotherapy techniques in paediatric radiotherapy..

The Cyberknife system (Accuray Inc., Sunnyvale, CA) enables radiotherapy using stereotactic ablative body radiotherapy (SABR) with a large number of non-coplanar beam orientations. Recently, a multileaf collimator has also been available to allow flexibility in field shaping. This work aims to evaluate the quality of treatment plans obtainable with the multileaf collimator. Specifically, the aim is to find a subset of beam orientations from a predetermined set of candidate directions, such that the treatment quality is maintained but the treatment time is reduced. An evolutionary algorithm is used to successively refine a randomly selected starting set of beam orientations. By using an efficient computational framework, clinically useful solutions can be found in several hours. It is found that 15 beam orientations are able to provide treatment quality which approaches that of the candidate beam set of 110 beam orientations, but with approximately half of the estimated treatment time. Choice of an efficient subset of beam orientations offers the possibility to improve the patient experience and maximise the number of patients treated..

OBJECTIVES: To determine the potential for dose escalation to a biological equivalent dose BED10 ≅ 100 Gy in hypofractionated radiotherapy for locally advanced pancreatic cancer (LAPC). MATERIALS AND METHODS: Ten unselected LAPC patients were retrospectively included in the study. Two fractionation regimens were compared (5 and 15 fractions). The aim was to cover 95% of the Planning Target Volume (PTV) with a BED10 = 54 Gy (base dose = 33 Gy in 5 fractions, 42.5 Gy in 15 fractions) whilst respecting organs-at-risk (OAR) constraints. Once the highest PTV coverage was achieved dose escalation to a BED10 ≅ 100 Gy (escalated dose = 50 Gy in 5 fractions, 67.5 Gy in 15 fractions) was attempted, limiting the PTV maximum dose to 130% of the escalated dose. RESULTS: In 5 fractions, 95% PTV coverage by both base and escalated doses could be achieved for one patient with PTV more than 1 cm away from OAR. 95% and 90% PTV coverage by the base dose was achieved in one and two patients respectively. In all other patients, coverage even by the base dose had to be compromised to comply with OAR constraints. In 15 fractions, 95% PTV coverage by the base dose was feasible for all patients except one. Dose escalation allowed improvement in target coverage by the base dose in both fractionation regimen whilst covering a sub-volume of the PTV with a BED10 ≅ 100 Gy. Both fractionation schemes were equivalent in terms of dose escalation potential. CONCLUSION: LAPC patients with OAR close to the PTV are generally not eligible for hypofractionation with dose escalation. However, this planning study shows that it is possible to cover PTV sub-volumes with a BED10 ≅ 100 Gy in addition to delivering a BED10 = 54 Gy to 90-95% of the PTV as commonly prescribed to this population. Combined with an adaptive approach, this may maximize PTV coverage by a high BED on days with favourable anatomy..

CT-based radiotherapy workflow is limited by poor soft tissue definition in the pelvis and reliance on rigid registration methods. Current image-guided radiotherapy and adaptive radiotherapy models therefore have limited ability to improve clinical outcomes. The advent of MRI-guided radiotherapy solutions provides the opportunity to overcome these limitations with the potential to deliver online real-time MRI-based plan adaptation on a daily basis, a true "plan of the day." This review describes the application of MRI guided radiotherapy in two pelvic tumour sites likely to benefit from this approach..

Radiotherapy (RT) aims to deliver a spatially conformal dose of radiation to tumours while maximizing the dose sparing to healthy tissues. However, the internal patient anatomy is constantly moving due to respiratory, cardiac, gastrointestinal and urinary activity. The long term goal of the RT community to 'see what we treat, as we treat' and to act on this information instantaneously has resulted in rapid technological innovation. Specialized treatment machines, such as robotic or gimbal-steered linear accelerators (linac) with in-room imaging suites, have been developed specifically for real-time treatment adaptation. Additional equipment, such as stereoscopic kilovoltage (kV) imaging, ultrasound transducers and electromagnetic transponders, has been developed for intrafraction motion monitoring on conventional linacs. Magnetic resonance imaging (MRI) has been integrated with cobalt treatment units and more recently with linacs. In addition to hardware innovation, software development has played a substantial role in the development of motion monitoring methods based on respiratory motion surrogates and planar kV or Megavoltage (MV) imaging that is available on standard equipped linacs. In this paper, we review and compare the different intrafraction motion monitoring methods proposed in the literature and demonstrated in real-time on clinical data as well as their possible future developments. We then discuss general considerations on validation and quality assurance for clinical implementation. Besides photon RT, particle therapy is increasingly used to treat moving targets. However, transferring motion monitoring technologies from linacs to particle beam lines presents substantial challenges. Lessons learned from the implementation of real-time intrafraction monitoring for photon RT will be used as a basis to discuss the implementation of these methods for particle RT..

MOTIVATION: With interlaced microbeam radiation therapy (MRT) a first kilovoltage radiotherapy (RT) concept combining spatially fractionated entrance beams and homogeneous dose distribution at the target exists. However, this technique suffers from its high sensitivity to positioning errors of the target relative to the radiation source. With spiral microbeam radiation therapy (spiralMRT), this publication introduces a new irradiation geometry, offering similar spatial fractionation properties as interlaced MRT, while being less vulnerable to target positioning uncertainties. METHODS: The dose distributions achievable with spiralMRT in a simplified human head geometry were calculated with Monte Carlo simulations based on Geant4 and the dependence of the result on the microbeam pitch, total field size, and photon energy were analysed. A comparison with interlaced MRT and conventional megavoltage tomotherapy was carried out. RESULTS: SpiralMRT can deliver homogeneous dose distributions, while using spatially fractionated entrance beams. The valley dose of spiralMRT entrance beams is by up to 40% lower than the corresponding tomotherapy dose, thus indicating a better normal tissue sparing. The optimum photon energy is found to be around [Formula: see text]. CONCLUSIONS: SpiralMRT is a promising approach to delivering homogeneous dose distributions with spatially fractionated entrance beams, possibly decreasing normal tissue side effects in hypofractionated RT..

PURPOSE: The use of dynamic arcs for delivery of stereotactic body radiation therapy (SBRT) on Cyberknife is investigated, with a view to improving treatment times. This study investigates the required modeling of robot and multileaf collimator (MLC) motion between control points in the trajectory and then uses this to develop an optimization method for treatment planning of a dynamic arc with Cyberknife. The resulting plans are compared in terms of dose-volume histograms and estimated treatment times with those produced by a conventional beam arrangement. METHODS: Five SBRT patient cases (prostate A - conventional, prostate B - brachytherapy-type, lung, liver, and partial left breast) were retrospectively studied. A suitable arc trajectory with control points spaced at 5° was proposed and treatment plans were produced for typical clinical protocols. The optimization consisted of a fluence optimization, segmentation, and direct aperture optimization using a gradient descent method. Dose delivered by the moving MLC was either taken to be the dose delivered discretely at the control points or modeled using effective fluence delivered between control points. The accuracy of calculated dose was assessed by recalculating after optimization using five interpolated beams and 100 interpolated apertures between each optimization control point. The resulting plans were compared using dose-volume histograms and estimated treatment times with those for a conventional Cyberknife beam arrangement. RESULTS: If optimization is performed based on discrete doses delivered at the arc control points, large differences of up to 40% of the prescribed dose are seen when recalculating with interpolation. When the effective fluence between control points is taken into account during optimization, dosimetric differences are <2% for most structures when the plans are recalculated using intermediate nodes, but there are differences of up to 15% peripherally. Treatment plan quality is comparable between the arc trajectory and conventional body path. All plans meet the relevant clinical goals, with the exception of specific structures which overlap with the planning target volume. Median estimated treatment time is 355 s (range 235-672 s) for arc delivery and 675 s (range 554-1025 s) for conventional delivery. CONCLUSIONS: The method of using effective fluence to model MLC motion between control points is sufficiently accurate to provide for accurate inverse planning of dynamic arcs with Cyberknife. The proposed arcing method produces treatment plans with comparable quality to the body path, with reduced estimated treatment delivery time..

OBJECTIVES: To systematically identify the preferred magnetic resonance imaging (MRI) sequences following volunteer imaging on a 1.5 Tesla (T) MR-Linear Accelerator (MR Linac) for future protocol development. METHODS: Non-patient volunteers were recruited to a Research and Ethics committee approved prospective MR-only imaging study on a 1.5T MR Linac system. Volunteers attended 1-3 imaging sessions that included a combination of mDixon, T1w, T2w sequences using 2-dimensional (2D) and 3-dimensional (3D) acquisitions. Each sequence was acquired over 2-7 minutes and reviewed by a panel of 3 observers to evaluate image quality using a visual grading analysis based on a 4-point Likert scale. Sequences were acquired and modified iteratively until deemed fit for purpose (online image matching or re-planning) and all observers agreed they were suitable in 3 volunteers. RESULTS: 26 volunteers underwent 31 imaging sessions of six general anatomical regions. Images were acquired in one or two of six general anatomical regions: male pelvis (n = 9), female pelvis (n = 4), chestwall/breast (n = 5), lung/oesophagus (n = 5), abdomen (n = 3) and head and neck (n = 5). Images were acquired using a pre-defined exam-card that on average, included six sequences (range 2-10), with a maximum scan time of approximately one hour. The majority of observers preferred T2-weighted sequences. The thorax teams were the only groups to prefer T1-weighted imaging. CONCLUSIONS: An iterative process identified sequence agreement in all anatomical regions. These sequences will now be evaluated in patient volunteers. ADVANCES IN KNOWLEDGE: This manuscript is the first publication sharing the results of the first systematic selection of MRI sequences for use in on-board MRI-guided radiotherapy by end-users (therapeutic radiographers and clinical oncologists) in healthy volunteers..

With the advent of MRI-guided radiotherapy, the suitability of commercially available radiation dose detectors needs to be assessed. The aim of this study was to investigate the effect of the magnetic field (B-field) on the response of the Gafchromic EBT-3 films. Moreover, as an independent study, we contribute to clarifying the inconsistency of the results of recent published studies, on the effect of B-field on the sensitivity of Gafchromic films. A 60Co beam was used to irradiate film samples in an electromagnet. An in-house PMMA phantom was designed to fit in the 5 cm gap between the two poles of the magnet. The phantom consists of two symmetrical plates where a film can be inserted. The absorbed dose rate of the 60Co beam for zero B-field was measured using alanine pellets in a Farmer-type holder. A 12-point response curve was created, representing [Formula: see text] as a function of dose, for each of five different B-field strengths (0 T to 2 T). This study finds that there is at most a small effect of the magnetic field on the response of EBT-3 film. In terms of netOD (red channel) the change in response varied from ‒0.0011 at 0.5 T to 0.0045 at 2.0 T, with a standard uncertainty of 0.0030. If uncorrected, this would lead to an error in film-measured dose, for the red channel, of 2.4% at 2 T, with a standard uncertainty on dose of 1.4%. Results are also presented for B-field strengths of 0.5 T, 1 T and 1.5 T, which are all zero within the measurement uncertainty. Comparison between other studies is also presented. Considering the small change on dose determined with EBT-3 when irradiated under the presence of B-field and taking into account the overall uncertainty in dosimetry using EBT-3 film achieved in this work, EBT-3 is assessed to be a suitable detector for relative and absolute dosimetry, with appropriate corrections, in MRI-guided radiotherapy. The results of the current work also elucidate the inconsistency on the reports from previous studies and demonstrate the necessity of similar investigations by independent teams, especially if the existing results may be in conflict..

Despite the utility of tumour characterisation using quantitative parameter maps from multi-b-value diffusion-weighted MRI (DWI), clinicians often prefer the use of the image with highest diffusion-weighting (b-value), for instance for defining regions of interest (ROIs). However, these images are typically degraded by noise, as they do not utilize the information from the full acquisition. We present a principal component analysis (PCA) approach for model-free denoising of DWI data. PCA-denoising was compared to synthetic MRI, where a diffusion model is fitted for each voxel and a denoised image at a given b-value is generated from the model fit. A quantitative comparison of systematic and random errors was performed on data simulated using several diffusion models (mono-exponential, bi-exponential, stretched-exponential and kurtosis). A qualitative visual comparison was also performed for in vivo images in six healthy volunteers and three pancreatic cancer patients. In simulations, the reduction in random errors from PCA-denoising was substantial (up to 55%) and similar to synthetic MRI (up to 53%). Model-based synthetic MRI denoising resulted in substantial (up to 29% of signal) systematic errors, whereas PCA-denoising was able to denoise without introducing systematic errors (less than 2%). In vivo, the signal-to-noise ratio (SNR) and sharpness of PCA-denoised images were superior to synthetic MRI, resulting in clearer tumour boundaries. In the presence of motion, PCA-denoising did not cause image blurring, unlike image averaging or synthetic MRI. Multi-b-value MRI can be denoised model-free with our PCA-denoising strategy that reduces noise to a level similar to synthetic MRI, but without introducing systematic errors associated with the synthetic MRI method..

Significant improvements in radiotherapy are likely to come from biological rather than technical optimization, for example increasing tumour radiosensitivity via combination with targeted therapies. Such paradigms must first be evaluated in preclinical models for efficacy, and recent advances in small animal radiotherapy research platforms allow advanced irradiation protocols, similar to those used clinically, to be carried out in orthotopic models. Dose assessment in such systems is complex however, and a lack of established tools and methodologies for traceable and accurate dosimetry is currently limiting the capabilities of such platforms and slowing the clinical uptake of new approaches. Here we report the creation of an anatomically correct phantom, fabricated from materials with tissue-equivalent electron density, into which dosimetry detectors can be incorporated for measurement as part of quality control (QC). The phantom also allows training in preclinical radiotherapy planning and cross-institution validation of dose delivery protocols for small animal radiotherapy platforms without the need to sacrifice animals, with high reproducibility. Mouse CT data was acquired and segmented into soft tissue, bone and lung. The skeleton was fabricated using 3D printing, whilst lung was created using computer numerical control (CNC) milling. Skeleton and lung were then set into a surface-rendered mould and soft tissue material added to create a whole-body phantom. Materials for fabrication were characterized for atomic composition and attenuation for x-ray energies typically found in small animal irradiators. Finally cores were CNC milled to allow intracranial incorporation of bespoke detectors (alanine pellets) for dosimetry measurement..

MR-guided radiotherapy treatment planning utilises the high soft-tissue contrast of MRI to reduce uncertainty in delineation of the target and organs at risk. Replacing 4D-CT with MRI-derived synthetic 4D-CT would support treatment plan adaptation on hybrid MR-guided radiotherapy systems for inter- and intrafractional differences in anatomy and respiration, whilst mitigating the risk of CT to MRI registration errors. Three methods were devised to calculate synthetic 4D and midposition (time-weighted mean position of the respiratory cycle) CT from 4D-T1w and Dixon MRI. The first approach employed intensity-based segmentation of Dixon MRI for bulk-density assignment (sCTD). The second step added spine density information using an atlas of CT and Dixon MRI (sCTDS). The third iteration used a polynomial function relating Hounsfield units and normalised T1w image intensity to account for variable lung density (sCTDSL). Motion information in 4D-T1w MRI was applied to generate synthetic CT in midposition and in twenty respiratory phases. For six lung cancer patients, synthetic 4D-CT was validated against 4D-CT in midposition by comparison of Hounsfield units and dose-volume metrics. Dosimetric differences found by comparing sCTD,DS,DSL and CT were evaluated using a Wilcoxon signed-rank test (p   =  0.05). Compared to sCTD and sCTDS, planning on sCTDSL significantly reduced absolute dosimetric differences in the planning target volume metrics to less than 98 cGy (1.7% of the prescribed dose) on average. When comparing sCTDSL and CT, average radiodensity differences were within 97 Hounsfield units and dosimetric differences were significant only for the planning target volume D99% metric. All methods produced clinically acceptable results for the organs at risk in accordance with the UK SABR consensus guidelines and the LungTech EORTC phase II trial. The overall good agreement between sCTDSL and CT demonstrates the feasibility of employing synthetic 4D-CT for plan adaptation on hybrid MR-guided radiotherapy systems..

Thermo-radiosensitisation is a promising approach for treatment of radio-resistant tumours such as those containing hypoxic subregions. Response prediction and treatment planning should account for tumour response heterogeneity, e.g. due to microenvironmental factors, and quantification of the biological effects induced. 3D tumour spheroids provide a physiological in vitro model of tumour response and a systems oncology framework for simulating spheroid response to radiation and hyperthermia is presented. Using a cellular automaton model, 3D oxygen diffusion, delivery of radiation and/or hyperthermia were simulated for many ([Formula: see text]) individual cells forming a spheroid. The iterative oxygen diffusion model was compared to an analytical oxygenation model and simulations were calibrated and validated against experimental data for irradiated (0-10 Gy) and/or heated (0-240 CEM43) HCT116 spheroids. Despite comparable clonogenic survival, spheroid growth differed significantly following radiation or hyperthermia. This dynamic response was described well by the simulation ([Formula: see text] > 0.85). Heat-induced cell death was implemented as a fast, proliferation-independent process, allowing reoxygenation and repopulation, whereas radiation was modelled as proliferation-dependent mitotic catastrophe. This framework stands out both through its experimental validation and its novel ability to predict spheroid response to multimodality treatment. It provides a good description of response where biological dose-weighting based on clonogenic survival alone was insufficient..

PURPOSE: To determine the 3-dimensional (3D) intrafractional motion of head and neck squamous cell carcinoma (HNSCC). METHODS AND MATERIALS: Dynamic contrast-enhanced magnetic resonance images from 56 patients with HNSCC in the treatment position were analyzed. Dynamic contrast-enhanced magnetic resonance imaging consisted of 3D images acquired every 2.9 seconds for 4 minutes 50 seconds. Intrafractional tumor motion was studied in the 3 minutes 43 seconds of images obtained after initial contrast enhancement. To assess tumor motion, rigid registration (translations only) was performed using a region of interest (ROI) mask around the tumor. The results were compared with bulk body motion from registration to all voxels. Motion was split into systematic motion and random motion. Correlations between the tumor site and random motion were tested. The within-subject coefficient of variation was determined from 8 patients with repeated baseline measures. Random motion was also assessed at the end of the first week (38 patients) and second week (25 patients) of radiation therapy to investigate trends of motion. RESULTS: Tumors showed irregular occasional rapid motion (eg, swallowing or coughing), periodic intermediate motion (respiration), and slower systematic drifts throughout treatment. For 95% of the patients, displacements due to systematic and random motion were <1.4 mm and <2.1 mm, respectively, 95% of the time. The motion without an ROI mask was significantly (P<.0001, Wilcoxon signed rank test) less than the motion with an ROI mask, indicating that tumors can move independently from the bony anatomy. Tumor motion was significantly (P=.005, Mann-Whitney U test) larger in the hypopharynx and larynx than in the oropharynx. The within-subject coefficient of variation for random motion was 0.33. The average random tumor motion did not increase notably during the first 2 weeks of treatment. CONCLUSIONS: The 3D intrafractional tumor motion of HNSCC is small, with systematic motion <1.4 mm and random motion <2.1 mm 95% of the time..

Combined radiotherapy and hyperthermia offer great potential for the successful treatment of radio-resistant tumours through thermo-radiosensitization. Tumour response heterogeneity, due to intrinsic, or micro-environmentally induced factors, may greatly influence treatment outcome, but is difficult to account for using traditional treatment planning approaches. Systems oncology simulation, using mathematical models designed to predict tumour growth and treatment response, provides a powerful tool for analysis and optimization of combined treatments. We present a framework that simulates such combination treatments on a cellular level. This multiscale hybrid cellular automaton simulates large cell populations (up to 107 cells) in vitro, while allowing individual cell-cycle progression, and treatment response by modelling radiation-induced mitotic cell death, and immediate cell kill in response to heating. Based on a calibration using a number of experimental growth, cell cycle and survival datasets for HCT116 cells, model predictions agreed well (R2 > 0.95) with experimental data within the range of (thermal and radiation) doses tested (0-40 CEM43, 0-5 Gy). The proposed framework offers flexibility for modelling multimodality treatment combinations in different scenarios. It may therefore provide an important step towards the modelling of personalized therapies using a virtual patient tumour..

Radiation therapy to the prostate involves increasingly sophisticated delivery techniques and changing fractionation schedules. With a low estimated α/β ratio, a larger dose per fraction would be beneficial, with moderate fractionation schedules rapidly becoming a standard of care. The integration of a magnetic resonance imaging (MRI) scanner and linear accelerator allows for accurate soft tissue tracking with the capacity to replan for the anatomy of the day. Extreme hypofractionation schedules become a possibility using the potentially automated steps of autosegmentation, MRI-only workflow, and real-time adaptive planning. The present report reviews the steps involved in hypofractionated adaptive MRI-guided prostate radiation therapy and addresses the challenges for implementation..

PURPOSE: Tumor hypoxia correlates with treatment failure in patients undergoing conventional radiation therapy. However, no published studies have investigated tumor hypoxia in patients undergoing stereotactic body radiation therapy (SBRT). We aimed to noninvasively quantify the tumor hypoxic volume (HV) in non-small cell lung cancer (NSCLC) tumors to elucidate the potential role of tumor vascular response and reoxygenation at high single doses. METHODS AND MATERIALS: Six SBRT-eligible patients with NSCLC tumors >1 cm were prospectively enrolled in an institutional review board-approved study. Dynamic positron emission tomography images were acquired at 0 to 120 minutes, 150 to 180 minutes, and 210 to 240 minutes after injection of 18F-fluoromisonidazole. Serial imaging was performed prior to delivery of 18 Gy and at approximately 48 hours and approximately 96 hours after SBRT. Tumor HVs were quantified using the tumor-to-blood ratio (>1.2) and rate of tracer influx (>0.0015 mL·min·cm-3). RESULTS: An elevated and in some cases persistent level of tumor hypoxia was observed in 3 of 6 patients. Two patients exhibited no detectable baseline tumor hypoxia, and 1 patient with high baseline hypoxia only completed 1 imaging session. On the basis of the tumor-to-blood ratio, in the remaining 3 patients, tumor HVs increased on day 2 after 18 Gy and then showed variable responses on day 4. In the 3 of 6 patients with detectable hypoxia at baseline, baseline tumor HVs ranged between 17% and 24% (mean, 21%), and HVs on days 2 and 4 ranged between 33% and 45% (mean, 40%) and between 18% and 42% (mean, 28%), respectively. CONCLUSIONS: High single doses of radiation delivered as part of SBRT may induce an elevated and in some cases persistent state of tumor hypoxia in NSCLC tumors. Hypoxia imaging with 18F-fluoromisonidazole positron emission tomography should be used in a larger cohort of NSCLC patients to determine whether elevated tumor hypoxia is predictive of treatment failure in SBRT..

Combined radiotherapy (RT) and hyperthermia (HT) treatments may improve treatment outcome by heat induced radio-sensitisation. We propose an empirical cell survival model (AlphaR model) to describe this multimodality therapy. The model is motivated by the observation that heat induced radio-sensitisation may be explained by a reduction in the DNA damage repair capacity of heated cells. We assume that this repair is only possible up to a threshold level above which survival will decrease exponentially with dose. Experimental cell survival data from two cell lines (HCT116, Cal27) were considered along with that taken from the literature (baby hamster kidney [BHK] and Chinese hamster ovary cells [CHO]) for HT and combined RT-HT. The AlphaR model was used to study the dependence of clonogenic survival on treatment temperature, and thermal dose R2 ≥ 0.95 for all fits). For HT survival curves (0-80 CEM43 at 43.5-57 °C), the number of free fit AlphaR model parameters could be reduced to two. Both parameters increased exponentially with temperature. We derived the relative biological effectiveness (RBE) or HT treatments at different temperatures, to provide an alternative description of thermal dose, based on our AlphaR model. For combined RT-HT, our analysis is restricted to the linear quadratic arm of the model. We show that, for the range used (20-80 CEM43, 0-12 Gy), thermal dose is a valid indicator of heat induced radio-sensitisation, and that the model parameters can be described as a function thereof. Overall, the proposed model provides a flexible framework for describing cell survival curves, and may contribute to better quantification of heat induced radio-sensitisation, and thermal dose in general..

.

Owing to its excellent soft-tissue contrast, magnetic resonance (MR) imaging has found an increased application in radiation therapy (RT). By harnessing these properties for treatment planning, automated segmentation methods can alleviate the manual workload burden to the clinical workflow. We investigated atlas-based segmentation methods of organs at risk (OARs) in the head and neck (H&N) region using one approach that selected the most similar atlas from a library of segmented images and two multi-atlas approaches. The latter were based on weighted majority voting and an iterative atlas-fusion approach called STEPS. We built the atlas library from pre-treatment T1-weighted MR images of 12 patients with manual contours of the parotids, spinal cord and mandible, delineated by a clinician. Following a leave-one-out cross-validation strategy, we measured the geometric accuracy by calculating Dice similarity coefficients (DSC), standard and 95% Hausdorff distances (HD and HD95), and the mean surface distance (MSD), whereby the manual contours served as the gold standard. To benchmark the algorithm, we determined the inter-observer variability (IOV) between three observers. To investigate the dosimetric effect of segmentation inaccuracies, we implemented an auto-planning strategy within the treatment planning system Monaco (Elekta AB, Stockholm, Sweden). For each set of auto-segmented OARs, we generated a plan for a 9-beam step and shoot intensity modulated RT treatment, designed according to our institution's clinical H&N protocol. Superimposing the dose distributions on the gold standard OARs, we calculated dose differences to OARs caused by delineation differences between auto-segmented and gold standard OARs. We investigated the correlations between geometric and dosimetric differences. The mean DSC was larger than 0.8 and the mean MSD smaller than 2 mm for the multi-atlas approaches, resulting in a geometric accuracy comparable to previously published results and within the range of the IOV. While dosimetric differences could be as large as 23% of the clinical goal, treatment plans fulfilled all imposed clinical goals for the gold standard OARs. Correlations between geometric and dosimetric measures were low with R2  <  0.5. The geometric accuracy and the ability to achieve clinically acceptable treatment plans indicate the suitability of using atlas-based contours for RT treatment planning purposes. The low correlations between geometric and dosimetric measures suggest that geometric measures alone are not sufficient to predict the dosimetric impact of segmentation inaccuracies on treatment planning for the data utilised in this study..

Adaptive radiotherapy (ART) strategies systematically monitor variations in target and neighbouring structures to inform treatment-plan modification during radiotherapy. This is necessary because a single plan designed before treatment is insufficient to capture the actual dose delivered to the target and adjacent critical structures during the course of radiotherapy. Magnetic resonance imaging (MRI) provides superior soft-tissue image contrast over current standard X-ray-based technologies without additional radiation exposure. With integrated MRI and radiotherapy platforms permitting motion monitoring during treatment delivery, it is possible that adaption can be informed by real-time anatomical imaging. This allows greater treatment accuracy in terms of dose delivered to target with smaller, individualised treatment margins. The use of functional MRI sequences would permit ART to be informed by imaging biomarkers, so allowing both personalised geometric and biological adaption. In this review, we discuss ART solutions enabled by MRI guidance and its potential gains for our patients across tumour types..

BACKGROUND: Current standard radiotherapy for oropharynx cancer (OPC) is associated with high rates of severe toxicities, shown to adversely impact patients' quality of life. Given excellent outcomes of human papilloma virus (HPV)-associated OPC and long-term survival of these typically young patients, treatment de-intensification aimed at improving survivorship while maintaining excellent disease control is now a central concern. The recent implementation of magnetic resonance image - guided radiotherapy (MRgRT) systems allows for individual tumor response assessment during treatment and offers possibility of personalized dose-reduction. In this 2-stage Bayesian phase II study, we propose to examine weekly radiotherapy dose-adaptation based on magnetic resonance imaging (MRI) evaluated tumor response. Individual patient's plan will be designed to optimize dose reduction to organs at risk and minimize locoregional failure probability based on serial MRI during RT. Our primary aim is to assess the non-inferiority of MRgRT dose adaptation for patients with low risk HPV-associated OPC compared to historical control, as measured by Bayesian posterior probability of locoregional control (LRC). METHODS: Patients with T1-2 N0-2b (as per AJCC 7th Edition) HPV-positive OPC, with lymph node <3 cm and <10 pack-year smoking history planned for curative radiotherapy alone to a dose of 70 Gy in 33 fractions will be eligible. All patients will undergo pre-treatment MRI and at least weekly intra-treatment MRI. Patients undergoing MRgRT will have weekly adaptation of high dose planning target volume based on gross tumor volume response. The stage 1 of this study will enroll 15 patients to MRgRT dose adaptation. If LRC at 6 months with MRgRT dose adaptation is found sufficiently safe as per the Bayesian model, stage 2 of the protocol will expand enrollment to an additional 60 patients, randomized to either MRgRT or standard IMRT. DISCUSSION: Multiple methods for safe treatment de-escalation in patients with HPV-positive OPC are currently being studied. By leveraging the ability of advanced MRI techniques to visualize tumor and soft tissues through the course of treatment, this protocol proposes a workflow for safe personalized radiation dose-reduction in good responders with radiosensitive tumors, while ensuring tumoricidal dose to more radioresistant tumors. MRgRT dose adaptation could translate in reduced long term radiation toxicities and improved survivorship while maintaining excellent LRC outcomes in favorable OPC. TRIAL REGISTRATION: ClinicalTrials.gov ID: NCT03224000; Registration date: 07/21/2017..

BACKGROUND AND PURPOSE: The superior soft-tissue contrast of 4D-T2w MRI motivates its use for delineation in radiotherapy treatment planning. We address current limitations of slice-selective implementations, including thick slices and artefacts originating from data incompleteness and variable breathing. MATERIALS AND METHODS: A method was developed to calculate midposition and 4D-T2w images of the whole thorax from continuously acquired axial and sagittal 2D-T2w MRI (1.5 × 1.5 × 5.0 mm3). The method employed image-derived respiratory surrogates, deformable image registration and super-resolution reconstruction. Volunteer imaging and a respiratory motion phantom were used for validation. The minimum number of dynamic acquisitions needed to calculate a representative midposition image was investigated by retrospectively subsampling the data (10-30 dynamic acquisitions). RESULTS: Super-resolution 4D-T2w MRI (1.0 × 1.0 × 1.0 mm3, 8 respiratory phases) did not suffer from data incompleteness and exhibited reduced stitching artefacts compared to sorted multi-slice MRI. Experiments using a respiratory motion phantom and colour-intensity projection images demonstrated a minor underestimation of the motion range. Midposition diaphragm differences in retrospectively subsampled acquisitions were <1.1 mm compared to the full dataset. 10 dynamic acquisitions were found sufficient to generate midposition MRI. CONCLUSIONS: A motion-modelling and super-resolution method was developed to calculate high quality 4D/midposition T2w MRI from orthogonal 2D-T2w MRI..

Microbeam radiation therapy (MRT) is still a preclinical approach in radiation oncology that uses planar micrometre wide beamlets with extremely high peak doses, separated by a few hundred micrometre wide low dose regions. Abundant preclinical evidence demonstrates that MRT spares normal tissue more effectively than conventional radiation therapy, at equivalent tumour control. In order to launch first clinical trials, accurate and efficient dose calculation methods are an inevitable prerequisite. In this work a hybrid dose calculation approach is presented that is based on a combination of Monte Carlo and kernel based dose calculation. In various examples the performance of the algorithm is compared to purely Monte Carlo and purely kernel based dose calculations. The accuracy of the developed algorithm is comparable to conventional pure Monte Carlo calculations. In particular for inhomogeneous materials the hybrid dose calculation algorithm out-performs purely convolution based dose calculation approaches. It is demonstrated that the hybrid algorithm can efficiently calculate even complicated pencil beam and cross firing beam geometries. The required calculation times are substantially lower than for pure Monte Carlo calculations..

2D cine MR imaging may be utilized to monitor rapidly moving tumors and organs-at-risk for real-time adaptive radiotherapy. This study systematically investigates the impact of geometric imaging parameters on the ability of 2D cine MR imaging to guide template-matching-driven autocontouring of lung tumors and abdominal organs. Abdominal 4D MR images were acquired of six healthy volunteers and thoracic 4D MR images were obtained of eight lung cancer patients. At each breathing phase of the images, the left kidney and gallbladder or lung tumor, respectively, were outlined as volumes of interest. These images and contours were used to create artificial 2D cine MR images, while simultaneously serving as 3D ground truth. We explored the impact of five different imaging parameters (pixel size, slice thickness, imaging plane orientation, number and relative alignment of images as well as strategies to create training images). For each possible combination of imaging parameters, we generated artificial 2D cine MR images as training and test images. A template-matching algorithm used the training images to determine the tumor or organ position in the test images. Subsequently, a 3D base contour was shifted to the determined position and compared to the ground truth via centroid distance and Dice similarity coefficient. The median centroid distance between adapted and ground truth contours was 1.56 mm for the kidney, 3.81 mm for the gallbladder and 1.03 mm for the lung tumor (median Dice similarity coefficient: 0.95, 0.72 and 0.93). We observed that a decrease in image resolution led to a modest decrease in localization accuracy, especially for the small gallbladder. However, for all volumes of interest localization accuracy varied substantially more between subjects than due to the different imaging parameters. Automated tumor and organ localization using 2D cine MR imaging and template-matching-based autocontouring is robust against variation of geometric imaging parameters. Future work and optimization efforts of 2D cine MR imaging for real-time adaptive radiotherapy is needed to characterize the influence of sequence- and anatomical site-specific imaging contrast..

Dosimetric quality assurance (QA) of the new Elekta Unity (MR-linac) will differ from the QA performed of a conventional linac due to the constant magnetic field, which creates an electron return effect (ERE). In this work we aim to validate PRESAGE® dosimetry in a transverse magnetic field, and assess its use to validate the research version of the Monaco TPS of the MR-linac. Cylindrical samples of PRESAGE® 3D dosimeter separated by an air gap were irradiated with a cobalt-60 unit, while placed between the poles of an electromagnet at 0.5 T and 1.5 T. This set-up was simulated in EGSnrc/Cavity Monte Carlo (MC) code and relative dose distributions were compared with measurements using 1D and 2D gamma criteria of 3% and 1.5 mm. The irradiation conditions were adapted for the MR-linac and compared with Monaco TPS simulations. Measured and EGSnrc/Cavity simulated profiles showed good agreement with a gamma passing rate of 99.9% for 0.5 T and 99.8% for 1.5 T. Measurements on the MR-linac also compared well with Monaco TPS simulations, with a gamma passing rate of 98.4% at 1.5 T. Results demonstrated that PRESAGE® can accurately measure dose and detect the ERE, encouraging its use as a QA tool to validate the Monaco TPS of the MR-linac for clinically relevant dose distributions at tissue-air boundaries..

In this work we describe an ultra-fast, low-latency implementation of the energy/mass transfer (EMT) mapping method to accumulate dose on deforming geometries such as lung using the central processing unit (CPU). It enables the computation of the actually delivered dose for intensity-modulated radiation therapy on 4D image data in real-time at 25 Hz. In order to accumulate the delivered dose onto a reference phase a pre-calculated deformable vector field is used. The aim of this study is to present an online dose accumulation technique that can be carried out in less than 40 ms to accommodate the machine log update rate of our research linac. Three speed optimisation strategies for the CPU are discussed: single-core optimisation, parallelisation for multiple cores and vectorisation. The single-core implementation accumulates dose in about 1.1 s on a typical high resolution grid for a lung stereotactic body radiation therapy case. Adding parallelisation decreased the runtime to about 50 ms while adding vectorisation satisfied our real-time constraint by further reducing the dose accumulation time to 15 ms without compromising on resolution or accuracy. The presented method allows real-time dose accumulation on deforming patient geometries and has the potential to enable online dose evaluation and re-planning scenarios..

Margins are employed in radiotherapy treatment planning to mitigate the dosimetric effects of geometric uncertainties for the clinical target volume (CTV). Unfortunately, whilst the use of margins can increase the probability that sufficient dose is delivered to the CTV, it can also result in delivering high dose of radiation to surrounding organs at risk (OARs). We expand on our previous work on beam-dependent margins and propose a novel adaptive margin concept, where margins are moulded away from selected OARs for better OAR-high-dose sparing, whilst maintaining similar dose coverage probability to the CTV. This, however, comes at a cost of a larger irradiation volume, and thus can negatively impact other structures. We investigate the impact of the adaptive margin concept when applied to prostate radiotherapy treatments, and compare treatment plans generated using our beam-dependent margins without adaptation, with adaption from the rectum and with adaptation from both the rectum and bladder. Five prostate patients were used in this planning study. All plans achieved similar dose coverage probability, and were able to ensure at least 90% population coverage with the target receiving at least 95% of the prescribed dose to [Formula: see text]. We observed overall better high-dose sparing to OARs that were considered when using the adapted beam-dependent PTVs, with the degree of sparing dependent on both the number of OARs under consideration as well as the relative position between the CTV and the OARs..

OBJECTIVE: Mimicking state-of-the-art patient radiotherapy with high-precision irradiators for small animals is expected to advance the understanding of dose-effect relationships and radiobiology in general. We work on the implementation of intensity-modulated radiotherapy-like irradiation schemes for small animals. As a first step, we present a fast analytical dose calculation algorithm for keV photon beams. METHODS: We follow a superposition-convolution approach adapted to kV X-rays, based on previous work for microbeam therapy. We assume local energy deposition at the photon interaction point due to the short electron ranges in tissue. This allows us to separate the dose calculation into locally absorbed primary dose and the scatter contribution, calculated in a point kernel approach. We validate our dose model against Geant4 Monte Carlo (MC) simulations and compare the results to Muriplan (XStrahl Ltd, Camberley, UK). RESULTS: For field sizes of (1 mm)2 to (1 cm)2 in water, the depth dose curves show a mean disagreement of 1.7% to MC simulations, with the largest deviations in the entrance region (4%) and at large depths (5% at 7 cm). Larger discrepancies are observed at water-to-bone boundaries, in bone and at the beam edges in slab phantoms and a mouse brain. Calculation times are in the order of 5 s for a single beam. CONCLUSION: The algorithm shows good agreement with MC simulations in an initial validation. It has the potential to become an alternative to full MC dose calculation. Advances in knowledge: The presented algorithm demonstrates the potential of kernel-based dose calculation for kV photon beams. It will be valuable in intensity-modulated radiotherapy and inverse treatment planning for high precision small-animal radiotherapy..

BACKGROUND AND PURPOSE: Computed tomography (CT) imaging is the current gold standard for radiotherapy treatment planning (RTP). The establishment of a magnetic resonance imaging (MRI) only RTP workflow requires the generation of a synthetic CT (sCT) for dose calculation. This study evaluates the feasibility of using a multi-atlas sCT synthesis approach (sCTa) for head and neck and prostate patients. MATERIAL AND METHODS: The multi-atlas method was based on pairs of non-rigidly aligned MR and CT images. The sCTa was obtained by registering the MRI atlases to the patient's MRI and by fusing the mapped atlases according to morphological similarity to the patient. For comparison, a bulk density assignment approach (sCTbda) was also evaluated. The sCTbda was obtained by assigning density values to MRI tissue classes (air, bone and soft-tissue). After evaluating the synthesis accuracy of the sCTs (mean absolute error), sCT-based delineations were geometrically compared to the CT-based delineations. Clinical plans were re-calculated on both sCTs and a dose-volume histogram and a gamma analysis was performed using the CT dose as ground truth. RESULTS: Results showed that both sCTs were suitable to perform clinical dose calculations with mean dose differences less than 1% for both the planning target volume and the organs at risk. However, only the sCTa provided an accurate and automatic delineation of bone. CONCLUSIONS: Combining MR delineations with our multi-atlas CT synthesis method could enable MRI-only treatment planning and thus improve the dosimetric and geometric accuracy of the treatment, and reduce the number of imaging procedures..

In radiotherapy, abdominal and thoracic sites are candidates for performing motion tracking. With real-time control it is possible to adjust the multileaf collimator (MLC) position to the target position. However, positions are not perfectly matched and position errors arise from system delays and complicated response of the electromechanic MLC system. Although, it is possible to compensate parts of these errors by using predictors, residual errors remain and need to be compensated to retain target coverage. This work presents a method to statistically describe tracking errors and to automatically derive a patient-specific, per-segment margin to compensate the arising underdosage on-line, i.e. during plan delivery. The statistics of the geometric error between intended and actual machine position are derived using kernel density estimators. Subsequently a margin is calculated on-line according to a selected coverage parameter, which determines the amount of accepted underdosage. The margin is then applied onto the actual segment to accommodate the positioning errors in the enlarged segment. The proof-of-concept was tested in an on-line tracking experiment and showed the ability to recover underdosages for two test cases, increasing [Formula: see text] in the underdosed area about [Formula: see text] and [Formula: see text], respectively. The used dose model was able to predict the loss of dose due to tracking errors and could be used to infer the necessary margins. The implementation had a running time of 23 ms which is compatible with real-time requirements of MLC tracking systems. The auto-adaptivity to machine and patient characteristics makes the technique a generic yet intuitive candidate to avoid underdosages due to MLC tracking errors..

Near real-time application of Monte Carlo (MC) dose calculation in clinic and research is hindered by the long computational runtimes of established software. Currently, fast MC software solutions are available utilising accelerators such as graphical processing units (GPUs) or clusters based on central processing units (CPUs). Both platforms are expensive in terms of purchase costs and maintenance and, in case of the GPU, provide only limited scalability. In this work we propose a cloud-based MC solution, which offers high scalability of accurate photon dose calculations. The MC simulations run on a private virtual supercomputer that is formed in the cloud. Computational resources can be provisioned dynamically at low cost without upfront investment in expensive hardware. A client-server software solution has been developed which controls the simulations and transports data to and from the cloud efficiently and securely. The client application integrates seamlessly into a treatment planning system. It runs the MC simulation workflow automatically and securely exchanges simulation data with the server side application that controls the virtual supercomputer. Advanced encryption standards were used to add an additional security layer, which encrypts and decrypts patient data on-the-fly at the processor register level. We could show that our cloud-based MC framework enables near real-time dose computation. It delivers excellent linear scaling for high-resolution datasets with absolute runtimes of 1.1 seconds to 10.9 seconds for simulating a clinical prostate and liver case up to 1% statistical uncertainty. The computation runtimes include the transportation of data to and from the cloud as well as process scheduling and synchronisation overhead. Cloud-based MC simulations offer a fast, affordable and easily accessible alternative for near real-time accurate dose calculations to currently used GPU or cluster solutions..

To tackle the problem of magnetic resonance imaging (MRI)-only radiotherapy treatment planning (RTP), we propose a multi-atlas information propagation scheme that jointly segments organs and generates pseudo x-ray computed tomography (CT) data from structural MR images (T1-weighted and T2-weighted). As the performance of the method strongly depends on the quality of the atlas database composed of multiple sets of aligned MR, CT and segmented images, we also propose a robust way of registering atlas MR and CT images, which combines structure-guided registration, and CT and MR image synthesis. We first evaluated the proposed framework in terms of segmentation and CT synthesis accuracy on 15 subjects with prostate cancer. The segmentations obtained with the proposed method were compared using the Dice score coefficient (DSC) to the manual segmentations. Mean DSCs of 0.73, 0.90, 0.77 and 0.90 were obtained for the prostate, bladder, rectum and femur heads, respectively. The mean absolute error (MAE) and the mean error (ME) were computed between the reference CTs (non-rigidly aligned to the MRs) and the pseudo CTs generated with the proposed method. The MAE was on average [Formula: see text] HU and the ME [Formula: see text] HU. We then performed a dosimetric evaluation by re-calculating plans on the pseudo CTs and comparing them to the plans optimised on the reference CTs. We compared the cumulative dose volume histograms (DVH) obtained for the pseudo CTs to the DVH obtained for the reference CTs in the planning target volume (PTV) located in the prostate, and in the organs at risk at different DVH points. We obtained average differences of [Formula: see text] in the PTV for [Formula: see text], and between [Formula: see text] and 0.05% in the PTV, bladder, rectum and femur heads for D mean and [Formula: see text]. Overall, we demonstrate that the proposed framework is able to automatically generate accurate pseudo CT images and segmentations in the pelvic region, potentially bypassing the need for CT scan for accurate RTP..

.

OBJECTIVES: The aim of this study was to develop and verify a method to obtain good temporal resolution T2-weighted 4-dimensional (4D-T2w) magnetic resonance imaging (MRI) by using motion information from T1-weighted 4D (4D-T1w) MRI, to support treatment planning in MR-guided radiotherapy. MATERIALS AND METHODS: Ten patients with primary non-small cell lung cancer were scanned at 1.5 T axially with a volumetric T2-weighted turbo spin echo sequence gated to exhalation and a volumetric T1-weighted stack-of-stars spoiled gradient echo sequence with golden angle spacing acquired in free breathing. From the latter, 20 respiratory phases were reconstructed using the recently developed 4D joint MoCo-HDTV algorithm based on the self-gating signal obtained from the k-space center. Motion vector fields describing the respiratory cycle were obtained by deformable image registration between the respiratory phases and projected onto the T2-weighted image volume. The resulting 4D-T2w volumes were verified against the 4D-T1w volumes: an edge-detection method was used to measure the diaphragm positions; the locations of anatomical landmarks delineated by a radiation oncologist were compared and normalized mutual information was calculated to evaluate volumetric image similarity. RESULTS: High-resolution 4D-T2w MRI was obtained. Respiratory motion was preserved on calculated 4D-T2w MRI, with median diaphragm positions being consistent with less than 6.6 mm (2 voxels) for all patients and less than 3.3 mm (1 voxel) for 9 of 10 patients. Geometrical positions were coherent between 4D-T1w and 4D-T2w MRI as Euclidean distances between all corresponding anatomical landmarks agreed to within 7.6 mm (Euclidean distance of 2 voxels) and were below 3.8 mm (Euclidean distance of 1 voxel) for 355 of 470 pairs of anatomical landmarks. Volumetric image similarity was commensurate between 4D-T1w and 4D-T2w MRI, as mean percentage differences in normalized mutual information (calculated over all respiratory phases and patients), between corresponding respiratory phases of 4D-T1w and 4D-T2w MRI and the tie-phase of 4D-T1w and 3-dimensional T2w MRI, were consistent to 0.41% ± 0.37%. Four-dimensional T2w MRI displayed tumor extent, structure, and position more clearly than corresponding 4D-T1w MRI, especially when mobile tumor sites were adjacent to organs at risk. CONCLUSIONS: A methodology to obtain 4D-T2w MRI that retrospectively applies the motion information from 4D-T1w MRI to 3-dimensional T2w MRI was developed and verified. Four-dimensional T2w MRI can assist clinicians in delineating mobile lesions that are difficult to define on 4D-T1w MRI, because of poor tumor-tissue contrast..

PURPOSE: Anchored electromagnetic transponders for tumor motion monitoring during lung radiotherapy were clinically evaluated. First, intrafractional motion patterns were analyzed as well as their interfractional variations. Second, intra- and interfractional changes of the geometric transponder positions were investigated. MATERIALS AND METHODS: Intrafractional motion data from 7 patients with an upper or middle lobe tumor and three implanted transponders each was used to calculate breathing amplitudes, overall motion amount and motion midlines in three mutual perpendicular directions and three-dimensionally (3D) for 162 fractions. For 6 patients intra- and interfractional variations in transponder distances and in the size of the triangle defined by the transponder locations over the treatment course were determined. RESULTS: Mean 3D values of all fractions were up to 4.0, 4.6 and 3.4 mm per patient for amplitude, overall motion amount and midline deviation, respectively. Intrafractional transponder distances varied with standard deviations up to 3.2 mm, while a maximal triangle shrinkage of 36.5% over 39 days was observed. CONCLUSIONS: Electromagnetic real-time motion monitoring was feasible for all patients. Detected respiratory motion was on average modest in this small cohort without lower lobe tumors, but changes in motion midline were of the same size as the amplitudes and greater midline motion can be observed in some fractions. Intra- and interfractional variations of the geometric transponder positions can be large, so for reliable motion management correlation between transponder and tumor motion needs to be evaluated per patient..

PURPOSE: Firstly, this study provides a real-time implementation of online dose reconstruction for tracked volumetric arc therapy (VMAT). Secondly, this study describes a novel offline quality assurance tool, based on commercial dose calculation algorithms. METHODS: Online dose reconstruction for VMAT is a computationally challenging task in terms of computer memory usage and calculation speed. To potentially reduce the amount of memory used, we analyzed the impact of beam angle sampling for dose calculation on the accuracy of the dose distribution. To establish the performance of the method, we planned two single-arc VMAT prostate stereotactic body radiation therapy cases for delivery with dynamic MLC tracking. For quality assurance of our online dose reconstruction method we have also developed a stand-alone offline dose reconstruction tool, which utilizes the RayStation treatment planning system to calculate dose. RESULTS: For the online reconstructed dose distributions of the tracked deliveries, we could establish strong resemblance for 72 and 36 beam co-planar equidistant beam samples with less than 1.2% deviation for the assessed dose-volume indicators (clinical target volume D98 and D2, and rectum D2). We could achieve average runtimes of 28-31 ms per reported MLC aperture for both dose computation and accumulation, meeting our real-time requirement. To cross-validate the offline tool, we have compared the planned dose to the offline reconstructed dose for static deliveries and found excellent agreement (3%/3 mm global gamma passing rates of 99.8%-100%). CONCLUSION: Being able to reconstruct dose during delivery enables online quality assurance and online replanning strategies for VMAT. The offline quality assurance tool provides the means to validate novel online dose reconstruction applications using a commercial dose calculation engine..

PURPOSE: This study investigates the feasibility and potential benefits of radiotherapy with a 1.5T MR-Linac for locally advanced non-small cell lung cancer (LA NSCLC) patients. MATERIAL AND METHODS: Ten patients with LA NSCLC were retrospectively re-planned six times: three treatment plans were created according to a protocol for conventionally fractionated radiotherapy and three treatment plans following guidelines for isotoxic target dose escalation. In each case, two plans were designed for the MR-Linac, either with standard (∼7mm) or reduced (∼3mm) planning target volume (PTV) margins, while one conventional linac plan was created with standard margins. Treatment plan quality was evaluated using dose-volume metrics or by quantifying dose escalation potential. RESULTS: All generated treatment plans fulfilled their respective planning constraints. For conventionally fractionated treatments, MR-Linac plans with standard margins had slightly increased skin dose when compared to conventional linac plans. Using reduced margins alleviated this issue and decreased exposure of several other organs-at-risk (OAR). Reduced margins also enabled increased isotoxic target dose escalation. CONCLUSION: It is feasible to generate treatment plans for LA NSCLC patients on a 1.5T MR-Linac. Margin reduction, facilitated by an envisioned MRI-guided workflow, enables increased OAR sparing and isotoxic target dose escalation for the respective treatment approaches..

PURPOSE: The purpose of this study is to use dynamic [18F]fluoromisonidazole ([18F]FMISO) positron emission tomography (PET) to compare estimates of tumor hypoxic fractions (HFs) derived by tracer kinetic modeling, tissue-to-blood ratios (TBR), and independent oxygen (pO2) measurements. PROCEDURES: BALB/c mice with EMT6 subcutaneous tumors were selected for PET imaging and invasive pO2 measurements. Data from 120-min dynamic [18F]FMISO scans were fit to two-compartment irreversible three rate constant (K 1, k 2, k 3) and Patlak models (K i). Tumor HFs were calculated and compared using K i, k 3, TBR, and pO2 values. The clinical impact of each method was evaluated on [18F]FMISO scans for three non-small cell lung cancer (NSCLC) radiotherapy patients. RESULTS: HFs defined by TBR (≥1.2, ≥1.3, and ≥1.4) ranged from 2 to 85 % of absolute tumor volume. HFs defined by K i (>0.004 ml min cm-3) and k 3 (>0.008 min-1) varied from 9 to 85 %. HF quantification was highly dependent on metric (TBR, k 3, or K i) and threshold. HFs quantified on human [18F]FMISO scans varied from 38 to 67, 0 to 14, and 0.1 to 27 %, for each patient, respectively, using TBR, k 3, and K i metrics. CONCLUSIONS: [18F]FMISO PET imaging metric choice and threshold impacts hypoxia quantification reliability. Our results suggest that tracer kinetic modeling has the potential to improve hypoxia quantification clinically as it may provide a stronger correlation with direct pO2 measurements..

BACKGROUND AND PURPOSE: Radiotherapy guidance based on magnetic resonance imaging (MRI) is currently becoming a clinical reality. Fast 2d cine MRI sequences are expected to increase the precision of radiation delivery by facilitating tumour delineation during treatment. This study compares four auto-contouring algorithms for the task of delineating the primary tumour in six locally advanced (LA) lung cancer patients. MATERIAL AND METHODS: Twenty-two cine MRI sequences were acquired using either a balanced steady-state free precession or a spoiled gradient echo imaging technique. Contours derived by the auto-contouring algorithms were compared against manual reference contours. A selection of eight image data sets was also used to assess the inter-observer delineation uncertainty. RESULTS: Algorithmically derived contours agreed well with the manual reference contours (median Dice similarity index: ⩾0.91). Multi-template matching and deformable image registration performed significantly better than feature-driven registration and the pulse-coupled neural network (PCNN). Neither MRI sequence nor image orientation was a conclusive predictor for algorithmic performance. Motion significantly degraded the performance of the PCNN. The inter-observer variability was of the same order of magnitude as the algorithmic performance. CONCLUSION: Auto-contouring of tumours on cine MRI is feasible in LA lung cancer patients. Despite large variations in implementation complexity, the different algorithms all have relatively similar performance..

Microbeam radiation therapy (MRT) is a treatment approach in radiation therapy where the treatment field is spatially fractionated into arrays of a few tens of micrometre wide planar beams of unusually high peak doses separated by low dose regions of several hundred micrometre width. In preclinical studies, this treatment approach has proven to spare normal tissue more effectively than conventional radiation therapy, while being equally efficient in tumour control. So far dose calculations in MRT, a prerequisite for future clinical applications are based on Monte Carlo simulations. However, they are computationally expensive, since scoring volumes have to be small. In this article a kernel based dose calculation algorithm is presented that splits the calculation into photon and electron mediated energy transport, and performs the calculation of peak and valley doses in typical MRT treatment fields within a few minutes. Kernels are analytically calculated depending on the energy spectrum and material composition. In various homogeneous materials peak, valley doses and microbeam profiles are calculated and compared to Monte Carlo simulations. For a microbeam exposure of an anthropomorphic head phantom calculated dose values are compared to measurements and Monte Carlo calculations. Except for regions close to material interfaces calculated peak dose values match Monte Carlo results within 4% and valley dose values within 8% deviation. No significant differences are observed between profiles calculated by the kernel algorithm and Monte Carlo simulations. Measurements in the head phantom agree within 4% in the peak and within 10% in the valley region. The presented algorithm is attached to the treatment planning platform VIRTUOS. It was and is used for dose calculations in preclinical and pet-clinical trials at the biomedical beamline ID17 of the European synchrotron radiation facility in Grenoble, France..

Microbeam Radiation Therapy is an innovative pre-clinical strategy which uses arrays of parallel, tens of micrometres wide kilo-voltage photon beams to treat tumours. These x-ray beams are typically generated on a synchrotron source. It was shown that these beam geometries allow exceptional normal tissue sparing from radiation damage while still being effective in tumour ablation. A final biological explanation for this enhanced therapeutic ratio has still not been found, some experimental data support an important role of the vasculature. In this work, the effect of microbeams on a normal microvascular network of the cerebral cortex was assessed in computer simulations and compared to the effect of homogeneous, seamless exposures at equal energy absorption. The anatomy of a cerebral microvascular network and the inflicted radiation damage were simulated to closely mimic experimental data using a novel probabilistic model of radiation damage to blood vessels. It was found that the spatial dose fractionation by microbeam arrays significantly decreased the vascular damage. The higher the peak-to-valley dose ratio, the more pronounced the sparing effect. Simulations of the radiation damage as a function of morphological parameters of the vascular network demonstrated that the distribution of blood vessel radii is a key parameter determining both the overall radiation damage of the vasculature and the dose-dependent differential effect of microbeam irradiation..

Radiotherapy treatment planning for use with high-energy photon beams currently employs a binary approach in defining the planning target volume (PTV). We propose a margin concept that takes the beam directions into account, generating beam-dependent PTVs (bdPTVs) on a beam-by-beam basis. The resulting degree of overlaps between the bdPTVs are used within the optimisation process; the optimiser effectively considers the same voxel to be both target and organ at risk (OAR) with fractional contributions. We investigate the impact of this novel approach when applied to prostate radiotherapy treatments, and compare treatment plans generated using beam dependent margins to conventional margins. Five prostate patients were used in this planning study, and plans using beam dependent margins improved the sparing of high doses to target-surrounding OARs, though a trade-off in delivering additional low dose to the OARs can be observed. Plans using beam dependent margins are observed to have a slightly reduced target coverage. Nevertheless, all plans are able to satisfy 90% population coverage with the target receiving at least 95% of the prescribed dose to [Formula: see text]..

Currently hard coherent x-ray radiation at high photon fluxes can only be produced with large and expensive radiation sources, such as 3[Formula: see text] generation synchrotrons. Especially in medicine, this limitation prevents various promising developments in imaging and therapy from being translated into clinical practice. Here we present a new concept of highly brilliant x-ray sources, line focus x-ray tubes (LFXTs), which may serve as a powerful and cheap alternative to synchrotrons and a range of other existing technologies. LFXTs employ an extremely thin focal spot and a rapidly rotating target for the electron beam which causes a change in the physical mechanism of target heating, allowing higher electron beam intensities at the focal spot. Monte Carlo simulations and numeric solutions of the heat equation are used to predict the characteristics of the LFXT. In terms of photon flux and coherence length, the performance of the line focus x-ray tube compares with inverse Compton scattering sources. Dose rates of up to 180 Gy [Formula: see text] can be reached in 50 cm distance from the focal spot. The results demonstrate that the line focus tube can serve as a powerful compact source for phase contrast imaging and microbeam radiation therapy. The production of a prototype seems technically feasible..

BACKGROUND AND PURPOSE: Adaptive radiotherapy (ART) using plan selection is being introduced clinically for bladder cancer, but the challenge of how to compensate for intra-fractional motion remains. The purpose of this study was to assess target coverage with respect to intra-fractional motion and the potential for normal tissue sparing in MRI-guided ART (MRIGART) using isotropic (MRIGARTiso), an-isotropic (MRIGARTanIso) and population-based margins (MRIGARTpop). MATERIALS AND METHODS: Nine bladder cancer patients treated in a phase II trial of plan selection underwent 6-7 weekly repeat MRI series, each with volumetric scans acquired over a 10 min period. Adaptive re-planning on the 0 min MRI scans was performed using density override, simulating a hypo-fractionated schedule. Target coverage was evaluated on the 10 min scan to quantify the impact of intra-fractional motion. RESULTS: MRIGARTanIso reduced the course-averaged PTV by median 304 cc compared to plan selection. Bladder shifts affected target coverage in individual fractions for all strategies. Two patients had a v95% of the bladder below 98% for MRIGARTiso. MRIGARTiso decreased the bowel V25 with 15-46 cc compared to MRIGARTpop. CONCLUSION: Online re-optimised ART has a considerable normal tissue sparing potential. MRIGART with online corrections for target shift during a treatment fraction should be considered in ART for bladder cancer..

PURPOSE: A study of real-time adaptive radiotherapy systems was performed to test the hypothesis that, across delivery systems and institutions, the dosimetric accuracy is improved with adaptive treatments over non-adaptive radiotherapy in the presence of patient-measured tumor motion. METHODS AND MATERIALS: Ten institutions with robotic(2), gimbaled(2), MLC(4) or couch tracking(2) used common materials including CT and structure sets, motion traces and planning protocols to create a lung and a prostate plan. For each motion trace, the plan was delivered twice to a moving dosimeter; with and without real-time adaptation. Each measurement was compared to a static measurement and the percentage of failed points for γ-tests recorded. RESULTS: For all lung traces all measurement sets show improved dose accuracy with a mean 2%/2mm γ-fail rate of 1.6% with adaptation and 15.2% without adaptation (p<0.001). For all prostate the mean 2%/2mm γ-fail rate was 1.4% with adaptation and 17.3% without adaptation (p<0.001). The difference between the four systems was small with an average 2%/2mm γ-fail rate of <3% for all systems with adaptation for lung and prostate. CONCLUSIONS: The investigated systems all accounted for realistic tumor motion accurately and performed to a similar high standard, with real-time adaptation significantly outperforming non-adaptive delivery methods..

BACKGROUND AND PURPOSE: There are concerns that radiotherapy doses delivered in a magnetic field might be distorted due to the Lorentz force deflecting secondary electrons. This study investigates this effect on lung stereotactic body radiotherapy (SBRT) treatments, conducted either with or without multileaf collimator (MLC) tumor tracking. MATERIAL AND METHODS: Lung SBRT treatments with an MR-linac were simulated for nine patients. Two different treatment techniques were compared: conventional, non-tracked deliveries and deliveries with real-time MLC tumor tracking, each conducted either with or without a 1.5T magnetic field. RESULTS: Slight dose distortions at air-tissue-interfaces were observed in the presence of the magnetic field. Most prominently, the dose to 2% of the skin increased by 1.4Gy on average. Regardless of the presence of the magnetic field, MLC tracking was able to spare healthy tissue, for example by decreasing the mean lung dose by 0.3Gy on average, while maintaining the target dose. CONCLUSIONS: Accounting for the magnetic field during treatment plan optimization allowed for design and delivery of clinically acceptable lung SBRT treatments with an MR-linac. Furthermore, the ability of MLC tumor tracking to decrease dose exposure of healthy tissue, was not inhibited by the magnetic field..

PURPOSE: Upcoming veterinary trials in microbeam radiation therapy (MRT) demand for more advanced irradiation techniques than in preclinical research with small animals. The treatment of deep-seated tumors in cats and dogs with MRT requires sophisticated irradiation geometries from multiple ports, which impose further efforts to spare the normal tissue surrounding the target. METHODS: This work presents the development and benchmarking of a precise patient alignment protocol for MRT at the biomedical beamline ID17 of the European Synchrotron Radiation Facility (ESRF). The positioning of the patient prior to irradiation is verified by taking x-ray projection images from different angles. RESULTS: Using four external fiducial markers of 1.7  mm diameter and computed tomography-based treatment planning, a target alignment error of less than 2  mm can be achieved with an angular deviation of less than 2(∘). Minor improvements on the protocol and the use of smaller markers indicate that even a precision better than 1  mm is technically feasible. Detailed investigations concerning the imaging dose lead to the conclusion that doses for skull radiographs lie in the same range as dose reference levels for human head radiographs. A currently used online dose monitor for MRT has been proven to give reliable results for the imaging beam. CONCLUSIONS: The ESRF biomedical beamline ID17 is technically ready to apply conformal image-guided MRT from multiple ports to large animals during future veterinary trials..

Conventional treatment planning in intensity-modulated radiation therapy (IMRT) is a trial-and-error process that usually involves tedious tweaking of optimization parameters. Here, we present an algorithm that automates part of this process, in particular the adaptation of voxel-based penalties within normal tissue. Thereby, the proposed algorithm explicitly considers a priori known physical limitations of photon irradiation. The efficacy of the developed algorithm is assessed during treatment planning studies comprising 16 prostate and 5 head and neck cases. We study the eradication of hot spots in the normal tissue, effects on target coverage and target conformity, as well as selected dose volume points for organs at risk. The potential of the proposed method to generate class solutions for the two indications is investigated. Run-times of the algorithms are reported. Physically constrained voxel-based penalty adaptation is an adequate means to automatically detect and eradicate hot-spots during IMRT planning while maintaining target coverage and conformity. Negative effects on organs at risk are comparably small and restricted to lower doses. Using physically constrained voxel-based penalty adaptation, it was possible to improve the generation of class solutions for both indications. Considering the reported run-times of less than 20 s, physically constrained voxel-based penalty adaptation has the potential to reduce the clinical workload during planning and automated treatment plan generation in the long run, facilitating adaptive radiation treatments. PACS number(s): 87.55.de..

INTRODUCTION: The efficacy of radiation therapy treatments for pancreatic cancer is compromised by abdominal motion which limits the spatial accuracy for dose delivery - especially for particles. In this work we investigate the potential of worst case optimization for interfractional offline motion mitigation in carbon ion treatments of pancreatic cancer. METHODS: We implement a worst case optimization algorithm that explicitly models the relative biological effectiveness of carbon ions during inverse planning. We perform a comparative treatment planning study for seven pancreatic cancer patients. Treatment plans that have been generated using worst case optimization are compared against (1) conventional intensity-modulated carbon ion therapy, (2) single field uniform dose carbon ion therapy, and (3) an ideal yet impractical scenario relying on daily re-planning. The dosimetric quality and robustness of the resulting treatment plans is evaluated using reconstructions of the daily delivered dose distributions on fractional control CTs. RESULTS: Idealized daily re-planning consistently gives the best dosimetric results with regard to both target coverage and organ at risk sparing. The absolute reduction of D 95 within the gross tumor volume during fractional dose reconstruction is most pronounced for conventional intensity-modulated carbon ion therapy. Single field uniform dose optimization exhibits no substantial reduction for six of seven patients and values for D 95 for worst case optimization fall in between. The treated volume (D>95 % prescription dose) outside of the gross tumor volume is reduced by a factor of two by worst case optimization compared to conventional optimization and single field uniform dose optimization. Single field uniform dose optimization comes at an increased radiation exposure of normal tissues, e.g. ≈2 Gy (RBE) in the mean dose in the kidneys compared to conventional and worst case optimization and ≈4 Gy (RBE) in D 1 in the spinal cord compared to worst case optimization. CONCLUSION: Interfractional motion substantially deteriorates dose distributions for carbon ion treatments of pancreatic cancer patients. Single field uniform dose optimization mitigates the negative influence of motion on target coverage at an increased radiation exposure of normal tissue. Worst case optimization enables an exploration of the trade-off between robust target coverage and organ at risk sparing during inverse treatment planning beyond margin concepts..

PURPOSE: To quantify the performance of the Clarity ultrasound (US) imaging system (Elekta AB, Stockholm, Sweden) for real-time dynamic multileaf collimator (MLC) tracking. METHODS: The Clarity calibration and quality assurance phantom was mounted on a motion platform moving with a periodic sine wave trajectory. The detected position of a 30 mm hypoechogenic sphere within the phantom was continuously reported via Clarity's real-time streaming interface to an in-house tracking and delivery software and subsequently used to adapt the MLC aperture. A portal imager measured MV treatment field/MLC apertures and motion platform positions throughout each experiment to independently quantify system latency and geometric error. Based on the measured range of latency values, a prostate stereotactic body radiation therapy (SBRT) delivery was performed with three realistic motion trajectories. The dosimetric impact of system latency on MLC tracking was directly measured using a 3D dosimeter mounted on the motion platform. RESULTS: For 2D US imaging, the overall system latency, including all delay times from the imaging and delivery chain, ranged from 392 to 424 ms depending on the lateral sector size. For 3D US imaging, the latency ranged from 566 to 1031 ms depending on the elevational sweep. The latency-corrected geometric root-mean squared error was below 0.75 mm (2D US) and below 1.75 mm (3D US). For the prostate SBRT delivery, the impact of a range of system latencies (400-1000 ms) on the MLC tracking performance was minimal in terms of gamma failure rate. CONCLUSIONS: Real-time MLC tracking based on a noninvasive US input is technologically feasible. Current system latencies are higher than those for x-ray imaging systems, but US can provide full volumetric image data and the impact of system latency was measured to be small for a prostate SBRT case when using a US-like motion input..

PURPOSE: This study provides a proof of concept for real-time 4D dose reconstruction for lung stereotactic body radiation therapy (SBRT) with multileaf collimator (MLC) tracking and assesses the impact of tumor tracking on the size of target margins. METHODS: The authors have implemented real-time 4D dose reconstruction by connecting their tracking and delivery software to an Agility MLC at an Elekta Synergy linac and to their in-house treatment planning software (TPS). Actual MLC apertures and (simulated) target positions are reported to the TPS every 40 ms. The dose is calculated in real-time from 4DCT data directly after each reported aperture by utilization of precalculated dose-influence data based on a Monte Carlo algorithm. The dose is accumulated onto the peak-exhale (reference) phase using energy-mass transfer mapping. To investigate the impact of a potentially reducible safety margin, the authors have created and delivered treatment plans designed for a conventional internal target volume (ITV) + 5 mm, a midventilation approach, and three tracking scenarios for four lung SBRT patients. For the tracking plans, a moving target volume (MTV) was established by delineating the gross target volume (GTV) on every 4DCT phase. These were rigidly aligned to the reference phase, resulting in a unified maximum GTV to which a 1, 3, or 5 mm isotropic margin was added. All scenarios were planned for 9-beam step-and-shoot IMRT to meet the criteria of RTOG 1021 (3 × 18 Gy). The GTV 3D center-of-volume shift varied from 6 to 14 mm. RESULTS: Real-time dose reconstruction at 25 Hz could be realized on a single workstation due to the highly efficient implementation of dose calculation and dose accumulation. Decreased PTV margins resulted in inadequate target coverage during untracked deliveries for patients with substantial tumor motion. MLC tracking could ensure the GTV target dose for these patients. Organ-at-risk (OAR) doses were consistently reduced by decreased PTV margins. The tracked MTV + 1 mm deliveries resulted in the following OAR dose reductions: lung V20 up to 3.5%, spinal cord D2 up to 0.9 Gy/Fx, and proximal airways D2 up to 1.4 Gy/Fx. CONCLUSIONS: The authors could show that for patient data at clinical resolution and realistic motion conditions, the delivered dose could be reconstructed in 4D for the whole lung volume in real-time. The dose distributions show that reduced margins yield lower doses to healthy tissue, whilst target dose can be maintained using dynamic MLC tracking..

PURPOSE: Microbeam radiation therapy is an innovative treatment approach in radiation therapy that uses arrays of a few tens of micrometer wide and a few hundreds of micrometer spaced planar x-ray beams as treatment fields. In preclinical studies these fields efficiently eradicated tumors while normal tissue could effectively be spared. However, development and clinical application of microbeam radiation therapy is impeded by a lack of suitable small scale sources. Until now, only large synchrotrons provide appropriate beam properties for the production of microbeams. METHODS: In this work, a conventional x-ray tube with a small focal spot and a specially designed collimator are used to produce microbeams for preclinical research. The applicability of the developed source is demonstrated in a pilot in vitro experiment. The properties of the produced radiation field are characterized by radiochromic film dosimetry. RESULTS: 50 μm wide and 400 μm spaced microbeams were produced in a 20 × 20 mm2 sized microbeam field. The peak to valley dose ratio ranged from 15.5 to 30, which is comparable to values obtained at synchrotrons. A dose rate of up to 300 mGy/s was achieved in the microbeam peaks. Analysis of DNA double strand repair and cell cycle distribution after in vitro exposures of pancreatic cancer cells (Panc1) at the x-ray tube and the European Synchrotron leads to similar results. In particular, a reduced G2 cell cycle arrest is observed in cells in the microbeam peak region. CONCLUSIONS: At its current stage, the source is restricted to in vitro applications. However, moderate modifications of the setup may soon allow in vivo research in mice and rats..

By adapting to the actual patient anatomy during treatment, tracked multi-leaf collimator (MLC) treatment deliveries offer an opportunity for margin reduction and healthy tissue sparing. This is assumed to be especially relevant for hypofractionated protocols in which intrafractional motion does not easily average out. In order to confidently deliver tracked treatments with potentially reduced margins, it is necessary to monitor not only the patient anatomy but also the actually delivered dose during irradiation. In this study, we present a novel real-time online dose reconstruction tool which calculates actually delivered dose based on pre-calculated dose influence data in less than 10 ms at a rate of 25 Hz. Using this tool we investigate the impact of clinical target volume (CTV) to planning target volume (PTV) margins on CTV coverage and organ-at-risk dose. On our research linear accelerator, a set of four different CTV-to-PTV margins were tested for three patient cases subject to four different motion conditions. Based on this data, we can conclude that tracking eliminates dose cold spots which can occur in the CTV during conventional deliveries even for the smallest CTV-to-PTV margin of 1 mm. Changes of organ-at-risk dose do occur frequently during MLC tracking and are not negligible in some cases. Intrafractional dose reconstruction is expected to become an important element in any attempt of re-planning the treatment plan during the delivery based on the observed anatomy of the day..

Recently we introduced interactive dose shaping (IDS) as a new IMRT planning strategy. This planning concept is based on a hierarchical sequence of local dose modification and recovery operations. The purpose of this work is to provide a feasibility study for the IDS planning strategy based on a small set of six prostate patients. The IDS planning paradigm aims to perform interactive local dose adaptations of an IMRT plan without compromising already established valuable dose features in real-time. Various IDS tools were developed in our in-house treatment planning software Dynaplan and were utilized to create IMRT treatment plans for six patients with an adeno-carcinoma of the prostate. The sequenced IDS treatment plans were compared to conventionally optimized clinically approved plans (9 beams, co-planar). For each patient, several IDS plans were created, with different trade-offs between organ sparing and target coverage. The reference dose distributions were imported into Dynaplan. For each patient, the IDS treatment plan with a similar or better trade-off between target coverage and OAR sparing was selected for plan evaluation, guided by a physician. For this initial study we were able to generate treatment plans for prostate geometries in 15-45 min. Individual local dose adaptations could be performed in less than one second. The average differences compared to the reference plans were for the mean dose: 0.0 Gy (boost) and 1.2 Gy (PTV), for D98% : -1.1 Gy and for D2% : 1.1 Gy (both target volumes). The dose-volume quality indicators were well below the Quantec constraints. However, we also observed limitations of our currently implemented approach. Most prominent was an increase of the non-tumor integral dose by 16.4% on average, demonstrating that further developments of our planning strategy are required..

In this work we present a novel treatment planning technique called interactive dose shaping (IDS) to be employed for the optimization of intensity modulated radiation therapy (IMRT). IDS does not rely on a Newton-based optimization algorithm which is driven by an objective function formed of dose volume constraints on pre-segmented volumes of interest (VOIs). Our new planning technique allows for direct, interactive adaptation of localized planning features. This is realized by a dose modification and recovery (DMR) planning engine which implements a two-step approach: firstly, the desired localized plan adaptation is imposed on the current plan (modification) while secondly inevitable, undesired disturbances of the dose pattern elsewhere are compensated for automatically by the recovery module. Together with an ultra-fast dose update calculation method the DMR engine has been implemented in a newly designed 3D therapy planning system Dynaplan enabling true real-time interactive therapy planning. Here we present the underlying strategy and algorithms of the DMR based planning concept. The functionality of the IDS planning approach is demonstrated for a phantom geometry of clinical resolution and size..

PURPOSE: To investigate the influence of interfractional changes on the delivered dose of intensity modulated proton (IMPT) and photon plans (IMXT). METHODS AND MATERIALS: Five postoperative head and neck cancer patients, previously treated with tomotherapy at our institute, were analyzed. The planning study is based on megavoltage (MV) control images. For each patient one IMPT plan and one IMXT plan were generated on the first MV-CT and recalculated on weekly control MV-CTs in the actual treatment position. Dose criteria for evaluation were coverage and conformity of the planning target volume (PTV), as well as mean dose to parotids and maximum dose to spinal cord. RESULTS: Considerable dosimetric changes were observed for IMPT and IMXT plans. Proton plans showed a more pronounced increase of maximum dose and decrease of minimum dose with local underdosage occurring even in the center of the PTV (worst IMPT vs. IMXT coverage: 66.7% vs. 85.0%). The doses to organs at risk (OARs) increased during the treatment period. However, the OAR doses of IMPT stayed below corresponding IMXT values at any time. For both modalities treatment plans did not necessarily worsen monotonically throughout the treatment. CONCLUSIONS: Although absolute differences between planned and reconstructed doses were larger in IMPT plans, doses to OARs were higher in IMXT plans. Tumor coverage was more stable in IMXT plans; IMPT dose distributions indicated a high risk for local underdosage during the treatment course..

PURPOSE: The authors investigated the potential of optimized noncoplanar irradiation trajectories for volumetric modulated arc therapy (VMAT) treatments of nasopharyngeal patients and studied the trade-off between treatment plan quality and delivery time in radiation therapy. METHODS: For three nasopharyngeal patients, the authors generated treatment plans for nine different delivery scenarios using dedicated optimization methods. They compared these scenarios according to dose characteristics, number of beam directions, and estimated delivery times. In particular, the authors generated the following treatment plans: (1) a 4π plan, which is a not sequenced, fluence optimized plan that uses beam directions from approximately 1400 noncoplanar directions and marks a theoretical upper limit of the treatment plan quality, (2) a coplanar 2π plan with 72 coplanar beam directions as pendant to the noncoplanar 4π plan, (3) a coplanar VMAT plan, (4) a coplanar step and shoot (SnS) plan, (5) a beam angle optimized (BAO) coplanar SnS IMRT plan, (6) a noncoplanar BAO SnS plan, (7) a VMAT plan with rotated treatment couch, (8) a noncoplanar VMAT plan with an optimized great circle around the patient, and (9) a noncoplanar BAO VMAT plan with an arbitrary trajectory around the patient. RESULTS: VMAT using optimized noncoplanar irradiation trajectories reduced the mean and maximum doses in organs at risk compared to coplanar VMAT plans by 19% on average while the target coverage remains constant. A coplanar BAO SnS plan was superior to coplanar SnS or VMAT; however, noncoplanar plans like a noncoplanar BAO SnS plan or noncoplanar VMAT yielded a better plan quality than the best coplanar 2π plan. The treatment plan quality of VMAT plans depended on the length of the trajectory. The delivery times of noncoplanar VMAT plans were estimated to be 6.5 min in average; 1.6 min longer than a coplanar plan but on average 2.8 min faster than a noncoplanar SnS plan with comparable treatment plan quality. CONCLUSIONS: The authors' study reconfirms the dosimetric benefits of noncoplanar irradiation of nasopharyngeal tumors. Both SnS using optimized noncoplanar beam ensembles and VMAT using an optimized, arbitrary, noncoplanar trajectory enabled dose reductions in organs at risk compared to coplanar SnS and VMAT. Using great circles or simple couch rotations to implement noncoplanar VMAT, however, was not sufficient to yield meaningful improvements in treatment plan quality. The authors estimate that noncoplanar VMAT using arbitrary optimized irradiation trajectories comes at an increased delivery time compared to coplanar VMAT yet at a decreased delivery time compared to noncoplanar SnS IMRT..

This study demonstrates the feasibility of automated marker tracking for the real-time detection of intrafractional target motion using noisy kilovoltage (kV) X-ray images degraded by the megavoltage (MV) treatment beam. The authors previously introduced the in-line imaging geometry, in which the flat-panel detector (FPD) is mounted directly underneath the treatment head of the linear accelerator. They found that the 121 kVp image quality was severely compromised by the 6 MV beam passing through the FPD at the same time. Specific MV-induced artefacts present a considerable challenge for automated marker detection algorithms. For this study, the authors developed a new imaging geometry by re-positioning the FPD and the X-ray tube. This improved the contrast-to-noise-ratio between 40% and 72% at the 1.2 mAs/image exposure setting. The increase in image quality clearly facilitates the quick and stable detection of motion with the aid of a template matching algorithm. The setup was tested with an anthropomorphic lung phantom (including an artificial lung tumour). In the tumour one or three Calypso beacons were embedded to achieve better contrast during MV radiation. For a single beacon, image acquisition and automated marker detection typically took around 76 ± 6 ms. The success rate was found to be highly dependent on imaging dose and gantry angle. To eliminate possible false detections, the authors implemented a training phase prior to treatment beam irradiation and also introduced speed limits for motion between subsequent images..

PURPOSE: Microbeam radiation therapy (MRT) is a still preclinical tumor therapy approach that uses arrays of a few tens of micrometer wide parallel beams separated by a few 100 μm. The production, measurement, and planning of such radiation fields are a challenge up to now. Here, the authors investigate the feasibility of radiochromic film dosimetry in combination with a microscopic readout as a tool to validate peak and valley doses in MRT, which is an important requirement for a future clinical application of the therapy. METHODS: Gafchromic(®) HD-810 and HD-V2 films are exposed to MRT fields at the biomedical beamline ID17 of the European Synchrotron Radiation Facility (ESRF) and are afterward scanned with a microscope. The measured dose is compared with Monte Carlo calculations. Image analysis tools and film handling protocols are developed that allow accurate and reproducible dosimetry. The performance of HD-810 and HD-V2 films is compared and a detailed analysis of the resolution, noise, and energy dependence is carried out. Measurement uncertainties are identified and analyzed. RESULTS: The dose was measured with a resolution of 5 × 1000 μm(2) and an accuracy of 5% in the peak and between 10% and 15% in the valley region. As main causes for dosimetry uncertainties, statistical noise, film inhomogeneities, and calibration errors were identified. Calibration errors strongly increase at low doses and exceeded 3% for doses below 50 and 70 Gy for HD-V2 and HD-810 films, respectively. While the grain size of both film types is approximately 2 μm, the statistical noise in HD-V2 is much higher than in HD-810 films. However, HD-810 films show a higher energy dependence at low photon energies. CONCLUSIONS: Both film types are appropriate for dosimetry in MRT and the microscope is superior to the microdensitometer used before at the ESRF with respect to resolution and reproducibility. However, a very careful analysis of the image data is required. Dosimetry at low photon energies should be performed with great caution due to the energy sensitivity of the films. In this respect, HD-V2 films showed to have an advantage over HD-810 films. However, HD-810 films have a lower statistical noise level. When a higher resolution is required, e.g., for the dosimetry of pencil beam irradiations, noise may render HD-V2 films inapplicable..

Monte-Carlo (MC) simulations are considered to be the most accurate method for calculating dose distributions in radiotherapy. Its clinical application, however, still is limited by the long runtimes conventional implementations of MC algorithms require to deliver sufficiently accurate results on high resolution imaging data. In order to overcome this obstacle we developed the software-package PhiMC, which is capable of computing precise dose distributions in a sub-minute time-frame by leveraging the potential of modern many- and multi-core CPU-based computers. PhiMC is based on the well verified dose planning method (DPM). We could demonstrate that PhiMC delivers dose distributions which are in excellent agreement to DPM. The multi-core implementation of PhiMC scales well between different computer architectures and achieves a speed-up of up to 37[Formula: see text] compared to the original DPM code executed on a modern system. Furthermore, we could show that our CPU-based implementation on a modern workstation is between 1.25[Formula: see text] and 1.95[Formula: see text] faster than a well-known GPU implementation of the same simulation method on a NVIDIA Tesla C2050. Since CPUs work on several hundreds of GB RAM the typical GPU memory limitation does not apply for our implementation and high resolution clinical plans can be calculated..

The risk of developing normal tissue injuries often limits the radiation dose that can be applied to the tumour in radiation therapy. Microbeam Radiation Therapy (MRT), a spatially fractionated photon radiotherapy is currently tested at the European Synchrotron Radiation Facility (ESRF) to improve normal tissue protection. MRT utilizes an array of microscopically thin and nearly parallel X-ray beams that are generated by a synchrotron. At the ion microprobe SNAKE in Munich focused proton microbeams ("proton microchannels") are studied to improve normal tissue protection. Here, we comparatively investigate microbeam/microchannel irradiations with sub-millimetre X-ray versus proton beams to minimize the risk of normal tissue damage in a human skin model, in vitro. Skin tissues were irradiated with a mean dose of 2 Gy over the irradiated area either with parallel synchrotron-generated X-ray beams at the ESRF or with 20 MeV protons at SNAKE using four different irradiation modes: homogeneous field, parallel lines and microchannel applications using two different channel sizes. Normal tissue viability as determined in an MTT test was significantly higher after proton or X-ray microchannel irradiation compared to a homogeneous field irradiation. In line with these findings genetic damage, as determined by the measurement of micronuclei in keratinocytes, was significantly reduced after proton or X-ray microchannel compared to a homogeneous field irradiation. Our data show that skin irradiation using either X-ray or proton microchannels maintain a higher cell viability and DNA integrity compared to a homogeneous irradiation, and thus might improve normal tissue protection after radiation therapy..

.

Stereotactic Synchrotron Radiotherapy (SSRT) and Microbeam Radiation Therapy (MRT) are both novel approaches to treat brain tumor and potentially other tumors using synchrotron radiation. Although the techniques differ by their principles, SSRT and MRT share certain common aspects with the possibility of combining their advantages in the future. For MRT, the technique uses highly collimated, quasi-parallel arrays of X-ray microbeams between 50 and 600 keV. Important features of highly brilliant Synchrotron sources are a very small beam divergence and an extremely high dose rate. The minimal beam divergence allows the insertion of so called Multi Slit Collimators (MSC) to produce spatially fractionated beams of typically ∼25-75 micron-wide microplanar beams separated by wider (100-400 microns center-to-center(ctc)) spaces with a very sharp penumbra. Peak entrance doses of several hundreds of Gy are extremely well tolerated by normal tissues and at the same time provide a higher therapeutic index for various tumor models in rodents. The hypothesis of a selective radio-vulnerability of the tumor vasculature versus normal blood vessels by MRT was recently more solidified. SSRT (Synchrotron Stereotactic Radiotherapy) is based on a local drug uptake of high-Z elements in tumors followed by stereotactic irradiation with 80 keV photons to enhance the dose deposition only within the tumor. With SSRT already in its clinical trial stage at the ESRF, most medical physics problems are already solved and the implemented solutions are briefly described, while the medical physics aspects in MRT will be discussed in more detail in this paper..

PURPOSE: The purpose of the work reported here was to investigate the influence of sub-millimeter size heterogeneities on the degradation of the distal edges of proton beams and to validate Monte Carlo (MC) methods' ability to correctly predict such degradation. METHODS: A custom-designed high-resolution plastic phantom approximating highly heterogeneous, lung-like structures was employed in measurements and in Monte Carlo simulations to evaluate the degradation of proton Bragg curves penetrating heterogeneous media. RESULTS: Significant differences in distal falloff widths and in peak dose values were observed in the measured and the Monte Carlo simulated curves compared to pristine proton Bragg curves. Furthermore, differences between simulations of beams penetrating CT images of the phantom did not agree well with the corresponding experimental differences. The distal falloff widths in CT image-based geometries were underestimated by up to 0.2 cm in water (corresponding to 0.8-1.4 cm in lung tissue), and the peak dose values of pristine proton beams were overestimated by as much as ˜35% compared to measured curves or depth-dose curves simulated on the basis of true geometry. The authors demonstrate that these discrepancies were caused by the limited spatial resolution of CT images that served as a basis for dose calculations and lead to underestimation of the impact of the fine structure of tissue heterogeneities. A convolution model was successfully applied to mitigate the underestimation. CONCLUSIONS: The results of this study justify further development of models to better represent heterogeneity effects in soft-tissue geometries, such as lung, and to correct systematic underestimation of the degradation of the distal edge of proton doses..

BACKGROUND AND PURPOSE: The latency of a multileaf collimator (MLC) tracking system used to overcome respiratory motion causes misalignment of the treatment beam with respect to the gross tumour volume, which may result in reduced target coverage. This study investigates the magnitude of this effect. MATERIAL AND METHODS: Simulated superior-inferior breathing motion was used to construct histograms of isocentre offset with respect to the gross tumour volume (GTV) for a variety of tracking latencies. Dose distributions for conformal volumetric modulated arc therapy (VMAT) arcs were then calculated at a range of offsets and summed according to these displacement histograms. The results were verified by delivering the plans to a Delta(4) phantom on a motion platform. RESULTS: In the absence of an internal target margin, a tracking latency of 150 ms reduces the GTV D95% by approximately 2%. With a margin of 2 mm, the same drop in dose occurs for a tracking latency of 450 ms. Lung V(13Gy) is unaffected by a range of latencies. These results are supported by the phantom measurements. CONCLUSIONS: Assuming that internal motion can be modelled by a rigid translation of the patient, MLC tracking of conformal VMAT can be effectively accomplished in the absence of an internal target margin for substantial breathing motion (4 s period and 20 mm peak-peak amplitude) so long as the system latency is less than 150 ms..

PURPOSE: The aim of this work was to compare and validate various computed tomography (CT) number calibration techniques with respect to cone beam CT (CBCT) dose calculation accuracy. METHODS: CBCT dose calculation accuracy was assessed for pelvic, lung, and head and neck (H&N) treatment sites for two approaches: (1) physics-based scatter correction methods (CBCTr); (2) density override approaches including assigning water density to the entire CBCT (W), assignment of either water or bone density (WB), and assignment of either water or lung density (WL). Methods for CBCT density assignment within a commercially available treatment planning system (RSauto), where CBCT voxels are binned into six density levels, were assessed and validated. Dose-difference maps and dose-volume statistics were used to compare the CBCT dose distributions with the ground truth of a planning CT acquired the same day as the CBCT. RESULTS: For pelvic cases, all CTN calibration methods resulted in average dose-volume deviations below 1.5 %. RSauto provided larger than average errors for pelvic treatments for patients with large amounts of adipose tissue. For H&N cases, all CTN calibration methods resulted in average dose-volume differences below 1.0 % with CBCTr (0.5 %) and RSauto (0.6 %) performing best. For lung cases, WL and RSauto methods generated dose distributions most similar to the ground truth. CONCLUSION: The RSauto density override approach is an attractive option for CTN adjustments for a variety of anatomical sites. RSauto methods were validated, resulting in dose calculations that were consistent with those calculated on diagnostic-quality CT images, for CBCT images acquired of the lung, for patients receiving pelvic RT in cases without excess adipose tissue, and for H&N cases..

PURPOSE: Real-time, markerless localization of lung tumors with kV imaging is often inhibited by ribs obscuring the tumor and poor soft-tissue contrast. This study investigates the use of dual-energy imaging, which can generate radiographs with reduced bone visibility, to enhance automated lung tumor tracking for real-time adaptive radiotherapy. METHODS: kV images of an anthropomorphic breathing chest phantom were experimentally acquired and radiographs of actual lung cancer patients were Monte-Carlo-simulated at three imaging settings: low-energy (70 kVp, 1.5 mAs), high-energy (140 kVp, 2.5 mAs, 1 mm additional tin filtration), and clinical (120 kVp, 0.25 mAs). Regular dual-energy images were calculated by weighted logarithmic subtraction of high- and low-energy images and filter-free dual-energy images were generated from clinical and low-energy radiographs. The weighting factor to calculate the dual-energy images was determined by means of a novel objective score. The usefulness of dual-energy imaging for real-time tracking with an automated template matching algorithm was investigated. RESULTS: Regular dual-energy imaging was able to increase tracking accuracy in left-right images of the anthropomorphic phantom as well as in 7 out of 24 investigated patient cases. Tracking accuracy remained comparable in three cases and decreased in five cases. Filter-free dual-energy imaging was only able to increase accuracy in 2 out of 24 cases. In four cases no change in accuracy was observed and tracking accuracy worsened in nine cases. In 9 out of 24 cases, it was not possible to define a tracking template due to poor soft-tissue contrast regardless of input images. The mean localization errors using clinical, regular dual-energy, and filter-free dual-energy radiographs were 3.85, 3.32, and 5.24 mm, respectively. Tracking success was dependent on tumor position, tumor size, imaging beam angle, and patient size. CONCLUSIONS: This study has highlighted the influence of patient anatomy on the success rate of real-time markerless tumor tracking using dual-energy imaging. Additionally, the importance of the spectral separation of the imaging beams used to generate the dual-energy images has been shown..

BACKGROUND: With the development of more conformal and precise radiation techniques such as Intensity-Modulated Radiotherapy (IMRT), Stereotactic Body Radiotherapy (SBRT) and Image-Guided Radiotherapy (IGRT), patients with hepatic tumors could be treated with high local doses by sparing normal liver tissue. However, frequently occurring large HCC tumors are still a dosimetric challenge in spite of modern high sophisticated RT modalities. This interventional clinical study has been set up to evaluate the value of different fiducial markers, and to use the modern imaging methods for further treatment optimization using physical and informatics approaches. METHODS AND DESIGN: Surgically implanted radioopaque or electromagnetic markers are used to detect tumor local-ization during radiotherapy. The required markers for targeting and observation during RT can be implanted in a previously defined optimal position during the oncologically indicated operation. If there is no indication for a surgical resection or open biopsy, markers may be inserted into the liver or tumor tissue by using ultrasound-guidance. Primary study aim is the detection of the patients' anatomy at the time of RT by observation of the marker position during the indicated irradiation (IGRT). Secondary study aims comprise detection and recording of 3D liver and tumor motion during RT. Furthermore, the study will help to develop technical strategies and mechanisms based on the recorded information on organ motion to avoid inaccurate dose application resulting from fast organ motion and deformation. DISCUSSION: This is an open monocentric non-randomized, prospective study for the evaluation of organ motion using interstitial markers or implantable radiotransmitter. The trial will evaluate the full potential of different fiducial markers to further optimize treatment of moving targets, with a special focus on liver lesions..

PURPOSE: Microbeam Radiation Therapy (MRT), an alternative preclinical treatment strategy using spatially modulated synchrotron radiation on a micrometer scale, has the great potential to cure malignant tumors (e.g., brain tumors) while having low side effects on normal tissue. Dose measurement and calculation in MRT is challenging because of the spatial accuracy required and the arising high dose differences. Dose calculation with Monte Carlo simulations is time consuming and their accuracy is still a matter of debate. In particular, the influence of photon polarization has been discussed in the literature. Moreover, it is controversial whether a complete knowledge of phase space trajectories, i.e., the simulation of the machine from the wiggler to the collimator, is necessary in order to accurately calculate the dose. METHODS: With Monte Carlo simulations in the Geant4 toolkit, the authors investigate the influence of polarization on the dose distribution and the therapeutically important peak to valley dose ratios (PVDRs). Furthermore, the authors analyze in detail phase space information provided by Martínez-Rovira et al. ["Development and commissioning of a Monte Carlo photon model for the forthcoming clinical trials in microbeam radiation therapy," Med. Phys. 39(1), 119-131 (2012)] and examine its influence on peak and valley doses. A simple source model is developed using parallel beams and its applicability is shown in a semiadjoint Monte Carlo simulation. Results are compared to measurements and previously published data. RESULTS: Polarization has a significant influence on the scattered dose outside the microbeam field. In the radiation field, however, dose and PVDRs deduced from calculations without polarization and with polarization differ by less than 3%. The authors show that the key consequences from the phase space information for dose calculations are inhomogeneous primary photon flux, partial absorption due to inclined beam incidence outside the field center, increased beam width and center to center distance due to the beam propagation from the collimator to the phantom surface and imperfect absorption in the absorber material of the Multislit Collimator. These corrections have an effect of approximately 10% on the valley dose and suffice to describe doses in MRT within the measurement uncertainties of currently available dosimetry techniques. CONCLUSIONS: The source for the first clinical pet trials in MRT is characterized with respect to its phase space and the photon polarization. The results suggest the use of a presented simplified phase space model in dose calculations and hence pave the way for alternative and fast dose calculation algorithms. They also show that the polarization is of minor importance for the clinical important peak and valley doses inside the microbeam field..

PURPOSE: To provide a proof of concept study for IMRT treatment planning through interactive dose shaping (IDS) by utilising the respective tools to create IMRT treatment plans for six prostate patients. METHODS: The IDS planning paradigm aims to perform interactive local dose adaptations of an IMRT plan without compromising already established valuable dose features in real-time. Various IDS tools are available in our in-house treatment planning software Dynaplan and were utilised to create IMRT treatment plans for six patients with an adeno-carcinoma of the prostate. The sequenced IDS treatment plans were compared to conventionally optimised clinically approved plans (9 beams, co-planar). The starting point consisted of open fields. The IDS tools were utilised to sculpt dose out of the rectum and bladder. For each patient, several IDS plans were created, with different trade-offs between organ sparing and target coverage. The reference dose distributions were imported into Dynaplan. For each patient, the IDS treatment plan with a similar or better trade-off between target coverage and OAR sparing was selected for plan evaluation, guided by a physician. Pencil beam dose calculation was performed on a grid with a voxel size of 1.95×1.95×2.0 mm(3) . D98%, D2%, mean dose and dose-volume indicators as specified by Quantec were calculated for plan evaluation. RESULTS: It was possible to utilise the software prototype to generate treatment plans for prostate patient geometries in 15-45 minutes. Individual local dose adaptations could be performed in less than one second. The average differences compared to the reference plans were for the mean dose: 0.0 Gy (boost) and 1.2 Gy (CTV), for D98%: -1.1 Gy and for D2%: 1.1 Gy (both target volumes). The dose-volume quality indicators were well below the Quantec constraints. CONCLUSION: Real-time treatment planning utilising IDS is feasible and has the potential to be implemented clinically. Research at The Institute of Cancer Research is supported by Cancer Research UK under Programme C46/A10588..

PURPOSE: Tumor hypoxia is correlated with treatment failure. To date, there are no published studies investigating hypoxia in non-small cell lung cancer (NSCLC) patients undergoing SBRT. We aim to use 18F-fluoromisonidazole (18F-FMISO) positron emission tomography (PET) imaging to non-invasively quantify the tumor hypoxic volume (HV), to elucidate potential roles of reoxygenation and tumor vascular response at high doses, and to identify an optimal prognostic imaging time-point. METHODS: SBRT-eligible patients with NSCLC tumors >1cm were prospectively enrolled in an IRB-approved study. Computed Tomography and dynamic PET images (0-120min, 150-180min, and 210-240min post-injection) were acquired using a Siemens BiographmCT PET/CT scanner. 18F-FMISO PET was performed on a single patient at 3 different time points around a single SBRT delivery of 18 Gy and HVs were compared using a tumor-to-blood ratio (TBR)>1.2 and rate of influx (Ki)>0.0015 (Patlak). RESULTS: Results from our first patient showed substantial temporal changes in HV following SBRT. Using a TBR threshold >1.2 and summed images 210-240min, the HVs were 19%, 31% and 13% of total tumor volume on day 0, 2 (48 hours post-SBRT), and 4 (96 hours post-SBRT). The absolute volume of hypoxia increased by nearly a factor of 2 after 18 Gy and then decreased almost to baseline 96 hours later. Selected imaging timepoints resulted in temporal changes in HV quantification obtained with TBR. Ki, calculated using 4-hour dynamic data, evaluated HVs as 22%, 75% and 21%, respectively. CONCLUSIONS: ith the results of only one patient, this novel pilot study highlights the potential benefit of 18F-FMISO PET imaging as results indicate substantial temporal changes in tumor HV post-SBRT. Analysis suggests that TBR is not a robust parameter for accurate HV quantification and heavily influenced by imaging timepoint selection. Kinetic modeling parameters are more sensitive and may aid in future treatment individualization based on patient-specific biological information..

Cone beam CT IGRT provides on-line anatomical data of the patient at the treatment couch while new hybrid MRI radiotherapy systems promise to provide this data during the actual radiation delivery itself. To exploit these data, and account for tissue rotations, deformations and tumor regression during radiation therapy delivery, fast, on-line IMRT (re)planning is required.On-line IMRT re-planning needs to re-generate a similar patient specific dose distribution as done in the pre-treatment planning but then for the new state of the anatomy. For this the pre-treatment information has to be propagated to the actual treatment. A fast dose engine is required, which may be, in the case of MRI based IGRT, Monte Carlo based in order to account for the magnetic field induced dose effects. Then robust re-planning methods or class solutions should be available that account efficiently for both rigid and non-rigid anatomical changes while preserving the patient specific pre-treatment dose considerations. Preferably this is all fully automatic, but also fast interactive re-planning is an option, especially in hypo-fractionated (boost) radiotherapy. LEARNING OBJECTIVES: 1. Preservation of pre-treatment prescriptions in the on-line IMRT re-planning 2. Fast re-planning techniques and speed limitations for dose engines 3. Pros and cons of interactive versus automatic re-planning..

PURPOSE: The dosimetric verification of treatment plans in helical tomotherapy usually is carried out via verification measurements. In this study, a method for independent dose calculation of tomotherapy treatment plans is presented, that uses a conventional treatment planning system with a pencil kernel dose calculation algorithm for generation of verification dose distributions based on patient CT data. METHODS: A pencil beam algorithm that directly uses measured beam data was configured for dose calculation for a tomotherapy machine. Tomotherapy treatment plans were converted into a format readable by an in-house treatment planning system by assigning each projection to one static treatment field and shifting the calculation isocenter for each field in order to account for the couch movement. The modulation of the fluence for each projection is read out of the delivery sinogram, and with the kernel-based dose calculation, this information can directly be used for dose calculation without the need for decomposition of the sinogram. The sinogram values are only corrected for leaf output and leaf latency. Using the converted treatment plans, dose was recalculated with the independent treatment planning system. Multiple treatment plans ranging from simple static fields to real patient treatment plans were calculated using the new approach and either compared to actual measurements or the 3D dose distribution calculated by the tomotherapy treatment planning system. In addition, dose-volume histograms were calculated for the patient plans. RESULTS: Except for minor deviations at the maximum field size, the pencil beam dose calculation for static beams agreed with measurements in a water tank within 2%/2 mm. A mean deviation to point dose measurements in the cheese phantom of 0.89% ± 0.81% was found for unmodulated helical plans. A mean voxel-based deviation of -0.67% ± 1.11% for all voxels in the respective high dose region (dose values >80%), and a mean local voxel-based deviation of -2.41% ± 0.75% for all voxels with dose values >20% were found for 11 modulated plans in the cheese phantom. Averaged over nine patient plans, the deviations amounted to -0.14% ± 1.97% (voxels >80%) and -0.95% ± 2.27% (>20%, local deviations). For a lung case, mean voxel-based deviations of more than 4% were found, while for all other patient plans, all mean voxel-based deviations were within ± 2.4%. CONCLUSIONS: The presented method is suitable for independent dose calculation for helical tomotherapy within the known limitations of the pencil beam algorithm. It can serve as verification of the primary dose calculation and thereby reduce the need for time-consuming measurements. By using the patient anatomy and generating full 3D dose data, and combined with measurements of additional machine parameters, it can substantially contribute to overall patient safety..

PURPOSE: To evaluate the performance of the Elekta Agility multileaf collimator (MLC) for dynamic real-time tumor tracking. METHODS: The authors have developed a new control software which interfaces to the Agility MLC to dynamically program the movement of individual leaves, the dynamic leaf guides (DLGs), and the Y collimators ("jaws") based on the actual target trajectory. A motion platform was used to perform dynamic tracking experiments with sinusoidal trajectories. The actual target positions reported by the motion platform at 20, 30, or 40 Hz were used as shift vectors for the MLC in beams-eye-view. The system latency of the MLC (i.e., the average latency comprising target device reporting latencies and MLC adjustment latency) and the geometric tracking accuracy were extracted from a sequence of MV portal images acquired during irradiation for the following treatment scenarios: leaf-only motion, jaw + leaf motion, and DLG + leaf motion. RESULTS: The portal imager measurements indicated a clear dependence of the system latency on the target position reporting frequency. Deducting the effect of the target frequency, the leaf adjustment latency was measured to be 38 ± 3 ms for a maximum target speed v of 13 mm/s. The jaw + leaf adjustment latency was 53 ± 3 at a similar speed. The system latency at a target position frequency of 30 Hz was in the range of 56-61 ms for the leaves (v ≤ 31 mm/s), 71-78 ms for the jaw + leaf motion (v ≤ 25 mm/s), and 58-72 ms for the DLG + leaf motion (v ≤ 59 mm/s). The tracking accuracy showed a similar dependency on the target position frequency and the maximum target speed. For the leaves, the root-mean-squared error (RMSE) was between 0.6-1.5 mm depending on the maximum target speed. For the jaw + leaf (DLG + leaf) motion, the RMSE was between 0.7-1.5 mm (1.9-3.4 mm). CONCLUSIONS: The authors have measured the latency and geometric accuracy of the Agility MLC, facilitating its future use for clinical tracking applications..

Electromagnetic (EM) tracking allows localization of small EM sensors in a magnetic field of known geometry without line-of-sight. However, this technique requires a cable connection to the tracked object. A wireless alternative based on magnetic fields, referred to as transponder tracking, has been proposed by several authors. Although most of the transponder tracking systems are still in an early stage of development and not ready for clinical use yet, Varian Medical Systems Inc. (Palo Alto, California, USA) presented the Calypso system for tumor tracking in radiation therapy which includes transponder technology. But it has not been used for computer-assisted interventions (CAI) in general or been assessed for accuracy in a standardized manner, so far. In this study, we apply a standardized assessment protocol presented by Hummel et al (2005 Med. Phys. 32 2371-9) to the Calypso system for the first time. The results show that transponder tracking with the Calypso system provides a precision and accuracy below 1 mm in ideal clinical environments, which is comparable with other EM tracking systems. Similar to other systems the tracking accuracy was affected by metallic distortion, which led to errors of up to 3.2 mm. The potential of the wireless transponder tracking technology for use in many future CAI applications can be regarded as extremely high..

PURPOSE: Various strategies to select beneficial beam ensembles for intensity-modulated radiation therapy (IMRT) have been suggested over the years. These beam angle selection (BAS) strategies are usually evaluated against reference configurations applying equispaced coplanar beams but they are not compared to one another. Here, the authors present a meta analysis of four BAS strategies that incorporates fluence optimization (FO) into BAS by combinatorial optimization (CO) and one BAS strategy that decouples FO from BAS, i.e., spherical cluster analysis (SCA). The underlying parameters of the BAS process are investigated and the dosimetric benefits of the BAS strategies are quantified. METHODS: For three intracranial lesions in proximity to organs at risk (OARs) the authors compare treatment plans applying equispaced coplanar beam ensembles with treatment plans using five different BAS strategies, i.e., four CO techniques and SCA, to establish coplanar and noncoplanar beam ensembles. Treatment plans applying 5, 7, 9, and 11 beams are investigated. For the CO strategies the authors perform BAS runs with a 5°, 10°, 15°, and 20° angular resolution, which corresponds to a minimum of 18 coplanar and a maximum of 1400 noncoplanar candidate beams. In total 272 treatment plans with different BAS settings are generated for every patient. The quality of the treatment plans is compared based on the protection of OARs yet integral dose, target homogeneity, and target conformity are also considered. RESULTS: It is possible to reduce the average mean and maximum doses in OARs by more than 4 Gy (1 Gy) with optimized noncoplanar (coplanar) beam ensembles found with BAS by CO or SCA. For BAS including FO by CO, the individual algorithm used and the angular resolution in the space of candidate beams does not have a crucial impact on the quality of the resulting treatment plans. All CO algorithms yield similar target conformity and slightly improved target homogeneity in comparison to equispaced coplanar setups. Furthermore, optimized coplanar (noncoplanar) beam ensembles enabled more than a 6% (5%) reduction of the integral dose. For SCA, however, integral dose was increased and target conformity was decreased in comparison to equispaced coplanar setups-especially for a small number of beams. CONCLUSION: Both BAS strategies incorporating FO by CO and independent BAS strategies excluding FO provide dose savings in OARs for optimized coplanar and especially noncoplanar beam ensembles; they should not be neglected in the clinic..

PURPOSE: Tumor hypoxia is correlated with treatment failure in radiotherapy. The hypoxic fraction of a tumor can be difficult to define using non-invasive imaging methods such as Positron Emission Tomography (PET). The purpose of this study is to use tracer kinetic modeling techniques to increase the accuracy and precision of hypoxic volume (HV) quantification with 18F-fluoromisonidazole (18F-FMISO) PET imaging. METHODS: Ten male BALB/c mice with EMT-6 tumors were selected for imaging. Anaesthetized animals were injected with 18F-FMISO and imaged using whole-body 120-min dynamic PET and Computed Tomography (CT). Data from dynamic 18F-FMISO scans for 3 mice were fit to a 2-compartment irreversible 3-rate constant model (K1, k2, k3) and a PATLAK model (Ki). Different thresholds were applied to the resultant Ki and k3 rate constants generated on a voxel-by-voxel basis to calculate the HV for each tumor. These kinetic-derived HVs were compared to HVs generated using conventional tumor-to-blood (TBR) ratios. RESULTS: Kinetically-derived HV (using the k3 tracer binding constant and the rate of tracer influx, Ki), calculated on a voxel-by-voxel basis, showed an increase in HV of up to 38% over the TBR method depending on the assumed threshold. Moreover, the PATLAK kinetic model using a Ki threshold may provide better delineation and calculation of the tumor HV than a 2-compartment pharmacokinetic model. CONCLUSION: For PET imaging with 18F-FMISO, pharmacokinetic modeling can improve the accuracy and precision of tumor hypoxia quantification over that obtained with tumor-to-blood ratios. In particular, the PATLAK method to define a tumor HV may be superior due to a reduction in noise and to better HV localization. Ongoing studies will correlate tumor HVs with Eppendorf electrode measurements and carbonic anhydrase IX staining of histologic sections. These studies could lead to improved techniques for identifying patients likely to benefit from therapies designed to overcome hypoxic radioresistance..

We have previously investigated the use of a conventional amorphous-silicon flat-panel detector (FPD) for intrafractional image guidance in the in-line geometry. In this configuration, the FPD is mounted between the patient and the treatment head, with the front of the FPD facing towards the patient. By geometrically separating signals from the diagnostic (kV) and treatment (MV) beams, it is possible to monitor the patient and treatment beam at the same time. In this study, we propose an FPD design based on existing technology with a 70% reduced up-stream areal density that is more suited to this new application. We have investigated our FPD model by means of a validated Monte Carlo simulation. Experimentally, simple rectangular fields were used to irradiate through the detector and observe the impact of removing detector components such as the support structure or the phosphor screen on the measured signal. The proposed FPD performs better than the conventional FPD: (i) attenuation of the MV beam is decreased by 60%; (ii) the MV signal is reduced by 20% for the primary MV field region which can avoid saturation of the FPD; and (iii) long range scatter from the MV into the kV region of the detector is greatly reduced..

Intensity modulated treatment plan optimization is a computationally expensive task. The feasibility of advanced applications in intensity modulated radiation therapy as every day treatment planning, frequent re-planning for adaptive radiation therapy and large-scale planning research severely depends on the runtime of the plan optimization implementation. Modern computational systems are built as parallel architectures to yield high performance. The use of GPUs, as one class of parallel systems, has become very popular in the field of medical physics. In contrast we utilize the multi-core central processing unit (CPU), which is the heart of every modern computer and does not have to be purchased additionally. In this work we present an ultra-fast, high precision implementation of the inverse plan optimization problem using a quasi-Newton method on pre-calculated dose influence data sets. We redefined the classical optimization algorithm to achieve a minimal runtime and high scalability on CPUs. Using the proposed methods in this work, a total plan optimization process can be carried out in only a few seconds on a low-cost CPU-based desktop computer at clinical resolution and quality. We have shown that our implementation uses the CPU hardware resources efficiently with runtimes comparable to GPU implementations, at lower costs..

PURPOSE: 4d cone-beam computed tomography (CBCT) scans are usually reconstructed by extracting the motion information from the 2d projections or an external surrogate signal, and binning the individual projections into multiple respiratory phases. In this "after-the-fact" binning approach, however, projections are unevenly distributed over respiratory phases resulting in inefficient utilization of imaging dose. To avoid excess dose in certain respiratory phases, and poor image quality due to a lack of projections in others, the authors have developed a novel 4d CBCT acquisition framework which actively triggers 2d projections based on the forward-predicted position of the tumor. METHODS: The forward-prediction of the tumor position was independently established using either (i) an electromagnetic (EM) tracking system based on implanted EM-transponders which act as a surrogate for the tumor position, or (ii) an external motion sensor measuring the chest-wall displacement and correlating this external motion to the phase-shifted diaphragm motion derived from the acquired images. In order to avoid EM-induced artifacts in the imaging detector, the authors devised a simple but effective "Faraday" shielding cage. The authors demonstrated the feasibility of their acquisition strategy by scanning an anthropomorphic lung phantom moving on 1d or 2d sinusoidal trajectories. RESULTS: With both tumor position devices, the authors were able to acquire 4d CBCTs free of motion blurring. For scans based on the EM tracking system, reconstruction artifacts stemming from the presence of the EM-array and the EM-transponders were greatly reduced using newly developed correction algorithms. By tuning the imaging frequency independently for each respiratory phase prior to acquisition, it was possible to harmonize the number of projections over respiratory phases. Depending on the breathing period (3.5 or 5 s) and the gantry rotation time (4 or 5 min), between ∼90 and 145 projections were acquired per respiratory phase resulting in a dose of ∼1.7-2.6 mGy per respiratory phase. Further dose savings and decreases in the scanning time are possible by acquiring only a subset of all respiratory phases, for example, peak-exhale and peak-inhale only scans. CONCLUSIONS: This study is the first experimental demonstration of a new 4d CBCT acquisition paradigm in which imaging dose is efficiently utilized by actively triggering only those projections that are desired for the reconstruction process..

PURPOSE: Pulse pileup occurring at high x-ray fluxes can severely degrade the energy resolution provided by a photon counting detector, which can represent a problem in spectroscopic CT when performing quantitative material discrimination tasks. As the effects of pileup can be most easily seen as a degradation of a detector's count rate linearity at high fluxes, it has been proposed previously to quantify and correct these nonlinearities. While this strategy has been applied successfully to materials without K-edges, it is currently unknown if this still prevails when using medical contrast agents. The purpose of this study is to close this gap. METHODS: A Medipix2MXR Hexa detector was employed, featuring a pixel pitch of 165 [micro sign]m and a 1 mm thick CdTe sensor. A phantom containing various concentrations of iodine and gadolinium contrast agents was subject to energy selective CT acquisitions, using a pulsable x-ray source operated at 70 kVp. These acquisitions were obtained at low and high photon fluxes of 1.0 × 10(6) and 1.3 × 10(7) mm(-2) s(-1), respectively. Nonlinearity corrections were applied to the high-flux projections and for each pixel separately. The results were compared to the results at low photon fluxes. RESULTS: At high fluxes, a general reduction of the reconstructed attenuation coefficients was observed, which could be partially recovered using the correction strategy applied. The spectroscopic separation of iodine from the phantom material, however, degraded with increasing x-ray flux. In contrast to this, gadolinium could still be discriminated almost as well as in the low flux case. CONCLUSIONS: Nonlinearity corrections applied to high flux measurements can help to recover attenuation coefficients normally obtained at low fluxes for low-Z materials, which do not exhibit an absorption edge in the relevant energy range. However, as a result of a significant change of the x-ray spectrum, the spectroscopic contrast normally observed for iodine was found to vanish with increasing x-ray flux. In other words, the authors' results indicate that nonlinearity corrections may be feasible only when the K-edge of interest is sufficiently high compared to the mean photon energy, and that spectroscopic CT at high x-ray fluxes may suffer from less limitations when using high-Z materials as contrast agents. A future study should aim to confirm these findings under clinical conditions..

PURPOSE: The advent of widespread kV-cone beam computer tomography in image guided radiation therapy and special therapeutic application of keV photons, e.g., in microbeam radiation therapy (MRT) require accurate and fast dose calculations for photon beams with energies between 40 and 200 keV. Multiple photon scattering originating from Compton scattering and the strong dependence of the photoelectric cross section on the atomic number of the interacting tissue render these dose calculations by far more challenging than the ones established for corresponding MeV beams. That is why so far developed analytical models of kV photon dose calculations fail to provide the required accuracy and one has to rely on time consuming Monte Carlo simulation techniques. METHODS: In this paper, the authors introduce a novel analytical approach for kV photon dose calculations with an accuracy that is almost comparable to the one of Monte Carlo simulations. First, analytical point dose and pencil beam kernels are derived for homogeneous media and compared to Monte Carlo simulations performed with the Geant4 toolkit. The dose contributions are systematically separated into contributions from the relevant orders of multiple photon scattering. Moreover, approximate scaling laws for the extension of the algorithm to inhomogeneous media are derived. RESULTS: The comparison of the analytically derived dose kernels in water showed an excellent agreement with the Monte Carlo method. Calculated values deviate less than 5% from Monte Carlo derived dose values, for doses above 1% of the maximum dose. The analytical structure of the kernels allows adaption to arbitrary materials and photon spectra in the given energy range of 40-200 keV. CONCLUSIONS: The presented analytical methods can be employed in a fast treatment planning system for MRT. In convolution based algorithms dose calculation times can be reduced to a few minutes..

This paper introduces the concept of analytical probabilistic modeling (APM) to quantify uncertainties in quality indicators of radiation therapy treatment plans. Assuming Gaussian probability densities over the input parameters of the treatment plan quality indicators, APM enables the calculation of the moments of the induced probability density over the treatment plan quality indicators by analytical integration. This paper focuses on analytical probabilistic dose calculation algorithms and the implications of APM regarding treatment planning. We derive closed-form expressions for the expectation value and the (co)variance of (1) intensity-modulated photon and proton dose distributions based on a pencil beam algorithm and (2) the standard quadratic objective function used in inverse planning. Complex correlation models of high dimensional uncertain input parameters and the different nature of random and systematic uncertainties in fractionated radiation therapy are explicitly incorporated into APM. APM variance calculations on phantom data sets show that the correlation assumptions and the difference of random and systematic uncertainties of the input parameters have a crucial impact on the uncertainty of the resulting dose. The derivations regarding the quadratic objective function show that APM has the potential to enable robust planning at almost the same computational cost like conventional inverse planning after a single probabilistic dose calculation. Beneficial applications of APM in the context of radiation therapy treatment planning are feasible..

PURPOSE: The authors have developed a system that monitors intrafractional target motion perpendicular to the treatment beam with the aid of radioopaque markers by means of separating kV image and megavoltage (MV) treatment field on a single flat-panel detector. METHODS: They equipped a research Siemens Artiste linear accelerator (linac) with a 41 × 41 cm(2) a-Si flat-panel detector underneath the treatment head. The in-line geometry allows kV (imaging) and MV (treatment) beams to share closely aligned beam axes. The kV source, usually mounted directly across from the flat-panel imager, was retracted toward the gantry by 13 cm to intentionally misalign kV and MV beams, resulting in a geometric separation of MV treatment field and kV image on the detector. Two consecutive images acquired within 140 ms (the first with MV-only and the second with kV and MV signal) were subtracted to generate a kV-only image. The images were then analyzed "online" with an automated threshold-based marker detection algorithm. They employed a 3D and a 4D phantom equipped with either a single radioopaque marker or three Calypso beacons to mimic respiratory motion. Measured room positions were either cross-referenced with a phantom voltage signal (single marker) or the Calypso system. The accuracy of the back-projection (from detected marker positions into room coordinates) was verified by a simulation study. RESULTS: A phantom study has demonstrated that the imaging framework is capable of automatically detecting marker positions and sending this information to the tracking tool at an update rate of 7.14 Hz. The system latency is 86.9 ± 1.0 ms for single marker detection in the absence of MV radiation. In the presence of a circular MV field of 5 cm diameter, the latency is 87.1 ± 0.9 ms. The total RMS position detection accuracy is 0.20 mm (without MV radiation) and 0.23 mm (with MV). CONCLUSIONS: Based on the evaluated motion patterns and MV field size, the positional accuracy and system latency indicate that this system is suitable for real-time adaptive applications..

PURPOSE: The goal of this work was to compare different methods of incorporating the additional dose of mega-voltage cone-beam CT (MV-CBCT) for image-guided intensity modulated radiotherapy (IMRT) of different tumor entities. MATERIAL AND METHODS: The absolute dose delivered by the MV-CBCT was calculated and considered by creating a scaled IMRT plan (scIMRT) by renormalizing the clinically approved plan (orgIMRT) so that the sum with the MV-CBCT dose yields the same prescribed dose. In the other case, a newly optimized plan (optIMRT) was generated by including the dose distribution of the MV-CBCT as pre-irradiation. Both plans were compared with the orgIMRT plan and a plan where the last fraction was skipped. RESULTS: No significant changes were observed regarding the 95% conformity index of the target volume. The mean dose of the organs at risk (OAR) increased by approx. 7% for the scIMRT plan and 5% for the optIMRT plan. A significant increase of the mean dose to the outline contour was observed, ranging from 3.1 ± 1.3% (optIMRT) to 13.0 ± 6.1% (scIMRT) for both methods over all entities. If the dose of daily MV-CBCT would have been ignored, the additional dose accumulated to nearly a whole treatment fraction with a general increase of approx. 10% to the OARs and approx. 4% to the target volume. CONCLUSION: Both methods of incorporating the additional MV-CBCT dose into the treatment plan are suitable for clinical practice. The dose distribution of the target volume could be achieved as conformal as with the orgIMRT plan, while only a moderate increase of mean dose to OAR was observed..

In this work, the image quality of a novel megavoltage cone-beam-computed tomography (CBCT) scanner is compared to three other image-guided radiation therapy devices by analysing images of different-sized quality assurance phantoms. The following devices are compared in terms of image uniformity, signal-to-noise ratio, contrast-to-noise ratio (CNR), electron density to HU conversion, presampling modulation transfer function (MTF(pre)) and combined spatial resolution and noise (Q-factor): (i) the Siemens Artiste kilovoltage (kV) (121 kV) CBCT device, (ii) the Artiste treatment beam line (TBL), 6 MV, (iii) the Tomotherapy (3.5 MV) fan-beam CT and (iv) Siemens' novel approach using a carbon target for a dedicated imaging beam line (IBL), 4.2 MV. Machine settings were selected to produce the same imaging dose for all devices. For a head phantom, IBL scans display CNR values 2.6 ± 0.3 times higher than for the TBL at the same dose level (for a CT-number range of -200 to -60 HU). kV CBCT, on the other hand, displays CNR values 7.9 ± 0.3 times higher than the IBL. There was no significant deviation in spatial resolution between IBL, TBL and Tomotherapy in terms of 50% and 10% MTF(pre). For kV CBCT, the MTF(pre) was significantly higher than those for other devices. In our Q-factor analysis, the IBL (14.6) scores higher than the TBL (7.9) and Tomotherapy (9.7) due to its lower noise level. The linearity of electron density to HU conversion is demonstrated for different-sized phantoms. Employing the IBL instead of the TBL significantly reduces the imaging dose by up to a factor of 5 at a constant image quality level, providing an immediate benefit for the patient..

We present a novel technical concept of a two-dimensional binary multileaf collimator (2D-bMLC) especially designed for fast dose delivery in rotational IMRT. The 2D-bMLC consists of individually controlled absorber channels, which are arranged side by side forming a 2D collimator aperture. In each channel three separate tungsten modules are arranged behind each other. To open and close an element, the central module is shifted between two positions. The purpose of this work is the presentation of the 2D-bMLC concept and its dosimetric evaluation. To determine the dosimetric properties, we designed a Monte Carlo model of an exemplary 2D-bMLC, consisting of 30 × 30 elements. A virtual source model of a flattening filter-free 7 MV linac was used to characterize the linac phase space. A primary radiation efficiency factor of 43% was calculated for the open 2D-bMLC by dividing the integral dose scored for a 2D-bMLC field by the integral dose scored for an open field with the same dimensions. The leakage calculated for the closed collimator was below 0.5%. Following the primary photon fluence distribution, the bixel intensity decreases with the distance of the element to the central axis of the treatment machine. From the collimator field's center toward its borders, the geometric bixel widths increase in a symmetric and predictable manner by up to 4%. The increase is explained by the specific design of the 2D-bMLC. Abutting element beams exhibit a slight tongue-and-groove effect if opened sequentially. This effect as well as the primary radiation efficiency is basically affected by the source size and the dimensions of the collimator elements. We successfully established and evaluated a dosimetric model of the 2D-bMLC. The results are promising, and we will therefore investigate on real patient plans, if the concept could be advantageous for fast rotational IMRT treatments..

Spectroscopic x-ray imaging by means of photon counting detectors has received growing interest during the past years. Critical to the image quality of such devices is their pixel pitch and the sensor material employed. This paper describes the imaging properties of Medipix2 MXR multi-chip assemblies bump bonded to 1 mm thick CdTe sensors. Two systems were investigated with pixel pitches of 110 and 165 μm, which are in the order of the mean free path lengths of the characteristic x-rays produced in their sensors. Peak widths were found to be almost constant across the energy range of 10 to 60 keV, with values of 2.3 and 2.2 keV (FWHM) for the two pixel pitches. The average number of pixels responding to a single incoming photon are about 1.85 and 1.45 at 60 keV, amounting to detective quantum efficiencies of 0.77 and 0.84 at a spatial frequency of zero. Energy selective CT acquisitions are presented, and the two pixel pitches' abilities to discriminate between iodine and gadolinium contrast agents are examined. It is shown that the choice of the pixel pitch translates into a minimum contrast agent concentration for which material discrimination is still possible. We finally investigate saturation effects at high x-ray fluxes and conclude with the finding that higher maximum count rates come at the cost of a reduced energy resolution..

Pencil beam algorithms are still considered as standard photon dose calculation methods in Radiotherapy treatment planning for many clinical applications. Despite their established role in radiotherapy planning their performance and clinical applicability has to be continuously adapted to evolving complex treatment techniques such as adaptive radiation therapy (ART). We herewith report on a new highly efficient version of a well-established pencil beam convolution algorithm which relies purely on measured input data. A method was developed that improves raytracing efficiency by exploiting the capability of modern CPU architecture for a runtime reduction. Since most of the current desktop computers provide more than one calculation unit we used symmetric multiprocessing extensively to parallelize the workload and thus decreasing the algorithmic runtime. To maximize the advantage of code parallelization, we present two implementation strategies - one for the dose calculation in inverse planning software, and one for traditional forward planning. As a result, we could achieve on a 16-core personal computer with AMD processors a superlinear speedup factor of approx. 18 for calculating the dose distribution of typical forward IMRT treatment plans..

PURPOSE: One limitation of accurate dose delivery in radiotherapy is intrafractional movement of the tumor or the entire patient which may lead to an underdosage of the target tissue or an overdosage of adjacent organs at risk. In order to compensate for this movement, different techniques have been developed. In this study the tracking performances of a multileaf collimator (MLC) tracking system and a robotic treatment couch tracking system were compared under equal conditions. METHODS: MLC tracking was performed using a tracking system based on the Siemens 160 MLC. A HexaPOD robotic treatment couch tracking system was also installed at the same linac. A programmable 4D motion stage was used to reproduce motion trajectories with different target phantoms. Motion localization of the target was provided by the 4D tracking system of Calypso Medical Inc. The gained positional data served as input signal for the control systems of the MLC and HexaPOD tracking systems attempting to compensate for the target motion. The geometric and dosimetric accuracy for the tracking of eight different respiratory motion trajectories was investigated for both systems. The dosimetric accuracy of both systems was also evaluated for the tracking of five prostate motion trajectories. RESULTS: For the respiratory motion the average root mean square error of all trajectories in y direction was reduced from 4.1 to 2.0 mm for MLC tracking and to 2.2 mm for HexaPOD tracking. In x direction it was reduced from 1.9 to 0.9 mm (MLC) and to 1.0 mm (HexaPOD). The average 2%/2 mm gamma pass rate for the respiratory motion trajectories was increased from 76.4% for no tracking to 89.8% and 95.3% for the MLC and the HexaPOD tracking systems, respectively. For the prostate motion trajectories the average 2%/2 mm gamma pass rate was 60.1% when no tracking was applied and was improved to 85.0% for MLC tracking and 95.3% for the HexaPOD tracking system. CONCLUSIONS: Both systems clearly increased the geometric and dosimetric accuracy during tracking of respiratory motion trajectories. Thereby, the geometric accuracy was increased almost equally by both systems, whereas the dosimetric accuracy of the HexaPOD tracking system was slightly better for all considered respiratory motion trajectories. Substantial improvement of the dosimetric accuracy was also observed during tracking of prostate motion trajectories during an intensity-modulated radiotherapy plan. Thereby, the HexaPOD tracking system showed better results than the MLC tracking..

We have previously developed a tumour tracking system, which adapts the aperture of a Siemens 160 MLC to electromagnetically monitored target motion. In this study, we exploit the use of a novel linac-mounted kilovoltage x-ray imaging system for MLC tracking. The unique in-line geometry of the imaging system allows the detection of target motion perpendicular to the treatment beam (i.e. the directions usually featuring steep dose gradients). We utilized the imaging system either alone or in combination with an external surrogate monitoring system. We equipped a Siemens ARTISTE linac with two flat panel detectors, one directly underneath the linac head for motion monitoring and the other underneath the patient couch for geometric tracking accuracy assessments. A programmable phantom with an embedded metal marker reproduced three patient breathing traces. For MLC tracking based on x-ray imaging alone, marker position was detected at a frame rate of 7.1 Hz. For the combined external and internal motion monitoring system, a total of only 85 x-ray images were acquired prior to or in between the delivery of ten segments of an IMRT beam. External motion was monitored with a potentiometer. A correlation model between external and internal motion was established. The real-time component of the MLC tracking procedure then relied solely on the correlation model estimations of internal motion based on the external signal. Geometric tracking accuracies were 0.6 mm (1.1 mm) and 1.8 mm (1.6 mm) in directions perpendicular and parallel to the leaf travel direction for the x-ray-only (the combined external and internal) motion monitoring system in spite of a total system latency of ~0.62 s (~0.51 s). Dosimetric accuracy for a highly modulated IMRT beam--assessed through radiographic film dosimetry--improved substantially when tracking was applied, but depended strongly on the respective geometric tracking accuracy. In conclusion, we have for the first time integrated MLC tracking with x-ray imaging in the in-line geometry and demonstrated highly accurate respiratory motion tracking..

Megavoltage imaging in image-guided radiotherapy usually suffers from the relatively small fraction of photons present in the energy range providing good soft tissue contrast, which corresponds to photon energies below 50 keV. As a consequence, comparatively high imaging doses are required to form low-noise images. Single-crystal targets can help to alleviate this problem through the emission of so-called coherent bremsstrahlung, amounting to a net increase in low-energy photons if the electron beam impinging on a target is carefully aligned with a major symmetry axis of the underlying crystal lattice. In this work, we present an overview of crystal materials and directions that appeared particularly promising during our studies of this phenomenon, based on theoretical considerations. We find that, while diamond targets perform best in absolute terms, those transition metals that exhibit a body-centred cubic lattice appear as interesting alternatives..

The beam angle selection (BAS) problem in intensity-modulated radiation therapy is often interpreted as a combinatorial optimization problem, i.e. finding the best combination of η beams in a discrete set of candidate beams. It is well established that the combinatorial BAS problem may be solved efficiently with metaheuristics such as simulated annealing or genetic algorithms. However, the underlying parameters of the optimization process, such as the inclusion of non-coplanar candidate beams, the angular resolution in the space of candidate beams, and the number of evaluated beam ensembles as well as the relative performance of different metaheuristics have not yet been systematically investigated. We study these open questions in a meta-analysis of four strategies for combinatorial optimization in order to provide a reference for future research related to the BAS problem in intensity-modulated radiation therapy treatment planning. We introduce a high-performance inverse planning engine for BAS. It performs a full fluence optimization for ≈3600 treatment plans per hour while handling up to 50 GB of dose influence data (≈1400 candidate beams). For three head and neck patients, we compare the relative performance of a genetic, a cross-entropy, a simulated annealing and a naive iterative algorithm. The selection of ensembles with 5, 7, 9 and 11 beams considering either only coplanar or all feasible candidate beams is studied for an angular resolution of 5°, 10°, 15° and 20° in the space of candidate beams. The impact of different convergence criteria is investigated in comparison to a fixed termination after the evaluation of 10 000 beam ensembles. In total, our simulations comprise a full fluence optimization for about 3000 000 treatment plans. All four combinatorial BAS strategies yield significant improvements of the objective function value and of the corresponding dose distributions compared to standard beam configurations with equi-spaced coplanar beams. The genetic and the cross-entropy algorithms showed faster convergence in the very beginning of the optimization but the simulated annealing algorithm eventually arrived at almost the same objective function values. These three strategies typically yield clinically equivalent treatment plans. The iterative algorithm showed the worst convergence properties. The choice of the termination criterion had a stronger influence on the performance of the simulated annealing algorithm than on the performance of the genetic and the cross-entropy algorithms. We advocate to terminate the optimization process after the evaluation of 1000 beam combinations without objective function decrease. For our simulations, this resulted in an average deviation of the objective function from the reference value after 10 000 evaluated beam ensembles of 0.5% for all metaheuristics. On average, there was only a minor improvement when increasing the angular resolution in the space of candidate beam angles from 20° to 5°. However, we observed significant improvements when considering non-coplanar candidate beams for challenging head and neck cases..

PURPOSE: To investigate in a simulation study whether using a variable relative biological effectiveness (RBE) in calculation and optimization of intensity-modulated proton therapy (IMPT) instead of using an RBE of 1.1 would result in significant changes in the RBE-weighted dose (RWD) distributions. METHODS AND MATERIALS: For 4 patients with head-and-neck tumors, three IMPT plans were prepared respectively. The first plan was physically optimized (IMPT-PO plan), and the RWD was calculated with a constant RBE of 1.1. Then the plan's RWD was recalculated (IMPT-R plan) using a variable RBE model taking into account the linear energy transfer (LET) and tissue-specific radiobiological parameters. The third IMPT plan was optimized using a biological optimization routine (IMPT-BO plan). RESULTS: Comparing the IMPT-PO and IMPT-R plans, we observed that the RWD in radioresistant tissues was more sensitive to the LET than in radiosensitive tissues. The IMPT-R plans were in general more inhomogeneous than the IMPT-PO plans. The differences of RWD distributions for all volumes between IMPT-PO and IMPT-BO plans complied with predefined dose-volume constraints. The average LET was significantly lower in IMPT-BO plans than in IMPT-R plans. CONCLUSION: In radioresistant normal tissues caution has to be used regarding the LET distribution because these are most sensitive to changes in the LET. Biological optimization of IMPT plans based on the organ-specific biological parameters and LET distributions is feasible..

PURPOSE: Dynamic multileaf collimator tracking represents a promising method for high-precision radiotherapy to moving tumors. In the present study, we report on the integration of electromagnetic real-time tumor position monitoring into a multileaf collimator-based tracking system. METHODS AND MATERIALS: The integrated system was characterized in terms of its geometric and radiologic accuracy. The former was assessed from portal images acquired during radiation delivery to a phantom in tracking mode. The tracking errors were calculated from the positions of the tracking field and of the phantom as extracted from the portal images. Radiologic accuracy was evaluated from film dosimetry performed for conformal and intensity-modulated radiotherapy applied to different phantoms moving on sinusoidal trajectories. A static radiation delivery to the nonmoving target served as a reference for the delivery to the moving phantom with and without tracking applied. RESULTS: Submillimeter tracking accuracy was observed for two-dimensional target motion despite the relatively large system latency of 500 ms. Film dosimetry yielded almost complete recovery of a circular dose distribution with tracking in two dimensions applied: 2%/2 mm gamma-failure rates could be reduced from 59.7% to 3.3%. For single-beam intensity-modulated radiotherapy delivery, accuracy was limited by the finite leaf width. A 2%/2 mm gamma-failure rate of 15.6% remained with tracking applied. CONCLUSION: The integrated system we have presented marks a major step toward the clinical implementation of high-precision dynamic multileaf collimator tracking. However, several challenges such as irregular motion traces or a thorough quality assurance still need to be addressed..

.

A method for simulating spot-scanned delivery to a moving tumour was developed which uses patient-specific image and plan data. The magnitude of interplay effects was investigated for two patient cases under different fractionation and respiratory motion variation scenarios. The use of volumetric rescanning for motion mitigation was also investigated. For different beam arrangements, interplay effects lead to severely distorted dose distributions for a single fraction delivery. Baseline shift variations for single fraction delivery reduced the dose to the clinical target volume (CTV) by up to 14.1 Gy. Fractionated delivery significantly reduced interplay effects; however, local overdosage of 12.3% compared to the statically delivered dose remained for breathing period variations. Variations of the tumour baseline position and respiratory period were found to have the largest influence on target inhomogeneity; these effects were reduced with fractionation. Volumetric rescanning improved the dose homogeneity. For the CTV, underdosage was improved by up to 34% in the CTV and overdosage to the lung was reduced by 6%. Our results confirm that rescanning potentially increases the dose homogeneity; however, it might not sufficiently compensate motion-induced dose distortions. Other motion mitigation techniques may be required to additionally treat lung tumours with scanned proton beams..

Prediction of respiratory motion is essential for real-time tracking of lung or liver tumours in radiotherapy to compensate for system latencies. This study compares the performance of respiratory motion prediction based on linear regression (LR), neural networks (NN), kernel density estimation (KDE) and support vector regression (SVR) for various sampling rates and system latencies ranging from 0.2 to 0.6 s. Root-mean-squared prediction errors are evaluated on 12 3D lung tumour motion traces acquired at 30 Hz during radiotherapy treatments. The effect of stationary predictor training versus continuous predictor retraining as well as full 3D motion processing versus independent coordinate-wise motion processing is investigated. Model parameter optimization is performed through a grid search in the model parameter space for each predictor and all considered latencies, sampling rates, training schemes and 3D data-processing modes. Comparison of the predictors is performed in the clinically applicable setting of patient-independent model parameters. The considered predictors roughly halve the prediction errors compared to using no prediction. When averaging over all sampling rates and latencies, prediction errors normalized to errors of using no prediction of 0.44, 0.46, 0.49 and 0.55 for NN, SVR, LR and KDE are observed. The small differences between the predictors emphasize the relative importance of adequate model parameter optimization compared to the actual prediction model selection. Thorough model parameter tuning is therefore essential for fair predictor comparisons..

.

.

PURPOSE: Advanced high quality radiation therapy techniques such as IMRT require an accurate delivery of precisely modulated radiation fields to the target volume. Interfractional and intrafractional motion of the patient's anatomy, however, may considerably deteriorate the accuracy of the delivered dose to the planned dose distributions. In order to compensate for these potential errors, a dynamic real-time capable MLC control system was designed. METHODS: The newly developed adaptive MLC control system contains specialized algorithms which are capable of continuous optimization and correction of the aperture of the MLC according to the motion of the target volume during the dose delivery. The algorithms calculate the new leaf positions based on target information provided online to the system. The algorithms were implemented in a dynamic target tracking control system designed for a Siemens 160 MLC. To assess the quality of the new target tracking system in terms of dosimetric accuracy, experiments with various types of motion patterns using different phantom setups were performed. The phantoms were equipped with radiochromic films placed between solid water slabs. Dosimetric results of exemplary deliveries to moving targets with and without dynamic MLC tracking applied were compared in terms of the gamma criterion to the reference dose delivered to a static phantom. RESULTS: Our measurements indicated that dose errors for clinically relevant two-dimensional target motion can be compensated by the new control system during the dose delivery of open fields. For a clinical IMRT dose distribution, the gamma success rate was increased from 19% to 77% using the new tracking system. Similar improvements were achieved for the delivery of a complete IMRT treatment fraction to a moving lung phantom. However, dosimetric accuracy was limited by the system's latency of 400 ms and the finite leaf width of 5 mm in the isocenter plane. CONCLUSIONS: Different experimental setups representing different target tracking scenarios proved that the tracking concept, the new algorithms and the dynamic control system make it possible to effectively compensate for dose errors due to target motion in real-time. These early results indicate that the method is suited to increasing the accuracy and the quality of the treatment delivery for the irradiation of moving tumors..

Coherent bremsstrahlung denotes the process of bremsstrahlung emission by electrons traversing a single crystal, causing prominent peaks in the resulting photon spectra at low energies. While this phenomenon has been known for decades, little attempt has been made to exploit its potential for megavoltage imaging, where the quality of images is affected by low contrast due to the lack of sufficient photons at the energy range suited for diagnostic purposes. We provide a theoretical foundation of coherent bremsstrahlung in the first-order Born approximation without confinement to high energies. Based on this theory, first evidence is given that diamond crystals are capable of boosting the amount of diagnostic photons by about 10-20%. It is shown that this behaviour is largely conserved for polychromatic electrons hitting thick targets, where multiple scattering dominates the energy distribution of the emitted photons..

An intuitive heuristic to establish beam configurations for intensity-modulated radiation therapy is introduced as an extension of beam ensemble selection strategies applying scalar scoring functions. It is validated by treatment plan comparisons for three intra-cranial, pancreas, and prostate cases each. Based on a patient specific matrix listing the radiological quality of candidate beam directions individually for every target voxel, a set of locally ideal beam angles is generated. The spherical distribution of locally ideal beam angles is characteristic for every treatment site and patient: ideal beam angles typically cluster around distinct orientations. We interpret the cluster centroids, which are identified with a spherical K-means algorithm, as irradiation angles of an intensity-modulated radiation therapy treatment plan. The fluence profiles are subsequently optimized during a conventional inverse planning process. The average computation time for the pre-optimization of a beam ensemble is six minutes on a state-of-the-art work station. The treatment planning study demonstrates the potential benefit of the proposed beam angle optimization strategy. For the three prostate cases under investigation, the standard treatment plans applying nine coplanar equi-spaced beams and treatment plans applying an optimized non-coplanar nine-beam ensemble yield clinically comparable dose distributions. For symmetric patient geometries, the dose distribution formed by nine equi-spaced coplanar beams cannot be improved significantly. For the three pancreas and intra-cranial cases under investigation, the optimized non-coplanar beam ensembles enable better sparing of organs at risk while guaranteeing equivalent target coverage. Beam angle optimization by spherical cluster analysis shows the biggest impact for target volumes located asymmetrically within the patient and close to organs at risk..

As an approach towards more biology-oriented treatment planning for external beam radiation therapy, we present the incorporation of local radiation damage models into three dimensional treatment planning. This allows effect based instead of dose based plan optimization which could potentially better match the biologically relevant tradeoff between target and normal tissues. In particular, our approach facilitates an effective comparison of different fractionation schemes. It is based on the linear quadratic model to describe the biological radiation effect. Effect based optimization was integrated into our inverse treatment planning software KonRad, and we demonstrate the resulting differences between conventional and biological treatment planning. Radiation damage can be analyzed both qualitatively and quantitatively in dependence of the fractionation scheme and tissue specific parameters in a three dimensional voxel based system. As an example the potential advantages as well as the associated risks of hypofractionation for prostate cancer are analyzed and visualized with the help of effective dose volume histograms. Our results suggest a very conservative view regarding alternative fractionation schemes since uncertainties in biological parameters are still too big to make reliable clinical predictions..

PURPOSE: Four-dimensional (4D) imaging is a key to motion-adapted radiotherapy of lung tumors. We evaluated in a ventilated ex vivo system how size and displacement of artificial pulmonary nodules are reproduced with helical 4D-CT, 4D-MRI, and linac-integrated cone beam CT (CBCT). METHODS AND MATERIALS: Four porcine lungs with 18 agarose nodules (mean diameters 1.3-1.9 cm), were ventilated inside a chest phantom at 8/min and subject to 4D-CT (collimation 24 x 1.2 mm, pitch 0.1, slice/increment 24 x 10(2)/1.5/0.8 mm, pitch 0.1, temporal resolution 0.5 s), 4D-MRI (echo-shared dynamic three-dimensional-flash; repetition/echo time 2.13/0.72 ms, voxel size 2.7 x 2.7 x 4.0 mm, temporal resolution 1.4 s) and linac-integrated 4D-CBCT (720 projections, 3-min rotation, temporal resolution approximately 1 s). Static CT without respiration served as control. Three observers recorded lesion size (RECIST-diameters x/y/z) and axial displacement. Interobserver- and interphase-variation coefficients (IO/IP VC) of measurements indicated reproducibility. RESULTS: Mean x/y/z lesion diameters in cm were equal on static and dynamic CT (1.88/1.87; 1.30/1.39; 1.71/1.73; p > 0.05), but appeared larger on MRI and CBCT (2.06/1.95 [p < 0.05 vs. CT]; 1.47/1.28 [MRI vs. CT/CBCT p < 0.05]; 1.86/1.83 [CT vs. CBCT p < 0.05]). Interobserver-VC for lesion sizes were 2.54-4.47% (CT), 2.29-4.48% (4D-CT); 5.44-6.22% (MRI) and 4.86-6.97% (CBCT). Interphase-VC for lesion sizes ranged from 2.28% (4D-CT) to 10.0% (CBCT). Mean displacement in cm decreased from static CT (1.65) to 4D-CT (1.40), CBCT (1.23) and MRI (1.16). CONCLUSIONS: Lesion sizes are exactly reproduced with 4D-CT but overestimated on 4D-MRI and CBCT with a larger variability due to limited temporal and spatial resolution. All 4D-modalities underestimate lesion displacement..

In recent experiments, quasi-monoenergetic and well-collimated very-high energy electron (VHEE) beams were obtained by laser-plasma accelerators. We investigate their potential use for radiation therapy. Monte Carlo simulations are used to study the influence of the experimental characteristics such as beam energy, energy spread and initial angular distribution on the dose distributions. It is found that magnetic focusing of the electron beam improves the lateral penumbra. The dosimetric properties of the laser-accelerated VHEE beams are implemented in our inverse treatment planning system for intensity-modulated treatments. The influence of the beam characteristics on the quality of a prostate treatment plan is evaluated. In comparison to a clinically approved 6 MV IMRT photon plan, a better target coverage is achieved. The quality of the sparing of organs at risk is found to be dependent on the depth. The bladder and rectum are better protected due to the sharp lateral penumbra at low depths, whereas the femoral heads receive a larger dose because of the large scattering amplitude at larger depths..

The purpose of this work is to investigate robust 4D optimization techniques which account for respiratory motion uncertainties. Two robust optimization techniques were applied to generate 4D optimized lung treatment plans. The probabilistic optimization approach minimizes the dose variance in the target volume while the worst case optimization minimizes a weighted combination of the nominal and worst case dose distributions which occur in the presence of respiratory motion variation. The two 4D optimization approaches were compared with a margin-based midventilation planning approach in five lung patients. Respiratory motion amplitude and baseline variations were quantified from tidal volume measurements during planning 4D CT acquisition. A similar target coverage was obtained for all three approaches, although the 4D optimization methods tended to be better at sparing the organs at risk. Both robust planning methods are suited for automatic determination of treatment plans which ensure target dose conformality under respiratory motion variations, while minimizing the dose burden of healthy lung tissue..

.

PURPOSE: Arc-modulated cone beam therapy (AMCBT) is a fast treatment technique deliverable in a single rotation with a conventional C-arm shaped linac. In this planning study, the authors assess the dosimetric properties of single-arc therapy in comparison to helical tomotherapy for three different tumor types. METHODS: Treatment plans for three patients with prostate carcinoma, three patients with anal cancer, and three patients with head and neck cancer were optimized for helical tomotherapy and AMCBT. The dosimetric comparison of the two techniques is based on physical quantities derived from dose-volume histograms. RESULTS: For prostate cancer, the quality of dose distributions calculated for AMCBT was of equal quality as that generated for tomotherapy with the additional benefits of a faster delivery and a lower integral dose. For highly complex geometries, the plan quality achievable with helical tomotherapy could not be achieved with arc-modulated cone beam therapy. CONCLUSIONS: Rotation therapy with a conventional linac in a single arc is capable to deliver a high and homogeneous dose to the target and spare organs at risk. Advantages of this technique are a fast treatment time and a lower integral dose in comparison to helical tomotherapy. For highly complex cases, e.g., with several target regions, the dose shaping capabilities of AMCBT are inferior to those of tomotherapy. However, treatment plans for AMCBT were also clinically acceptable..

It is still an unanswered question whether a relatively low dose of radiation to a large volume or a higher dose to a small volume produces the higher cancer incidence. This is of interest in view of modalities like IMRT or rotation therapy where high conformity to the target volume is achieved at the cost of a large volume of normal tissue exposed to radiation. Knowledge of the shape of the dose response for radiation-induced cancer is essential to answer the question of what risk of second cancer incidence is implied by which treatment modality. This study therefore models the dose response for radiation-induced second cancer after radiation therapy of which the exact mechanisms are still unknown. A second cancer risk estimation tool for treatment planning is presented which has the potential to be used for comparison of different treatment modalities, and risk is estimated on a voxel basis for different organs in two case studies. The presented phenomenological model summarises the impact of microscopic biological processes into effective parameters of mutation and cell sterilisation. In contrast to other models, the effective radiosensitivities of mutated and non-mutated cells are allowed to differ. Based on the number of mutated cells present after irradiation, the model is then linked to macroscopic incidence by summarising model parameters and modifying factors into natural cancer incidence and the dose response in the lower-dose region. It was found that all principal dose-response functions discussed in the literature can be derived from the model. However, from the investigation and due to scarcity of adequate data, rather vague statements about likelihood of dose-response functions can be made than a definite decision for one response. Based on the predicted model parameters, the linear response can probably be rejected using the dynamics described, but both a flattening response and a decrease appear likely, depending strongly on the effective cell sterilisation of the mutated cells. Thus insights could be gained into the impact of parameters describing the effective mutation or cell sterilisation of non-mutated as well as of mutated cells, which constitute precursors of cancer. The biggest drawbacks in the estimation of second cancer incidence remain the low statistical power of clinical studies on radiation induction of cancer and the inability to isolate the effect due to radiation alone - if the latter is possible at all. We conclude that at the present stage of knowledge, further investigations have to be carried out in order to really compare treatment modalities with respect to the second cancer risk they imply..

The problem of the enormous amount of scattered radiation in kV CBCT (kilo voltage cone beam computer tomography) is addressed. Scatter causes undesirable streak- and cup-artifacts and results in a quantitative inaccuracy of reconstructed CT numbers, so that an accurate dose calculation might be impossible. Image contrast is also significantly reduced. Therefore we checked whether an appropriate implementation of the fast iterative scatter correction algorithm we have developed for MV (mega voltage) CBCT reduces the scatter contribution in a kV CBCT as well. This scatter correction method is based on a superposition of pre-calculated Monte Carlo generated pencil beam scatter kernels. The algorithm requires only a system calibration by measuring homogeneous slab phantoms with known water-equivalent thicknesses. In this study we compare scatter corrected CBCT images of several phantoms to the fan beam CT images acquired with a reduced cone angle (a slice-thickness of 14 mm in the isocenter) at the same system. Additional measurements at a different CBCT system were made (different energy spectrum and phantom-to-detector distance) and a first order approach of a fast beam hardening correction will be introduced. The observed image quality of the scatter corrected CBCT images is comparable concerning resolution, noise and contrast-to-noise ratio to the images acquired in fan beam geometry. Compared to the CBCT without any corrections the contrast of the contrast-and-resolution phantom with scatter correction and additional beam hardening correction is improved by a factor of about 1.5. The reconstructed attenuation coefficients and the CT numbers of the scatter corrected CBCT images are close to the values of the images acquired in fan beam geometry for the most pronounced tissue types. Only for extreme dense tissue types like cortical bone we see a difference in CT numbers of 5.2%, which can be improved to 4.4% with the additional beam hardening correction. Cupping is reduced from 20% to 4% with scatter correction and 3% with an additional beam hardening correction. After 3 iterations (small phantoms) and 6 to 7 iterations (large phantoms) the algorithm converges. Therefore the algorithm is very fast, that means 1.3 seconds per projection for 3 iterations on a standard PC..

PURPOSE: In radiotherapy with hadrons, it is anticipated that carbon ions are superior to protons, mainly because of their biological properties: the relative biological effectiveness (RBE) for carbon ions is supposedly higher in the target than in the surrounding normal tissue, leading to a therapeutic advantage over protons. The purpose of this report is to investigate this effect by using biological model calculations. METHODS AND MATERIALS: We compared spread-out Bragg peaks for protons and carbon ions by using physical and biological optimization. The RBE for protons and carbon ions was calculated according to published biological models. These models predict increased RBE values in regions of high linear energy transfer (LET) and an inverse dependency of the RBE on dose. RESULTS: For pure physical optimization, protons yield a better dose distribution along the central axis. In biologically optimized plans, RBE variations for protons were relatively small. For carbon ions, high RBE values were found in the high-LET target region, as well as in the low-dose region outside the target. This means that the LET dependency and dose dependency of the RBE can cancel each other. We show this for radioresistant tissues treated with two opposing beams, for which the predicted carbon RBE within the target volume was lower than outside. CONCLUSIONS: For tissue parameters used in this study, the model used does not predict a biologic advantage of carbon ions. More reliable model parameters and clinical trials are necessary to explore the true potential of radiotherapy with carbon ions..

PURPOSE: A study to quantitatively compare the image quality of four different image guided radiotherapy (IGRT) devices based on phantom measurements with respect to the additional dose delivered to the patient. METHODS: Images of three different head-sized phantoms (diameter 16-18 cm) were acquired with the following four IGRT-CT solutions: (i) the Siemens Primatom single slice fan beam computed tomography (CT) scanner with an acceleration voltage of 130 kV, (ii) a Tomotherapy HI-ART II unit using a fan beam scanner with an energy of 3.5 MeV and (iii) the Siemens Artíste prototype, providing the possibility to perform kV (121 kV) and MV (6 MV) cone beam (CB) CTs. For each device three scan protocols (named low, normal, high) were selected to yield the same weighted computed tomography dose index (CTDI(w)). Based on the individual inserts of the different phantoms the image quality achieved with each device at a certain dose level was characterized in terms of homogeneity, spatial resolution, signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR) and electron density-to-CT-number conversion. RESULTS: Based on the current findings for head-sized phantoms all devices show an electron density-to-CT-number conversion almost independent of the imaging parameters and hence can be suited for treatment planning purposes. The evaluation of the image quality, however, points out clear differences due to the different energies and geometries. The Primatom standard CT scanner shows throughout the best performance, especially for soft tissue contrast and spatial resolution with low imaging doses. Reasonable soft tissue contrast can be obtained with slightly higher doses compared to the CT scanner with the kVCB and the Tomotherapy unit. In order to get similar results with the MVCB system a much higher dose needs to be applied to the patient. CONCLUSION: Considering the entire investigations, especially in terms of contrast and spatial resolution, a rough tendency for decreasing image quality can be given: Primatom, Artíste prototype kVCB, Tomotherapy, Artíste prototype MVCB..

New technical developments constantly aim at improving the outcome of radiation therapy. With the use of a computer-controlled multileaf collimator (MLC), the quality of the treatment and the efficiency in patient throughput is significantly increased. New MLC designs aim to further enhance the advantages. In this article, we present the first detailed experimental investigation of the new 160 MLC, Siemens Medical Solutions. The assessment included the experimental investigation of typical MLC characteristics such as leakage, tongue-and-groove effect, penumbra, leaf speed, and leaf positioning accuracy with a 6 MV treatment beam. The leakage is remarkably low with an average of 0.37% due to a new design principle of slightly tilted leaves instead of the common tongue-and-groove design. But due to the tilt, the triangular tongue-and-groove effect occurs. Its magnitude of approximately 19% is similar to the dose defect measured for MLCs with the common tongue-and-groove design. The average longitudinal penumbra measured at depth d(max) = 15 mm with standard 100 x 100 mm2 fields is 4.1 +/- 0.5 mm for the central range and increases to 4.9 +/- 1.3 mm for the entire field range of 400 x 400 mm2. The increase is partly due to the single-focusing design and the large distance between the MLC and the isocenter enabling a large patient clearance. Regarding the leaf speed, different velocity tests were performed. The positions of the moving leaves were continuously recorded with the kilovoltage-imaging panel. The maximum leaf velocities measured were 42.9 +/- 0.6 mm/s. In addition, several typical intensity-modulated radiation therapy treatments were performed and the delivery times compared to the Siemens OPTIFOCUS MLC. An average decrease of 11% in delivery time was observed. The experimental results presented in this article indicate that the dosimetric characteristics of the 160 MLC are capable of improving the quality of dose delivery with respect to precision and dose conformity..

In proton scanning systems that employ active energy variation for depth modulation, a switch of the particle energy might typically require 1-2 s. For plans comprising many energy slices, these seconds could sum up to a non-negligible fraction of the total treatment duration. We have applied the Nyquist-Shannon sampling theorem to determine an efficient spatial arrangement of Bragg peaks in a target volume. This pre-determined schedule of increasing energy spacing with higher energy allows us to reduce the number of used energy slices without compromising the physical dosimetric quality of a plan. Our results suggest that the advantage of such a simple implementation would be especially significant for larger, deep-seated tumors such as the prostate; the number of energy slices was cut by a factor of 2-6..

An optimal plan in modern treatment planning tools is found through the use of an iterative optimization algorithm, which deals with a high amount of patient-related data and number of treatment parameters to be optimized. Thus, calculating a good plan is a very time-consuming process which limits the application for patients in clinics and for research activities aiming for more accuracy. A common technique to handle the vast amount of radiation dose data is the concept of the influence matrix (DIJ), which stores the dose contribution of each bixel to the patient in the main memory of the computer. This study revealed that a bottleneck for the optimization time arises from the data transfer of the dose data between the memory and the CPU. In this note, we introduce a new method which speeds up the data transportation from stored dose data to the CPU. As an example we used the DIJ approach as is implemented in our treatment planning tool KonRad, developed at the German Cancer Research Center (DKFZ) in Heidelberg. A data cycle reordering method is proposed to take the advantage of modern memory hardware. This induces a minimal eviction policy which results in a memory behaviour exhibiting a 2.6 times faster algorithm compared to the naive implementation. Although our method is described for the DIJ approach implemented in KonRad, we believe that any other planning tool which uses a similar approach to store the dose data will also benefit from the described methods..

BACKGROUND: To determine whether proton radiotherapy has clinical advantages over photon radiotherapy, we modeled the dose characteristics of both to critical normal tissue volumes using data from patients with four types of childhood brain tumors. PROCEDURES: Three-dimensional imaging and treatment planning data, including targeted tumor and normal tissues contours, were acquired for 40 patients, 10 each with optic pathway glioma (OPG), craniopharyngioma (CR), infratentorial ependymoma (EP), or medulloblastoma (MB). Dose-volume data were collected for the entire brain, temporal lobes, cochlea, and hypothalamus from each patient. The data were averaged and compared based on treatment modality (protons vs. photons) using dose-cognitive effects models. Outcomes were estimated over 5 years. RESULTS: Relatively small critical normal tissue volumes such as the cochlea and hypothalamus may be spared from radiation exposure when not adjacent to the primary tumor volume. Larger normal tissue volumes such as the supratentorial brain or temporal lobes receive less of the low and intermediate doses. When applied to longitudinal models of radiation dose-cognitive effects, these differences resulted in clinically significant higher IQ scores for patients with MB and CR and academic reading scores in patients with OPG. Extreme differences between proton and photon dose distributions precluded meaningful comparison of protons and photons for patients with EP. CONCLUSIONS: Differences in the overall dose distributions, as indicated by modeling changes in cognitive function, showed that a reduction in the lower-dose volumes or mean dose would have long-term, clinical advantages for children with MB, CR, and OPG..

Intrafractional organ motion remains a source of error in conformal radiotherapy of dynamic targets such as tumours of the lung or of the prostate. The purpose of this work was to devise a method for the continuous and routine measurement of intrafractional organ motion. The method consists of a combination of an electromagnetic (EM), internal marker-based tracking system with the on-board kilovoltage x-ray imaging system of a modern treatment machine. The EM system continuously tracks the target, while x-ray images can be acquired simultaneously if demand arises. An image processing algorithm has been developed to automatically localize and track the EM markers in the x-ray images. We have demonstrated simultaneous target tracking using the EM system and x-ray imaging of a mobile target inside a programmable thorax phantom. The target motion was very well reproduced by both systems. The comparability of the target locations reported by both systems was established (better than 0.25 mm up to target velocities of 3 cm s(-1)). One immediate use of the synchronized system was shown: the generation of a 4D cone beam computed tomography data set using the EM system for the measurement of motion. In conclusion, we have developed a system for the routine measurement of intrafractional motion that continuously provides the 3D position of the target with the ability to acquire images of the treatment field only when needed, thereby eliminating avoidable imaging dose to the patient..

Laser-induced particle accelerators have been recognized as a potential proton source for radiotherapeutic applications in recent years. However, there are still major difficulties--especially regarding the resulting proton spectra--to overcome for a successful application in the clinic. Here we elaborate on the physics of double-layer targets to propose a tentative 'optical gantry' setup. The spectral requirements for a quality dose deposition of the fast protons are estimated. Plasma simulations of the one-dimensional expansion of microstructured targets are performed according to various target dimensions, rear proton densities and substrate masses. Subsequently, the dependence of the resulting proton spectra on these parameters is evaluated and compared to previously published analytical considerations. Quasi-monoenergetic proton beams, which would be suitable for high-quality dose delivery, could be achieved from pure proton targets if one were able to select out the rear layer of those targets. However, much more realistic heavy substrate layered targets are not able to preserve this high spectral standard, partly due to a second Coulomb-expansion in the center-of-mass frame of the fast protons. This expansion can be mitigated by a reduction of the total positive charge in the rear layer, resulting in a comparable spectral quality as the previous target types. In conclusion, the promising spectral results as well as an estimation of the total number of fast protons which can be expected from such a setup, suggest that the introduction of laser-based proton accelerators into the clinic might be possible in the future..

The sharp dose gradients which are possible in intensity modulated proton therapy (IMPT) not only offer the possibility of generating excellent target coverage while sparing neighbouring organs at risk, but can also lead to treatment plans which are very sensitive to uncertainties in treatment variables such as the range of individual Bragg peaks. We developed a method to account for uncertainties of treatment variables in the optimization based on a worst case dose distribution. The worst case dose distribution is calculated using several possible realizations of the uncertainties. This information is used by the objective function of the inverse treatment planning system to generate treatment plans which are acceptable under all considered realizations of the uncertainties. The worst case optimization method was implemented in our in-house treatment planning software KonRad in order to demonstrate the usefulness of this approach for clinical cases. In this paper, we investigated range uncertainties, setup uncertainties and a combination of both uncertainties. Using our method the sensitivity of the resulting treatment plans to these uncertainties is considerably reduced..

BACKGROUND: The aim of this treatment planning study was to investigate the potential advantages of intensity-modulated (IM) proton therapy (IMPT) compared with IM photon therapy (IMRT) in nasopharyngeal carcinoma (NPC). METHODS: Eight NPC patients were chosen. The dose prescriptions in cobalt Gray equivalent (GyE) for gross tumor volumes of the primary tumor (GTV-T), planning target volumes of GTV-T and metastatic (PTV-TN) and elective (PTV-N) lymph node stations were 72.6 GyE, 66 GyE, and 52.8 GyE, respectively. For each patient, nine coplanar fields IMRT with step-and-shoot technique and 3D spot-scanned three coplanar fields IMPT plans were prepared. Both modalities were planned in 33 fractions to be delivered with a simultaneous integrated boost technique. All plans were prepared and optimized by using the research version of the inverse treatment planning system KonRad (DKFZ, Heidelberg). RESULTS: Both treatment techniques were equal in terms of averaged mean dose to target volumes. IMPT plans significantly improved the tumor coverage and conformation (P < 0.05) and they reduced the averaged mean dose to several organs at risk (OARs) by a factor of 2-3. The low-to-medium dose volumes (0.33-13.2 GyE) were more than doubled by IMRT plans. CONCLUSION: In radiotherapy of NPC patients, three-field IMPT has greater potential than nine-field IMRT with respect to tumor coverage and reduction of the integral dose to OARs and non-specific normal tissues. The practicality of IMPT in NPC deserves further exploration when this technique becomes available on wider clinical scale..

In intensity modulated treatment techniques, the modulation of each treatment field is obtained using an optimization algorithm. Multiple optimization algorithms have been proposed in the literature, e.g. steepest descent, conjugate gradient, quasi-Newton methods to name a few. The standard optimization algorithm in our in-house inverse planning tool KonRad is a quasi-Newton algorithm. Although this algorithm yields good results, it also has some drawbacks. Thus we implemented an improved optimization algorithm based on the limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) routine. In this paper the improved optimization algorithm is described. To compare the two algorithms, several treatment plans are optimized using both algorithms. This included photon (IMRT) as well as proton (IMPT) intensity modulated therapy treatment plans. To present the results in a larger context the widely used conjugate gradient algorithm was also included into this comparison. On average, the improved optimization algorithm was six times faster to reach the same objective function value. However, it resulted not only in an acceleration of the optimization. Due to the faster convergence, the improved optimization algorithm usually terminates the optimization process at a lower objective function value. The average of the observed improvement in the objective function value was 37%. This improvement is clearly visible in the corresponding dose-volume-histograms. The benefit of the improved optimization algorithm is particularly pronounced in proton therapy plans. The conjugate gradient algorithm ranked in between the other two algorithms with an average speedup factor of two and an average improvement of the objective function value of 30%..

In radiotherapy with scanned particle beams, tissue heterogeneities lateral to the beam direction are problematic in two ways: they pose a challenge to dose calculation algorithms, and they lead to a high sensitivity to setup errors. In order to quantify and avoid these problems, a heterogeneity number H(i) as a method to quantify lateral tissue heterogeneities of single beam spot i is introduced. To evaluate this new concept, two kinds of potential errors were investigated for single beam spots: First, the dose calculation error has been obtained by comparing the dose distribution computed by a simple pencil beam algorithm to more accurate Monte Carlo simulations. The resulting error is clearly correlated with H(i). Second, the analysis of the sensitivity to setup errors of single beam spots also showed a dependence on H(i). From this data it is concluded that H(i) can be used as a criterion to assess the risks of a compromised delivered dose due to lateral tissue heterogeneities. Furthermore, a method how to incorporate this information into the inverse planning process for intensity modulated proton therapy is presented. By suppressing beam spots with a high value of H(i), the unfavorable impact of lateral tissue heterogeneities can be reduced, leading to treatment plans which are more robust to dose calculation errors of the pencil beam algorithm. Additional possibilities to use the information of H(i) are outlined in the discussion..

A three-dimensional (3D) intensity modulated proton therapy treatment plan to be delivered by magnetic scanning may comprise thousands of discrete beam positions. This research presents the minimization of the total scan path length by application of a fast simulated annealing (FSA) optimization algorithm. Treatment plans for clinical prostate and head and neck cases were sequenced for continuous raster scanning in two ways, and the resulting scan path lengths were compared: (1) A simple back-and-forth, top-to-bottom (zigzag) succession, and (2) an optimized path produced as a solution of the FSA algorithm. Using a first approximation of the scanning dynamics, the delivery times for the scan sequences before and after path optimization were calculated for comparison. In these clinical examples, the FSA optimization shortened the total scan path length for the 3D target volumes by approximately 13%-56%. The number of extraneous spilled particles was correspondingly reduced by about 13%-54% due to the more efficient scanning maps that eliminated multiple crossings through regions of zero fluence. The relative decrease in delivery time due to path length minimization was estimated to be less than 1%, due to both a high scanning speed and time requirements that could not be altered by optimization (e.g., time required to change the beam energy). In a preliminary consideration of application to rescanning techniques, the decrease in delivery time was estimated to be 4%-20%..

BACKGROUND AND PURPOSE: Currently, inverse planning for intensity-modulated radiotherapy (IMRT) can be a time-consuming trial and error process. This is because many planning objectives are inherently contradictory and cannot reach their individual optimum all at the same time. Therefore in clinical practice the potential of IMRT cannot be fully exploited for all patients. Multicriteria (multiobjective) optimization combined with interactive plan navigation is a promising approach to overcome these problems. PATIENTS AND METHODS: We developed a new inverse planning system called "Multicriteria Interactive Radiotherapy Assistant (MIRA)". The optimization result is a database of patient specific, Pareto-optimal plan proposals. The database is explored with an intuitive user interface that utilizes both a new interactive element for plan navigation and familiar dose visualizations in form of DVH and isodose projections. Two clinical test cases, one paraspinal meningioma case and one prostate case, were optimized using MIRA and compared with the clinically approved planning program KonRad. RESULTS: Generating the databases required no user interaction and took approx. 2-3h per case. The interactive exploration required only a few minutes until the best plan was identified, resulting in a significant reduction of human planning time. The achievable plan quality was comparable to KonRad with the additional benefit of having plan alternatives at hand to perform a sensitivity analysis or to decide for a different clinical compromise. CONCLUSIONS: The MIRA system provides a complete database and interactive exploration of the solution space in real time. Hence, it is ideally suited for the inherently multicriterial problem of inverse IMRT treatment planning..

In this paper, we propose an optimization concept for a rotation therapy technique which is referred to as arc-modulated cone beam therapy (AMCBT). The aim is a reduction of the treatment time while achieving a treatment plan quality equal to or better than that of IMRT. Therefore, the complete dose is delivered in one single gantry rotation and the beam is modulated by a multileaf collimator. The degrees of freedom are the field shapes and weights for a predefined number of beam directions. In the new optimization loop, the beam weights are determined by a gradient algorithm and the field shapes by a tabu search algorithm. We present treatment plans for AMCBT for two clinical cases. In comparison to step-and-shoot IMRT treatment plans, it was possible by AMCBT to achieve dose distributions with a better dose conformity to the target and a lower mean dose for the most relevant organ at risk. Furthermore, the number of applied monitor units was reduced for AMCBT in comparison to IMRT treatment plans..

Advanced radiotherapeutical techniques like intensity-modulated radiation therapy (IMRT) are based on an accurate knowledge of the location of the radiation target. An accurate dose delivery, therefore, requires a method to account for the inter- and intrafractional target motion and the target deformation occurring during the course of treatment. A method to compensate in real time for changes in the position and shape of the target is the use of a dynamic multileaf collimator (MLC) technique which can be devised to automatically arrange the treatment field according to real-time image information. So far, various approaches proposed for leaf sequencers have had to rely on a priori known target motion data and have aimed to optimize the overall treatment time. Since for a real-time dose delivery the target motion is not known a priori, the velocity range of the leading leaves is restricted by a safety margin to c x v(max) while the following leaves can travel with an additional maximum speed to compensate for the respective target movements. Another aspect to be considered is the tongue and groove effect. A uniform radiation field can only be achieved if the leaf movements are synchronized. The method presented in this note is the first to combine a synchronizing sequencer and real-time tracking with a dynamic MLC. The newly developed algorithm is capable of online optimizing the leaf velocities by minimizing the overall treatment time while at the same time it synchronizes the leaf trajectories in order to avoid the tongue and groove effect. The simultaneous synchronization is performed with the help of an online-calculated mid-time leaf trajectory which is common for all leaf pairs and which takes into account the real-time target motion and deformation information..

The most recent experimental results obtained with laser-plasma accelerators are applied to radio-therapy simulations. The narrow electron beam, produced during the interaction of the laser with the gas jet, has a high charge (0.5 nC) and is quasimonoenergetic (170 +/- 20 MeV). The dose deposition is calculated in a water phantom placed at different distances from the diverging electron source. We show that, using magnetic fields to refocus the electron beam inside the water phantom, the transverse penumbra is improved. This electron beam is well suited for delivering a high dose peaked on the propagation axis, a sharp and narrow tranverse penumbra combined with a deep penetration..

The aim of the work was to investigate in advance the dosimetric properties of a new multileaf collimator (MLC) concept with the help of Monte Carlo (MC) simulations prior to the production of a prototype. The geometrical design of the MLC was implemented in the MC code GEANT4. For the simulation of a 6 MV treatment beam, an experimentally validated phase space and a virtual spatial Gaussian-shaped model placed in the origin were used. For the simulation of the geometry in GEANT4, the jaws and the two leaf packages were implemented with the help of computer-aided design data. First, transmission values for different tungsten alloys were extracted using the simulation codes GEANT4 and BEAMnrc and compared to experimental measurements. In a second step, high-resolution simulations were performed to detect the leakage at depth of maximum dose. The 20%-80% penumbra along the travel direction of the leaves was determined using 10 x 10 cm2 fields shifted along the x- and y-axis. The simulated results were compared with measured data. The simulation of the transmission values for different tungsten alloys showed a good agreement with the experimental measurements (within 2.0%). This enabled an accurate estimation of the attenuation coefficient for the various leaf materials. Simulations with varying width of the spatial Gaussian distribution showed that the leakage and the penumbra depend very much on this parameter: for instance, for widths of 2 and 4 mm, the interleaf leakage is below 0.3% and 0.75%, respectively. The results for the leakage and the penumbra (4.7+/-0.5 mm) are in good agreement with the measurements. This study showed that GEANT4 is appropriate for the investigation of the dosimetric properties of a multileaf collimator. In particular, a quantification of the leakage, the penumbra, and the tongue-and-groove effect and an evaluation of the influence of the beam parameters such as the width of the Gaussian distribution was possible..

In this paper, we deal with the effects of interfractional organ motion during radiation therapy. We consider two problems: first, treatment plan evaluation in the presence of motion, and second, the incorporation of organ motion into IMRT optimization. Concerning treatment plan evaluation, we face the problem that the delivered dose cannot be predicted with certainty at the time of treatment planning but is associated with uncertainties. We present a method to simulate stochastic properties of the dose distribution. This provides the treatment planner with information about motion-related risks of different plans and may support the decision for or against a treatment plan. This information includes the display of probabilities of individual voxels to receive doses from a therapeutical interval or above critical levels, as well as a diagram that shows the variability of the dose volume histogram. Concerning the incorporation of organ motion into IMRT planning, we further analyse the approach of inverse planning based on probability distributions of possible patient geometries. We consider three different sources of uncertainty, namely uncertainty about the amplitude of motion, a systematic error and a random error. We analyse the impact of these sources of uncertainty on the optimized treatment plans for prostate cancer..

A new online imaging approach, linac-integrated cone beam CT (CBCT), has been developed over the past few years. It has the advantage that a patient can be examined in their treatment position directly before or during a radiotherapy treatment. Unfortunately, respiratory organ motion, one of the largest intrafractional organ motions, often leads to artefacts in the reconstructed 3D images. One way to take this into account is to register the breathing phase during image acquisition for a phase-correlated image reconstruction. Therefore, the main focus of this work is to present a system which has the potential to investigate the correlation between internal (movement of the diaphragm) and external (data of a respiratory gating system) information about breathing phase and amplitude using an inline CBCT scanner. This also includes a feasibility study about using the acquired information for a respiratory-correlated 4D CBCT reconstruction. First, a moving lung phantom was used to develop and to specify the required methods which are based on an image reconstruction using only projections belonging to a certain moving phase. For that purpose, the corresponding phase has to be detected for each projection. In the case of the phantom, an electrical signal allows one to track the movement in real time. The number of projections available for the image reconstruction depends on the breathing phase and the size of the position range from which projections should be used for the reconstruction. The narrower this range is, the better the inner structures can be located, but also the noise of the images increases due to the limited number of projections. This correlation has also been analysed. In a second step, the methods were clinically applied using data sets of patients with lung tumours. In this case, the breathing phase was detected by an external gating system (AZ-733V, Anzai Medical Co.) based on a pressure sensor attached to the patient's abdominal region with a fixation belt. The comparison of the reconstructed 4D CBCT images and the corresponding 4D CT images used for the treatment planning provides the required information for the calculation of possible setup errors. So, a repositioning of the patient is feasible even though the patient moves due to respiration. In addition to the external signal, the position of the diaphragm in the cranial-caudal direction could be extracted from each projection. Both independent sources of information show a very good agreement of the phase and even the amplitude of the movement and the external signal respectively. This suggests the usability of such a system for a gated dose delivery approach. However, more studies involving patients with different incidences have to be carried out to confirm these first results..

Radiotherapy treatment planning is associated with uncertainties. Examples are uncertainties in the tumour location due to organ movement or the inter/intra observer variability in target definition. Different approaches to incorporate uncertainties into IMRT optimization have been proposed. In this note, we point out a relation between two previously published methods: the coverage probability approach and the concept of optimizing the expectation value of an objective function that depends on a set of random variables. Both concepts are generally different, but turn out to be equivalent in special cases..

In this paper, we present a new technique for simultaneous multifield optimization of the biological effect (i.e. relative biological effectiveness times dose) for intensity modulated radiotherapy with ion beams. It offers complete inverse treatment planning by taking into account planning constraints for the target volume as well as for organs at risk. The approach is based on the mixed irradiation formalism of the linear-quadratic model from radiobiology. We employ a novel objective function to directly optimize the biological effect rather than the physical dose. The required biological input data are reduced to a minimum and are completely independent from the optimization itself. They can be derived from any radiobiological model or even from directly measured data. The new optimization method was fully integrated into our inverse treatment planning tool KonRad. Comparisons with the TRiP98 treatment planning code were done for simple spread-out Bragg peaks as well as for three-dimensional treatment plans, where all fields were optimized separately. While the agreement between both planning systems was very good, the calculation time was substantially reduced in KonRad. By enabling the multifield optimization, the quality of the treatment plans and the sparing of healthy tissues can be clearly improved..

BACKGROUND: The purpose of the study was the clinical implementation of a kV cone beam CT (CBCT) for setup correction in radiotherapy. PATIENTS AND METHODS: For evaluation of the setup correction workflow, six tumor patients (lung cancer, sacral chordoma, head-and-neck and paraspinal tumor, and two prostate cancer patients) were selected. All patients were treated with fractionated stereotactic radiotherapy, five of them with intensity modulated radiotherapy (IMRT). For patient fixation, a scotch cast body frame or a vacuum pillow, each in combination with a scotch cast head mask, were used. The imaging equipment, consisting of an x-ray tube and a flat panel imager (FPI), was attached to a Siemens linear accelerator according to the in-line approach, i.e. with the imaging beam mounted opposite to the treatment beam sharing the same isocenter. For dose delivery, the treatment beam has to traverse the FPI which is mounted in the accessory tray below the multi-leaf collimator. For each patient, a predefined number of imaging projections over a range of at least 200 degrees were acquired. The fast reconstruction of the 3D-CBCT dataset was done with an implementation of the Feldkamp-David-Kress (FDK) algorithm. For the registration of the treatment planning CT with the acquired CBCT, an automatic mutual information matcher and manual matching was used. RESULTS AND DISCUSSION: Bony landmarks were easily detected and the table shifts for correction of setup deviations could be automatically calculated in all cases. The image quality was sufficient for a visual comparison of the desired target point with the isocenter visible on the CBCT. Soft tissue contrast was problematic for the prostate of an obese patient, but good in the lung tumor case. The detected maximum setup deviation was 3 mm for patients fixated with the body frame, and 6 mm for patients positioned in the vacuum pillow. Using an action level of 2 mm translational error, a target point correction was carried out in 4 cases. The additional workload of the described workflow compared to a normal treatment fraction led to an extra time of about 10-12 minutes, which can be further reduced by streamlining the different steps. CONCLUSION: The cone beam CT attached to a LINAC allows the acquisition of a CT scan of the patient in treatment position directly before treatment. Its image quality is sufficient for determining target point correction vectors. With the presented workflow, a target point correction within a clinically reasonable time frame is possible. This increases the treatment precision, and potentially the complex patient fixation techniques will become dispensable..

One of the most prominent imaging techniques in image-guided radiotherapy (IGRT) is the acquisition of cone beam computed tomographies (CBCTs) at the linac with the patient in treatment position. CBCTs provide accurate 3-dimensional (3D) knowledge about the patient's anatomy for every treatment fraction and are therefore well suited for all adaptive corrections of errors related to interfractional uncertainties of the treatment process. In this paper, we first describe the technical development and implementation of this new imaging technique at our linac, i.e., the hardware components and their operating parameters are discussed in detail for a standard image acquisition of CBCTs. Then, an extension of this approach for the acquisition of complete images for extended field of views--the "shifted detector" technique--is presented followed by a first investigation of how CBCTs can be reliably used for adaptive dose calculations. Finally, a first clinical application, the process of automatic patient positioning based on CBCT images, is discussed. From our investigations, we conclude that the technical development of linac-integrated CBCTs bears an enormous potential for the correction of interfractional treatment errors. However, image quality and reconstruction speed of the images leave room for improvement. The development of clinical strategies for the optimal application of this new image modality in a clinical environment is one the major tasks for the future..

Today, inverse treatment planning for intensity modulated proton therapy (IMPT) usually employs a constant relative biological effectiveness (RBE). In this paper, the potential clinical relevance of RBE variations for scanning techniques in IMPT is investigated, and a new strategy to include the RBE into the inverse planning process is presented. Three-dimensional RBE distributions are calculated based on a phenomenological model that describes the RBE as a function of dose, linear energy transfer (LET) and tissue type in the framework of the linear-quadratic model. This RBE model is integrated into the optimization loop of inverse planning by using a modified version of the standard quadratic objective function, where the physical dose is replaced by the biological effect. This system for "biological optimization" was implemented into a research version of the inverse planning software KonRad and allows the direct optimization of the product of RBE and physical dose. Several treatment plans for a prostate case are presented, which compare the biological with the conventional physical dose optimization for IMPT scanning techniques, in particular distal edge tracking (DET) and the full three-dimensional (3D) modulation of beam spots. Mainly due to their different LET distributions, the RBE effects for these two techniques are quite different: while the RBE distribution was more or less homogeneous in the planning target volume (PTV) for 3D modulation, considerable RBE variations within the PTV were observed for DET. These unfavorable effects could be compensated for by employing the new biological objective function, which led to a more homogeneous distribution of the product of RBE and physical dose in the PTV. The computation time increased by a factor of 2 compared to the optimization of the physical dose. In conclusion, the proposed method allows the simultaneous multifield optimization of the biological effect in a reasonable time, and is therefore well suited for studying the influence of a variable RBE in IMPT as well as for minimizing potentially adverse effects..

PURPOSE: A planning study to analyze the impact of different leaf widths on the achievable dose distributions with intensity modulated radiation therapy (IMRT). METHODS: Five patients (3 intra- and 2 extra-cranial) with projected planning target volume (PTV) sizes smaller than 10 cm by 10 cm were re-planned with four different multileaf collimators (MLC). Two internal collimators with an isocentric leaf width of 4 and 10 mm and two add-on collimators with an isocentric leaf width of 2.75 and were evaluated. The inverse treatment planning system KonRad (Siemens Medical Solutions) was used to create IMRT 'step & shoot' plans. For each patient the same arrangement of beams and the same parameters for the optimization were used for all MLCs. The beamlet size for all treatment plans was chosen to coincide with the leaf width of the respective MLC. To evaluate the treatment plans 3D dose distributions and dose volume histograms were analyzed. As indicators for the quality of the PTV dose distribution the minimum dose, maximum dose and the standard deviation were used. For the organs at risk (OAR) the equivalent uniform dose (EUD) was calculated. To measure the dose coverage of the PTV the volume (V(90)) that received doses higher than 90% of the prescribed dose was calculated where for the conformity the dose conformity index given by Baltas et al. was determined. RESULTS: The MLC with the smallest leaf width yields the best mean value of all five patients for the PTV coverage and for the conformity. For the MLCs with the same leaf width, the add-on MLC leads to superior treatment plans than the internal MLC. This is due to the sharper penumbra of the add-on MLC. The number of IMRT field segments to deliver increased by approximately a factor of two if 2. MLC leafs are used instead of the standard 10 mm leafs. In case of the para-spinal patients the EUD value for the spinal cord is only reduced slightly by using MLCs with leaf widths smaller than 5 mm. For the intra-cranial the EUD value for some organs improved with reduced leaf widths while for some organs the 10 mm MLC leafs give comparable values. CONCLUSION: As expected the MLC with the smallest leaf width always yields the best PTV coverage. Reducing the leaf width from 4 to 2.75 mm results in a slight enhancement of the PTV coverage. With the selected organ parameters no significant improvement for most OAR was found. The disadvantage of the reduction of the leaf width is the increasing number of segments due to the more complex fluence patterns and therefore an increased delivery time..

.

We present a method to calculate dose uncertainties due to inter-fraction organ movements in fractionated radiotherapy, i.e. in addition to the expectation value of the dose distribution a variance distribution is calculated. To calculate the expectation value of the dose distribution in the presence of organ movements, one estimates a probability distribution of possible patient geometries. The respective variance of the expected dose distribution arises for two reasons: first, the patient is irradiated with a finite number of fractions only and second, the probability distribution of patient geometries has to be estimated from a small number of images and is therefore not exactly known. To quantify the total dose variance, we propose a method that is based on the principle of Bayesian inference. The method is of particular interest when organ motion is incorporated in inverse IMRT planning by means of inverse planning performed on a probability distribution of patient geometries. In order to make this a robust approach, it turns out that the dose variance should be considered (and minimized) in the optimization process. As an application of the presented concept of Bayesian inference, we compare three approaches to inverse planning based on probability distributions that account for an increasing degree of uncertainty. The Bayes theorem further provides a concept to interpolate between patient specific data and population-based knowledge on organ motion which is relevant since the number of CT images of a patient is typically small..

One aim of adaptive radiotherapy (ART) is the observation of organ motion followed by a subsequent adaptation of the treatment plan. One way of achieving this goal is a kV x-ray source mounted at a linear accelerator in combination with a flat-panel imager. Two imaging hardware configurations were evaluated for their potential for online tracking and the subsequent correction of organ motion by using fluoroscopic images: x-ray tube positioned with (A) 90 degrees and (B) 180 degrees offset to the MV beam. For one lung case two IMRT plans with five coplanar beams and the table positioned at 0 degrees were optimized for two multileaf collimators (MLCs) with 10 mm and 2.75 mm leaf width. Respiratory motion, modelled by rigid transformation in the lungs, was investigated for different amplitudes. The 3D dose distributions for different cases (no movement, uncorrected movement, correction for the movement perpendicular to the respective kV beam) were evaluated with the help of dose volume histograms (DVHs) and a modified conformity (Baltas et al 1998 Int. J. Radiat. Oncol. Biol. Phys. 40 515-24) and coverage index using the 90% isodose. For the corrected treatment plans the influence of the observed displacement vector caused by organ movement was accounted for by a respective displacement of the target point. For the simulated movement with a small amplitude (3 mm) in the anterior-posterior (AP) direction the dose distributions resulting from the correction of the displacement vector using imaging system A or B showed similar results for both systems and were in good agreement with the dose distribution of the static (not moving) patient. Increasing the amplitude in the AP direction to 6 mm or even 9 mm leads for both amplitudes and both MLCs to almost the same conformity and coverage index as the static dose distribution if imaging system B is used for the online correction. For the dose distribution obtained with correction based on imaging system A the deviation between the optimal and the corrected dose distribution is increasing with increasing amplitude. For the MLC with the smaller leaf width the difference between the optimal and the corrected dose distributions is always significantly larger than for the less conformal dose distributions created by the MLC with the 10 mm leaves. These results can be explained by the fact that system A cannot observe movement in the AP-LR plane perpendicular to the MV beam and therefore cannot correct for these movements whereas system B only fails to observe the motion in the beam direction which for photon irradiation has less impact on the dose distribution..

We investigate an off-line strategy to incorporate inter fraction organ movements in IMRT treatment planning. Nowadays, imaging modalities located in the treatment room allow for several CT scans of a patient during the course of treatment. These multiple CT scans can be used to estimate a probability distribution of possible patient geometries. This probability distribution can subsequently be used to calculate the expectation value of the delivered dose distribution. In order to incorporate organ movements into the treatment planning process, it was suggested that inverse planning could be based on that probability distribution of patient geometries instead of a single snapshot. However, it was shown that a straightforward optimization of the expectation value of the dose may be insufficient since the expected dose distribution is related to several uncertainties: first, this probability distribution has to be estimated from only a few images. And second, the distribution is only sparsely sampled over the treatment course due to a finite number of fractions. In order to obtain a robust treatment plan these uncertainties should be considered and minimized in the inverse planning process. In the current paper, we calculate a 3D variance distribution in addition to the expectation value of the dose distribution which are simultaneously optimized. The variance is used as a surrogate to quantify the associated risks of a treatment plan. The feasibility of this approach is demonstrated for clinical data of prostate patients. Different scenarios of dose expectation values and corresponding variances are discussed..

Respiratory organ motion is known to be one of the largest intrafractional organ motions. Therefore, it is important to investigate the potential benefit of gated dose delivery approaches which aim to account for the respective dose uncertainties. In this study respiration is simulated by a moving lung phantom; the movement is not restricted to a normal sinusoidal progression and simulates the one of the embedded lung tumour in the cranial-caudal direction. An IMRT plan with a total of 29 beam segments was designed for the treatment of this tumour. It was irradiated in its resting position-which is the position at exhalation-and during movement. Furthermore the irradiation was triggered using different amplitude thresholds, which means that the irradiation only proceeded if the deviation of the tumour's position from its resting position is smaller than the given threshold. We determined the gating-related increase of the treatment time for various gating procedures. We also measured the resulting dose distribution in specific slices of the phantom perpendicular to the direction of the movement using film dosimetry and compared it to the dose distribution of the static case. Since these film measurements cannot be done inside the whole tumour, additionally the movement and gating was simulated using the planning software to calculate the 3D dose distribution inside the tumour and to generate dose volume histograms for different treatment modalities. The total treatment time was observed to increase by 20%-100% depending on the individual gating threshold and can be calculated easily. The analysis of the films showed that irradiation without gating leads to significant underdosages up to 33%, especially at the edge of the tumour. With gating it is possible to considerably reduce this underdosage down to 9% depending on the trigger threshold. The calculation of the dose volume histograms makes it possible to find a reasonable compromise between the improvement of the dose distribution and the increase of the treatment time..

The aim of this treatment planning comparison study was to explore different spinal irradiation techniques with respect to the risk of late side-effects, particularly radiation-induced cancer. The radiotherapy techniques compared were conventional photon therapy, intensity modulated x-ray therapy (IMXT), conventional electron therapy, intensity/energy modulated electron therapy (IMET) and proton therapy (IMPT).CT images for radiotherapy use from five children, median age 8 and diagnosed with medulloblastoma, were selected for this study. Target volumes and organs at risk were defined in 3-D. Treatment plans using conventional photon therapy, IMXT, conventional electron therapy, IMET and IMPT were set up. The probability of normal tissue complication (NTCP) and the risk of cancer induction were calculated using models with parameters-sets taken from published data for the general population; dose data were taken from dose volume histograms (DVH). Similar dose distributions in the targets were achieved with all techniques but the absorbed doses in the organs-at-risk varied significantly between the different techniques. The NTCP models based on available data predicted very low probabilities for side-effects in all cases. However, the effective mean doses outside the target volumes, and thus the predicted risk of cancer induction, varied significantly between the techniques. The highest lifetime risk of secondary cancers was estimated for IMXT (30%). The lowest risk was found with IMPT (4%). The risks associated with conventional photon therapy, electron therapy and IMET were 20%, 21% and 15%, respectively. This model study shows that spinal irradiation of young children with photon and electron techniques results in a substantial risk of radiation-induced secondary cancers. Multiple beam IMXT seems to be associated with a particularly high risk of secondary cancer induction. To minimise this risk, IMPT should be the treatment of choice. If proton therapy is not available, advanced electron therapy may provide a better alternative..

In this paper, we investigate an off-line strategy to incorporate inter-fraction organ motion in IMRT treatment planning. It was suggested that inverse planning could be based on a probability distribution of patient geometries instead of a single snap shot. However, this concept is connected to two intrinsic problems: first, this probability distribution has to be estimated from only a few images; and second, the distribution is only sparsely sampled over the treatment course due to a finite number of fractions. In the current work, we develop new concepts of inverse planning which account for these two problems..

To study the effects of a variable relative biological effectiveness (RBE) in inverse treatment planning for proton therapy, fast methods for three-dimensional RBE calculations are required. We therefore propose a simple phenomenological model for the RBE in therapeutic proton beams. It describes the RBE as a function of the dose, the linear energy transfer (LET) and tissue specific parameters. Published experimental results for the dependence of the parameters alpha and beta from the linear-quadratic model on the dose averaged LET were evaluated. Using a linear function for alpha(LET) in the relevant LET region below 30 keV per micrometre and a constant beta, a simple formula for the RBE could be derived. The new model was able to reproduce the basic dependences of RBE on dose and LET, and the RBE values agreed well with experimental results. The model was also applied to spread-out Bragg peaks (SOBP), where the main effects of a variable RBE are an increase of the RBE along the SOBP plateau, and a shift in depth of the distal falloff. The new method allows fast RBE estimations and has therefore potential applications in iterative treatment planning for proton therapy..

BACKGROUND AND PURPOSE: In recent years photon intensity modulated radiation therapy (IMRT) has gained attention due to its ability to improve conformity of dose distributions. A potential advantage of electron-IMRT is that the dose fall off in the depth dose curve makes it possible to modulate the dose distribution in the direction of the beam by selecting different electron energies. This paper examines the use of a computer based energy selection in combination with the IMRT technique to optimise the electron dose distribution. MATERIALS AND METHODS: One centimetre square electron beamlets ranging from 2.5 to 50 MeV were pre-calculated in water using Monte Carlo methods. A modified IMRT optimisation tool was then used to find an optimum mix of electron energies and intensities. The main principles used are illustrated in some simple geometries and tested on two clinical cases of post-operated ca. mam. RESULTS: It is clearly illustrated that the energy optimisation procedure lowers the dose to lung and heart and makes the dose in the target more homogeneous. Increasing the energy at steep gradients compensates for lack of target coverage at beam edges and steep gradients. Comparison with a clinically acceptable four segment plan indicates the advantage of the used electron IMRT technique. CONCLUSIONS: Using an intensity optimised mix of computer selected electron energies has the potential to improve electron treatments for mastectomy patients with good target coverage and reduced dose to normal tissue such as lung and heart..

Complex dose-delivery techniques, as currently applied in intensity-modulated radiation therapy (IMRT), require a highly efficient treatment-verification process. The present paper deals with the problem of the scatter correction for therapy verification by use of portal images obtained by an electronic portal imaging device (EPID) based on amorphous silicon. It also presents an iterative method for the scatter correction of portal images based on Monte Carlo-generated scatter kernels. First applications of this iterative scatter-correction method for the verification of intensity-modulated treatments are discussed on the basis of MVCT- and dose reconstruction. Several experiments with homogeneous and anthropomorphic phantoms were performed in order to validate the scatter correction method and to investigate the precision and relevance in view of its clinical applicability. It is shown that the devised concept of scatter correction significantly improves the results of MVCT- and dose reconstruction models, which is in turn essential for an exact online IMRT verification..

A common requirement of radiation therapy is that treatment planning for different radiation modalities is devised on the basis of the same treatment planning system (TPS). The present study presents a novel multi-modal TPS with separate modules for the dose calculation, the optimization engine and the graphical user interface, which allows to integrate different treatment modalities. For heavy-charged particles, both most promising techniques, the distal edge tracking (DET) and the 3-dimensional scanning (3D) technique can be optimized. As a first application, the quality of optimized intensity-modulated treatment plans for photons (IMXT) and protons (IMPT) was analyzed in one clinical case on the basis of the achieved physical dose distributions. A comparison of the proton plans with the photon plans showed no significant improvement in terms of target volume dose, however there was an improvement in terms of organs at risk as well as a clear reduction of the total integral dose. For the DET technique, it is possible to create a treatment plan with almost the same quality of the 3D technique, however with a clearly reduced number (factor of 5) of beam spots as well as a reduced optimization time. Due to its modular design, the system can be easily expanded to more sophisticated dose-calculation algorithms or to modeling of biological effects..

The clinical impact of a variable relative biological effectiveness (RBE) in proton therapy is still under investigation. As the RBE depends on the linear energy transfer (LET), we present a fast method for three-dimensional calculations of the dose-averaged LET. The corresponding algorithm, based on an analytical expression for the LET and on Monte Carlo simulations, accounts for tissue inhomogeneities and allows LET calculations for realistic treatment plans and patient geometries given by computed tomography (CT) data sets. LET distributions were calculated for different proton spot scanning techniques, in particular for the full three-dimensional (3D) modulation and for the distal edge tracking technique (DET). In the case of 3D modulation, a more homogeneous distribution of LET in the planning target volume was observed than for the DET. Consequently, the proposed method allows to assess the impact of potential variations in RBE for various dose delivery techniques..

As the relative biological effectiveness of protons depends on the linear energy transfer (LET), simple methods for LET calculations are desired for the optimization of proton therapy. This work provides an analytical model for the LET on the central axis of broad proton beams in water, which can also be applied to spread-out Bragg peaks. For realistic treatment situations with polyenergetic beams, the LET is here defined as a local mean of the proton stopping power, weighted by the local energy spectrum. The proposed model considers only Coulomb interactions and neglects nonelastic nuclear interactions. By assuming a Gaussian shape for the energy spectrum and by using a suitable parametrization of the stopping power, analytical expressions for the track averaged and the dose averaged LET are derived, which account for range straggling as well as for the initial width of the energy spectrum. The analytical model was evaluated by Monte Carlo simulations with GEANT 3.21. Local energy spectra were simulated to obtain LET distributions for several cases, using clinical energies between 70 and 250 MeV and varying widths of the initial energy spectrum. Good agreement was found between the analytical model and the Monte Carlo simulations (with maximum deviations of 0.5 keV per micrometer), which justifies the assumptions used in the derivation of the analytical model..

For proton dose calculations in heterogeneous media, it was shown in a previous work that the conventional pencil beam approach based on pathlength scaling does not properly account for scattering effects in nonwater media (Szymanowski and Oelfke 2002 Phys. Med. Biol. 47 3313-30). A two-dimensional scaling method was therefore introduced, which is able to predict with high accuracy the propagation of proton pencil beams both along the depth and the lateral directions in inhomogeneous media. In order to integrate this improved pencil beam algorithm in a CT based treatment planning system, two CT calibration curves are needed. The first one relates the Hounsfield numbers to the relative stopping powers, as for the conventional pencil beam approach. The second curve is to relate the Hounsfield numbers to the material-specific lateral scaling factors. The purpose of this work is to provide the CT calibration curves needed for the integration of the pencil beam algorithm featuring the two-dimensional scaling method. Similarly to as suggested by Schneider et al (1996 Phys. Med. Biol. 41 111-24) for the calibration curve in terms of stopping powers, we follow a stoichiometric procedure to get the calibration curve in terms of material-specific lateral scaling factors. The calibration curves for a CT scanner of the type Siemens Somatom Plus 4 are obtained from the analytical calculation of the CT Hounsfield numbers, relative stopping powers and material-specific lateral scaling factors for human biological tissues..

To investigate the role of sophisticated dose calculation methods for treatment planning, we compared conventional pencil beam optimized 6 and 15 MV intensity-modulated treatment plans with optimizations based on the superposition technique. Five lung and five head and neck IMRT cases with spatial resolutions of bixels and dose voxels usually employed in clinical practice were considered for tumor volumes between 15 and 500 cm3. We investigated the systematic error of the pencil beam algorithm and the pencil beam induced error to the optimal solution of bixel weights. For the lung cases, the pencil beam overestimated the mean dose deposited inside the planning target volume (PTV) by about 8%, for small lung tumors even up to 20.6%. In the head and neck cases only a slight overestimation in mean PTV dose of 1.5% was observed. The optimization with the superposition method substantially improved the dose coverage of the considered radiation targets. Additionally, for the head and neck cases, the brainstem was significantly spared by about 4% mean PTV dose through the use of the superposition technique. Our studies showed that, in target regions with intricate tissue inhomogeneities, superposition or Monte Carlo techniques have to be used for the optimization and the final dose calculation of intensity-modulated treatment plans..

Intensity modulated radiotherapy with high enery photons (IMRT) and with charged particles (IMPT) refer to the most advanced development in conformal radiation therapy. Their general aim is to increase local tumor control rates while keeping the radiation induced complications below desired thresholds. IMRT is currently widely introduced in clinical practice. However, the more complicated IMPT is still under development. Especially, spot- scanning techniques integrated in rotating gantries that can deliver proton or light ion-beams to a radiation target from any direction will be available in the near future. We describe the basic concepts of intensity modulated particle therapy (IMPT). Starting from the potential advantages of hadron therapy inverse treatment planning strategies are discussed for various dose delivery techniques of IMPT. Of special interest are the techniques of distal edge tracking (DET) and 3D-scanning. After the introduction of these concepts a study of comparative inverse treatment planning is presented. The study aims to identify the potential advantages of achievable physical dose distributions with proton and carbon beams, if different dose delivery techniques are employed. Moreover, a comparison to standard photon IMRT is performed. The results of the study are summarized as: i) IMRT with photon beams is a strong competitor to intensity modulated radiotherapy with charged particles. The most obvious benefit observed for charged particles is the reduction of medium and low doses in organs at risk. ii) The 3D-scanning technique could not improve the dosimetric results achieved with DET, although 10-15 times more beam spots were employed for 3D-scanning than for DET. However, concerns may arise about the application of DET, if positioning errors of the patient or organ movements have to be accounted for. iii) Replacing protons with carbon ions leads to further improvements of the physical dose distributions. However, the additional degree of improvement due to carbon ions is modest. The main clinical potential of heavy ion beams is probably related to their radiobiological properties..

BACKGROUND AND PURPOSE: Dose calculation algorithms play a central role for the optimization and verification of treatment plans in radiation therapy. Complex treatment techniques like intensity modulated radiotherapy (IMRT) require accurate and fast dose algorithms especially for clinical cases which involve severe tissue inhomogeneities. For these cases the standard dose engine in current treatment planning systems--the convolution of photon pencil beams--usually fails to predict the dose with the required accuracy. The role of more accurate but time consuming dose calculations like superposition algorithms or Monte Carlo simulations in clinical practice is under investigation at several therapy centers. PATIENTS AND METHODS: The paper presents the design, implementation and the first application of a superposition algorithm in a clinical setting at the German Cancer Research Center (DKFZ). It first describes in detail how the superposition algorithm is adapted to the dose delivery system at DKFZ in terms of standard dosimetric data. Then details of the implementation of the algorithm are given with a focus on various methods for the reduction of dose computation times. Next, the algorithm is evaluated in various experiments with dosimetric phantoms. These studies are employed for the development of time efficient sampling strategies of the elemental dose kernels. Finally, the algorithm is applied to dose calculations of clinical cases with tumors adjacent to lung tissue. RESULTS: Severe differences in dose coverage of the tumors and dose burden of the surrounding tissues in comparison to standard pencil beam calculations are observed. A standard 4-7 beam plan in a convenient dose grid (approximately 3 mm in each direction) is calculated in about 30 min on a Pentium 4 (1.9 GHz) applying the superposition algorithm described here..

In order to provide automatic IMRT dose delivery with an add-on MMLC a technical integration of a MMLC system with a linear accelerator was realized. The principle of this integration and the changes and enhancements of the existing hard- and software are briefly described. The system was tested by the automatic delivery of an IMRT plan designed for a head and neck phantom. A verification of dose delivery was performed with film dosimetry. The plan consisting of 78 "step and shoot" segments could be delivered within 17 minutes. A high spatial accuracy of the fluence pattern at the isocentre was reached by a resolution of 2.75x2.75 mm(2). The measured dose profiles were within 3% of the maximum dose of the calculated profiles..

In inverse planning for intensity-modulated radiotherapy, the dose calculation is a crucial element limiting both the maximum achievable plan quality and the speed of the optimization process. One way to integrate accurate dose calculation algorithms into inverse planning is to precalculate the dose contribution of each beam element to each voxel for unit fluence. These precalculated values are stored in a big dose calculation matrix. Then the dose calculation during the iterative optimization process consists merely of matrix look-up and multiplication with the actual fluence values. However, because the dose calculation matrix can become very large, this ansatz requires a lot of computer memory and is still very time consuming, making it not practical for clinical routine without further modifications. In this work we present a new method to significantly reduce the number of entries in the dose calculation matrix. The method utilizes the fact that a photon pencil beam has a rapid radial dose falloff, and has very small dose values for the most part. In this low-dose part of the pencil beam, the dose contribution to a voxel is only integrated into the dose calculation matrix with a certain probability. Normalization with the reciprocal of this probability preserves the total energy, even though many matrix elements are omitted. Three probability distributions were tested to find the most accurate one for a given memory size. The sampling method is compared with the use of a fully filled matrix and with the well-known method of just cutting off the pencil beam at a certain lateral distance. A clinical example of a head and neck case is presented. It turns out that a sampled dose calculation matrix with only 1/3 of the entries of the fully filled matrix does not sacrifice the quality of the resulting plans, whereby the cutoff method results in a suboptimal treatment plan..

New dose delivery techniques with proton beams, such as beam spot scanning or raster scanning, require fast and accurate dose algorithms which can be applied for treatment plan optimization in clinically acceptable timescales. The clinically required accuracy is particularly difficult to achieve for the irradiation of complex, heterogeneous regions of the patient's anatomy. Currently applied fast pencil beam dose calculations based on the standard inhomogeneity correction of pathlength scaling often cannot provide the accuracy required for clinically acceptable dose distributions. This could be achieved with sophisticated Monte Carlo simulations which are still unacceptably time consuming for use as dose engines in optimization calculations. We therefore present a new algorithm for proton dose calculations which aims to resolve the inherent problem between calculation speed and required clinical accuracy. First, a detailed derivation of the new concept, which is based on an additional scaling of the lateral proton fluence is provided. Then, the newly devised two-dimensional (2D) scaling method is tested for various geometries of different phantom materials. These include standard biological tissues such as bone, muscle and fat as well as air. A detailed comparison of the new 2D pencil beam scaling with the current standard pencil beam approach and Monte Carlo simulations, performed with GEANT, is presented. It was found that the new concept proposed allows calculation of absorbed dose with an accuracy almost equal to that achievable with Monte Carlo simulations while requiring only modestly increased calculation times in comparison to the standard pencil beam approach. It is believed that this new proton dose algorithm has the potential to significantly improve the treatment planning outcome for many clinical cases encountered in highly conformal proton therapy..

The concept of inverse planning for intensity-modulated radiation therapy and its application for photon and charged particle beams is presented. Starting from theoretical solutions of the "inverse problem" in radiation therapy, a clinically applied optimization approach is discussed. A central topic is the mathematical formulation of clinical objectives in terms of physical parameters such as dose levels and irradiated volumes. Examples for practical inverse treatment planning and its clinical application for photon beams are provided. Inverse treatment planning of dose delivery techniques with charged particle beams is discussed by extending the conventional planning concept. A new multimodality inverse planning tool is described and applied to an example of comparative planning between photon and proton IMRT..

The optimization of intensity modulated radiotherapy (IMRT) for charged particle beams is a necessary prerequisite to evaluate the clinical potential of this treatment modality in comparison to IMRT with high energy photons. A theoretical study for IMRT with charged particle beams delivered by rotation therapy is presented. First, the inverse problem for two-dimensional rotation therapy with arbitrary depth dose curves is formulated. Then a numerical strategy is devised to calculate fluence profiles for the simplified case of arbitrary rotationally invariant dose distributions. This mathematical framework is applied to study various aspects of charged particle IMRT. A central topic of the investigation is the evaluation of dose delivery, based on distal edge tracking (DET) and intensity modulation. The potential of DET-IMRT with charged particle beams is studied in comparison to an optimal, conventional dose delivery technique, which employs the concept of a spread-out Bragg peak (SOBP). Moreover, a comparison to photon IMRT is provided for simple geometric dose patterns. The technique of DET-IMRT for the delivery of a homogeneous target dose is only feasible for targets up to a critical radius, depending on the individual shape of the employed Bragg peak. The irradiation of larger targets requires energy modulation in addition to the range modulation for DET. The accurate placement of the Bragg peak with respect to the target edge is found to be of potential importance. Comparing dose delivery via DET-IMRT with the optimal SOBP technique revealed a significant advantage of DET-IMRT, especially a saved integral dose in target-adjacent healthy tissues of up to 30%, and a reduction of the penumbra at the target edge by almost 50%. A saving in integral dose to healthy tissues by a factor of 2-3 was observed for DET-IMRT in comparison to IMRT with high energy photons..

UNLABELLED: The following question is investigated: How narrow do the leaves of a multileaf collimator have to be such that further reduction of the leaf width does not lead to physical improvements of the dose distribution. Because of the physical principles of interaction between radiation and matter, dose distributions in radiotherapy are generally relatively smooth. According to the theory of sampling, the dose distribution can therefore be represented by a set of evenly spaced samples. The distance between the samples is identified with the distance between the leaf centers of a multileaf collimator. The optimum sampling distance is derived from the 20% to 80% field edge penumbra through the concept of the dose deposition kernel, which is approximated by a Gaussian. The leaf width of the multileaf collimator is considered to be independent from the sampling distance. Two cases are studied in detail: (i) the leaf width equals the sampling distance, which is the regular case, and (ii) the leaf width is twice the sampling distance. The practical delivery of the latter treatment geometry requires a couch movement or a collimator rotation. The optimum sampling distance equals the 20%-80% penumbra divided by 1.7 and is on the order of 1.5-2 mm for a typical 6 MV beam in soft tissue. The optimum leaf width equals this sampling distance in the regular case. A relatively small deterioration results if the leaf width is doubled, while the sampling distance remains the same. The deterioration can be corrected for by deconvolving the fluence profile with an inverse filter. CONCLUSIONS: With the help of the sampling theory and, more generally, the theory of linear systems, one can find a general answer to the question about the optimum leaf width of a multileaf collimator from a physical point of view. It is important to distinguish between the sampling distance and the leaf width. The sampling distance is more critical than the leaf width. The leaf width can be up to twice as large as the sampling width. Furthermore, the derived sampling distance can be used to select the optimum resolution of both the fluence and the dose grid in dose calculation and inverse planning algorithms..

Rotation therapy with photons is currently under investigation for the delivery of intensity modulated radiotherapy (IMRT). An analytical approach for inverse treatment planning of this radiotherapy technique is described. The inverse problem for the delivery of arbitrary 2D dose profiles is first formulated and then solved analytically. In contrast to previously applied strategies for solving the inverse problem, it is shown that the most general solution for the fluence profiles consists of two independent solutions of different parity. A first analytical expression for both fluence profiles is derived. The mathematical derivation includes two different strategies, an elementary expansion of fluence and dose into polynomials and a more practical approach in terms of Fourier transforms. The obtained results are discussed in the context of previous work on this problem..

A new algorithm for the reconstruction of two-dimensional (2D) images from projections is described. The algorithm is based on the decomposition of the projections into Chebyshev polynomials of the second kind, which are the ideal basis functions for this application. The Chebyshev decomposition is done via the fast discrete sine transform. A discrete reconstruction filter is applied that corresponds to the ramp filter used in standard filtered backprojection (FBP) reconstruction. In contrast to FBP, the filter is applied to the Chebyshev coefficients and not to the Fourier coefficients of the projections. Then the reconstructed image is simply obtained by means of backprojection. Consequently, the method can be considered as a Chebyshev-domain filtered backprojection (CD-FBP). The total calculation time is dominated by the backprojection step only and is comparable to FBP. The merits of CD-FBP as compared with standard FBP are that: (a) The result is exact if the 2D function to be reconstructed can be decomposed into polynomials of finite degree, and if the sampling is adequate. Otherwise a polynomial approximation results. (b) The algorithm is inherently discrete. (c) It is particularly well suited for reconstructions from projections with non-equidistant samples that occur for instance in 2D PET (positron emission tomography) imaging and in a special form of fan beam scanning. Examples of applications comprise reconstructions of the Shepp and Logan head phantom in various sampling geometries, and a real PET test object. In the PET example an increased resolution is observed in comparison with standard FBP..

A fast optimization algorithm for range modulation in clinical and experimental applications of proton therapy is described. The method is versatile towards the number of parameters provided for range modulation, i.e. the trade-off between accuracy and simplicity of the latter can be chosen freely. The approach is, therefore, adaptable to most operating proton therapy facilities. It requires only a few basic measurements as input data and results in a depth dose uniformity of better than 2%. A typical calculation takes less than 90 s on a DEC VAXstation 3100, FORTRAN. The method has been extensively tested at the TRIUMF Proton Therapy Facility in Vancouver, BC..

.

The feasibility of using PET for proton dose monitoring is examined here in detail. First experimental studies in a Lucite phantom have been performed at the medical TRIUMF proton beamline for proton energies of 62 MeV and 110 MeV. The proton dose delivered to the phantom ranged from 16 Gy up to 317 Gy. The induced activity was analysed 20-40 min after the irradiation with a PET scanner. The obtained depth activity profiles were compared to our calculation based on a model using available isotope production cross-section data. Both the observed absolute count rates and the activity profiles were found to agree very well with this model. Effects such as proton range straggling, inelastic nuclear interactions and the energy spectrum of the emitted positrons were studied in detail and found to change the activities by 5-10%. The lateral deposition of dose in the phantom could be very well localized by the induced activity. However, the spatial correlation between dose depth profiles and depth activity profiles was found to be poor, hence the extraction of isodose profiles from activity profiles seems to be very difficult..

Measurements of relative biological effectiveness (RBE) have been made on the range-modulated 70 MeV proton beam at TRIUMF using a precise cell sorting survival assay. In this study, Chinese hamster V79-WNRE cells were suspended in medium containing liquid gelatin at 37 degrees C in irradiation tubes and the gel was allowed to solidify by cooling to 4 degrees C. Complete cell survival responses were measured at 11 positions with 2 mm spacing within a proton stopping peak width of approximately 2 cm. Survival responses after proton irradiation were compared with responses to 60Co gamma rays measured at the same time, and RBE values were determined as a function of both dose and depth. Above doses of 4 Gy, the average RBE for these cells throughout the modulated proton stopping distribution was 1.21 +/- 0.05, measured at a survival of 1%. However, we also observed that, within the spread-out Bragg peak, the RBE increased with increasing depth, from approximately 1.2 at the proximal part to > 1.3 at the distal part of the peak. At the distal edge of the stopping distribution, the RBE value increased significantly, to an extent that may be of concern when this region of the treatment volume is close to sensitive tissues. Below 4 Gy, the RBE value was also dependent on radiation dose, increasing significantly to values of approximately 1.37 and 1.56 at 2 and 1 Gy, respectively. Our results illustrate that the use of a single RBE value in different irradiation protocols can be an oversimplification, and argues for the use of "proton gray doses" rather than "gamma-ray equivalent grays.".

.

.

.

.

.

.

.

.