Global democratisation of proton radiotherapy.

Proton radiotherapy is an advanced treatment option compared with conventional x-ray treatment, delivering much lower doses of radiation to healthy tissues surrounding the tumour. However, proton therapy is currently not widely available. In this Review, we summarise the evolution of proton therapy to date, together with the benefits to patients and society. These developments have led to an exponential growth in the number of hospitals using proton radiotherapy worldwide. However, the gap between the number of patients who should be treated with proton radiotherapy and those who have access to it remains large. We summarise the ongoing research and development that is contributing to closing this gap, including the improvement of treatment efficiency and efficacy, and advances in fixed-beam treatments that do not require an enormously large, heavy, and costly gantry. The ultimate goal of decreasing the size of proton therapy machines to fit into standard treatment rooms appears to be within reach, and we discuss future research and development opportunities to achieve this goal.

Simultaneous T2 -weighted real-time MRI of two orthogonal slices.

PURPOSE: MR guidance is used during therapy to detect and compensate for lesion motion. T2 -weighted MRI often has a superior lesion contrast in comparison to T1 -weighted real-time imaging. The purpose of this work was to design a fast T2 -weighted sequence capable of simultaneously acquiring two orthogonal slices, enabling real-time tracking of lesions. METHODS: To generate a T2 contrast in two orthogonal slices simultaneously, a sequence (Ortho-SFFP-Echo) was designed that samples the T2 -weighted spin echo (S- ) signal in a TR-interleaved acquisition of two slices. Slice selection and phase-encoding directions are swapped between the slices, leading to a unique set of spin-echo signal conditions. To minimize motion-related signal dephasing, additional flow-compensation strategies are implemented. In both the abdominal breathing phantom and in vivo experiments, a time series was acquired using Ortho-SSFP-Echo. The centroid of the target was tracked in postprocessing steps. RESULTS: In the phantom, the lesion could be identified and delineated in the dynamic images. In the volunteer experiments, the kidney was visualized with a T2 contrast at a temporal resolution of 0.45 s under free-breathing conditions. A respiratory belt demonstrated a strong correlation with the time course of the kidney centroid in the head-foot direction. A hypointense saturation band at the slice overlap did not inhibit lesion tracking in the semi-automatic postprocessing steps. CONCLUSION: The Ortho-SFFP-Echo sequence delivers real-time images with a T2 -weighted contrast in two orthogonal slices. The sequence allows for simultaneous acquisition, which could be beneficial for real-time motion tracking in radiotherapy or interventional MRI.

Combined clinical and research training in medical physics in a multi-institutional setting: 13-year experience of Harvard Medical Physics Residency Program.

PURPOSE: This manuscript describes the structure, management and outcomes of a multi-institutional clinical and research medical physics residency program (Harvard Medical Physics Residency Program, or HMPRP) to provide potentially useful information to the centers considering a multi-institutional approach for their training programs. METHODS: Data from the program documents and public records was used to describe HMPRP and obtain statistics about participating faculty, enrolled residents, and graduates. Challenges associated with forming and managing a multi-institutional program and developed solutions for effective coordination between several clinical centers are described. RESULTS: HMPRP was formed in 2009 and was accredited by the Commission on Accreditation of Medical Physics Education Programs (CAMPEP) in 2011. It is a 3-year therapy program, with a dedicated year of research and the 2 years of clinical training at three academic hospitals. A CAMPEP-accredited Certificate Program is embedded in HMPRP to allow enrolled residents to complete a formal didactic training in medical physics if necessary. The clinical training covers the material required by CAMPEP. In addition, training in protons, CyberKnife, MR-linac, and at network locations is included. The clinical training and academic record of the residents is outstanding. All graduates have found employment within clinical medical physics, mostly at large academic centers and graduates had a 100% pass rate at the oral American Board of Radiology exams. On average, three manuscripts per resident are published during residency, and multiple abstracts are presented at conferences. CONCLUSIONS: A multi-institutional medical physics residency program can be successfully formed and managed. With a collaborative administrative structure, the program creates an environment for high-quality clinical training of the residents and high productivity in research. The main advantage of such program is access to a wide variety of resources. The main challenge is creating a structure for efficient management of multiple resources at different locations. This report may provide valuable information to centers considering starting a multi-institutional residency program.

