PET/CT imaging for treatment verification after proton therapy: a study with plastic phantoms and metallic implants

Evaluation of Delivered Doses in Proton Beam Therapy for Prostate Cancer Using Positron Emission Tomography/Computed Tomography Imaging

AIMS

Proton beams deposit energy along their paths and stop abruptly without penetrating the opposite side, making it difficult to detect their actual paths. However, confirming the path may lead to evaluating the actual doses to organs at risk in proton therapy for prostate cancer. As proton beams produce positron emitters through nuclear fragmentation reactions, theoretically, proton beam paths can be measured by positron emission tomography/computed tomography (PET/CT). Therefore, this study investigated whether conducting PET/CT examinations immediately after proton beam therapy helps to assess the doses delivered to the rectal and urinary bladder walls, which are the major sites of radiation-related toxicity.

MATERIALS AND METHODS

Between June 2022 and June 2023, 51 consecutive patients with prostate cancer who underwent proton beam therapy were enrolled and imaged with PET/CT to measure these radioactive particles and validate the actual dose delivered to the rectal and urinary bladder walls.

RESULTS

The delivered doses assessed using PET/CT after proton beam therapy strongly correlated with the planned volume for proton beam treatment.

CONCLUSIONS

PET/CT exhibited potential as a valuable tool for validating the irradiated dose to organs at risk.

A torus source and its application for non-primary radiation evaluation

Objective. Non-primary radiation doses to normal tissues from proton therapy may be associated with an increased risk of secondary malignancies, particularly in long-term survivors. Thus, a systematic method to evaluate if the dose level of non-primary radiation meets the IEC standard requirements is needed.Approach. Different from the traditional photon radiation therapy system, proton therapy systems are composed of several subsystems in a thick bunker. These subsystems are all possible sources of non-primary radiation threatening the patient. As a case study, 7 sources in the P-Cure synchrotron-based proton therapy system are modeled in Monte Carlo (MC) code: tandem injector, injection, synchrotron ring, extraction, beam transport line, scanning nozzle and concrete reflection/scattering. To accurately evaluate the synchrotron beam loss and non-primary dose, a new model called the torus source model is developed. Its parametric equations define the position and direction of the off-orbit particle bombardment on the torus pipe shell in the Cartesian coordinate system. Non-primary doses are finally calculated by several FLUKA simulations.Main results. The ratios of summarized non-primary doses from different sources to the planned dose of 2 Gy are all much smaller than the IEC requirements in both the 15-50 cm and 50-200 cm regions. Thus, the P-Cure synchrotron-based proton therapy system is clean and patient-friendly, and there is no need an inner shielding concrete between the accelerator and patient.Significance. Non-primary radiation dose level is a very important indicator to evaluate the quality of a PT system. This manuscript provides a feasible MC procedure for synchrotron-based proton therapy with new beam loss model. Which could help people figure out precisely whether this level complies with the IEC standard before the system put into clinical treatment. What' more, the torus source model could be widely used for bending magnets in gantries and synchrotrons to evaluate non-primary doses or other radiation doses.

Single pulse protoacoustic range verification using a clinical synchrocyclotron

Objective.Proton therapy as the next generation radiation-based cancer therapy offers dominant advantages over conventional radiation therapy due to the utilization of the Bragg peak; however, range uncertainty in beam delivery substantially mitigates the advantages of proton therapy. This work reports using protoacoustic measurements to determine the location of proton Bragg peak deposition within a water phantom in real time during beam delivery.Approach.In protoacoustics, proton beams have a definitive range, depositing a majority of the dose at the Bragg peak; this dose is then converted to heat. The resulting thermoelastic expansion generates a 3D acoustic wave, which can be detected by acoustic detectors to localize the Bragg peak.Main results.Protoacoustic measurements were performed with a synchrocyclotron proton machine over the exhaustive energy range from 45.5 to 227.15 MeV in clinic. It was found that the amplitude of the acoustic waves is proportional to proton dose deposition, and therefore encodes dosimetric information. With the guidance of protoacoustics, each individual proton beam (7 pC/pulse) can be directly visualized with sub-millimeter (<0.7 mm) resolution using single beam pulse for the first time.Significance.The ability to localize the Bragg peak in real-time and obtain acoustic signals proportional to dose within tumors could enable precision proton therapy and hope to progress towardsin vivomeasurements.

Quantitative assessment of radionuclide production yields in in-beam and offline PET measurements at different proton irradiation facilities

Objective. Reliable radionuclide production yield data are a prerequisite for positron-emission-tomography (PET) based in vivo proton treatment verification. In this context, activation data acquired at two different treatment facilities with different imaging systems were analyzed to provide experimentally determined radionuclide yields in thick targets and were compared with each other to investigate the impact of the respective imaging technique. Approach. Homogeneous thick targets (PMMA, gelatine, and graphite) were irradiated with mono-energetic proton pencil-beams at two distinct energies. Material activation was measured (i) in-beam during and after beam delivery with a double-head prototype PET camera and (ii) offline shortly after beam delivery with a commercial full-ring PET/CT scanner. Integral as well as depth-resolved β +-emitter yields were determined for the dominant positron-emitting radionuclides 11C, 15O, 13N and (in-beam only) 10C. In-beam data were used to investigate the qualitative impact of different monitoring time schemes on activity depth profiles and their quantitative impact on count rates and total activity. Main results. Production yields measured with the in-beam camera were comparable to or higher compared to respective offline results. Depth profiles of radionuclide-specific yields obtained from the double-head camera showed qualitative differences to data acquired with the full-ring camera with a more convex profile shape. Considerable impact of the imaging timing scheme on the activity profile was observed for gelatine only with a range variation of up to 3.5 mm. Evaluation of the coincidence rate and the total number of observed events in the considered workflows confirmed a strongly decreasing rate in targets with a large oxygen fraction. Significance. The observed quantitative and qualitative differences between the datasets underline the importance of a thorough system commissioning. Due to the lack of reliable cross-section data, in-house phantom measurements are still considered a gold standard for careful characterization of the system response and to ensure a reliable beam range verification.

Proton Therapy for Prostate Cancer: Challenges and Opportunities

Simple Summary

