Dose perturbations from implanted helical gold markers in proton therapy of prostate cancer

Proton Beam Therapy for Pancreatic Tumors: A Consensus Statement from the Particle Therapy Cooperative Group Gastrointestinal Subcommittee

Radiation therapy manages pancreatic cancer in various settings; however, the proximity of gastrointestinal (GI) luminal organs at risk (OARs) poses challenges to conventional radiation therapy. Proton beam therapy (PBT) may reduce toxicities compared to photon therapy. This consensus statement summarizes PBT's safe and optimal delivery for pancreatic tumors. Our group has specific expertise using PBT for GI indications and has developed expert recommendations for treating pancreatic tumors with PBT. Computed tomography (CT) simulation: Patients should be simulated supine (arms above head) with custom upper body immobilization. For stomach/duodenum filling consistency, patients should restrict oral intake within 3 hours before simulation/treatments. Fiducial markers may be implanted for image guidance; however, their design and composition require scrutiny. The reconstruction field-of-view should encompass all immobilization devices at the target level (CT slice thickness 2-3 mm). Four-dimensional CT should quantify respiratory motion and guide motion mitigation. Respiratory gating is recommended when motion affects OAR sparing or reduces target coverage. Treatment planning: Beam-angle selection factors include priority OAR-dose minimization, water-equivalent-thickness stability along the beam path, and enhanced relative biological effect consideration due to the increased linear energy transfer at the proton beam end-of-range. Posterior and right-lateral beam angles that avoid traversing GI luminal structures are preferred (minimizing dosimetric impacts of variable anatomies). Pencil beam scanning techniques should use robust optimization. Single-field optimization is preferable to increase robustness, but if OAR constraints cannot be met, multifield optimization may be used. Treatment delivery: Volumetric image guidance should be used daily. CT scans should be acquired ad hoc as necessary (at minimum every other week) to assess the dosimetric impacts of anatomy changes. Adaptive replanning should be performed as required. Our group has developed recommendations for delivering PBT to safely and effectively manage pancreatic tumors.

A Novel Polymer-Encapsulated Multi-Imaging Modality Fiducial Marker with Positive Signal Contrast for Image-Guided Radiation Therapy

BACKGROUND

Current fiducial markers (FMs) in external-beam radiotherapy (EBRT) for prostate cancer (PCa) cannot be positively visualized on magnetic resonance imaging (MRI) and create dose perturbation and significant imaging artifacts on computed tomography (CT) and MRI. We report our initial experience with clinical imaging of a novel multimodality FM, NOVA.

METHODS

We tested Gold Anchor [G-FM], BiomarC [carbon, C-FM], and NOVA FMs in phantoms imaged with kilovoltage (kV) X-rays, transrectal ultrasound (TRUS), CT, and MRI. Artifacts of the FMs on CT were quantified by the relative streak artifacts level (rSAL) metric. Proton dose perturbations (PDPs) were measured with Gafchromic EBT3 film, with FMs oriented either perpendicular to or parallel with the beam axis. We also tested the performance of NOVA-FMs in a patient.

RESULTS

NOVA-FMs were positively visualized on all 4 imaging modalities tested. The rSAL on CT was 0.750 ± 0.335 for 2-mm reconstructed slices. In F-tests, PDP was associated with marker type and depth of measurement (p < 10-6); at 5-mm depth, PDP was significantly greater for the G-FM (12.9%, p = 10-6) and C-FM (6.0%, p = 0.011) than NOVA (4.5%). EBRT planning with MRI/CT image co-registration and daily alignments using NOVA-FMs in a patient was feasible and reproducible.

CONCLUSIONS

NOVA-FMs were positively visible and produced less PDP than G-FMs or C-FMs. NOVA-FMs facilitated MRI/CT fusion and identification of regions of interest.

Radioactivation effects of titanium caused by clinical proton beam: a simulation study

Objective. In proton beam therapy (PBT), metals in the patient body perturb the dose distribution, and their radioactivation may affect the dose distribution around the metal; however, the radioactivation effect has been not clarified with PBT. In this study, we aimed to evaluate the radioactivation effect of metal depending on proton energies and secondary neutrons with a clinical proton beam using a Monte Carlo (MC) simulation.Approach.The radionuclides produced from a titanium alloy (Ti-6Al-4V) and their radioactivity were calculated using a 210-MeV passive scattering proton beam with a 60-mm Spread-out Bragg Peak, and the deposited doses caused by the radioactivation were computed using the MC simulation. The position of metal was changed according to the proton mean energy in water. To assess neutron effects on the radioactivation, we calculated the radioactivation in following three situations: (i) full MC simulation with neutrons, (ii) simulation without secondary neutrons generated from the beamline components, and (iii) simulation without any secondary neutrons.Main results.Immediately after the irradiation, the radionuclide with the largest activity was Sc-45 m (half-life of 318 ms) regardless of the proton energy and the presence of neutrons. Total radioactivity tended to increase according to the proton energy. The accumulated dose for 24 h caused by the metal activation showed an increasing trend with the proton energy, with a maximum increase rate of 0.045% to the prescribed dose. The accumulated dose at a distance of 10 mm from the metal was lower than 1/10 of that at a distance of 1 mm.Significance.The radioactivation effect of the titanium was comprehensively evaluated in the clinical passive scattering proton beam. We expect that radioactivation effects on the clinical dose distribution would be small. We consider that these results will help the clinical handling of high-Z metals in PBT.

Improved sub-milimeter range-verification method for proton therapy using a composite hadron tumour marker (HTM)