Diffusion tensor transformation for personalizing target volumes in radiation therapy.

Diffusion tensor imaging (DTI) is used in tumor growth models to provide information on the infiltration pathways of tumor cells into the surrounding brain tissue. When a patient-specific DTI is not available, a template image such as a DTI atlas can be transformed to the patient anatomy using image registration. This study investigates a model, the invariance under coordinate transform (ICT), that transforms diffusion tensors from a template image to the patient image, based on the principle that the tumor growth process can be mapped, at any point in time, between the images using the same transformation function that we use to map the anatomy. The ICT model allows the mapping of tumor cell densities and tumor fronts (as iso-levels of tumor cell density) from the template image to the patient image for inclusion in radiotherapy treatment planning. The proposed approach transforms the diffusion tensors to simulate tumor growth in locally deformed anatomy and outputs the tumor cell density distribution over time. The ICT model is validated in a cohort of ten brain tumor patients. Comparative analysis with the tumor cell density in the original template image shows that the ICT model accurately simulates tumor cell densities in the deformed image space. By creating radiotherapy target volumes as tumor fronts, this study provides a framework for more personalized radiotherapy treatment planning, without the use of additional imaging.

A reference-point-method-based online proton treatment plan re-optimization strategy and a novel solution to planning constraint infeasibility problem.

Objective. Propose a highly automated treatment plan re-optimization strategy suitable for online adaptive proton therapy. The strategy includes a rapid re-optimization method that generates quality replans and a novel solution that efficiently addresses the planning constraint infeasibility issue that can significantly prolong the re-optimization process.Approach. We propose a systematic reference point method (RPM) model that minimizes the l-infinity norm from the initial treatment plan in the daily objective space for online re-optimization. This model minimizes the largest objective value deviation among the objectives of the daily replan from their reference values, leading to a daily replan similar to the initial plan. Whether a set of planning constraints is feasible with respect to the daily anatomy cannot be known before solving the corresponding optimization problem. The conventional trial-and-error-based relaxation process can cost a significant amount of time. To that end, we propose an optimization problem that first estimates the magnitude of daily violation of each planning constraint. Guided by the violation magnitude and clinical importance of the constraints, the constraints are then iteratively converted into objectives based on their priority until the infeasibility issue is solved.Main results.The proposed RPM-based strategy generated replans similar to the offline manual replans within the online time requirement for six head and neck and four breast patients. The average targetD95and relevant organ at risk sparing parameter differences between the RPM replans and clinical offline replans were -0.23, -1.62 Gy for head and neck cases and 0.29, -0.39 Gy for breast cases. The proposed constraint relaxation solution made the RPM problem feasible after one round of relaxation for all four patients who encountered the infeasibility issue.Significance. We proposed a novel RPM-based re-optimization strategy and demonstrated its effectiveness on complex cases, regardless of whether constraint infeasibility is encountered.

Integrating muscle fiber orientation from visible human data into radiotherapy target volumes.

Objective.A major challenge in treatment of tumors near skeletal muscle is defining the target volume for suspected tumor invasion into the muscle. This study develops a framework that generates radiation target volumes with muscle fiber orientation directly integrated into their definition. The framework is applied to nineteen sacral tumor patients with suspected infiltration into surrounding muscles.Approach.To compensate for the poor soft-tissue contrast of CT images, muscle fiber orientation is derived from cryo-images of two cadavers from the human visible project (VHP). The approach consists of (a) detecting image gradients in the cadaver images representative of muscle fibers, (b) mapping this information onto the patient image, and (c) embedding the muscle fiber orientation into an expansion method to generate patient-specific clinical target volumes (CTV). The validation tested the consistency of image gradient orientation across VHP subjects for the piriformis, gluteus maximus, paraspinal, gluteus medius, and gluteus minimus muscles. The model robustness was analyzed by comparing CTVs generated using different VHP subjects. The difference in shape between the new CTVs and standard CTV was analyzed for clinical impact.Main results.Good agreement was found between the image gradient orientation across VHP subjects, as the voxel-wise median cosine similarity was at least 0.86 (for the gluteus minimus) and up to 0.98 for the piriformis. The volume and surface similarity between the CTVs generating from different VHP subjects was on average at least 0.95 and 5.13 mm for the Dice Similarity Coefficient and the Hausdorff 95% Percentile Index, showing excellent robustness. Finally, compared to the standard CTV with different margins in muscle and non-muscle tissue, the new CTV margins are reduced in muscle tissue depending on the chosen clinical margins.Significance.This study implements a method to integrate muscle fiber orientation into the target volume without the need for additional imaging.