Reported clinical outcomes of proton therapy (PT) for localized prostate cancer are similar to photon-based external beam radiotherapy. Apparently, the dosimetric advantages of PT have yet to be translated to clinical benefits. The suboptimal clinical outcomes of PT might be attributable to inadequate dose prescription, as indicated by the ASCENDE-RT trial. Moreover, uncertainties involved in the treatment planning and delivery processes, as well as technological limitations in PT treatment systems, may lead to discrepancies between planned doses and actual doses delivered to patients. In this article, we reviewed the current status of PT for prostate cancer and discussed different clinical implementations that could potentially improve the clinical outcome of PT for prostate cancer. Various technological advancements under which uncertainties in dose calculations can be minimized, including MRI-guided PT, dual-energy photon-counting CT and high-resolution Monte Carlo-based treatment planning systems, are highlighted.

Abstract

The dosimetric advantages of proton therapy (PT) treatment plans are demonstrably superior to photon-based external beam radiotherapy (EBRT) for localized prostate cancer, but the reported clinical outcomes are similar. This may be due to inadequate dose prescription, especially in high-risk disease, as indicated by the ASCENDE-RT trial. Alternatively, the lack of clinical benefits with PT may be attributable to improper dose delivery, mainly due to geometric and dosimetric uncertainties during treatment planning, as well as delivery procedures that compromise the dose conformity of treatments. Advanced high-precision PT technologies, and treatment planning and beam delivery techniques are being developed to address these uncertainties. For instance, external magnetic resonance imaging (MRI)-guided patient setup rooms are being developed to improve the accuracy of patient positioning for treatment. In-room MRI-guided patient positioning systems are also being investigated to improve the geometric accuracy of PT. Soon, high-dose rate beam delivery systems will shorten beam delivery time to within one breath hold, minimizing the effects of organ motion and patient movements. Dual-energy photon-counting computed tomography and high-resolution Monte Carlo-based treatment planning systems are available to minimize uncertainties in dose planning calculations. Advanced in-room treatment verification tools such as prompt gamma detector systems will be used to verify the depth of PT. Clinical implementation of these new technologies is expected to improve the accuracy and dose conformity of PT in the treatment of localized prostate cancers, and lead to better clinical outcomes. Improvement in dose conformity may also facilitate dose escalation, improving local control and implementation of hypofractionation treatment schemes to improve patient throughput and make PT more cost effective.

Radioactivity and Space Range of Ultra-Low-Activity for in vivo Off-line PET Verification of Proton and Carbon Ion Beam-A Phantom Study

Purpose: The radioactivity induced by proton and heavy ion beam belongs to the ultra-low-activity (ULA). Therefore, the radioactivity and space range of commercial off-line positron emission tomography (PET) acquisition based on ULA should be evaluated accurately to guarantee the reliability of clinical verification. The purpose of this study is to quantify the radioactivity and space range of off-line PET acquisition by simulating the ULA triggered by proton and heavy ion beam.

Methods: PET equipment validation phantom and low activity ¹⁸F-FDG were used to simulate the ULA with radioactivity of 11.1–1480 Bq/mL. The radioactivity of ULA was evaluated by comparing the radioactivity in the images with the values calculated from the decay function with a radioactivity error tolerance of 5%. The space range of ULA was evaluated by comparing the width of the R50 analyzed activity distribution curve with the actual width of the container with a space range error tolerance of 4 mm.

Results: When radioactivity of ULA was >148 Bq/mL, the radioactivity error was <5%. When radioactivity of ULA was >30 Bq/mL, the space range error was below 4 mm.

Conclusions: Off-line PET can be used to quantify the radioactivity of proton and heavy ion beam when the ULA exceeds 148 Bq/mL, both in radioactivity and in space range.

A proof-of-concept study of an in-situ partial-ring time-of-flight PET scanner for proton beam verification

Development of a PET system capable of in-situ imaging requires a design that can accommodate the proton treatment beam nozzle. Among the several PET instrumentation approaches developed thus far, the dual-panel PET scanner is often used as it is simpler to develop and integrate within the proton therapy gantry. Partial-angle coverage of these systems can however lead to limited-angle artefacts in the reconstructed PET image. We have previously demonstrated via simulations that time-of-flight (TOF) reconstruction reduces the artifacts accompanying limited-angle data, and permits proton range measurement with 1-2 mm accuracy and precision. In this work we show measured results from a small proof-of-concept dual-panel PET system that uses TOF information to reconstruct PET data acquired after proton irradiation. The PET scanner comprises of two detector modules, each comprised of an array of 4×4×30 mm³ lanthanum bromide scintillator. Measurements are performed with an oxygen-rich gel-water, an adipose tissue equivalent material, and in vitro tissue phantoms. For each phantom measurement, 2 Gy dose was deposited using 54 - 100 MeV proton beams. For each phantom, a Monte Carlo simulation generating the expected distribution of PET isotope from the corresponding proton irradiation was also performed. Proton range was calculated by drawing multiple depth-profiles over a central region encompassing the proton dose deposition. For each profile, proton range was calculated using two techniques (a) 50% pick-off from the distal edge of the profile, and (b) comparing the measured and Monte Carlo profile to minimize the absolute sum of differences over the entire profile. A 10 min PET acquisition acquired with minimal delay post proton-irradiation is compared with a 10 min PET scan acquired after a 20 min delay. Measurements show that PET acquisition with minimal delay is necessary to collect ¹⁵O signal, and maximize ¹¹C signal collection with a short PET acquisition. In comparison with the 50% pick-off technique, the shift technique is more robust and offers better precision in measuring the proton range for the different phantoms. Range measurements from PET images acquired with minimal delay, and the shift technique demonstrate the ability to achieve <1.5 mm accuracy and precision in estimating proton range.

Novel On-line PET Imaging for Intra-Beam Range Verification and Delivery Optimization: A Simulation Feasibility Study

On-line PET image-based method uses an initial particle beam to measure the particle beam range (BR) within the same fraction so that any measured range-shift with respect to the predicted BR can be compensated before the rest therapeutic beam deliveries. However, the method requires to use a low-dose initial beam to minimize the risk of beam overshooting, which leads to low image count and inaccurate BR measurement. In this in-silico study, we evaluated the feasibility of a new on-line PET imaging method that measures BR at the mid-plane of a target volume with part of the high-dose therapy beams to verify BR and guide adaptive treatment re-planning. Simulations included various processes of proton beam radiations to a tumor inside a human brain phantom, positron and PET image generation at the mid-plane with initial beams, activity range measurement, and range-shift compensated beam delivery. The results demonstrated that the new method, under the simulated conditions, can achieve ~1.1 mm mid-plane BR measurement accuracy and closely match the delivered range-shift compensated dose distribution with the planned one. Overall, it is promising that this new method may significantly improve particle therapy accuracy.