Objective.The results of a follow-up experiment investigating a novel method for sub-milimetre range verification (RV) in proton therapy (PT) are presented.Approach.The method consists of implanting a hadron tumour marker (HTM) near the planned treatment volume, and measuring theγ-ray signals emitted as a result of activation by the proton beam. These signals are highly correlated with the energy of the beam impinging on the HTM and can provide an absolute measurement of the range of the beam relative to the position of the HTM, which is independent of any uncertainties in beam delivery.Main results.Three candidate HTM materials were identified and combined into a single composite HTM, which makes use of the strongest reaction in each material. The setup of the previous experiment was improved on by using high-purity germanium detectors to measure theγ-ray signal with a higher resolution than was previously achieved. A PMMA phantom was also used to simulate theγ-ray background from tissue activation. HTM RV using the data collected in this study yielded range measurements whose average deviation from the expected value was 0.13(22)mm.Significance.Range uncertainty in PT limits the prescribed treatment plan for cancer patients with large safety margins and constrains the direction of the proton beam in relation to any organ at risk. The sub-milimetre range uncertainty achieved in this study using HTM RV, if implemented clinically, would allow for a reduction in the size of safety margins, increasing the therapeutic window for PT.

Experimental Comparison of Fiducial Markers Used in Proton Therapy: Study of Different Imaging Modalities and Proton Fluence Perturbations Measured With CMOS Pixel Sensors

Fiducial markers are used for image guidance to verify the correct positioning of the target for the case of tumors that can suffer interfractional motion during proton therapy. The markers should be visible on daily imaging, but at the same time, they should produce minimal streak artifacts in the CT scans for treatment planning and induce only slight dose perturbations during particle therapy. In this work, these three criteria were experimentally investigated at the Heidelberg Ion Beam Therapy Center. Several small fiducial markers with different geometries and materials (gold, platinum, and carbon-coated ZrO2) were evaluated. The streak artifacts on treatment planning CT were measured with and without iMAR correction, showing significantly smaller artifacts from markers lighter than 6 mg and a clear improvement with iMAR correction. Daily imaging as X-ray projections and in-room mobile CT were also performed. Markers heavier than 6 mg showed a better contrast in the X-ray projections, whereas on the images from the in-room mobile CT, all markers were clearly visible. In the other part of this work, fluence perturbations of proton beams were measured for the same markers by using a tracker system of several high spatial resolution CMOS pixel sensors. The measurements were performed for single-energy beams, as well as for a spread-out Bragg peak. Three-dimensional fluence distributions were computed after reconstructing all particle trajectories. These measurements clearly showed that the ZrO2 markers and the low-mass gold/platinum markers (0.35mm diameter) induce perturbations being 2–3 times lower than the heavier gold or platinum markers of 0.5mm diameter. Monte Carlo simulations, using the FLUKA code, were used to compute dose distributions and showed good agreement with the experimental data after adjusting the phase space of the simulated proton beam compared to the experimental beam.

Dosimetric Effect of Biozorb Markers for Accelerated Partial Breast Irradiation in Proton Therapy

Purpose

To investigate dosimetric implications of biodegradable Biozorb (BZ) markers for proton accelerated partial breast irradiation (APBI) plans.

Materials and Methods

Six different BZs were placed within in-house breast phantoms to acquire computed tomography (CT) images. A contour correction method with proper mass density overriding for BZ titanium clip and surrounding tissue was applied to minimize inaccuracies found in the CT images in the RayStation planning system. Each breast phantom was irradiated by a monoenergetic proton beam (103.23 MeV and 8×8 cm²) using a pencil-beam scanning proton therapy system. For a range perturbation study, doses were measured at 5 depths below the breast phantoms by using an ionization chamber and compared to the RayStation calculations with 3 scenarios for the clip density: the density correction method (S1: 1.6 g/cm³), raw CT (S2), and titanium density (S3: 4.54 g/cm³). For the local dose perturbation study, the radiographic EDR2 film was placed at 0 and 2 cm below the phantoms and compared to the RayStation calculations. Clinical effects of the perturbations were retrospectively examined with 10 APBI plans for the 3 scenarios (approved by our institutional review board).

Results

In the range perturbation study, the S1 simulation showed a good agreement with the chamber measurements, while excess pullbacks of 1∼2 mm were found in the S2 and S3 simulations. The film study showed local dose shadowing and perturbation by the clips that RayStation could not predict. In the plan study, no significant differences in the plan quality were found among the 3 scenarios. However, substantial range pullbacks were observed for S3.

Conclusion

The density correction method could minimize the dose and range difference between measurement and RayStation prediction. It should be avoided to simply override the known physical density of the BZ clips for treatment planning owing to overestimation of the range pullback.

Feasibility of transrectal and transperineal fiducial marker placement for prostate cancer before proton therapy

BACKGROUND

To compare the feasibility of transrectal and transperineal fiducial marker placement for prostate cancer before proton therapy.

MATERIALS AND METHODS

From 2013 to 2015, the first 40 prostate cancer patients that were scheduled for proton therapy underwent transrectal fiducial marker placement, and the next 40 patients underwent transperineal fiducial marker placement (the first series). Technical and clinical success and pain scores were evaluated. In the second series (n = 280), the transrectal or transperineal approach was selected depending on the presence/absence of comorbidities, such as blood coagulation abnormalities. Seven patients refused to undergo the procedure. Thus, the total number of patients across both series was 353 (262 and 91 underwent the transrectal and transperineal approach, respectively). Technical and clinical success, complications, marker migration and the distance between the two markers were evaluated.

RESULTS

In the first series, the technical and clinical success rates were 100% in both groups. The transrectal group exhibited lower pain scores than the transperineal group. The overall technical success rates of the transrectal and transperineal groups were 100% (262/262) and 99% (90/91), respectively (P > 0.05). The overall clinical success rate was 100% in both groups, and there were no major complications in either group. The migration rates of the two groups did not differ significantly. The mean distance between the two markers was 25.6 ± 7.1 mm (mean ± standard deviation) in the transrectal group and 31.9 ± 5.2 mm in the transperineal group (P < 0.05).

CONCLUSION

Both the transrectal and transperineal fiducial marker placement methods are feasible and safe.

A novel proton counting detector and method for the validation of tissue and implant material maps for Monte Carlo dose calculation