The influence of anisotropy on the clinical target volume of brain tumor patients.

Objective.Current radiotherapy guidelines for glioma target volume definition recommend a uniform margin expansion from the gross tumor volume (GTV) to the clinical target volume (CTV), assuming uniform infiltration in the invaded brain tissue. However, glioma cells migrate preferentially along white matter tracts, suggesting that white matter directionality should be considered in an anisotropic CTV expansion. We investigate two models of anisotropic CTV expansion and evaluate their clinical feasibility.Approach.To incorporate white matter directionality into the CTV, a diffusion tensor imaging (DTI) atlas is used. The DTI atlas consists of water diffusion tensors that are first spatially transformed into local tumor resistance tensors, also known as metric tensors, and secondly fed to a CTV expansion algorithm to generate anisotropic CTVs. Two models of spatial transformation are considered in the first step. The first model assumes that tumor cells experience reduced resistance parallel to the white matter fibers. The second model assumes that the anisotropy of tumor cell resistance is proportional to the anisotropy observed in DTI, with an 'anisotropy weighting parameter' controlling the proportionality. The models are evaluated in a cohort of ten brain tumor patients.Main results.To evaluate the sensitivity of the model, a library of model-generated CTVs was computed by varying the resistance and anisotropy parameters. Our results indicate that the resistance coefficient had the most significant effect on the global shape of the CTV expansion by redistributing the target volume from potentially less involved gray matter to white matter tissue. In addition, the anisotropy weighting parameter proved useful in locally increasing CTV expansion in regions characterized by strong tissue directionality, such as near the corpus callosum.Significance.By incorporating anisotropy into the CTV expansion, this study is a step toward an interactive CTV definition that can assist physicians in incorporating neuroanatomy into a clinically optimized CTV.

Embedding machine learning based toxicity models within radiotherapy treatment plan optimization.

Objective.This study addresses radiation-induced toxicity (RIT) challenges in radiotherapy (RT) by developing a personalized treatment planning framework. It leverages patient-specific data and dosimetric information to create an optimization model that limits adverse side effects using constraints learned from historical data.Approach.The study uses the optimization with constraint learning (OCL) framework, incorporating patient-specific factors into the optimization process. It consists of three steps: optimizing the baseline treatment plan using population-wide dosimetric constraints; training a machine learning (ML) model to estimate the patient's RIT for the baseline plan; and adapting the treatment plan to minimize RIT using ML-learned patient-specific constraints. Various predictive models, including classification trees, ensembles of trees, and neural networks, are applied to predict the probability of grade 2+ radiation pneumonitis (RP2+) for non-small cell lung (NSCLC) cancer patients three months post-RT. The methodology is assessed with four high RP2+ risk NSCLC patients, with the goal of optimizing the dose distribution to constrain the RP2+ outcome below a pre-specified threshold. Conventional and OCL-enhanced plans are compared based on dosimetric parameters and predicted RP2+ risk. Sensitivity analysis on risk thresholds and data uncertainty is performed using a toy NSCLC case.Main results.Experiments show the methodology's capacity to directly incorporate all predictive models into RT treatment planning. In the four patients studied, mean lung dose and V20 were reduced by an average of 1.78 Gy and 3.66%, resulting in an average RP2+ risk reduction from 95% to 42%. Notably, this reduction maintains tumor coverage, although in two cases, sparing the lung slightly increased spinal cord max-dose (0.23 and 0.79 Gy).Significance.By integrating patient-specific information into learned constraints, the study significantly reduces adverse side effects like RP2+ without compromising target coverage. This unified framework bridges the gap between predicting toxicities and optimizing treatment plans in personalized RT decision-making.