In Vivo 3-D Dose Verification Using PET/CT Images After Carbon-Ion Radiation Therapy

Objective

To investigate the usefulness of positron emission tomography (PET) images obtained after carbon-ion irradiation for dose verification in carbon-ion radiotherapy.

Methods and Materials

An anthropomorphic head phantom was used in this study. Three cubes with volumes of 1, 4, and 10 ml were contoured as targets in the phantom CT through a treatment planning system. Treatment plans with six prescriptions from 2.5 to 10 Gy (2.5, 3, 5, 6, 8, and 10 Gy effective dose) were designed and delivered by 90° fixed carbon-ion beams, respectively. After irradiation of the phantom, a PET/CT scan was performed to fuse the treatment-planning CT image with the PET/CT image. The relationship between target volume and the standard uptake value (SUV) in PET/CT was evaluated for corresponding plan prescription. The MIM Maestro software was used for the image fusion and data analysis.

Results

SUV in the target had an approximate linear relationship with the effective dose. The same effective dose could generate a roughly equal SUV for different target volumes (p < 0.05).

Conclusions

It is feasible to verify the actual 3-D dose distribution of carbon-ion radiotherapy by the approach in this study.

Denoising PET images for proton therapy using a residual U-net

The use of proton therapy has the advantage of high dose concentration as it is possible to concentrate the dose on the tumor while suppressing damage to the surrounding normal organs. However, the range uncertainty significantly affects the actual dose distribution in the vicinity of the proton range, limiting the benefit of proton therapy for reducing the dose to normal organs. By measuring the annihilation gamma rays from the produced positron emitters, it is possible to obtain a proton induced positron emission tomography (pPET) image according to the irradiation region of the proton beam. Smoothing with a Gaussian filter is generally used to denoise PET images; however, this approach lowers the spatial resolution. Furthermore, other conventional smoothing processing methods may deteriorate the steep region of the pPET images. In this study, we proposed a denoising method based on a Residual U-Net for pPET images. We conducted the Monte Carlo simulation and irradiation experiment on a human phantom to obtain pPET data. The accuracy of the range estimation and the image similarity were evaluated for pPET images using the Residual U-Net, a Gaussian filter, a median filter, the block-matching and 3D-filtering (BM3D), and a total variation (TV) filter. Usage of the Residual U-Net yielded effective results corresponding to the range estimation; however, the results of peak-signal-to-noise ratio were identical to those for the Gaussian filter, median filter, BM3D, and TV filter. The proposed method can contribute to improving the accuracy of treatment verification and shortening the PET measurement time.

Spatial Agreement of Brainstem Dose Distributions Depending on Biological Model in Proton Therapy for Pediatric Brain Tumors

Purpose

During radiation therapy for pediatric brain tumors, the brainstem is a critical organ at risk, possibly with different radio-sensitivity across its substructures. In proton therapy, treatment planning is currently performed using a constant relative biological effectiveness (RBE) of 1.1 (RBE_1.1), whereas preclinical studies point toward spatial variability of this factor. To shed light on this biological uncertainty, we investigated the spatial agreement between isodose maps produced by different RBE models, with emphasis on (smaller) substructures of the brainstem.

Methods and Materials

Proton plans were recalculated using Monte Carlo simulations in 3 anonymized pediatric patients with brain tumors (a craniopharyngioma, a low-grade glioma, and a posterior fossa ependymoma) to obtain dose and linear energy transfer distributions. Doses and volume metrics for the brainstem and its substructures were calculated using a constant RBE_1.1, 4 phenomenological RBE models with varying (α/β)_x parameters, and with a simpler linear energy transfer-dependent model. The spatial agreement between the dose distributions of constant RBE_1.1 versus the variable RBE models was compared using the Dice similarity coefficient.

Results

The spatial agreement between the variable RBE dose distributions and RBE_1.1 decreased with increasing isodose levels in all patient cases. The patient with ependymoma showed the greatest variation in dose and dose volumes, where V_50Gy(RBE) in the brainstem increased from 32% (RBE_1.1) to 35% to 49% depending on the applied model, corresponding to a spatial agreement (Dice similarity coefficient) between 0.79 and 0.95. The remaining patients showed similar trends, however, with lower absolute values due to lower brainstem doses.

Conclusions

All phenomenological RBE models fully enclosed the isodose volumes of the constant RBE_1.1, and the volumes based on variable RBE spatially agreed. The spatial agreement was dependent on the isodose level, where higher isodose levels showed larger expansions and less agreement between the variable RBE models and RBE_1.1.

Feasibility of quasi-prompt PET-based range verification in proton therapy

Compared to photon therapy, proton therapy allows a better conformation of the dose to the tumor volume with reduced radiation dose to co-irradiated tissues. In vivo verification techniques including positron emission tomography (PET) have been proposed as quality assurance tools to mitigate proton range uncertainties. Detection of differences between planned and actual dose delivery on a short timescale provides a fast trigger for corrective actions. Conventional PET-based imaging of ¹⁵O (T_1/2 = 2 min) and ¹¹C (T_1/2 = 20 min) distributions precludes such immediate feedback. We here present a demonstration of near real-time range verification by means of PET imaging of ¹²N (T_1/2 = 11 ms). PMMA and graphite targets were irradiated with a 150 MeV proton pencil beam consisting of a series of pulses of 10 ms beam-on and 90 ms beam-off. Two modules of a modified Siemens Biograph mCT PET scanner (21 × 21 cm² each), installed 25 cm apart, were used to image the beam-induced PET activity during the beam-off periods. The modifications enable the detectors to be switched off during the beam-on periods. ¹²N images were reconstructed using planar tomography. Using a 1D projection of the 2D reconstructed ¹²N image, the activity range was obtained from a fit of the activity profile with a sigmoid function. Range shifts due to modified target configurations were assessed for multiples of the clinically relevant 10⁸ protons per pulse (approximately equal to the highest intensity spots in the pencil beam scanning delivery of a dose of 1 Gy over a cubic 1 l volume). The standard deviation of the activity range, determined from 30 datasets obtained from three irradiations on PMMA and graphite targets, was found to be 2.5 and 2.6 mm (1σ) with 10⁸ protons per pulse and 0.9 and 0.8 mm (1σ) with 10⁹ protons per pulse. Analytical extrapolation of the results from this study shows that using a scanner with a solid angle coverage of 57%, with optimized detector switching and spot delivery times much smaller than the ¹²N half-life, an activity range measurement precision of 2.0 mm (1σ) and 1.3 mm (1σ) within 50 ms into an irradiation with 4 × 10⁷ and 10⁸ protons per pencil beam spot can be potentially realized. Aggregated imaging of neighboring spots or, if possible, increasing the number of protons for a few probe beam spots will enable the realization of higher precision range measurement.