The presence of artificial implants complicates the delivery of proton therapy due to inaccurate characterization of both the implant and the surrounding tissues. In this work, we describe a method to characterize implant and human tissue mimicking materials in terms of relative stopping power (RSP) using a novel proton counting detector. Each proton is tracked by directly measuring the deposited energy along the proton track using a fast, pixelated spectral detector AdvaPIX-TPX3 (TPX3). We considered three scenarios to characterize the RSPs. First, in-air measurements were made in the presence of metal rods (Al, Ti and CoCr) and bone. Then, measurements of energy perturbations in the presence of metal implants and bone in an anthropomorphic phantom were performed. Finally, sampling of cumulative stopping power (CSP) of the phantom were made at different locations of the anthropomorphic phantom. CSP and RSP information were extracted from energy spectra at each beam path. To quantify the RSP of metal rods we used the shift in the most probable energy (MPE) of CSP from the reference CSP without a rod. Overall, the RSPs were determined as 1.48, 2.06, 3.08, and 5.53 from in-air measurements; 1.44, 1.97, 2.98, and 5.44 from in-phantom measurements, for bone, Al, Ti and CoCr, respectively. Additionally, we sampled CSP for multiple paths of the anthropomorphic phantom ranging from 18.63 to 25.23 cm deriving RSP of soft tissues and bones in agreement within 1.6% of TOPAS simulations. Using minimum error of these multiple CSP, optimal mass densities were derived for soft tissue and bone and they are within 1% of vendor-provided nominal densities. The preliminary data obtained indicates the proposed novel method can be used for the validation of material and density maps, required by proton Monte Carlo Dose calculation, provided by competing multi-energy computed tomography and metal artifact reduction techniques.

Biological impact of dosimetric perturbations of a fiducial marker and the daily number of fields in proton therapy for prostate cancer

The purpose of this study was to estimate the biological impact of dosimetric perturbations of a fiducial marker and the daily number of fields in proton therapy for prostate cancer. Using a linear-quadratic model, normalized total doses (NTDs) of points where deposited dose was reduced from the prescribed dose by dosimetric perturbation of a fiducial marker were calculated in two hypothetical prostate cancer treatment schedules: a) irradiation of both parallel-opposed lateral fields and b) irradiation of alternate field in each daily treatment. The impact of hypofractionation and sublethal damage repair between irradiation on NTD was also estimated. The NTD of two fields/day schedule becomes lower than that of one field/day schedule. The difference becomes larger as dose reduction from one of two fields becomes more enhanced. The NTD reduction from the total dose in the two fields/day schedule is largest (30% of total dose) where the dose from one beam is completely lost by a fiducial marker. In contrast, the NTD reduction from the total dose in the one field/day schedule is largest (9% of total dose) where the half dose from one beam is decreased by a fiducial marker. In addition, the NTD reduction becomes larger as the fractional dose increases in a hypofractionated regimen, and when the effect of sublethal damage repair was incorporated. These influences become significant in prostate cancer since the radiobiological sensitivity α/β of prostate cancer is lower than other cancer types and normal tissues late complication. Treating with one alternate field in a daily treatment can improve a deteriorating treatment effect by dosimetric distortion of a fiducial marker in prostate cancer treatment. However, the choice of the number of beams in a fraction must also be determined by considering the sparing of normal tissues and patient-specific status.

Proton therapy range verification method via delayed γ-ray spectroscopy of a molybdenum tumour marker

In this work, a new method of range verification for proton therapy (PT) is experimentally demonstrated for the first time. If a metal marker is implanted near the tumour site, its response to proton activation will result in the emission of characteristic γ rays. The relative intensity of γ rays originating from competing fusion-evaporation reaction channels provides a unique signature of the average proton energy at the marker, and by extension the beam's range, in vivo and in real time. The clinical feasibility of this method was investigated at the PT facility at TRIUMF with a proof-of-principle experiment which irradiated a naturally-abundant molybdenum foil at various proton beam energies. Delayed characteristic γ rays were measured with two Compton-shielded LaBr₃ scintillators. The technique was successfully demonstrated by relating the relative intensity of two γ-ray peaks to the energy of the beam at the Mo target, opening the door to future clinical applications where the range of the beam can be verified in real time.

GEANT4 simulation of a range verification method using delayed γ spectroscopy of a 92Mo marker

In this work, we propose a novel technique for in-vivo proton therapy range verification. This technique makes use of a molybdenum hadron tumour marker, implanted at a short distance from the clinical treatment volume. Signals emitted from the marker during treatment can provide a direct measurement of the proton beam energy at the marker's position. Fusion-evaporation reactions between the proton beam and marker nucleus result in the emission of delayed characteristic γ rays, which are detected off-beam for an improved signal-to-noise ratio. In order to determine the viability of this technique and to establish an experimental setup for future work, the Monte Carlo package GEANT4 was used in combination with ROOT to simulate a treatment scenario with the new method outlined in this work. These simulations show that the intensity of delayed γ rays produced from competing reactions yields a precise measurement of the range of the proton beam relative to the marker, with sub-millimetre uncertainty.

The use of tumour markers in oesophageal cancer to quantify setup errors and baseline shifts during treatment

•

Implantation of solid gold markers safe.

•

Inter-fractional motion for markers in distal oesophagus largest cranio-caudally.

•

Reduced radiotherapy treatment margins with soft-tissue vs. bony-anatomy matching.

•

Impact of intra-fractional baseline shifts on margin calculation rather small.

Purpose

To prospectively evaluate the feasibility of solid gold marker placement in oesophageal cancer patients and to quantify inter-fractional and intra-fractional (baseline shift) marker motion during radiation treatment. Radiotherapy target margins and matching strategies were investigated.

Materials/methods

Thirty-four markers were implanted by echo-endoscopy in 10 patients. Patients received a planning 4D CT, daily pre-treatment cone-beam CT (CBCT) and a post-treatment CBCT for at least five fractions. For fractions with both pre- and post-treatment CBCT, marker displacement between planning CT and pre-treatment CBCT (inter-fractional) and between pre-treatment and post-treatment CBCT (intra-fractional; only for fractions without rotational treatment couch correction) were calculated in left–right (LR), cranio-caudal (CC) and anterior-posterior (AP) direction after bony-anatomy and soft-tissue matching. Systematic/random setup errors were estimated; treatment margins were calculated.

Results