Validation of dual-energy CT-based composition analysis using fresh animal tissues and composition-optimized tissue equivalent samples.

Objective Proton therapy allows for highly conformal dose deposition, but is sensitive to range uncertainties. Several approaches currently under development measure composition-dependent secondary radiation to monitor the delivered proton range in-vivo. To fully utilize these methods, an estimate of the elemental composition of the patient's tissue is often needed. Approach A published dual-energy computed tomography (DECT)-based composition-extraction algorithm was validated against reference compositions obtained with two independent methods. For this purpose, a set of phantoms containing either fresh porcine tissue or tissue-mimicking samples with known, realistic compositions were imaged with a CT scanner at two different energies. Then, the prompt gamma-ray (PG) signal during proton irradiation was measured with a PG detector prototype. The PG workflow used pre-calculated Monte Carlo simulations to obtain an optimized estimate of the sample's carbon and oxygen contents. The compositions were also assessed with combustion analysis (CCA), and the stopping-power ratio (SPR) was measured with a multi-layer ionization chamber. The DECT images were used to calculate SPR-, density- and elemental composition maps, and to assign voxel-wise compositions from a selection of human tissues. For a more comprehensive set of reference compositions, the original selection was extended by 135 additional tissues, corresponding to spongiosa, high-density bones and low-density tissues. Results The root-mean-square error for the soft tissue carbon and oxygen content was 8.5wt% and 9.5wt% relative to the CCA result and 2.1wt% and 10.3wt% relative to the PG result. The phosphorous and calcium content were predicted within 0.4wt% and 1.1wt% of the CCA results, respectively. The largest discrepancies were encountered in samples whose composition deviated the most from tabulated compositions or that were more inhomogeneous. Significance Overall, DECT-based composition estimations of relevant elements were in equal or better agreement with the ground truth than the established SECT-approach and could contribute to in-vivo dose verification measurements.

Breast imaging with an ultra-low field MRI scanner: a pilot study.

Breast cancer screening is necessary to reduce mortality due to undetected breast cancer. Current methods have limitations, and as a result many women forego regular screening. Magnetic resonance imaging (MRI) can overcome most of these limitations, but access to conventional MRI is not widely available for routine annual screening. Here, we used an MRI scanner operating at ultra-low field (ULF) to image the left breasts of 11 women (mean age, 35 years ±13 years) in the prone position. Three breast radiologists reviewed the imaging and were able to discern the breast outline and distinguish fibroglandular tissue (FGT) from intramammary adipose tissue. Additionally, the expert readers agreed on their assessment of the breast tissue pattern including fatty, scattered FGT, heterogeneous FGT, and extreme FGT. This preliminary work demonstrates that ULF breast MRI is feasible and may be a potential option for comfortable, widely deployable, and low-cost breast cancer diagnosis and screening.

Feasibility of Using Distal Endpoints for In-room PET Range Verification of Proton Therapy.

In an effort to verify the dose delivery in proton therapy, Positron Emission Tomography (PET) scans have been employed to measure the distribution of ß+ radioactivity produced from nuclear reactions of the protons with native nuclei. Because the dose and PET distributions are difficult to compare directly, the range verification is currently carried out by comparing measured and Monte Carlo (MC) simulation predicted PET distributions. In order to reduce the reliance on MC, simulated PET (simPET) and dose distal endpoints were compared to explore the feasibility of using distal endpoints for in-room PET range verification. MC simulations were generated for six head and neck patients with corrections for radiological decay, biological washout, and PET resolution. One-dimensional profiles of the dose and simPET were examined along the direction of the beam and covering the cross section of the beam. The chosen endpoints of the simPET (x-intercept of the linear fit to the distal falloff) and planned dose (20-50% of maximum dose) correspond to where most of the protons are below the threshold energy for the nuclear reactions. The difference in endpoint range between the distal surfaces of the dose and MC-PET were compared and the spread of range differences were assessed. Among the six patients, the mean difference between MC-PET and dose depth was found to be -1.6 mm to +0.5 mm between patients, with a standard deviation of 1.1 to 4.0 mm across the individual beams. In clinical practice, regions with deviations beyond the safety margin need to be examined more closely and can potentially lead to adjustments to the treatment plan.