Implementation of a double scattering nozzle for Monte Carlo recalculation of proton plans with variable relative biological effectiveness

A constant relative biological effectiveness (RBE) of 1.1 is currently used in clinical proton therapy. However, theRBEvaries with factors such as dose level, linear energy transfer (LET) and tissue type. MultipleRBEmodels have been developed to account for this biological variation. To enable recalculation of patients treated with double scattering (DS) proton therapy, includingLETand variableRBE, we implemented and commissioned a Monte Carlo (MC) model of a DS treatment nozzle. The main components from the IBA nozzle were implemented in the FLUKA MC code. We calibrated and verified the following entities to experimental measurements: range of pristine Bragg peaks (PBPs) and spread-out Bragg peaks (SOBPs), energy spread, lateral profiles, compensator range degradation, and absolute dose. We recalculated two patients with different field setups, comparing FLUKA vs. treatment planning system (TPS) dose, also obtainingLETand variableRBEdoses. We achieved good agreement between FLUKA and measurements. The range differences between FLUKA and measurements were for the PBPs within ±0.9 mm (83% ⩽ 0.5 mm), and for SOBPs ±1.6 mm (82% ⩽ 0.5 mm). The differences in modulation widths were below 5 mm (79% ⩽ 2 mm). The differences in the distal dose fall off (D80%-D20%) were below 0.5 mm for all PBPs and the lateral penumbras diverged from measurements by less than 1 mm. The mean dose difference (RBE= 1.1) in the target between the TPS and FLUKA were below 0.4% in a three-field plan and below 1.4% in a four-field plan. A dose increase of 9.9% and 7.2% occurred when using variableRBEfor the two patients, respectively. We presented a method to recalculate DS proton plans in the FLUKA MC code. The implementation was used to obtainLETand variableRBEdose and can be used for investigating variableRBEfor previously treated patients.

Particle beam microstructure reconstruction and coincidence discrimination in PET monitoring for hadron therapy

Positron emission tomography is one of the most mature techniques for monitoring the particles range in hadron therapy, aiming to reduce treatment uncertainties and therefore the extent of safety margins in the treatment plan. In-beam PET monitoring has been already performed using inter-spill and post-irradiation data, i.e. while the particle beam is off or paused. The full beam acquisition procedure is commonly discarded because the particle spills abruptly increase the random coincidence rates and therefore the image noise. This is because random coincidences cannot be separated by annihilation photons originating from radioactive decays and cannot be corrected with standard random coincidence techniques due to the time correlation of the beam-induced background with the ion beam microstructure. The aim of this paper is to provide a new method to recover in-spill data to improve the images obtained with full-beam PET acquisitions. This is done by estimating the temporal microstructure of the beam and thus selecting input PET events that are less likely to be random ones. The PET detector we used was the one developed within the INSIDE project and tested at the CNAO synchrotron-based facility. The data were taken on a PMMA phantom irradiated with 72 MeV proton pencil beams. The obtained results confirm the possibility of improving the acquired PET data without any external signal coming from the synchrotron or ad hoc detectors.

Evaluation of Proton Therapy Accuracy Using a PMMA Phantom and PET Prediction Module

Purpose: Positron emission tomography (PET) scanning is a widely used method of proton therapy verification. In this study, a proton radiotherapy accuracy verification process was developed by comparing predicted and measured PET data to verify the correctness of PET prediction and was tested at the Shanghai Proton and Heavy Ion Center.

Method: Irradiation was performed on a polymethyl methacrylate (PMMA) phantom. There were two dose groups, to which 2 and 4 Gy doses were delivered, and each dose group had different designed dose depths ranging from 5 to 20 cm. The predicted PET results were obtained using a PET prediction calculation module. The measured data were collected with a PET/computed tomography device. The predicted and measured PET data were normalized to similar PET amplitude values before comparison and were compared using depth and lateral profiles for the position error. The error was evaluated at the position corresponding to 50% of the maximum on the PET curves. The mean and standard deviation were calculated based on the data sampled in the scoring area. Gamma index analysis is also applied in the comparison.

Results: In the depth comparison, the 2 and 4 Gy dose cases yielded similar mean depth errors between 1 and −1 mm, and the deviation was <2 mm. In the lateral comparison, the 2 Gy cases had a mean lateral error around 1 mm, and the 4 Gy cases had a mean lateral error <1 mm, with a standard deviation <1 mm for both the 2 and 4 Gy cases. All the cases have a gamma passing rate over 95%.

Conclusion: The comparison of these PMMA phantom cases revealed good agreement between the predicted and measured PET data, with depth and lateral position errors <2 mm in total, considering the uncertainty. The comparison results demonstrate that the PET predictions obtained in PMMA phantom tests for single proton beam therapy verification are reliable and that the research can be extended to verification in human body treatment with further investigation.

Gel dosimetry for three dimensional proton range measurements in anthropomorphic geometries