No serious adverse events occurred. Twenty-three (67.6%) markers were visible during radiotherapy (n = 3 middle oesophagus, n = 16 distal oesophagus, n = 4 proximal stomach). Margins for inter-fractional displacement after bony-anatomy match depended on the localisation of the primary tumour and were 11.2 mm (LR), 16.4 mm (CC) and 8.2 mm (AP) for distal markers. Soft-tissue matching reduced the CC margin for these markers (16.4 mm to 10.5 mm). The mean intra-fractional shift of 12 distal markers was 0.4 mm (LR), 2.3 mm (CC) and 0.7 mm (AP). Inclusion of this shift resulted in treatment margins for distal markers of 12.8 mm (LR), 17.3 mm (CC) and 10.4 mm (AP) after bony-anatomy matching and 12.4 mm (LR), 11.4 mm (CC) and 9.7 mm (AP) after soft-tissue matching.

Conclusion

This study demonstrated that the implantation of gold markers was safe, albeit less stable compared to other marker types. Inter-fractional motion was largest cranio-caudally for markers in the distal oesophagus, which was reduced after soft-tissue compared to bony-anatomy matching. The impact of intra-fractional baseline shifts on margin calculation was rather small.

On-line range verification for proton beam therapy using spherical ionoacoustic waves with resonant frequency

In contrast to conventional X-ray therapy, proton beam therapy (PBT) can confine radiation doses to tumours because of the presence of the Bragg peak. However, the precision of the treatment is currently limited by the uncertainty in the beam range. Recently, a unique range verification methodology has been proposed based on simulation studies that exploit spherical ionoacoustic waves with resonant frequency (SPIREs). SPIREs are emitted from spherical gold markers in tumours initially introduced for accurate patient positioning when the proton beam is injected. These waves have a remarkable property: their amplitude is linearly correlated with the residual beam range at the marker position. Here, we present proof-of-principle experiments using short-pulsed proton beams at the clinical dose to demonstrate the feasibility of using SPIREs for beam-range verification with submillimetre accuracy. These results should substantially contribute to reducing the range uncertainty in future PBT applications.

Dose perturbation caused by metallic port in breast tissue expander in proton beam therapy

Proton beam treatment is being looked favourably now in breast treatment. Tissue expanders are often placed after mastectomy that contains metallic port for saline injection which produces dose perturbations in proton beam therapy with uncertain dosimetry. Dose perturbation for a stainless-steel injection port from a breast implant is investigated in this study. Measurements, Monte-Carlo simulation, and calculated dose distribution of plans based on kVCT and MVCT images are compared. Treatment plans are performed on kVCT and MVCT images to observe the effect of metal artifact from the breast implant. The kVCT based plan underestimates the beam range due to the overestimated water equivalent thickness of the metal ports as a result of image degradation. Compared to the measurement with metal port in the proton beam, the MVCT-based treatment planning provides more accurate dose calculation than the kVCT-based results. The dose perturbation factor calculated from MVCT planning is within 10% of the measurement results while HU corrected kVCT plan still shows dose difference as large as 100% due to the incorrect range pull back calculation caused by the misrepresentation of the volume between the plastic cap and the stainless-steel base. The dose enhancement observed at the metal and solid water interface is as large as 15%, which needs to be accounted for in the planning process if there is a clinical concern. Dose reduction as large as 16% is observed with depth from 1 cm to 4 cm underneath the thickest part of the metallic port whereas lateral dose perturbation is also seen up to 7 mm. The measurement data are supported by the Monte-Carlo simulated results with a maximum dose difference of 6%. It is concluded that if proton beam is used with metallic port, MVCT imaging data is recommended. In lieu of MVCT, DECT, CT scanner with metal artifact reduction software or in the very least, extended HU range should be used to reduce the streaking artifact as well as to produce a more accurate image of the metallic port.

Monte Carlo calculations of radiotherapy dose in "homogeneous" anatomy

Given the substantial literature on the use of Monte Carlo (MC) simulations to verify treatment planning system (TPS) calculations of radiotherapy dose in heterogeneous regions, such as head and neck and lung, this study investigated the potential value of running MC simulations of radiotherapy treatments of nominally homogeneous pelvic anatomy. A pre-existing in-house MC job submission and analysis system, built around BEAMnrc and DOSXYZnrc, was used to evaluate the dosimetric accuracy of a sample of 12 pelvic volumetric arc therapy (VMAT) treatments, planned using the Varian Eclipse TPS, where dose was calculated with both the Analytical Anisotropic Algorithm (AAA) and the Acuros (AXB) algorithm. In-house TADA (Treatment And Dose Assessor) software was used to evaluate treatment plan complexity, in terms of the small aperture score (SAS), modulation index (MI) and a novel exposed leaf score (ELS/ELA). Results showed that the TPS generally achieved closer agreement with the MC dose distribution when treatments were planned for smaller (single-organ) targets rather than larger targets that included nodes or metastases. Analysis of these MC results with reference to the complexity metrics indicated that while AXB was useful for reducing dosimetric uncertainties associated with density heterogeneity, the residual TPS dose calculation uncertainties resulted from treatment plan complexity and TPS model simplicity. The results of this study demonstrate the value of using MC methods to recalculate and check the dose calculations provided by commercial radiotherapy TPSs, even when the treated anatomy is assumed to be comparatively homogeneous, such as in the pelvic region.

An avoidance method to minimize dose perturbation effects in proton pencil beam scanning treatment of patients with small high-Z implants