Our journeys through the fascinating world of proton radiation therapy.

The purpose of this article is to share the excitement of the science of proton therapy, told by two physicists, who started their career in this area at different times. The authors' journey spans the evolution of proton therapy over the past 30 years, taking the reader from the time when it was an extremely exotic treatment modality until its more common use today. Over this time period, the authors' research and development aimed at an improved understanding of the physical benefits of intensity-modulated proton therapy and arc therapy, treatment planning and optimization to take proton-specific uncertainties into account, and imaging to measure the proton range in the patient. The final section focuses on emerging themes to democratize proton therapy by substantially reducing its size and price, for much greater affordability and global availability of this treatment modality.

Exploring trade-offs in treatment planning for brain tumor cases with a probabilistic definition of the clinical target volume.

PURPOSE: This study demonstrates how a novel probabilistic clinical target volume (CTV) concept-the clinical target distribution (CTD)-can be used to navigate the trade-off between target coverage and organ sparing with a semi-interactive treatment planning approach. METHODS: Two probabilistic treatment planning methods are presented that use tumor probabilities to balance tumor control with organ-at-risk (OAR) sparing. The first method explores OAR dose reduction by systematically discarding x % $x\%$ of CTD voxels with an unfavorable dose-to-probability ratio from the minimum dose coverage objective. The second method sequentially expands the target volume from the GTV edge, calculating the CTD coverage versus OAR sparing trade-off after dosing each expansion. Each planning method leads to estimated levels of tumor control under specific statistical models of tumor infiltration: an independent tumor islets model and contiguous circumferential tumor growth model. The methods are illustrated by creating proton therapy treatment plans for two glioblastoma patients with the clinical goal of sparing the hippocampus and brainstem. For probabilistic plan evaluation, the concept of a dose-expected-volume histogram is introduced, which plots the dose to the expected tumor volume ⟨ v ⟩ $\langle v \rangle$ considering tumor probabilities. RESULTS: Both probabilistic planning approaches generate a library of treatment plans to interactively navigate the planning trade-offs. In the first probabilistic approach, a significant reduction of hippocampus dose could be achieved by excluding merely 1% of CTD voxels without compromising expected tumor control probability (TCP) or CTD coverage: the hippocampus D 2 % $D_{2\%}$ dose reduces with 9.5 and 5.3 Gy for Patient 1 and 2, while the TCP loss remains below 1%. Moreover, discarding up to 10% of the CTD voxels does not significantly diminish the expected CTD dose, even though evaluation with a binary volume suggests poor CTD coverage. In the second probabilistic approach, the expected CTD D ⟨ 98 % ⟩ $D_{\langle 98\%\rangle }$ and TCP depend more strongly on the extent of the high-dose region: the target volume margin cannot be reduced by more than 2 mm if one aims at keeping the expected CTD D ⟨ 98 % ⟩ $D_{\langle 98\%\rangle }$ loss and TCP loss under 1 Gy and 2%, respectively. Therefore, there is less potential for improved OAR sparing without compromising TCP or expected CTD coverage. CONCLUSIONS: This study proposes and implements treatment planning strategies to explore trade-offs using tumor probabilities.

Recurrence analysis of glioblastoma cases based on distance and dose information.

BACKGROUND AND PURPOSE: Radiation therapy for glioblastoma (GBM) typically involves large target volumes. The aim of this study was to examine the recurrence pattern of GBM following modern radiochemotherapy according to EORTC guidelines and provide dose and distance information for the choice of optimal target volume margins. MATERIALS AND METHODS: In this study, the recurrences of 97 GBM patients, treated with radiochemotherapy from 2013 to 2017 at the Medical Center- University of Freiburg, Germany were analysed. Dose and distance based metrices were used to derive recurrence patterns. RESULTS: The majority of recurrences (75%) occurred locally within the primary tumor area. Smaller GTVs had a higher rate of distant recurrences. Larger treated volumes did not show a clinical benefit regarding progression free and overall survival. CONCLUSION: The identified recurrence pattern suggests that adjustments or reductions in target volume margins are feasible and could result in similar survival rates, potentially combined with a lower risk of side effects.