Proton beams used for radiotherapy have potential for superior sparing of normal tissue, although range uncertainties are among the main limiting factors in the accuracy of dose delivery. The aim of this study was to benchmark an N-vinylpyrrolidone based polymer gel to perform three-dimensional measurement of geometric proton beam characteristics and especially to test its suitability as a range probe in combination with an anthropomorphic phantom. For single proton pencil beams as well as for 3×3cm² mono-energy layers depth dose profiles, lateral dose distribution at different depths and proton range were evaluated in simple cubic gel phantoms at different energies from 75 to 115MeV and different dose levels. In addition, a 90MeV mono-energetic beam was delivered to an anthropomorphic 3D printed head phantom, which was filled with gel. Subsequently, all phantoms underwent magnetic resonance imaging using an axial pixel size of 0.68-0.98mm and with slice thicknesses of 2 or 3mm to derive a 3-dimensional distribution of the T₂ relaxation time, which correlates with radiation dose. Indices describing lateral dose distribution and proton range were compared against predictions from a treatment planning system (TPS, for cubic and head phantoms) and Monte Carlo simulations (MC, for the head phantom) after manual rigid co-registration with the T₂ relaxation time datasets. For all pencil beams, the FWHM agreement with TPS was better than 1mm or 7%. For the mono-energetic layer, the agreement with TPS in this respect was even better than 0.3mm in each case. With respect to range, results from gel measurements differed no more than 0.9mm (1.6%) from values predicted by TPS. In case of the anthropomorphic phantom, deviations with respect to a nominal range of about 61mm as well as in FWHM were slightly higher, namely within 1.0mm and 1.1mm respectively. Average deviations between gel and TPS/MC were similar (-0.3mm±0.4mm/-0.2±0.5mm). In conclusion, polymer gel dosimetry was found to be a valuable tool to determine geometric proton beam properties three-dimensionally and with high spatial resolution in simple cubic as well as in a more complex anthropomorphic phantom. Post registration range errors of the order of 1mm could be achieved. The additional registration uncertainty (95%) was 1mm.

Carbon ions beam therapy monitoring with the INSIDE in-beam PET

In vivo range monitoring techniques are necessary in order to fully take advantage of the high dose gradients deliverable in hadrontherapy treatments. Positron emission tomography (PET) scanners can be used to monitor beam-induced activation in tissues and hence measure the range. The INSIDE (Innovative Solutions for In-beam DosimEtry in Hadrontherapy) in-beam PET scanner, installed at the Italian National Center of Oncological Hadrontherapy (CNAO, Pavia, Italy) synchrotron facility, has already been successfully tested in vivo during a proton therapy treatment. We discuss here the system performance evaluation with carbon ion beams, in view of future in vivo tests. The work is focused on the analysis of activity images obtained with therapeutic treatments delivered to polymethyl methacrylate (PMMA) phantoms, as well as on the test of an innovative and robust Monte Carlo simulation technique for the production of reliable prior activity maps. Images are reconstructed using different integration intervals, so as to monitor the activity evolution during and after the treatment. Three procedures to compare activity images are presented, namely Pearson correlation coefficient, Beam's eye view and overall view. Images of repeated irradiations of the same treatments are compared to assess the integration time necessary to provide reproducible images. The range agreement between simulated and experimental images is also evaluated, so as to validate the simulation capability to provide sound prior information. The results indicate that at treatment end, or at most 20 s afterwards, the range measurement is reliable within 1-2 mm, when comparing both different experimental sessions and data with simulations. In conclusion, this work shows that the INSIDE in-beam PET scanner performance is promising towards its in vivo test with carbon ions.

Design of a new tracking device for on-line beam range monitor in carbon therapy

Charged particle therapy is a technique for cancer treatment that exploits hadron beams, mostly protons and carbon ions. A critical issue is the monitoring of the beam range so to check the correct dose deposition to the tumor and surrounding tissues. The design of a new tracking device for beam range real-time monitoring in pencil beam carbon ion therapy is presented. The proposed device tracks secondary charged particles produced by beam interactions in the patient tissue and exploits the correlation of the charged particle emission profile with the spatial dose deposition and the Bragg peak position. The detector, currently under construction, uses the information provided by 12 layers of scintillating fibers followed by a plastic scintillator and a pixelated Lutetium Fine Silicate (LFS) crystal calorimeter. An algorithm to account and correct for emission profile distortion due to charged secondaries absorption inside the patient tissue is also proposed. Finally detector reconstruction efficiency for charged particle emission profile is evaluated using a Monte Carlo simulation considering a quasi-realistic case of a non-homogenous phantom.

Nuclear physics in particle therapy: a review

Charged particle therapy has been largely driven and influenced by nuclear physics. The increase in energy deposition density along the ion path in the body allows reducing the dose to normal tissues during radiotherapy compared to photons. Clinical results of particle therapy support the physical rationale for this treatment, but the method remains controversial because of the high cost and of the lack of comparative clinical trials proving the benefit compared to x-rays. Research in applied nuclear physics, including nuclear interactions, dosimetry, image guidance, range verification, novel accelerators and beam delivery technologies, can significantly improve the clinical outcome in particle therapy. Measurements of fragmentation cross-sections, including those for the production of positron-emitting fragments, and attenuation curves are needed for tuning Monte Carlo codes, whose use in clinical environments is rapidly increasing thanks to fast calculation methods. Existing cross sections and codes are indeed not very accurate in the energy and target regions of interest for particle therapy. These measurements are especially urgent for new ions to be used in therapy, such as helium. Furthermore, nuclear physics hardware developments are frequently finding applications in ion therapy due to similar requirements concerning sensors and real-time data processing. In this review we will briefly describe the physics bases, and concentrate on the open issues.

The FLUKA Code: An Accurate Simulation Tool for Particle Therapy

Monte Carlo (MC) codes are increasingly spreading in the hadrontherapy community due to their detailed description of radiation transport and interaction with matter. The suitability of a MC code for application to hadrontherapy demands accurate and reliable physical models capable of handling all components of the expected radiation field. This becomes extremely important for correctly performing not only physical but also biologically based dose calculations, especially in cases where ions heavier than protons are involved. In addition, accurate prediction of emerging secondary radiation is of utmost importance in innovative areas of research aiming at in vivo treatment verification. This contribution will address the recent developments of the FLUKA MC code and its practical applications in this field. Refinements of the FLUKA nuclear models in the therapeutic energy interval lead to an improved description of the mixed radiation field as shown in the presented benchmarks against experimental data with both (4)He and (12)C ion beams. Accurate description of ionization energy losses and of particle scattering and interactions lead to the excellent agreement of calculated depth-dose profiles with those measured at leading European hadron therapy centers, both with proton and ion beams. In order to support the application of FLUKA in hospital-based environments, Flair, the FLUKA graphical interface, has been enhanced with the capability of translating CT DICOM images into voxel-based computational phantoms in a fast and well-structured way. The interface is capable of importing also radiotherapy treatment data described in DICOM RT standard. In addition, the interface is equipped with an intuitive PET scanner geometry generator and automatic recording of coincidence events. Clinically, similar cases will be presented both in terms of absorbed dose and biological dose calculations describing the various available features.