To implement a multi-field-optimization (MFO) technique for treating patients with high-Z implants in pencil beam scanning proton-therapy and generate treatment plans that avoids small implants. Two main issues were addressed: (i) the assessment of the optimal CT acquisition and segmentation technique to define the dimension of the implant and (ii) the distance of pencil beams from the implant (avoidance margin) to assure that it does not affect dose distribution. Different CT reconstruction protocols (by O-MAR or standard reconstruction and by 12 bit or 16 bit dynamic range) followed by thresholding segmentation were tested on a phantom with lead spheres of different sizes. The proper avoidance margin was assessed on a dedicated phantoms of different materials (copper/tantalum and lead), shape (square slabs and spheres) and detectors (two-dimensional array chamber and radio-chromic films). The method was then demonstrated on a head-and-neck carcinoma patient, who underwent carotid artery embolization with a platinum coil close to the target. Regardless the application of O-MAR reconstruction, the CT protocol with a full 16 bit dynamic range allowed better estimation of the sphere volumes, with maximal error around -5% in the greater sphere only. Except the configuration with a shallow target (which required a pre-absorber), particularly with a retracted snout, an avoidance margin of around 0.9-1.3 cm allowed to keep the difference between planned and measured dose below 5-10%. The patient plan analysis showed adequate plan quality and confirmed effective implant avoidance. Potential target under-dosage can be produced by patient misalignment, which could be minimized by daily alignment on the implant, identifiable on orthogonal kilovolt images. By implant avoidance MFO it was possible to minimize potential dose perturbation effects produced by small high-Z implants. An advantage of such approach lies in its potential applicability for any type of implant, regardless the precise knowledge of its composition.

Fluence perturbation from fiducial markers due to edge-scattering measured with pixel sensors for 12C ion beams

Fiducial markers are nowadays a common tool for patient positioning verification before radiotherapy treatment. These markers should be visible on x-ray projection imaging, produce low streak artifacts on CTs and induce small dose perturbations due to edge-scattering effects during the ion-beam therapy treatment. In this study, the latter effect was investigated and the perturbations created by the markers were evaluated with a new measurement method using a tracker system composed of six CMOS pixel sensors. The present method enables the determination of the particle trajectory before and after the target. The experiments have been conducted at the Marburg Ion Beam Therapy Center with carbon ion beams and the measurement concept was validated by comparison with radiochromic films. This work shows that the new method is very efficient and precise to measure the perturbations due to fiducial markers with a tracker system. Three dimensional fluence distributions of all particle trajectories were reconstructed and the maximum cold spots due to the markers and their position along the beam axis were quantified. In this study, four small commercial markers with different geometries and materials (gold and carbon-coated ZrO₂) were evaluated. The gold markers showed stronger perturbations than the lower density ones. However, it is important to consider that low density and low atomic number fiducial markers are not always visible on x-ray projections.

Dosimetric evaluation of MR-derived synthetic-CTs for MR-only proton treatment planning

PURPOSE

To evaluate proton dose calculation accuracy of optimized pencil beam scanning (PBS) plans on MR-derived synthetic-CTs for prostate patients.

MATERIAL AND METHODS

Ten patient datasets with both a CT and an MRI were planned with opposed lateral proton beams optimized to single field uniform dose under an IRB-approved study. The proton plans were created on CT datasets generated by a commercial synthetic CT-based software called MRCAT (MR for Calculating ATtenuation) routinely used in our clinic for photon-based MR-only planning. A standard prescription of 79.2 Gy (RBE) and 68.4 Gy (RBE) was used for intact prostate and prostate bed cases, respectively. Proton plans were first generated and optimized using the MRCAT synthetic-CT (syn-CT), and then recalculated on the planning CT rigidly aligned with the syn-CT (aligned-CT) and a deformed planning CT (deformed-CT), which was deformed to match outer contour between syn-CT and aligned-CT. The same beam arrangement, total MUs, MUs/spot, spot positions were used to recalculate dose on deformed-CT and aligned-CT without renormalization. DVH analysis was performed on aligned-CT, deformed-CT, and syn-CT to compare D_98%, V_100%, V_95% for PTV, PTVeval, and GTV as well as V_70Gy, V_50Gy for OARs.

RESULTS

The relative percentage dose difference between syn-CT and deformed-CT, were (0.17 ± 0.33 %) for PTVeval D_98% and (0.07 ± 0.1 %) for CTV D_98%. Rectum V_70Gy, V_50Gy, and Bladder V_70Gy were (2.76 ± 4.01 %), (11.6 ± 11.2 %), and (3.41 ± 2.86 %), respectively for the syn-CT, and (3.23 ± 3.63 %), (11.3 ± 8.18 %), and (3.29 ± 2.76 %), respectively for the deformed-CT, and (1.37 ± 1.84 %), (8.48 ± 6.67 %), and (4.91 ± 3.65 %), respectively for aligned-CT.

CONCLUSION

Dosimetric analysis shows that MR-only proton planning is feasible using syn-CT based on current clinical margins that account for a range uncertainty.

Detectability and structural stability of a liquid fiducial marker in fresh ex vivo pancreas tumour resection specimens on CT and 3T MRI

PURPOSE

To test the detectability of a liquid fiducial marker injected into ex vivo pancreas tumour tissue on magnetic resonance imaging (MRI) and computed tomography (CT). Furthermore, its injection performance using different needle sizes and its structural stability after fixation in formaldehyde were investigated.

METHODS

Liquid fiducial markers with a volume of 20-100 µl were injected into freshly resected pancreas specimens of three patients with suspected adenocarcinoma. X‑ray guided injection was performed using different needle sizes (18 G, 22 G, 25 G). The specimens were scanned on MRI and CT with clinical protocols. The markers were segmented on CT by signal thresholding. Marker detectability in MRI was assessed in the registered segmentations. Marker volume on CT was compared to the injected volume as a measure of backflow.

RESULTS

Markers with a volume ≥20 µl were detected as hyperintensity on X‑ray and CT. On T₁- and T₂-weighted 3T MRI, marker sizes ranging from 20-100 µl were visible as hypointensity. Since most markers were non-spherical, MRI detectability was poor and their differentiation from hypointensities caused by air cavities or surgical clips was only feasible with a reference CT. Marker backflow was only observed when using an 18-G needle. A volume decrease of 6.6 ± 13.0% was observed after 24 h in formaldehyde and, with the exception of one instance, no wash-out occurred.

CONCLUSIONS

The liquid fiducial marker injected in ex vivo pancreatic resection specimen was visible as hyperintensity on kV X‑ray and CT and as hypointensity on MRI. The marker's size was stable in formaldehyde. A marker volume of ≥50 µL is recommended in clinically used MRI sequences. In vivo injection is expected to improve the markers sphericity due to persisting metabolism and thereby enhance detectability on MRI.