Online adaptive planning methods for intensity-modulated radiotherapy.

Online adaptive radiation therapy aims at adapting a patient's treatment plan to their current anatomy to account for inter-fraction variations before daily treatment delivery. As this process needs to be accomplished while the patient is immobilized on the treatment couch, it requires time-efficient adaptive planning methods to generate a quality daily treatment plan rapidly. The conventional planning methods do not meet the time requirement of online adaptive radiation therapy because they often involve excessive human intervention, significantly prolonging the planning phase. This article reviews the planning strategies employed by current commercial online adaptive radiation therapy systems, research on online adaptive planning, and artificial intelligence's potential application to online adaptive planning.

Automated Segmentation of Sacral Chordoma and Surrounding Muscles Using Deep Learning Ensemble.

PURPOSE: The manual segmentation of organ structures in radiation oncology treatment planning is a time-consuming and highly skilled task, particularly when treating rare tumors like sacral chordomas. This study evaluates the performance of automated deep learning (DL) models in accurately segmenting the gross tumor volume (GTV) and surrounding muscle structures of sacral chordomas. METHODS AND MATERIALS: An expert radiation oncologist contoured 5 muscle structures (gluteus maximus, gluteus medius, gluteus minimus, paraspinal, piriformis) and sacral chordoma GTV on computed tomography images from 48 patients. We trained 6 DL auto-segmentation models based on 3-dimensional U-Net and residual 3-dimensional U-Net architectures. We then implemented an average and an optimally weighted average ensemble to improve prediction performance. We evaluated algorithms with the average and standard deviation of the volumetric Dice similarity coefficient, surface Dice similarity coefficient with 2- and 3-mm thresholds, and average symmetric surface distance. One independent expert radiation oncologist assessed the clinical viability of the DL contours and determined the necessary amount of editing before they could be used in clinical practice. RESULTS: Quantitatively, the ensembles performed the best across all structures. The optimal ensemble (volumetric Dice similarity coefficient, average symmetric surface distance) was (85.5 ± 6.4, 2.6 ± 0.8; GTV), (94.4 ± 1.5, 1.0 ± 0.4; gluteus maximus), (92.6 ± 0.9, 0.9 ± 0.1; gluteus medius), (85.0 ± 2.7, 1.1 ± 0.3; gluteus minimus), (92.1 ± 1.5, 0.8 ± 0.2; paraspinal), and (78.3 ± 5.7, 1.5 ± 0.6; piriformis). The qualitative evaluation suggested that the best model could reduce the total muscle and tumor delineation time to a 19-minute average. CONCLUSIONS: Our methodology produces expert-level muscle and sacral chordoma tumor segmentation using DL and ensemble modeling. It can substantially augment the streamlining and accuracy of treatment planning and represents a critical step toward automated delineation of the clinical target volume in sarcoma and other disease sites.

Validation of an MR-based multimodal method for molecular composition and proton stopping power ratio determination using ex vivo animal tissues and tissue-mimicking phantoms.