Investigating the limits of PET/CT imaging at very low true count rates and high random fractions in ion-beam therapy monitoring

PURPOSE

External beam radiotherapy with protons and heavier ions enables a tighter conformation of the applied dose to arbitrarily shaped tumor volumes with respect to photons, but is more sensitive to uncertainties in the radiotherapeutic treatment chain. Consequently, an independent verification of the applied treatment is highly desirable. For this purpose, the irradiation-induced β(+)-emitter distribution within the patient is detected shortly after irradiation by a commercial full-ring positron emission tomography/x-ray computed tomography (PET/CT) scanner installed next to the treatment rooms at the Heidelberg Ion-Beam Therapy Center (HIT). A major challenge to this approach is posed by the small number of detected coincidences. This contribution aims at characterizing the performance of the used PET/CT device and identifying the best-performing reconstruction algorithm under the particular statistical conditions of PET-based treatment monitoring. Moreover, this study addresses the impact of radiation background from the intrinsically radioactive lutetium-oxyorthosilicate (LSO)-based detectors at low counts.

METHODS

The authors have acquired 30 subsequent PET scans of a cylindrical phantom emulating a patientlike activity pattern and spanning the entire patient counting regime in terms of true coincidences and random fractions (RFs). Accuracy and precision of activity quantification, image noise, and geometrical fidelity of the scanner have been investigated for various reconstruction algorithms and settings in order to identify a practical, well-suited reconstruction scheme for PET-based treatment verification. Truncated listmode data have been utilized for separating the effects of small true count numbers and high RFs on the reconstructed images. A corresponding simulation study enabled extending the results to an even wider range of counting statistics and to additionally investigate the impact of scatter coincidences. Eventually, the recommended reconstruction scheme has been applied to exemplary postirradiation patient data-sets.

RESULTS

Among the investigated reconstruction options, the overall best results in terms of image noise, activity quantification, and accurate geometrical recovery were achieved using the ordered subset expectation maximization reconstruction algorithm with time-of-flight (TOF) and point-spread function (PSF) information. For this algorithm, reasonably accurate (better than 5%) and precise (uncertainty of the mean activity below 10%) imaging can be provided down to 80,000 true coincidences at 96% RF. Image noise and geometrical fidelity are generally improved for fewer iterations. The main limitation for PET-based treatment monitoring has been identified in the small number of true coincidences, rather than the high intrinsic random background. Application of the optimized reconstruction scheme to patient data-sets results in a 25% - 50% reduced image noise at a comparable activity quantification accuracy and an improved geometrical performance with respect to the formerly used reconstruction scheme at HIT, adopted from nuclear medicine applications.

CONCLUSIONS

Under the poor statistical conditions in PET-based treatment monitoring, improved results can be achieved by considering PSF and TOF information during image reconstruction and by applying less iterations than in conventional nuclear medicine imaging. Geometrical fidelity and image noise are mainly limited by the low number of true coincidences, not the high LSO-related random background. The retrieved results might also impact other emerging PET applications at low counting statistics.

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.

Mapping (15)O production rate for proton therapy verification

PURPOSE

This work was a proof-of-principle study for the evaluation of oxygen-15 ((15)O) production as an imaging target through the use of positron emission tomography (PET), to improve verification of proton treatment plans and to study the effects of perfusion.

METHODS AND MATERIALS

Dynamic PET measurements of irradiation-produced isotopes were made for a phantom and rabbit thigh muscles. The rabbit muscle was irradiated and imaged under both live and dead conditions. A differential equation was fitted to phantom and in vivo data, yielding estimates of (15)O production and clearance rates, which were compared to live versus dead rates for the rabbit and to Monte Carlo predictions.

RESULTS

PET clearance rates agreed with decay constants of the dominant radionuclide species in 3 different phantom materials. In 2 oxygen-rich materials, the ratio of (15)O production rates agreed with the expected ratio. In the dead rabbit thighs, the dynamic PET concentration histories were accurately described using (15)O decay constant, whereas the live thigh activity decayed faster. Most importantly, the (15)O production rates agreed within 2% (P>.5) between conditions.

CONCLUSIONS

We developed a new method for quantitative measurement of (15)O production and clearance rates in the period immediately following proton therapy. Measurements in the phantom and rabbits were well described in terms of (15)O production and clearance rates, plus a correction for other isotopes. These proof-of-principle results support the feasibility of detailed verification of proton therapy treatment delivery. In addition, (15)O clearance rates may be useful in monitoring permeability changes due to therapy.

Dosimetry of very high energy electrons (VHEE) for radiotherapy applications: using radiochromic film measurements and Monte Carlo simulations

Very high energy electrons (VHEE) in the range from 100-250 MeV have the potential of becoming an alternative modality in radiotherapy because of their improved dosimetry properties compared with MV photons from contemporary medical linear accelerators. Due to the need for accurate dosimetry of small field size VHEE beams we have performed dose measurements using EBT2 Gafchromic® film. Calibration of the film has been carried out for beams of two different energy ranges: 20 MeV and 165 MeV from conventional radio frequency linear accelerators. In addition, EBT2 film has been used for dose measurements with 135 MeV electron beams produced by a laser-plasma wakefield accelerator. The dose response measurements and percentage depth dose profiles have been compared with calculations carried out using the general-purpose FLUKA Monte Carlo (MC) radiation transport code. The impact of induced radioactivity on film response for VHEEs has been evaluated using the MC simulations. A neutron yield of the order of 10(-5) neutrons cm(-2) per incident electron has been estimated and induced activity due to radionuclide production is found to have a negligible effect on total dose deposition and film response. Neutron and proton contribution to the equivalent doses are negligible for VHEE. The study demonstrates that EBT2 Gafchromic film is a reliable dosimeter that can be used for dosimetry of VHEE. The results indicate an energy-independent response of the dosimeter for 20 MeV and 165 MeV electron beams and has been found to be suitable for dosimetry of VHEE.

Automation and uncertainty analysis of a method for in-vivo range verification in particle therapy