A novel range-verification method using ionoacoustic wave generated from spherical gold markers for particle-beam therapy: a simulation study

This study proposes a novel alternative range-verification method for proton beam with acoustic waves generated from spherical metal markers. When proton beam is incident on metal markers, most of the resulting pressure waves are confined in the markers because of the large difference in acoustic impedance between the metal and tissue. However, acoustic waves with frequency equal to marker’s resonant frequency escape this confinement; the marker briefly acts as an acoustic transmitter. Herein, this phenomenon is exploited to measure the range of the proton beam. We test the proposed strategy in 3-D simulations, combining the dose calculations with modelling of acoustic-wave propagation. A spherical gold marker of 2.0 mm diameter was placed in water with a 60 MeV proton beam incident on it. We investigated the dependence of pressure waves on the width of beam pulse and marker position. At short beam pulse, specific high-frequency acoustic waves of 1.62 MHz originating from the marker were observed in wave simulations, whose amplitude correlated with the distance between the marker and Bragg peak. Results indicate that the Bragg peak position can be estimated by measuring the acoustic wave amplitudes from the marker, using a single detector properly designed for the resonance frequency.

Liquid fiducial marker applicability in proton therapy of locally advanced lung cancer

BACKGROUND AND PURPOSE

We investigated the clinical applicability of a novel liquid fiducial marker (LFM) for image-guided pencil beam scanned (PBS) proton therapy (PBSPT) of locally advanced lung cancer (LALC).

MATERIALS AND METHODS

The relative proton stopping power (RSP) of the LFM was calculated and measured. Dose perturbations of the LFM and three solid markers, in a phantom, were measured. PBSPT treatment planning on computer tomography scans of five patients with LALC with the LFM implanted was performed with 1-3 fields.

RESULTS

The RSP was experimentally determined to be 1.164 for the LFM. Phantom measurements revealed a maximum relative deviation in dose of 4.8% for the LFM in the spread-out Bragg Peak, compared to 12-67% for the solid markers. Using the experimentally determined RSP, the maximum proton range error introduced by the LFM is about 1mm. If the marker was displaced at PBSPT, the maximum dosimetric error was limited to 2 percentage points for 3-field plans.

CONCLUSION

The dose perturbations introduced by the LFM were considerably smaller than the solid markers investigated. The RSP of the fiducial marker should be corrected in the treatment planning system to avoid errors. The investigated LFM introduced clinically acceptable dose perturbations for image-guided PBSPT of LALC.

A new scheme for real-time high-contrast imaging in lung cancer radiotherapy: a proof-of-concept study

Visualization of anatomy in real time is of critical importance for motion management in lung cancer radiotherapy. To achieve real-time, and high-contrast in-treatment imaging, we propose a novel scheme based on the measurement of Compton scatter photons. In our method, a slit x-ray beam along the superior-inferior direction is directed to the patient, (intersecting the lung region at a 2D plane) containing most of the tumor motion trajectory. X-ray photons are scattered off this plane primarily due to the Compton interaction. An imager with a pinhole or a slat collimator is placed at one side of the plane to capture the scattered photons. The resulting image, after correcting for incoming fluence inhomogeneity, x-ray attenuation, scatter angle variation, and outgoing beam geometry, represents the linear attenuation coefficient of Compton scattering. This allows the visualization of the anatomy on this plane. We performed Monte Carlo simulation studies both on a phantom and a patient for proof-of-principle purposes. In the phantom case, a small tumor-like structure could be clearly visualized. The contrast-resolution calculated using tumor/lung as foreground/background for kV fluoroscopy, cone beam computed tomography (CBCT), and scattering image were 0.037, 0.70, and 0.54, respectively. In the patient case, tumor motion could be clearly observed in the scatter images. Imaging dose to the voxels directly exposed by the slit beam was ~0.4 times of that under a single CBCT projection. These studies demonstrated the potential feasibility of the proposed imaging scheme to capture the instantaneous anatomy of a patient on a 2D plane with a high image contrast. Clear visualization of the tumor motion may facilitate marker-less tumor tracking.

Tumour Movement in Proton Therapy: Solutions and Remaining Questions: A Review

Movement of tumours, mostly by respiration, has been a major problem for treating lung cancer, liver tumours and other locations in the abdomen and thorax. Organ motion is indeed one component of geometrical uncertainties that includes delineation and target definition uncertainties, microscopic disease and setup errors. At present, minimising motion seems to be the easiest to implement in clinical practice. If combined with adaptive approaches to correct for gradual anatomical variations, it may be a practical strategy. Other approaches such as repainting and tracking could increase the accuracy of proton therapy delivery, but advanced 4D solutions are needed. Moreover, there is a need to perform clinical studies to investigate which approach is the best in a given clinical situation. The good news is that existing and emerging technology and treatment planning systems as will without doubt lead in the forthcoming future to practical solutions to tackle intra-fraction motion in proton therapy. These developments may also improve motion management in photon therapy as well.

The physics of proton therapy

The physics of proton therapy has advanced considerably since it was proposed in 1946. Today analytical equations and numerical simulation methods are available to predict and characterize many aspects of proton therapy. This article reviews the basic aspects of the physics of proton therapy, including proton interaction mechanisms, proton transport calculations, the determination of dose from therapeutic and stray radiations, and shielding design. The article discusses underlying processes as well as selected practical experimental and theoretical methods. We conclude by briefly speculating on possible future areas of research of relevance to the physics of proton therapy.

Depth dose perturbation by a hydrogel fiducial marker in a proton beam

The purpose of this study was to evaluate proton depth dose perturbation caused by a radio‐opaque hydrogel fiducial marker. Electronic proton stopping powers in the hydrogel were calculated for energies 0.5–250 MeV, and Monte Carlo simulations were generated of hydrogel vs. gold markers placed at various water phantom depths in a generic proton beam. Across the studied energy range, the gel/water stopping power ratio was 1.0146 to 1.0160. In the Monte Carlo simulation, the hydrogel marker caused no discernible perturbation of the proton percent depth‐dose (PDD) curve. In contrast, the gold marker caused dose reductions of as much as 20% and dose shadowing regions as long as 6.5 cm. In contrast to gold markers, the radio‐opaque hydrogel marker causes negligible proton depth dose perturbation. This factor may be taken into consideration for image‐guided proton therapy at facilities with suitable imaging modalities.