Objective. Range uncertainty in proton therapy is an important factor limiting clinical effectiveness. Magnetic resonance imaging (MRI) can measure voxel-wise molecular composition and, when combined with kilovoltage CT (kVCT), accurately determine mean ionization potential (Im), electron density, and stopping power ratio (SPR). We aimed to develop a novel MR-based multimodal method to accurately determine SPR and molecular compositions. This method was evaluated in tissue-mimicking andex vivoporcine phantoms, and in a brain radiotherapy patient.Approach. Four tissue-mimicking phantoms with known compositions, two porcine tissue phantoms, and a brain cancer patient were imaged with kVCT and MRI. Three imaging-based values were determined: SPRCM(CT-based Multimodal), SPRMM(MR-based Multimodal), and SPRstoich(stoichiometric calibration). MRI was used to determine two tissue-specific quantities of the Bethe Bloch equation (Im, electron density) to compute SPRCMand SPRMM. Imaging-based SPRs were compared to measurements for phantoms in a proton beam using a multilayer ionization chamber (SPRMLIC).Main results. Root mean square errors relative to SPRMLICwere 0.0104(0.86%), 0.0046(0.45%), and 0.0142(1.31%) for SPRCM, SPRMM, and SPRstoich, respectively. The largest errors were in bony phantoms, while soft tissue and porcine tissue phantoms had <1% errors across all SPR values. Relative to known physical molecular compositions, imaging-determined compositions differed by approximately ≤10%. In the brain case, the largest differences between SPRstoichand SPRMMwere in bone and high lipids/fat tissue. The magnitudes and trends of these differences matched phantom results.Significance. Our MR-based multimodal method determined molecular compositions and SPR in various tissue-mimicking phantoms with high accuracy, as confirmed with proton beam measurements. This method also revealed significant SPR differences compared to stoichiometric kVCT-only calculation in a clinical case, with the largest differences in bone. These findings support that including MRI in proton therapy treatment planning can improve the accuracy of calculated SPR values and reduce range uncertainties.

Consensus guide on CT-based prediction of stopping-power ratio using a Hounsfield look-up table for proton therapy.

BACKGROUND AND PURPOSE: Studies have shown large variations in stopping-power ratio (SPR) prediction from computed tomography (CT) across European proton centres. To standardise this process, a step-by-step guide on specifying a Hounsfield look-up table (HLUT) is presented here. MATERIALS AND METHODS: The HLUT specification process is divided into six steps: Phantom setup, CT acquisition, CT number extraction, SPR determination, HLUT specification, and HLUT validation. Appropriate CT phantoms have a head- and body-sized part, with tissue-equivalent inserts in regard to X-ray and proton interactions. CT numbers are extracted from a region-of-interest covering the inner 70% of each insert in-plane and several axial CT slices in scan direction. For optimal HLUT specification, the SPR of phantom inserts is measured in a proton beam and the SPR of tabulated human tissues is computed stoichiometrically at 100 MeV. Including both phantom inserts and tabulated human tissues increases HLUT stability. Piecewise linear regressions are performed between CT numbers and SPRs for four tissue groups (lung, adipose, soft tissue, and bone) and then connected with straight lines. Finally, a thorough but simple validation is performed. RESULTS: The best practices and individual challenges are explained comprehensively for each step. A well-defined strategy for specifying the connection points between the individual line segments of the HLUT is presented. The guide was tested exemplarily on three CT scanners from different vendors, proving its feasibility. CONCLUSION: The presented step-by-step guide for CT-based HLUT specification with recommendations and examples can contribute to reduce inter-centre variations in SPR prediction.

Multicriteria VMAT optimization.

PURPOSE: To make the planning of volumetric modulated arc therapy (VMAT) faster and to explore the tradeoffs between planning objectives and delivery efficiency. METHODS: A convex multicriteria dose optimization problem is solved for an angular grid of 180 equi-spaced beams. This allows the planner to navigate the ideal dose distribution Pareto surface and select a plan of desired target coverage versus organ at risk sparing. The selected plan is then made VMAT deliverable by a fluence map merging and sequencing algorithm, which combines neighboring fluence maps based on a similarity score and then delivers the merged maps together, simplifying delivery. Successive merges are made as long as the dose distribution quality is maintained. The complete algorithm is called VMERGE. RESULTS: VMERGE is applied to three cases: a prostate, a pancreas, and a brain. In each case, the selected Pareto-optimal plan is matched almost exactly with the VMAT merging routine, resulting in a high quality plan delivered with a single arc in less than 5 min on average. CONCLUSIONS: VMERGE offers significant improvements over existing VMAT algorithms. The first is the multicriteria planning aspect, which greatly speeds up planning time and allows the user to select the plan, which represents the most desirable compromise between target coverage and organ at risk sparing. The second is the user-chosen epsilon-optimality guarantee of the final VMAT plan. Finally, the user can explore the tradeoff between delivery time and plan quality, which is a fundamental aspect of VMAT that cannot be easily investigated with current commercial planning systems.