We introduce the automation of the range difference calculation deduced from particle-irradiation induced β(+)-activity distributions with the so-called most-likely-shift approach, and evaluate its reliability via the monitoring of algorithm- and patient-specific uncertainty factors. The calculation of the range deviation is based on the minimization of the absolute profile differences in the distal part of two activity depth profiles shifted against each other. Depending on the workflow of positron emission tomography (PET)-based range verification, the two profiles under evaluation can correspond to measured and simulated distributions, or only measured data from different treatment sessions. In comparison to previous work, the proposed approach includes an automated identification of the distal region of interest for each pair of PET depth profiles and under consideration of the planned dose distribution, resulting in the optimal shift distance. Moreover, it introduces an estimate of uncertainty associated to the identified shift, which is then used as weighting factor to 'red flag' problematic large range differences. Furthermore, additional patient-specific uncertainty factors are calculated using available computed tomography (CT) data to support the range analysis. The performance of the new method for in-vivo treatment verification in the clinical routine is investigated with in-room PET images for proton therapy as well as with offline PET images for proton and carbon ion therapy. The comparison between measured PET activity distributions and predictions obtained by Monte Carlo simulations or measurements from previous treatment fractions is performed. For this purpose, a total of 15 patient datasets were analyzed, which were acquired at Massachusetts General Hospital and Heidelberg Ion-Beam Therapy Center with in-room PET and offline PET/CT scanners, respectively. Calculated range differences between the compared activity distributions are reported in a 2D map in beam-eye-view. In comparison to previously proposed approaches, the new most-likely-shift method shows more robust results for assessing in-vivo the range from strongly varying PET distributions caused by differing patient geometry, ion beam species, beam delivery techniques, PET imaging concepts and counting statistics. The additional visualization of the uncertainties and the dedicated weighting strategy contribute to the understanding of the reliability of observed range differences and the complexity in the prediction of activity distributions. The proposed method promises to offer a feasible technique for clinical routine of PET-based range verification.

A fast Monte Carlo code for proton transport in radiation therapy based on MCNPX

An important requirement for proton therapy is a software for dose calculation. Monte Carlo is the most accurate method for dose calculation, but it is very slow. In this work, a method is developed to improve the speed of dose calculation. The method is based on pre-generated tracks for particle transport. The MCNPX code has been used for generation of tracks. A set of data including the track of the particle was produced in each particular material (water, air, lung tissue, bone, and soft tissue). This code can transport protons in wide range of energies (up to 200 MeV for proton). The validity of the fast Monte Carlo (MC) code is evaluated with data MCNPX as a reference code. While analytical pencil beam algorithm transport shows great errors (up to 10%) near small high density heterogeneities, there was less than 2% deviation of MCNPX results in our dose calculation and isodose distribution. In terms of speed, the code runs 200 times faster than MCNPX. In the Fast MC code which is developed in this work, it takes the system less than 2 minutes to calculate dose for 10(6) particles in an Intel Core 2 Duo 2.66 GHZ desktop computer.

Application of activity pencil beam algorithm using measured distribution data of positron emitter nuclei for therapeutic SOBP proton beam

PURPOSE

Recently, much research on imaging the clinical proton-irradiated volume using positron emitter nuclei based on target nuclear fragment reaction has been carried out. The purpose of this study is to develop an activity pencil beam (APB) algorithm for a simulation system for proton-activated positron-emitting imaging in clinical proton therapy using spread-out Bragg peak (SOBP) beams.

METHODS

The target nuclei of activity distribution calculations are (12)C nuclei, (16)O nuclei, and (40)Ca nuclei, which are the main elements in a human body. Depth activity distributions with SOBP beam irradiations were obtained from the material information of ridge filter (RF) and depth activity distributions of compounds of the three target nuclei measured by BOLPs-RGp (beam ON-LINE PET system mounted on a rotating gantry port) with mono-energetic Bragg peak (MONO) beam irradiations. The calculated data of depth activity distributions with SOBP beam irradiations were sorted in terms of kind of nucleus, energy of proton beam, SOBP width, and thickness of fine degrader (FD), which were verified. The calculated depth activity distributions with SOBP beam irradiations were compared with the measured ones. APB kernels were made from the calculated depth activity distributions with SOBP beam irradiations to construct a simulation system using the APB algorithm for SOBP beams.

RESULTS

The depth activity distributions were prepared using the material information of RF and the measured depth activity distributions with MONO beam irradiations for clinical therapy using SOBP beams. With the SOBP width widening, the distal fall-offs of depth activity distributions and the difference from the depth dose distributions were large. The shapes of the calculated depth activity distributions nearly agreed with those of the measured ones upon comparison between the two. The APB kernels of SOBP beams were prepared by making use of the data on depth activity distributions with SOBP beam irradiations that were made from the depth activity distributions with MONO beam irradiations and sorted in terms of energy, SOBP width, and thickness of FD. The data on APB kernels of SOBP beams were determined as installment data for the simulation system using the APB algorithm for SOBP beam irradiations.

CONCLUSIONS

A method of obtaining the depth activity distributions and the APB algorithm for clinical use of SOBP beams have been developed. It is suggested that the simulation system for imaging the clinical irradiated volume with the APB algorithm can be used in clinical proton therapy using SOBP beams by preparing and investigating the data on APB kernels of SOBP beams.

An algorithm to assess the need for clinical Monte Carlo dose calculation for small proton therapy fields based on quantification of tissue heterogeneity

PURPOSE

In proton therapy, complex density heterogeneities within the beam path constitute a challenge to dose calculation algorithms. This might question the reliability of dose distributions predicted by treatment planning systems based on analytical dose calculation. For cases in which substantial dose errors are expected, resorting to Monte Carlo dose calculations might be essential to ensure a successful treatment outcome and therefore the benefit is worth a presumably long computation time. The aim of this study was to define an indicator for the accuracy of dose delivery based on analytical dose calculations in treatment planning systems for small proton therapy fields to identify those patients for which Monte Carlo dose calculation is warranted.

METHODS

Fourteen patients treated at our facility with small passively scattered proton beams (apertures diameters below 7 cm) were selected. Plans were generated in the XiO treatment planning system in combination with a pencil beam algorithm developed at the Massachusetts General Hospital and compared to Monte Carlo dose calculations. Differences in the dose to the 50% of the gross tumor volume (D50, GTV) were assessed in a field-by-field basis. A simple and fast methodology was developed to quantify the inhomogeneity of the tissue traversed by a single small proton beam using a heterogeneity index (HI)-a concept presented by Plugfelder et al. [Med. Phys. 34, 1506-1513 (2007)] for scanned proton beams. Finally, the potential correlation between the error made by the pencil beam based treatment planning algorithm for each field and the level of tissue heterogeneity traversed by the proton beam given by the HI was evaluated.