PACS number: 87.55.Qr

Dosimetric impact of gold markers implanted closely to lung tumors: a Monte Carlo simulation

We are developing an innovative dynamic tumor tracking irradiation technique using gold markers implanted around a tumor as a surrogate signal, a real‐time marker detection system, and a gimbaled X‐ray head in the Vero4DRT. The gold markers implanted in a normal organ will produce uncertainty in the dose calculation during treatment planning because the photon mass attenuation coefficient of a gold marker is much larger than that of normal tissue. The purpose of this study was to simulate the dose variation near the gold markers in a lung irradiated by a photon beam using the Monte Carlo method. First, the single‐beam and the opposing‐beam geometries were simulated using both water and lung phantoms. Subsequently, the relative dose profiles were calculated using a stereotactic body radiotherapy (SBRT) treatment plan for a lung cancer patient having gold markers along the anteriorposterior (AP) and right‐left (RL) directions. For the single beam, the dose at the gold marker‐phantom interface laterally along the perpendicular to the beam axis increased by a factor of 1.35 in the water phantom and 1.58 in the lung phantom, respectively. Furthermore, the entrance dose at the interface along the beam axis increased by a factor of 1.63 in the water phantom and 1.91 in the lung phantom, while the exit dose increased by a factor of 1.00 in the water phantom and 1.12 in the lung phantom, respectively. On the other hand, both dose escalations and dose de‐escalations were canceled by each beam for opposing portal beams with the same beam weight. For SBRT patient data, the dose at the gold marker edge located in the tumor increased by a factor of 1.30 in both AP and RL directions. In clinical cases, dose escalations were observed at the small area where the distance between a gold marker and the lung tumor was ≤ 5 mm, and it would be clinically negligible in multibeam treatments, although further investigation may be required.

PACS number: 87.10.Rt

Model-guided respiratory organ motion prediction of the liver from 2D ultrasound

With the availability of new and more accurate tumour treatment modalities such as high-intensity focused ultrasound or proton therapy, accurate target location prediction has become a key issue. Various approaches for diverse application scenarios have been proposed over the last decade. Whereas external surrogate markers such as a breathing belt work to some extent, knowledge about the internal motion of the organs inherently provides more accurate results. In this paper, we combine a population-based statistical motion model and information from 2d ultrasound sequences in order to predict the respiratory motion of the right liver lobe. For this, the motion model is fitted to a 3d exhalation breath-hold scan of the liver acquired before prediction. Anatomical landmarks tracked in the ultrasound images together with the model are then used to reconstruct the complete organ position over time. The prediction is both spatial and temporal, can be computed in real-time and is evaluated on ground truth over long time scales (5.5 min). The method is quantitatively validated on eight volunteers where the ultrasound images are synchronously acquired with 4D-MRI, which provides ground-truth motion. With an average spatial prediction accuracy of 2.4 mm, we can predict tumour locations within clinically acceptable margins.

Benchmark measurements and simulations of dose perturbations due to metallic spheres in proton beams

Monte Carlo simulations are increasingly used for dose calculations in proton therapy due to its inherent accuracy. However, dosimetric deviations have been found using Monte Carlo code when high density materials are present in the proton beam line. The purpose of this work was to quantify the magnitude of dose perturbation caused by metal objects. We did this by comparing measurements and Monte Carlo predictions of dose perturbations caused by the presence of small metal spheres in several clinical proton therapy beams as functions of proton beam range, spread-out Bragg peak width and drift space. Monte Carlo codes MCNPX, GEANT4 and Fast Dose Calculator (FDC) were used. Generally good agreement was found between measurements and Monte Carlo predictions, with the average difference within 5% and maximum difference within 17%. The modification of multiple Coulomb scattering model in MCNPX code yielded improvement in accuracy and provided the best overall agreement with measurements. Our results confirmed that Monte Carlo codes are well suited for predicting multiple Coulomb scattering in proton therapy beams when short drift spaces are involved.

Dose perturbations by electromagnetic transponders in the proton environment

Surgically implanted electromagnetic transponders have been used in external beam radiotherapy for target localization and position monitoring in real time. The effect of transponders on proton therapy dose distributions has not been reported. A Monte Carlo implementation of the transponder geometry is validated against film measurements in a proton SOBP and subsequently used to generate dose distributions for transponders at different positions and orientations in the proton SOBP. The maximum dose deficit is extracted in each case. Dose shadows of up to 60% occur for transponders positioned very near the end of range of the Bragg peak. However, if transponders are positioned further than 5 mm from the end of range, and are not oriented parallel to the beam direction, then the dose deficit can be kept below 10%.

Determination of optimal fiducial marker across image-guided radiation therapy (IGRT) modalities: visibility and artifact analysis of gold, carbon, and polymer fiducial markers

The purpose of this study was to evaluate the visibility and artifact created by gold, carbon, and polymer fiducial markers in a simple phantom across computed tomography (CT), kilovoltage (kV), and megavoltage (MV) linear accelerator imaging and MV tomotherapy imaging. Three types of fiducial markers (gold, carbon, and polymer) were investigated for their visibility and artifacts in images acquired with various modalities and with different imaging parameters (kV, mAs, slice thickness). The imaging modalities include kV CT, 2D linac‐based kilovoltage and megavoltage X‐ray imaging systems, kV cone‐beam CT, and normal and fine tomotherapy imaging. The images were acquired on a phantom constructed using Superflab bolus in which markers of each type were inserted into the center layer. The visibility and artifacts produced by each marker were assessed qualitatively and quantitatively. All tested markers could be identified clearly on the acquired CT and linac‐based kV images; gold markers demonstrated the highest contrast. On the CT images, gold markers produced a significant artifact, while no artifacts were observed with polymer markers. Only gold markers were visible when using linac‐based MV and tomotherapy imaging. For linac‐based kV images, the contrast increased with kV and mAs values for all the markers, with the gold being the most pronounced. On CT images, the contrast increased with kV for the gold markers, while decreasing for the polymer and carbon marker. With the bolus phantom used, we found that when kV imaging‐based treatment verification equipment is available, polymer and carbon markers may be the preferred choice for target localization and patient treatment positioning verification due to less image artifacts. If MV imaging will be the sole modality for positioning verification, it may be necessary to use gold markers despite the artifacts they create on the simulation CT images.