RESULTS

Discrepancies up to 5.4% were found in D50 for single fields, although dose differences were within clinical tolerance levels (<3%) when combining all of the fields involved in the treatment. The discrepancies found for each field exhibited a strong correlation to their associated HI-values (Spearman's ρ=0.8, p<0.0001); the higher the level of tissue inhomogeneities for a particular field, the larger the error by the analytical algorithm. With the established correlation a threshold for HI can be set by choosing a tolerance level of 2-3%-commonly accepted in radiotherapy.

CONCLUSIONS

The HI is a good indicator for the accuracy of proton field delivery in terms of GTV prescription dose coverage when small fields are delivered. Each HI-value was obtained from the CT image in less than 3 min on a computer with 2 GHz CPU allowing implementation of this methodology in clinical routine. For HI-values exceeding the threshold, either a change in beam direction (if feasible) or a recalculation of the dose with Monte Carlo would be highly recommended.

Dosimetric accuracy of proton therapy for chordoma patients with titanium implants

PURPOSE

To investigate dosimetric errors in proton therapy treatment planning due to titanium implants, and to determine how these affect postoperative passively scattered proton therapy for chordoma patients with orthopedic hardware.

METHODS

The presence of titanium hardware near the tumor may affect the dosimetric accuracy of proton therapy. Artifacts in the computed tomography (CT) scan can cause errors in the proton stopping powers used for dose calculation, which are derived from CT numbers. Also, clinical dose calculation algorithms may not accurately simulate proton beam transport through the implants, which have very different properties as compared to human tissue. The authors first evaluated the impact of these two main issues. Dose errors introduced by metal artifacts were studied using phantoms with and without titanium inserts, and patient scans on which a metal artifact reduction method was applied. Pencil-beam dose calculations were compared to models of nuclear interactions in titanium and Monte Carlo simulations. Then, to assess the overall impact on treatment plans for chordoma, the authors compared the original clinical treatment plans to recalculated dose distributions employing both metal artifact reduction and Monte Carlo methods.

RESULTS

Dose recalculations of clinical proton fields showed that metal artifacts cause range errors up to 6 mm distal to regions affected by CT artifacts. Monte Carlo simulations revealed dose differences >10% in the high-dose area, and range differences up to 10 mm. Since these errors are mostly local in nature, the large number of fields limits the impact on target coverage in the chordoma treatment plans to a small decrease of dose homogeneity.

CONCLUSIONS

In the presence of titanium implants, CT metal artifacts and the approximations of pencil-beam dose calculations cause considerable errors in proton dose calculation. The spatial distribution of the errors however limits the overall impact on passively scattered proton therapy for chordoma.

Feasibility of proton-activated implantable markers for proton range verification using PET

Proton beam range verification using positron emission tomography (PET) currently relies on proton activation of tissue, the products of which decay with a short half-life and necessitate an on-site PET scanner. Tissue activation is, however, negligible near the distal dose fall-off region of the proton beam range due to their high interaction energy thresholds. Therefore Monte Carlo simulation is often supplemented for comparison with measurement; however, this also may be associated with systematic and statistical uncertainties. Therefore, we sought to test the feasibility of using long-lived proton-activated external materials that are inserted or infused into the target volume for more accurate proton beam range verification that could be performed at an off-site PET scanner. We irradiated samples of ≥98% (18)O-enriched water, natural Cu foils, and >97% (68)Zn-enriched foils as candidate materials, along with samples of tissue-equivalent materials including (16)O water, heptane (C7H16), and polycarbonate (C16H14O3)n, at four depths (ranging from 100% to 3% of center of modulation (COM) dose) along the distal fall-off of a modulated 160 MeV proton beam. Samples were irradiated either directly or after being embedded in Plastic Water® or balsa wood. We then measured the activity of the samples using PET imaging for 20 or 30 min after various delay times. Measured activities of candidate materials were up to 100 times greater than those of the tissue-equivalent materials at the four distal dose fall-off depths. The differences between candidate materials and tissue-equivalent materials became more apparent after longer delays between irradiation and PET imaging, due to the longer half-lives of the candidate materials. Furthermore, the activation of the candidate materials closely mimicked the distal dose fall-off with offsets of 1 to 2 mm. Also, signals from the foils were clearly visible compared to the background from the activated Plastic Water® and balsa wood phantoms. These results indicate that markers made from these candidate materials could be used for in vivo proton range verification using an off-site PET scanner.

Proton therapy verification with PET imaging

Proton therapy is very sensitive to uncertainties introduced during treatment planning and dose delivery. PET imaging of proton induced positron emitter distributions is the only practical approach for in vivo, in situ verification of proton therapy. This article reviews the current status of proton therapy verification with PET imaging. The different data detecting systems (in-beam, in-room and off-line PET), calculation methods for the prediction of proton induced PET activity distributions, and approaches for data evaluation are discussed.

Validation of density-elasticity relationships for finite element modeling of human pelvic bone by modal analysis

In total hip arthroplasty and particularly in revision surgery, computer assisted pre-operative prediction of the best possible anchorage strategy for implant fixation would be a great help to the surgeon. Computer simulation relies on validated numerical models. In the current study, three density-elasticity relationships (No. 1-3) from the literature for inhomogeneous material parameter assignment from CT data in automated finite element (FE) modeling of long bones were evaluated for their suitability for FE modeling of human pelvic bone. Numerical modal analysis was conducted on 10 FE models of hemipelvic bone specimens and compared to the gold standard provided by experimental modal analysis results from a previous in-vitro study on the same specimens. Overall, calculated resonance frequencies came out lower than measured values. Magnitude of mean relative deviation of numerical resonance frequencies with regard to measured values is lowest for the density-elasticity relationship No. 3 (-15.9%) and considerably higher for both density-elasticity relationships No. 1 (-41.1%) and No. 2 (-45.0%). Mean MAC values over all specimens amount to 77.8% (No. 1), 78.5% (No. 2), and 83.0% (No. 3). MAC results show, that mode shapes are only slightly influenced by material distribution. Calculated resonance frequencies are generally lower than measured values, which indicates, that numerical models lack stiffness. Even when using the best suited (No. 3) out of three investigated density-elasticity relationships, in FE modeling of pelvic bone a considerable underestimation of model stiffness has to be taken into account.