PACS number: 87

Proton dose perturbations caused by high-voltage leads from implanted cardioverter defibrillators

An increasing number of patients undergoing proton radiotherapy have cardiac implantable electrical devices (CIEDs). We recently encountered a situation in which a high‐voltage coil on a lead from an implanted cardiac defibrillator was located within the clinical treatment volume for a patient receiving proton radiotherapy for esophageal cancer. To study the effects of the lead on the dose delivery, we placed a high‐Z CIED lead at both the center and the distal edge of a clinical spread‐out Bragg peak (SOBP) in a water phantom, in both a stationary position and with the lead moving in a periodic pattern to simulate cardiorespiratory movement. We then calculated planned doses using a commercial proton treatment planning system (TPS), and compared them with the doses delivered in the phantom, measured using radiographic film. Dose profiles from TPS‐calculated and measured dose distributions showed large pertubrations in the delivered proton dose in the vicinity of the CIED lead when it was not moving. The TPS predicted perturbations up to 20% and measurements revealed perturbations up to 35%. However, the perturbations were less than 3% when the lead was moving. Greater dose perturbations were seen when the lead was placed at the distal edge of the SOBP than when it was placed in the center of the SOBP. We conclude that although cardiorespiratory motion of the lead mitigates some of the perturbations, the effects of the leads should be considered and steps taken to reduce these effects during the treatment planning process.

PACS numbers: 87.55.D‐,87.55.ne, 87.85.M

Investigation of dose perturbations and the radiographic visibility of potential fiducials for proton radiation therapy of the prostate

Image guidance using implanted fiducial markers is commonly used to ensure accurate and reproducible target positioning in radiation therapy for prostate cancer. The ideal fiducial marker is clearly visible in kV imaging, does not perturb the therapeutic dose in the target volume and does not cause any artifacts on the CT images used for treatment planning. As yet, ideal markers that fully meet all three of these criteria have not been reported. In this study, 12 fiducial markers were evaluated for their potential clinical utility in proton radiation therapy for prostate cancer. In order to identify the good candidates, each fiducial was imaged using a CT scanner as well as a kV imaging system. Additionally, the dose perturbation caused by each fiducial was quantified using radiochromic film and a clinical proton beam. Based on the results, three fiducials were identified as good candidates for use in proton radiotherapy of prostate cancer.

Dose perturbations and image artifacts caused by carbon-coated ceramic and stainless steel fiducials used in proton therapy for prostate cancer

Image-guided radiation therapy using implanted fiducial markers is a common solution for prostate localization to improve targeting accuracy. However, fiducials that are typically used for conventional photon radiotherapy cause large dose perturbations in patients who receive proton radiotherapy. A proposed solution has been to use fiducials of lower atomic number (Z) materials to minimize this effect in tissue, but the effects of these fiducials on dose distributions have not been quantified. The objective of this study was to analyze the magnitude of the dose perturbations caused by select lower-Z fiducials (a carbon-coated zirconium dioxide fiducial and a plastic-coated stainless steel fiducial) and compare them to perturbations caused by conventional gold fiducials. Sets of phantoms were used to assess select components of the effects on dose. First, the fiducials were assessed for radiographic visibility using both conventional computed tomography (CT) and an on-board kilovoltage imaging device at our proton therapy center. CT streak artifacts from the fiducials were also measured in a separate phantom. Second, dose perturbations were measured downstream of the fiducials using radiochromic film. The magnitude of dose perturbation was characterized as a function of marker material, implantation depth and orientation with respect to the beam axis. The radiographic visibility of the markers was deemed to be acceptable for clinical use. The dose measurements showed that the perpendicularly oriented zirconium dioxide and stainless steel fiducials located near the center of modulation of the proton beam perturbed the dose by less than 10%, but that the same fiducials in a parallel orientation near the end of the range of the beam could perturb the dose by as much as 38%. This suggests that carbon-coated and stainless steel fiducials could be used in proton therapy if they are located far from the end of the range of the beam and if they are oriented perpendicular to the beam axis.

Monte Carlo simulation of the neutron spectral fluence and dose equivalent for use in shielding a proton therapy vault

Neutron production is of principal concern when designing proton therapy vault shielding. Conventionally, neutron calculations are based on analytical methods, which do not accurately consider beam shaping components and nozzle shielding. The goal of this study was to calculate, using Monte Carlo modeling, the neutron spectral fluence and neutron dose equivalent generated by a realistic proton therapy nozzle and evaluate how these data could be used in shielding calculations. We modeled a contemporary passive scattering proton therapy nozzle in detail with the MCNPX simulation code. The neutron spectral fluence and dose equivalent at various locations in the treatment room were calculated and compared to those obtained from a thick iron target bombarded by parallel proton beams, the simplified geometry on which analytical methods are based. The neutron spectral fluence distributions were similar for both methods, with deeply penetrating high-energy neutrons (E > 10 MeV) being most prevalent along the beam central axis, and low-energy neutrons predominating the neutron spectral fluence in the lateral region. However, unlike the inverse square falloff used in conventional analytical methods, this study shows that the neutron dose equivalent per therapeutic dose in the treatment room decreased with distance approximately following a power law, with an exponent of about -1.63 in the lateral region and -1.73 in the downstream region. Based on the simulated data according to the detailed nozzle modeling, we developed an empirical equation to estimate the neutron dose equivalent at any location and distance in the treatment vault, e.g. for cases in which detailed Monte Carlo modeling is not feasible. We applied the simulated neutron spectral fluence and dose equivalent to a shielding calculation as an example.