TOPAS simulation of photoneutrons in radiotherapy: accuracy and speed with variance reduction

We provide optimal particle split numbers for speeding up TOPAS Monte Carlo simulations of linear accelerator (linac) treatment heads while maintaining accuracy. In addition, we provide a new TOPAS physics module for simulating photoneutron production and transport. TOPAS simulation of a Siemens Oncor linac was used to determine the optimal number of splits for directional bremsstrahlung splitting as a function of the field size for 6MV and 18MV x-ray beams. The linac simulation was validated against published data of lateral dose profiles and percentage depth-dose curves (PDD) for the largest square field (40cm side). In separate simulations, neutron particle split and the custom TOPAS physics module was used to generate and transport photoneutrons, called "TsPhotoNeutron". Verification of accuracy was performed by comparing simulations with published measurements of: 1) neutron yields as a function of beam energy for thick targets of Al, Cu, Ta, W, Pb and concrete; and 2) photoneutron energy spectrum at 40cm laterally from the isocenter of the linac from an 18MV beam with closed jaws and MLC. The optimal number of splits obtained for directional bremsstrahlung splitting enhanced the computational efficiency by two orders of magnitude. The efficiency decreased with increasing beam energy and field size. Calculated lateral profiles in the central region agreed within 1mm/2% from measured data, PDD curves within 1 mm/1%. For the TOPAS physics module, at a split number of 146, the efficiency of computing photoneutron yields was enhanced by a 27.6 factor, whereas it improved the accuracy over existing Geant4 physics modules. This work provides simulation parameters and a new TOPAS physics module to improve the efficiency and accuracy of TOPAS simulations that involve photonuclear processes occurring in high-Z materials found in linac components, patient devices, and treatment rooms, as well as to explore new therapeutic modalities such as very-high-energy electron therapy.

A simulation study on the radiosensitization properties of gold nanorods

Objective. Gold nanorods (GNRs) have emerged as versatile nanoparticles with unique properties, holding promise in various modalities of cancer treatment through drug delivery and photothermal therapy. In the rapidly evolving field of nanoparticle radiosensitization (NPRS) for cancer therapy, this study assessed the potential of gold nanorods as radiosensitizing agents by quantifying the key features of NPRS, such as secondary electron emission and dose enhancement, using Monte Carlo simulations.Approach. Employing the TOPAS track structure code, we conducted a comprehensive evaluation of the radiosensitization behavior of spherical gold nanoparticles and gold nanorods. We systematically explored the impact of nanorod geometry (in particular size and aspect ratio) and orientation on secondary electron emission and deposited energy ratio, providing validated results against previously published simulations.Main results. Our findings demonstrate that gold nanorods exhibit comparable secondary electron emission to their spherical counterparts. Notably, nanorods with smaller surface-area-to-volume ratios (SA:V) and alignment with the incident photon beam proved to be more efficient radiosensitizing agents, showing superiority in emitted electron fluence. However, in the microscale, the deposited energy ratio (DER) was not markedly influenced by the SA:V of the nanorod. Additionally, our findings revealed that the geometry of gold nanoparticles has a more significant impact on the emission of M-shell Auger electrons (with energies below 3.5 keV) than on higher-energy electrons.Significance. This research investigated the radiosensitization properties of gold nanorods, positioning them as promising alternatives to the more conventionally studied spherical gold nanoparticles in the context of cancer research. With increasing interest in multimodal cancer therapy, our findings have the potential to contribute valuable insights into the perspective of gold nanorods as effective multipurpose agents for synergistic photothermal therapy and radiotherapy. Future directions may involve exploring alternative metallic nanorods as well as further optimizing the geometry and coating materials, opening new possibilities for more effective cancer treatments.

The influence of different versions of FLUKA and GEANT4 on the calculation of response functions of ionization chambers in clinical proton beams

Objective.To investigate the influence of different versions of the Monte Carlo codesgeant4 andflukaon the calculation of overall response functionsfQof air-filled ionization chambers in clinical proton beams.Approach. fQfactors were calculated for six plane-parallel and four cylindrical ionization chambers withgeant4 andfluka. These factors were compared to already published values that were derived using older versions of these codes.Main results.Differences infQfactors calculated with different versions of the same Monte Carlo code can be up to ∼1%. Especially forgeant4, the updated version leads to a more pronounced dependence offQon proton energy and to smallerfQfactors for high energies.Significance.Different versions of the same Monte Carlo code can lead to differences in the calculation offQfactors of up to ∼1% without changing the simulation setup, transport parameters, ionization chamber geometry modeling, or employed physics lists. These findings support the statement that the dominant contributor to the overall uncertainty of Monte Carlo calculatedfQfactors are type-B uncertainties.

Microdosimetry-based investigation of biological effectiveness of252Cf brachytherapy source: TOPAS Monte Carlo study

Objective.To investigate biological effectiveness of252Cf brachytherapy source using Monte Carlo-calculated microdosimetric distributions.Approach.252Cf source capsule was placed at the center of the spherical water phantom and phase-space data were scored as a function of radial distance in water (R= 1-5 cm) using TOPAS Monte Carlo code. The phase-space data were used to calculate microdosimetric distributions at 1μm site size. Using these distributions, Relative Biological Effectiveness (RBE), mean quality factor (Q̅) and Oxygen Enhancement Ratio (OER) were calculated as a function ofR.Main results.The overall shapes of the microdosimetric distributions are comparable at all the radial distances in water. However, slight variation in the bin-wise yield is observed withR. RBE,Q̅and OER are insensitive to R over the range 1-5 cm. Microdosimetric kinetic model based RBE values are about 2.3 and 2.8 for HSG tumour cells and V79 cells, respectively, whereas biological weighting function-based RBE is about 2.8. ICRP 60 and ICRU 40 recommendation-basedQ̅values are about 14.5 and 16, respectively. Theory of dual radiation action based RBE is 11.4. The calculated value of OER is 1.6.Significance.This study demonstrates the relative insensitivity of RBE,Q̅and OER radially away from the252Cf source along the distances of 1-5 cm in water.

Extended in-field and out-of-field validation of a compact Monte Carlo model of an IBA PROTEUS®ONE proton beam in TOPAS/GEANT4

Objective:This study evaluates a compact Monte Carlo (MC) model of a pencil beam scanning clinical proton beam using TOPAS to estimate the dose out-of-field (OOF). Compact modelling means that the model starts from a pristine proton beam at the nozzle exit, customised based on acceptance and commissioning data, instead of modelling the full treatment head and room.Approach: First, in-field validation tests were performed. Then, the OOF dose was validated in an RW3 phantom with bubble detectors for personal neutron dosimetry (measuring the neutron dose equivalent) and thermoluminiescent detectors (measuring the absorbed dose by protons and gammas). Measurements were performed at 15 and 35 cm from the distal edge of the field for five different irradiation plans, covering different beam orientations, proton energies and a 40 mm range shifter. TOPAS simulations were performed with QGSP Binary Cascade HP (BIC) and QGSP Bertini HP (Bertini) hadron physics lists.Main results: In-field validation shows that MC simulations agree with point dose measurements within -2.5 % and +1.5 % at locations on- and off-axis and before, in and after the Bragg peak or plateau. The gamma passing rate 2%/3mm of four simulated treatment plans compared to the dose distribution calculated by the TPS exceeds 97 % agreement score. OOF dose simulations showed an average overestimation of 27 % of the neutron dose equivalent for the BIC hadron physics list and an average underestimation of 20 % for the Bertini hadron physics list. The simulated absorbed dose of protons and gammas showed a systematic underestimation which was on average 21 % and 51 % for BIC and Bertini respectively.Significance: Our study demonstrates that a compact MC model can reliably produce in-field data, while out-of-field dose data are within the uncertainties of the detector systems and MC simulations nuclear models, and do so with shorter modelling and faster calculation time.

Monte Carlo simulation for proton dosimetry in magnetic fields: Fano test and magnetic field correction factorskBfor Farmer-type ionization chambers

Objective. In this contribution we present a special Fano test for charged particles in presence of magnetic fields in the MC code TOol for PArticle Simulation (TOPAS), as well as the determination of magnetic field correction factorskBfor Farmer-type ionization chambers using proton beams.Approach. Customized C++ extensions for TOPAS were implemented to model the special Fano tests in presence of magnetic fields for electrons and protons. The Geant4-specific transport parameters,DRoverRandfinalRange,were investigated to optimize passing rate and computation time. ThekBwas determined for the Farmer-type PTW 30013 ionization chamber, and 5 custom built ionization chambers with same geometry but varying inner radius, testing magnetic flux density ranging from 0 to 1.0 T and two proton beam energies of 157.43 and 221.05 MeV.Main results. Using the investigated parameters, TOPAS passed the Fano test within 0.39 ± 0.15% and 0.82 ± 0.42%, respectively for electrons and protons. The chamber response (kB,M,Q) gives a maximum at different magnetic flux densities depending of the chamber size, 1.0043 at 1.0 T for the smallest chamber and 1.0051 at 0.2 T for the largest chamber. The local dose differencecBremained ≤ 0.1% for both tested energies. The magnetic field correction factorkB, for the chamber PTW 30013, varied from 0.9946 to 1.0036 for both tested energies.Significance. The developed extension for the special Fano test in TOPAS MC code with the adjusted transport parameters, can accurately transport electron and proton particles in magnetic field. This makes TOPAS a valuable tool for the determination ofkB. The ionization chambers we tested showed thatkBremains small (≤0.72%). To the best of our knowledge, this is the first calculations ofkBfor proton beams. This work represents a significant step forward in the development of MRgPT and protocols for proton dosimetry in presence of magnetic field.

Fast DM,Mcalculation in LDR brachytherapy using deep learning methods

\textbf{Objective:} The Monte Carlo (MC) method provides a complete solution to the tissue heterogeneity effects in low-energy low-dose rate (LDR) brachytherapy. However, long computation times limit the clinical implementation of MC-based treatment planning solutions. This work aims to apply deep learning (DL) methods, specifically a model trained with MC simulations, to predict accurate dose to medium in medium (D\textsubscript{M,M}) distributions in LDR prostate brachytherapy.\\

\textbf{Approach:} To train the DL model, 2369 single-seed configurations, corresponding to 44 prostate patient plans, were used. These patients underwent LDR brachytherapy treatments in which $^{125}$I SelectSeed sources were implanted. For each seed configuration, the patient geometry, the MC dose volume and the single-seed plan volume were used to train a 3D Unet convolutional neural network. Previous knowledge was included in the network as an r\textsuperscript{2} kernel related to the first-order dose dependency in brachytherapy. MC and DL dose distributions were compared through the dose maps, isodose lines, and dose-volume histograms. Features enclosed in the model were visualized.\\

\textbf{Main results:} Model features started from the symmetrical kernel and finalized with an anisotropic representation that considered the patient organs and their interfaces, the source position, and the low- and high-dose regions. For a full prostate patient, small differences were seen below the 20\% isodose line. When comparing DL-based and MC-based calculations, the predicted CTV D\textsubscript{90} metric had an average difference of -0.1\%. Average differences for OARs were -1.3\%, 0.07\%, and 4.9\% for the rectum D\textsubscript{2cc}, the bladder D\textsubscript{2cc}, and the urethra D\textsubscript{0.1cc}. The model took 1.8 ms to predict a complete 3D D\textsubscript{M,M} volume (1.18 M voxels).\\

\textbf{Significance:} The proposed DL model stands for a simple and fast engine which includes prior physics knowledge of the problem. Such an engine considers the anisotropy of a brachytherapy source and the patient tissue composition.

Validation of the TOPAS Monte Carlo toolkit for LDR brachytherapy simulations

PURPOSE

This work has the purpose of validating the Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for low-dose-rate (LDR) brachytherapy uses.

METHODS AND MATERIALS

Simulations of 12 LDR sources and 2 COMS eye plaques (10 mm and 20 mm in diameter) and comparisons with published reference data from the Carleton Laboratory for Radiotherapy Physics (CLRP), the TG-43 consensus data and the TG-129 consensus data were performed. Sources from the IROC Houston Source Registry were modeled. The OncoSeed 6711 and the SelectSeed 130.002 were also modeled for historical reasons. For each source, the dose rate constant, the radial dose function and the anisotropy functions at 0.5, 1 and 5 cm were extracted. For the eye plaques (loaded with ¹²⁵I sources), dose distribution maps, dose profiles along the central axis and transverse axis were calculated.

RESULTS

Dose rate constants for 11 of the 12 sources are within 4% of the consensus data and within 2% of the CLRP data. The radial dose functions and anisotropy functions are mostly within 2% of the CLRP data. In average, 92% of all voxels are within 1% of the CLRP data for the eye plaques dose distributions. The dose profiles are within 0.5% (central axis) and 1% (transverse axis) of the reference data.

CONCLUSION

The TOPAS MC toolkit was validated for LDR brachytherapy applications. Single-seed and multi-seed results agree with the published reference data. TOPAS has several benefits such as a simplified approach to MC simulations and an accessible brachytherapy package including comprehensive learning resources.

Powder X-ray diffraction as a powerful tool to exploit in organic electronics: shedding light on the first N,N',N''-trialkyldiindolocarbazole

The first crystal structure of a fully N-alkylated diindolocarbazole derivative, namely, 5,8,14-tributyldiindolo[3,2-b;2',3'-h]carbazole (1, C₃₆H₃₉N₃), has been determined from laboratory powder X-ray diffraction (PXRD) data. A complex trigonal structure with a high-volume unit cell of 12987 ų was found, with a very long a(=b) [52.8790 (14) Å] and a very short c [5.36308 (13) Å] unit-cell parameter (hexagonal setting). The detailed analysis of the intermolecular interactions observed in the crystal structure of 1 highlights its potential towards the implementation of this core as a semiconductor in organic thin-film transistor (OTFT) devices. Since the molecule has a flat configuration reflecting its π-conjugated system, neighbouring molecules are found to stack atop each other in a slipped parallel fashion via π-π stacking interactions between planes of ca 3.30 Å, with a centroid-centroid distance between the aromatic rings corresponding to the shortest axis of the unit cell (i.e. c). The alkylation of the three N atoms proves to be a decisive feature since it favours the presence of C-H...π interactions in all directions, which strengthens the crystal packing. As a whole, PXRD proves to be a valuable option for the resolution of otherwise inaccessible organic crystal structures of interest in different areas.

Structure modeling and quantitative X-ray diffraction of C-(A)-S-H

Quantitative X-ray diffraction of nanocrystalline calcium silicate hydrate (C-S-H) and its aluminium-substituted variants (C-A-S-H) has so far been limited by a lack of appropriate structure models. In this study, atomistic structure models derived from tobermorite were combined with a supercell approach using TOPAS. By accounting for nanostructural features such as isolated layers, turbostratic disorder and, for the first time, fibrils, characteristic reflections and asymmetric bands were more accurately simulated than before, providing the means for phase quantification and refinement of structural sites. This improved methodology is applied to autoclaved aerated concrete and the experimental study of related hydrothermal reactions. Scanning electron microscopy indicated a fibrillar morphology for intermediate C-(A)-S-H, and energy-dispersive X-ray spectroscopy constrained its Ca/Si ratio to 1.31-1.35. As a first step, the direct quantification of C-(A)-S-H via structure models was assessed by a series of X-ray diffraction measurements using corundum as an internal standard. Secondly, the verified structure model was applied to evaluate in situ X-ray diffraction experiments at 457, 466 and 473 K (1.1, 1.35 and 1.55 MPa, respectively). Finally, a quantitative study of industrially produced autoclaved aerated concrete was conducted, determining 20-30 wt% C-(A)-S-H at Ca/Si ratios < 1.0. In general, the developed structure models advance the study of Portland cement concrete and related materials, including autoclaved aerated concrete, and the supercell approach may be universally applicable to other nanocrystalline materials.

Monte Carlo simulation of Cherenkov imaging for Total Skin Electron Treatment with CT DICOM realistic patient geometry

This Monte Carlo (MC) simulation study provides an evaluation of dose uniformity in a patient and the difference between dose and Cherenkov distributions, which is invaluable in developing conversion factors to relate observed Cherenkov images to actual dose distributions for TSET patients. This MC simulations with TOPAS is performed using realistic patient geometries obtained with a 3D scanner during total skin electron treatments (TSET) at UPenn. For each treatment posture in the Stanford technique, the differences between Cherenkov photon distributions and dose distributions produced in MC are consistent with the differences observed between a Cherenkov imaging camera and in-vivo dose measurement with OSLD on patient skin. According to MC studies of a flat rectangular PVC board, the difference between Cherenkov and dose is mostly due to the spoiler. This is confirmed by observing consistent dose and Cherenkov distributions in clinical measurements on a PVC board without the spoiler. The accumulated dose and Cherenkov distributions for each patient are obtained by projecting the MC output of the 6 postures of the TSET treatment together onto a finite element model of the patient.

Biological Impact of Target Fragments on Proton Treatment Plans: An Analysis Based on the Current Cross-Section Data and a Full Mixed Field Approach

Simple Summary

Proton therapy is now an established external radiotherapy modality for cancer treatment. Clinical routine currently neglects the radiobiological impact of nuclear target fragments even if experimental evidences show a significant enhancement in cell-killing effect due to secondary particles. This paper quantifies the contribution of proton target fragments of different charge in different irradiation scenarios and compares the computationally predicted corrections to the overall biological dose with experimental data.

Abstract

Clinical routine in proton therapy currently neglects the radiobiological impact of nuclear target fragments generated by proton beams. This is partially due to the difficult characterization of the irradiation field. The detection of low energetic fragments, secondary protons and fragments, is in fact challenging due to their very short range. However, considering their low residual energy and therefore high LET, the possible contribution of such heavy particles to the overall biological effect could be not negligible. In this context, we performed a systematic analysis aimed at an explicit assessment of the RBE (relative biological effectiveness, i.e., the ratio of photon to proton physical dose needed to achieve the same biological effect) contribution of target fragments in the biological dose calculations of proton fields. The TOPAS Monte Carlo code has been used to characterize the radiation field, i.e., for the scoring of primary protons and fragments in an exemplary water target. TRiP98, in combination with LEM IV RBE tables, was then employed to evaluate the RBE with a mixed field approach accounting for fragments’ contributions. The results were compared with that obtained by considering only primary protons for the pristine beam and spread out Bragg peak (SOBP) irradiations, in order to estimate the relative weight of target fragments to the overall RBE. A sensitivity analysis of the secondary particles production cross-sections to the biological dose has been also carried out in this study. Finally, our modeling approach was applied to the analysis of a selection of cell survival and RBE data extracted from published in vitro studies. Our results indicate that, for high energy proton beams, the main contribution to the biological effect due to the secondary particles can be attributed to secondary protons, while the contribution of heavier fragments is mainly due to helium. The impact of target fragments on the biological dose is maximized in the entrance channels and for small α/β values. When applied to the description of survival data, model predictions including all fragments allowed better agreement to experimental data at high energies, while a minor effect was observed in the peak region. An improved description was also obtained when including the fragments’ contribution to describe RBE data. Overall, this analysis indicates that a minor contribution can be expected to the overall RBE resulting from target fragments. However, considering the fragmentation effects can improve the agreement with experimental data for high energy proton beams.

Dual-Frequency Polarized Reconfigurable Terahertz Antenna Based on Graphene Metasurface and TOPAS

A hybrid dual-frequency polarized reconfigurable terahertz antenna is designed and studied. Graphene and TOPAS are employed as the polarization conversion metasurface and dielectric substrate, respectively, enabling tunable polarization conversion and circular polarization. TOPAS is a good substrate material for broadband THz components due to its low absorption. By adjusting the chemical potential of graphene between 0 eV and 0.5 eV, the polarization state in the band of 1 THz (0.76-1.02 THz) and 2.5 THz (2.43-2.6 THz) can be reconstructed. Thanks to the multilayer graphene structure and low absorption TOPAS, the graphene metasurface exhibits a broad bandwidth of 0.26 and 0.17 THz, respectively, in the band of 1 THz and 2.5 THz. The working state of the circularly polarized antenna and linearly polarized antenna can be switched in the bands around 1 THz (0.7-0.75 THz, 0.96-1.04 THz) and 2.5 THz (2.42-2.52 THz), respectively, without changing the physical geometry. Moreover, the graphene antenna, metasurface, and hybrid structure are tested, respectively, to verify that the components do not interfere with each other in performance. The hybrid antenna shows great potential in tunable terahertz devices and related applications.

Investigating the feasibility of TOPAS-nBio for Monte Carlo track structure simulations by adapting GEANT4-DNA examples application

Purpose.The purpose of this work is to investigate the feasibility of TOPAS-nBio for track structure simulations using tuple scoring and ROOT/Python-based post-processing.Materials and methods.There are several example applications implemented in GEANT4-DNA demonstrating track structure simulations. These examples are not implemented by default in TOPAS-nBio. In this study, the tuple scorer was used to re-simulate these examples. The simulations contained investigations of different physics lists, calculation of energy-dependent range, stopping power, mean free path andW-value. Additionally, further applications of the TOPAS-nBio tool were investigated, focusing on physical interactions and deposited energies of electrons with initial energies in the range of 10-60 eV, not covered in the recently published GEANT4-DNA simulations. Low-energetic electrons are currently of great interest in the radiobiology research community due to their high effectiveness towards the induction of biological damage.Results.The quantities calculated with TOPAS-nBio show a good agreement with the simulations of GEANT4-DNA with deviations of 5% at maximum. Thus, we have presented a feasible way to implement the example applications included in GEANT4-DNA in TOPAS-nBio. With the extended simulations, an insight could be given, which further tracking information can be gained with the track structure code and how cross sections and physics models influence a particle's fate.Conclusion.With our results, we could show the potentials of applying the tuple scorer in TOPAS-nBio Monte Carlo track structure simulations. Using this scorer, a large amount of information about the track structure can be accessed, which can be analyzed as preferred after the simulation.

The relation between microdosimetry and induction of direct damage to DNA by alpha particles

In radiopharmaceutical treatmentsα-particles are employed to treat tumor cells. However, the mechanism that drives the biological effect induced is not well known. Being ionizing radiation,α-particles can affect biological organisms by producing damage to the DNA, either directly or indirectly. Following the principle that microdosimetry theory accounts for the stochastic way in which radiation deposits energy in sub-cellular sized volumes via physical collisions, we postulate that microdosimetry represents a reasonable framework to characterize the statistical nature of direct damage induction byα-particles to DNA. We used the TOPAS-nBio Monte Carlo package to simulate direct damage produced by monoenergetic alpha particles to different DNA structures. In separate simulations, we obtained the frequency-mean lineal energy (yF) and dose-mean lineal energy (yD) of microdosimetric distributions sampled with spherical sites of different sizes. The total number of DNA strand breaks, double strand breaks (DSBs) and complex strand breaks per track were quantified and presented as a function of eitheryForyD.The probability of interaction between a track and the DNA depends on how the base pairs are compacted. To characterize this variability on compactness, spherical sites of different size were used to match these probabilities of interaction, correlating the size-dependent specific energy (z) with the damage induced. The total number of DNA strand breaks per track was found to linearly correlate withyFandzFwhen using what we defined an effective volume as microdosimetric site, while the yield of DSB per unit dose linearly correlated withyDorzD,being larger for compacted than for unfolded DNA structures. The yield of complex breaks per unit dose exhibited a quadratic behavior with respect toyDand a greater difference among DNA compactness levels. Microdosimetric quantities correlate with the direct damage imparted on DNA.

Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations

PURPOSE

The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications.

METHODS AND MATERIALS

Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. The four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated.

RESULTS

The air-kerma strength was (9.772 ± 0.001) × 10-8 U Bq-1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h-1 U-1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data.

CONCLUSIONS

TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.

Radiobiological impact of gadolinium neutron capture from proton therapy and alternative neutron sources using TOPAS-nBio

PURPOSE

A multi-scale investigation of the biological properties of gadolinium neutron capture (GdNC) therapy with applications in particle therapy is conducted using the TOPAS Monte Carlo (MC) simulation code. The simulation results are used to quantify the amount of gadolinium dose enhancement produced as a result of the secondary neutron production from proton therapy scaled by measured data.

MATERIALS AND METHODS

MC modeling was performed using the radiobiology extension TOol for PArticle Simulation TOPAS-nBio MC simulation code to study the radiobiological effects produced from GdNC on a segment of DNA, a spherical cellular model, and from the modeling of previous experimental measurements. The average RBE values were calculated from two methods, microdosimetric kinematic (MK) and biological weighting r(y) within a 2 nm DNA segment for GdNC. The single-strand breaks (SSBs) and double-strand breaks (DSBs) were calculated from within the nucleus of a 20 µm diameter, spherical cell model. From a previous experimental proton therapy measurement using a spread-out Bragg peak (SOBP) of 4.5-9.5 cm and a delivered absorbed dose of 10.4 Gy, the amount of Gd neutron captures was calculated and used to quantify the amount of GdNC absolute dose from particle therapy.

RESULTS

The average RBE from microdosimetric kinematic and biological weighting was 1.35, and 1.70 for a 10% cell survival on HSG cell-line and weighting function data from early intestinal tolerance of mice. From a central isotropic GdNC source, the energy deposition is found to decrease from roughly 2.7 eV per capture down to approximately 0.01 eV per capture, a drop of two orders of magnitude within 50 nm. This result suggests that Gd needs to be close to the DNA (within 10-20 nm) in order for neutron capture to induce a significant dose enhancement due to the short-range electrons emitted after Gd neutron capture. Within a spherical cell model, the SSBs, and DSBs were determined to be 39 and 1.5 per neutron capture, respectively. From the total neutron captures produced from an experimental proton therapy measurement on a 3000 PPM Gd solution, an insignificant absolute Gd dose enhancement was quantified to be 5.4 × 10^-6  Gy per Gy of administered proton dose.

CONCLUSION

From this study and literature review, the production of secondary thermal neutrons from proton therapy is determined to be a limiting factor and unlikely to produce a clinically useful dose enhancement for secondary neutron capture therapy. Moreover, alternative neutron sources, such as, a compact deuterium-tritium (D-T) neutron generator, a "high yield" deuterium-deuterium (D-D) generator, or an industrial strength (100 mg) ²⁵² Cf source were investigated, with the ²⁵² Cf source the most likely to be capable of producing enough neutrons for 1 Gy of localized GdNC absolute dose within a reasonable treatment time.

Experiments and Monte Carlo simulations on multiple Coulomb scattering of protons

BACKGROUND AND PURPOSE

Monte Carlo simulations as well as analytical computations of proton transport in material media require accurate values of multiple Coulomb scattering (MCS) angles. High-quality experimental data on MCS angles in the energy range for proton therapy are, however, sparse. In this work, MCS modeling in proton transport was evaluated employing an experimental method to measure these angles on a medical proton beamline in clinically relevant materials. Results are compared to Monte Carlo simulations and analytical models.

MATERIALS AND METHODS

Aluminum, brass, and lucite (PMMA) scatterers of clinically relevant thicknesses were irradiated with protons at 100, 160, and 220 MeV. Resulting spatial distributions of individual pencil beams were measured with a scintillating screen. The MCS angles were determined by deconvolution and a virtual point source approach. Results were compared to those obtained with the Monte Carlo codes PENH, TOPAS, and RayStation Monte Carlo, as well as the analytical models RayStation Pencil Beam Algorithm and the Molière/Fano/Hanson variant of the Molière theory.

RESULTS

Experimental data obtained with the presented methodology agree with previously published results within 6%, with an average deviation of 3%. The combined average uncertainty of the experimental data yielded 1.8%, while the combined maximum uncertainty was below 4%. The obtained Monte Carlo results for PENH, TOPAS, and RayStation deviate on average for all considered energies, materials and thicknesses, by 2.5%, 3.4%, and 2.8% from the experimental data, respectively. For the analytical models, the average deviations amount to 4.5% and 2.9% for the RayStation Pencil Beam Algorithm and the Molière/Fano/Hanson model, respectively.

CONCLUSION

The experimental method developed for the present work allowed to measure MCS angles in clinical proton facilities with good accuracy. The presented method permits to extend the database on experimental MCS angles which is rather limited. This work further provides benchmark data for lucite in thicknesses relevant for clinical applications. The data may serve to validate dose engines of treatment planning systems and secondary dose check software. The Monte Carlo and analytical algorithms studied are capable of reproducing MCS data within the required accuracy for clinical applications.

Development and validation of the Dynamic Collimation Monte Carlo simulation package for pencil beam scanning proton therapy

PURPOSE

The aim of this work was to develop and experimentally validate a Dynamic Collimation Monte Carlo (DCMC) simulation package specifically designed for the simulation of collimators in pencil beam scanning proton therapy (PBS-PT). The DCMC package was developed using the TOPAS Monte Carlo platform and consists of a generalized PBS source model and collimator component extensions.

METHODS

A divergent point-source model of the IBA dedicated nozzle (DN) at the Miami Cancer Institute (MCI) was created and validated against on-axis commissioning measurements taken at MCI. The beamline optics were mathematically incorporated into the source to model beamlet deflections in the X and Y directions at the respective magnet planes. Off-axis measurements taken at multiple planes in air were used to validate both the off-axis spot size and divergence of the source model. The DCS trimmers were modeled and incorporated as TOPAS geometry extensions that linearly translate and rotate about the bending magnets. To validate the collimator model, a series of integral depth dose (IDD) and lateral profile measurements were acquired at MCI and used to benchmark the DCMC performance for modeling both pristine and range shifted beamlets. The water equivalent thickness (WET) of the range shifter was determined by quantifying the shift in the depth of the 80% dose point distal to the Bragg peak between the range shifted and pristine uncollimated beams.

RESULTS

A source model of the IBA DN system was successfully commissioned against on- and off-axis IDD and lateral profile measurements performed at MCI. The divergence of the source model was matched through an optimization of the source-to-axis distance and comparison against in-air spot profiles. The DCS model was then benchmarked against collimated IDD and in-air and in-phantom lateral profile measurements. Gamma analysis was used to evaluate the agreement between measured and simulated lateral profiles and IDDs with 1%/1 mm criteria and a 1% dose threshold. For the pristine collimated beams, the average 1%/1 mm gamma pass rates across all collimator configurations investigated were 99.8% for IDDs and 97.6% and 95.2% for in-air and in-phantom lateral profiles. All range shifted collimated IDDs passed at 100% while in-air and in-phantom lateral profiles had average pass rates of 99.1% and 99.8%, respectively. The measured and simulated WET of the polyethylene range shifter was determined to be 40.9 and 41.0 mm, respectively.

CONCLUSIONS

We have developed a TOPAS-based Monte Carlo package for modeling collimators in PBS-PT. This package was then commissioned to model the IBA DN system and DCS located at MCI using both uncollimated and collimated measurements. Validation results demonstrate that the DCMC package can be used to accurately model other aspects of a DCS implementation via simulation.

Monte Carlo calculation of perturbation correction factors for air-filled ionization chambers in clinical proton beams using TOPAS/GEANT

INTRODUCTION

Current dosimetry protocols for clinical protons using air-filled ionization chambers assume that the perturbation correction factor is equal to unity for all ionization chambers and proton energies. Since previous Monte Carlo based studies suggest that perturbation correction factors might be significantly different from unity this study aims to determine perturbation correction factors for six plane-parallel and four cylindrical ionization chambers in proton beams at clinical energies.

MATERIALS AND METHODS

The dose deposited in the air cavity of the ionization chambers was calculated with the help of the Monte Carlo code TOPAS/Geant4 while specific constructive details of the chambers were removed step by step. By comparing these dose values the individual perturbation correction factors p_cel, p_stem, p_sleeve, p_wall, p_cav⋅p_dis as well as the total perturbation correction factor p_Q were derived for typical clinical proton energies between 80 and 250MeV.

RESULTS

The total perturbation correction factor p_Q was smaller than unity for almost every ionization chamber and proton energy and in some cases significantly different from unity (deviation larger than 1%). The maximum deviation from unity was 2.0% for cylindrical and 1.5% for plane-parallel ionization chambers. Especially the factor p_wall was found to differ significantly from unity. It was shown that this is due to the fact that secondary particles, especially alpha particles and fragments, are scattered from the chamber wall into the air cavity resulting in an overresponse of the chamber.

CONCLUSION

Perturbation correction factors for ionization chambers in proton beams were calculated using Monte Carlo simulations. In contrast to the assumption of current dosimetry protocols the total perturbation correction factor p_Q can be significantly different from unity. Hence, beam quality correction factors [Formula: see text] that are calculated with the help of perturbation correction factors that are assumed to be unity come with a corresponding additional uncertainty.

Monte Carlo (MC) study of dose distribution and Cherenkov imaging in total skin electron therapy (TSET) with TOPAS

Malignant tissues can be effectively treated by Total Skin Electron Therapy (TSET) over the entire body surface using 6 MeV electron beams. During the radiation treatment, Cherenkov photons are emitted from the patient's skin, and can potentially be used for in-vivo imaging of the radiation dose distribution. A Monte Carlo (MC) simulation toolkit TOPAS is used to study the generation and propagation of Cherenkov photons that are generated from the interaction of electron radiation with human tissues, and to understand the relationship between the dose distributions and the Cherenkov photon distributions. Validation of MC simulations with experiments are performed at 100 SSD and 500 SSD, and simulations of a patient phantom in realistic clinical treatment setups have been done. These simulations with TOPAS show that the emitted Cherenkov distributions at phantom surfaces closely follow their corresponding dose distributions.

Ab initio structure solution of proteins at atomic resolution using charge-flipping techniques and cloud computing

Large protein structures at atomic resolution can be solved in minutes using charge-flipping techniques operating on hundreds of virtual machines (computers) on the Amazon Web Services cloud-computing platform driven by the computer programs TOPAS or TOPAS-Academic at a small financial cost. The speed of operation has allowed charge-flipping techniques to be investigated and modified, leading to two strategies that can solve a large range of difficult protein structures at atomic resolution. Techniques include the use of space-group symmetry restraints on the electron density as well as increasing the intensity of a randomly chosen high-intensity electron-density peak. It is also shown that the use of symmetry restraints increases the chance of finding a solution for low-resolution data. Finally, a flipping strategy that negates `uranium atom solutions' has been developed for structures that exhibit such solutions during charge flipping.

Impact of air gap, range shifter, and delivery technique on skin dose in proton therapy

OBJECTIVE

Side effects of radiation therapy may include skin damage. The surface dose is of great interest and contains the buildup effect. In particular, the proton therapy community requires further experimental data to quantify doses in the surface region. This specification includes the skin dose, which is defined according to ICRU Report No. 39 at 70 μm water equivalent depth. The aim of this study is to gather more knowledge of the skin dose by varying key parameters defined by the patient treatment plan. This consists of clinical aspects such as the influence of the air gap, the application of a range shifter (RS), or the proton delivery technique.

MATERIAL/METHODS

Skin doses were determined with a PTW 23391 extrapolation chamber with three thin Kapton® entrance windows operated as a conventional ionization chamber. The impact on the skin dose for quasi-monoenergetic pencil beam scanning (PBS) proton beams was evaluated for clinical air gaps between 3.5 and 51.1 cm. The differences in skin dose were assessed by irradiating equivalent fields with an RS of 51 mm water equivalent thickness (RS51) and without. Furthermore, the delivery techniques PBS, uniform scanning (US), and double scattering (DS) were compared by defining a spread-out Bragg peak (SOBP). TOPAS (V.3.1.2) was used to model an IBA nozzle with PBS and to score dose to water at the surface of a water phantom.

RESULTS

For the monoenergetic fields without the application of the RS the skin dose was constant down to an air gap of 6.2 cm. A lower air gap of 3.5 cm showed a variation in skin dose by up to 2.4% compared to the results obtained with larger air gaps. With the inserted RS51 an increase in the skin dose was found for air gaps smaller than 11.3 cm. Experimentally, a dose difference of 1.4% was recorded for an air gap of 6.2 cm by inserting an RS and none. With the Monte Carlo calculations the largest dose increase was observed at the air gap of 3.5 cm with 1.7% and 4.0% relative to the skin dose results without the RS and to the largest evaluated air gap of 51.1 cm, respectively. The SOBP comparison of the beam modalities at the measuring plane at the isocenter revealed higher skin doses without RS (including RS) by up to +1.9% (+1.5%) for DS and +1.3% (+1.1%) for US compared to PBS. For all three techniques an approx. 2% rise in skin dose was observed for the largest evaluated air gap of 37.7 cm to an air gap of 6.2 cm when using an RS51.

CONCLUSION

The study investigated aspects of skin dose of a water equivalent phantom by varying key parameters of a proton treatment plan. Parameters like the RS, the air gap, and the delivery modality have an impact on the order of 4.0% for the skin dose at the depth of 70 μm. The increases in skin dose are the effects of the contribution of the increased electron fluence at small air gaps and the emitted hadronic particles produced by the RS.

Single pencil beam benchmark of a module for Monte Carlo simulation of proton transport in the PENELOPE code

BACKGROUND AND PURPOSE

PENH is a recently coded module for simulation of proton transport in conjunction with the Monte Carlo code PENELOPE. PENELOPE applies class II simulation to all type of interactions, in particular, to elastic collisions. PENH uses calculated differential cross sections for proton elastic collisions that include electron screening effects as well as nuclear structure effects. Proton-induced nuclear reactions are simulated from information in the ENDF-6 database or from alternative nuclear databases in ENDF format. The purpose of this work is to benchmark this module by simulating absorbed dose distributions from a single finite spot size proton pencil beam in water.

MATERIALS AND METHODS

Monte Carlo simulations with PENH are compared with simulation results from TOPAS Monte Carlo (v3.1p2) and RayStation Monte Carlo (v6). Different beam models are examined in terms of mean energy and energy spread to match the measured profiles. The phase-space file is derived from experimental measurements. Simulated absorbed dose distributions are compared to experimental data obtained with the ionization chamber array MatriXX 2D detector (IBA Dosimetry) in a water tank. The experiments were conducted with a clinical IBA pencil beam scanning dedicated nozzle. In all simulations a Fermi-Eyges phase-space representation of a single finite spot size proton pencil beam is used.

RESULTS

In general, there is a good agreement between simulated results and experimental data up to a distance of 3 cm from the central axis. In the core region (region where the dose is more than 10% of the maximum dose) PENH shows, overall, the smallest deviations from experimental data, with the largest radial rms (root mean square) smaller than 0.2. The results achieved by TOPAS and RayStation in that region are very close to those of PENH. For the halo region, that is the area of the dose distribution outside the core region reaching down to 0.01% of the maximum intensity, the largest rms achieved by TOPAS is always smaller than 0.5, yielding better results than the rest of the codes.

CONCLUSION

The physics modeling of the PENELOPE/PENH code yields results consistent with measurements in the dose range relevant for proton therapy. The discrepancies between PENH appearing at distances larger than 3 cm from the central-beam axis are accountable to the lack of neutron simulation in this code. In contradistinction, TOPAS has a better agreement with experimental data at large distances from the central-beam axis because of the simulation of neutrons.

Producing a Beam Model of the Varian ProBeam Proton Therapy System using TOPAS Monte Carlo Toolkit

PURPOSE

A Geant4-based TOPAS Monte Carlo toolkit was utilized to model a Varian ProBeam proton therapy system, with the aim of providing an independent computational platform for validating advanced dosimetric methods.

MATERIALS AND METHODS

The model was tested for accuracy of dose and linear energy transfer (LET) prediction relative to the commissioning data, which included integral depth dose (IDD) in water and spot profiles in air measured at varying depths (for energies of 70 to 240 MeV in increments of 10 MeV, and 242 MeV), and absolute dose calibration. Emittance was defined based on depth-dependent spot profiles and Courant-Snyder's particle transport theory, which provided spot size and angular divergence along the inline and crossline plane. Energy spectra were defined as Gaussian distributions that best matched the range and maximum dose of the IDD. The validity of the model was assessed based on measurements of range, dose to peak difference, mean point to point difference, spot sizes at different depths, and spread-out Bragg peak (SOBP) IDD and was compared to the current treatment planning software (TPS).

RESULTS

Simulated and commissioned spot sizes agreed within 2.5%. The single spot IDD range, maximum dose, and mean point to point difference of each commissioned energy agreed with the simulated profiles generally within 0.07 mm, 0.4%, and 0.6%, respectively. A simulated SOBP plan agreed with the measured dose within 2% for the plateau region. The protons/MU and absolute dose agreed with the current TPS to within 1.6% and exhibited the greatest discrepancy at higher energies.

CONCLUSIONS

The TOPAS model agreed well with the commissioning data and included inline and crossline asymmetry of the beam profiles. The discrepancy between the measured and TOPAS-simulated SOBP plan may be due to beam modeling simplifications of the current TPS and the nuclear halo effect. The model can compute LET, and motivates future studies in understanding equivalent dose prediction in treatment planning, and investigating scintillation quenching.

Investigation of the effect of air gap size on the spatial resolution in proton- and helium radio- and tomography

PURPOSE

Proton computed (transmission) tomography (pCT) refers to the process of imaging an object by letting protons pass through it, while measuring their energy after, and their position and (optionally) direction both before and after their traversal through that object. The so far experimental technique has potential to improve treatment planning of proton therapy by enabling the direct acquisition of a proton stopping power map of tissue, thus removing the need to obtain it by converting X-ray CT attenuation data and thereby eliminating uncertainties which arise in the mentioned conversion process. The image reconstruction in pCT requires accurate estimates of the proton trajectories. In experimental pCT detector setups where the direction of the protons is not measured, the air gap between the detector planes and the imaged object worsens the spatial resolution of the image obtained. In this work we determined the mean proton paths and the corresponding spatial uncertainty, taking into account the presence of the air gap.

METHODS

We used Monte Carlo simulations of radiation transport to systematically investigate the effect of the air gap size between detector and patient on the spatial resolution of proton (ion) computed tomography for protons with an energy of 200MeV and 250MeV as well as for helium ions (He-4) with an energy of 798MeV. For the simulations we used TOPAS which itself is based on Geant4.

RESULTS

For all particles, which are detected at the same entrance and exit coordinate, the average ion path and the corresponding standard deviation was computed. From this information, the dependence of the spatial resolution on the air gap size and the angular confusion of the particle beam was inferred.

CONCLUSION

The presence of the airgap does not pose a problem for perfect fan beams. In realistic scenarios, where the initial angular confusion is around 5mrad and for typical air gap sizes up to 10cm, using an energy of 200MeV a spatial resolution of about 1.6mm can be achieved. Using protons with E=250MeV a spatial resolution of about 1.1mm and using helium ions (He-4) with E=798MeV even a spatial resolution below 0.7mm respectively is attainable.

Validation of a Monte Carlo model for multi leaf collimator based electron delivery

PURPOSE

To develop and validate a Monte Carlo model of the Varian TrueBeam to study electron collimation using the existing photon multi-leaf collimators (pMLC), instead of conventional electron applicators and apertures.

MATERIALS AND METHODS

A complete Monte Carlo model of the Varian TrueBeam was developed using Tool for particle simulation (TOPAS) (version 3.1.p3). Vendor-supplied information was used to model the treatment head components and the source parameters. A phase space plane was setup above the collimating jaws and captured particles were reused until a statistical uncertainty of 1% was achieved in the central axis. Electron energies 6, 9, 12, 16, and 20 MeV with a jaw-defined field of 20 × 20 cm² at iso-center, pMLC-defined fields of 6.8 × 6.8 cm² and 11.4 × 11.4 cm² at 80 cm source-to-surface distance (SSD) and an applicator-defined field of 10 × 10 cm² at iso-center were evaluated. All the measurements except the applicator-defined fields were measured using an ionization chamber in a water tank using 80 cm SSD. The dose difference, distance-to-agreement and gamma index were used to evaluate the agreement between the Monte Carlo calculations and measurements. Contributions of electron scattering off pMLC leaves and inter-leaf leakage on dose profiles were evaluated and compared with Monte Carlo calculations. Electron transport through a heterogeneous phantom was simulated and the resulting dose distributions were compared with film measurements. The validated Monte Carlo model was used to simulate several clinically motivated cases to demonstrate the benefit of pMLC-based electron delivery compared to applicator-based electron delivery.

RESULTS

Calculated and measured percentage depth-dose (PDD) curves agree within 2% after normalization. The agreement between normalized percentage depth dose curves were evaluated using one-dimensional gamma analysis with a local tolerance of 2%/1 mm and the %points passing gamma criteria was 100% for all energies. For jaw-defined fields, calculated profiles agree with measurements with pass rates of >97% for 2%/2 mm gamma criteria. Calculated FWHM and penumbra width agree with measurements within 0.4 cm. For fields with tertiary collimation using an pMLC or applicator, the average gamma pass rate of compared profiles was 98% with 2%/2 mm gamma criteria. The profiles measured to evaluate the pMLC leaf scattering agreed with Monte Carlo calculations with an average gamma pass rate of 96.5% with 3%/2 mm gamma criteria. Measured dose profiles below the heterogenous phantom agreed well with calculated profiles and matched within 2.5% for most points. The calculated clinically applicable cases using TOPAS MC and Eclipse TPS for single enface electron beam, electron-photon mixed beam and a matched electron-electron beam exhibited a reasonable agreement in PDDs, profiles and dose volume histograms.

CONCLUSION

We present a validation of a Monte Carlo model of Varian TrueBeam for pMLC-based electron delivery. Monte Carlo calculations agreed with measurements satisfying gamma criterion of 1%/1 mm for depth dose curves and 2%/1 mm for dose profiles. The simulation of clinically applicable cases demonstrated the clinical utility of pMLC-based electrons and the use of MC simulations for development of advanced radiation therapy techniques.

Design of a focused collimator for proton therapy spot scanning using Monte Carlo methods

PURPOSE

When designing a collimation system for pencil beam spot scanning proton therapy, a decision must be made whether or not to rotate, or focus, the collimator to match beamlet deflection as a function of off-axis distance. If the collimator is not focused, the beamlet shape and fluence will vary as a function of off-axis distance due to partial transmission through the collimator. In this work, we quantify the magnitude of these effects and propose a focused dynamic collimation system (DCS) for use in proton therapy spot scanning.

METHODS

This study was done in silico using a model of the Miami Cancer Institute's (MCI) IBA Proteus Plus system created in Geant4-based TOPAS. The DCS utilizes rectangular nickel trimmers mounted on rotating sliders that move in synchrony with the pencil beam to provide focused collimation at the edge of the target. Using a simplified setup of the DCS, simulations were performed at various off-axis locations corresponding to beam deflection angles ranging from 0° to 2.5°. At each off-axis location, focused (trimmer rotated) and unfocused (trimmer not rotated) simulations were performed. In all simulations, a 4 cm water equivalent thickness range shifter was placed upstream of the collimator, and a voxelized water phantom that scored dose was placed downstream, each with 4 cm airgaps.

RESULTS

Increasing the beam deflection angle for an unfocused trimmer caused the collimated edge of the beamlet profile to shift 0.08-0.61 mm from the baseline 0° simulation. There was also an increase in low-dose regions on the collimated edge ranging from 14.6% to 192.4%. Lastly, the maximum dose, D max , was 0-5% higher for the unfocused simulations. With a focused trimmer design, the profile shift and dose increases were all eliminated.

CONCLUSIONS

We have shown that focusing a collimator in spot scanning proton therapy reduces dose at the collimated edge compared to conventional, unfocused collimation devices and presented a simple, mechanical design for achieving focusing for a range of source-to-collimator distances.

Determination of surface dose in pencil beam scanning proton therapy

PURPOSE/OBJECTIVE

Quantification of surface dose within the first few hundred water equivalent µm is challenging. Nevertheless, it is of large interest for the proton therapy community to study dose effects in the skin. The experimental determination is affected by the detector properties, such as the detector volume and material. The International Commission on Radiation Units and Measurements in its report 39 recommends assessing the skin dose at a depth of 0.07 mm. The aim of this study is the estimation of the absorbed dose at and around a depth of 70 µm. We used various dosimetric approaches in conjunction with proton pencil beam scanning delivery to determine the skin dose in a clinical setting.

MATERIAL/METHODS

Five different detectors were tested for determining the surface dose in water: EBT3 and HD-V2 GAFCHROMIC™ radiochromic film, LiF:Mg,Ti thermoluminescent dosimeter, IBA PPC05 plane-parallel ionization chamber, and PTW 23391 extrapolation chamber. The irradiation setup consisted of quasi-monoenergetic scanned proton pencil beams with kinetic energies of 100, 150, and 226.7 MeV, respectively. Radiochromic films were placed within a vertical stack and in wedge geometry and were analyzed with FilmQA Pro™ adopting triple channel dosimetry. The extrapolation chamber PTW 23391, which served as a reference in the current work, was used in a conventional ionization chamber setup with a fixed electrode gap of 2 mm. Three Kapton® entrance windows with thicknesses of 25, 50, and 75 µm were employed. Thermoluminescent dosimeters were provided as powder and were pressed onto a sheet of aluminum. Furthermore, the Monte Carlo code TOol for PArticle Simulation (TOPAS) in version 3.1.p2 was used to model an IBA pencil beam scanning nozzle and score dose to water in a water phantom.

RESULTS

The resulting depth dose curves were normalized to their 100% dose at the reference depth of 3 cm. We obtained the skin doses with the extrapolation chamber and with TOPAS. For the experimental approach this resulted in 79.7 ± 0.3%, 86.0 ± 0.6%, and 87.1 ± 0.1% for the proton energies 100, 150, and 226.7 MeV, respectively. The results for TOPAS were 80.1 ± 0.2% (100 MeV), 87.1 ± 0.5% (150 MeV), and 86.9 ± 0.4% (226.7 MeV), respectively. Based on the experimental results of the skin dose, we provided a clinically relevant surface extrapolation factor for the common measurement methods. This allows the result of the first measurement depth of a detector to be scaled to the dose at the skin depth. Most practical would be the use of the surface extrapolation factor for the PPC05 chamber, due to its direct reading, the wide availability in clinics and the low uncertainties. The calculated factors were 0.986 ± 0.004 for 100 MeV, 0.961 ± 0.008 for 150 MeV, and 0.963 ± 0.003 for 226.7 MeV.

CONCLUSIONS

In this study, dissimilar experimental approaches were evaluated with respect to measurements at depths close to the surface. The experimental depth dose curves are in good agreement with the simulation with TOPAS Monte Carlo. To the author's knowledge this was the first experimental determination of the skin dose according to the International Commission on Radiation Units and Measurements 39 definition in proton pencil beam scanning.

A routine for the determination of the microstructure of stacking-faulted nickel cobalt aluminium hydroxide precursors for lithium nickel cobalt aluminium oxide battery materials

Stacking faults and interstratification faults in a cobalt- and aluminium-bearing nickel layered double hydroxide used as a precursor for Li(Ni_1−x−yCo_xAl_y)O₂ battery materials were quantified by a combination of a grid-search approach and a recursive routine for generating and averaging supercells of stacking-faulted layered substances implemented in the TOPAS software.

The microstructures of six stacking-faulted industrially produced cobalt- and aluminium-bearing nickel layered double hydroxide (LDH) samples that are used as precursors for Li(Ni_1−x−yCo_xAl_y)O₂ battery materials were investigated. Shifts from the brucite-type (AγB)□(AγB)□ stacking pattern to the CdCl₂-type (AγB)□(CβA)□(BαC)□ and the CrOOH-type (BγA)□(AβC)□(CαB)□ stacking order, as well as random intercalation of water molecules and carbonate ions, were found to be the main features of the microstructures. A recursive routine for generating and averaging supercells of stacking-faulted layered substances implemented in the TOPAS software was used to calculate diffraction patterns of the LDH phases as a function of the degree of faulting and to refine them against the measured diffraction data. The microstructures of the precursor materials were described by a model containing three parameters: transition probabilities for generating CdCl₂-type and CrOOH-type faults and a transition probability for the random intercalation of water/carbonate layers. Automated series of simulations and refinements were performed, in which the transition probabilities were modified incrementally and thus the microstructures optimized by a grid search. All samples were found to exhibit the same fraction of CdCl₂-type and CrOOH-type stacking faults, which indicates that they have identical Ni, Co and Al contents. Different degrees of interstratification faulting were determined, which could be correlated to different heights of intercalation-water-related mass-loss steps in the thermal analyses.

The thriving of older people assessment scale: Psychometric evaluation and short-form development

AIM

To evaluate the psychometric properties and performance of the 32-item Thriving of Older People Assessment Scale (TOPAS) and to explore reduction into a short-form.

BACKGROUND

The 32-item TOPAS has been used in studies of place-related well-being as a positive measure in long-term care to assess nursing home resident thriving; however, item redundancy has not previously been explored.

DESIGN

Cross-sectional study.

METHOD

Staff members completed the 32-item TOPAS as proxy raters for a random sample of Swedish nursing home residents (N = 4,831) between November 2013 - September 2014. Reliability analysis, exploratory factor analysis and item response theory-based analysis were undertaken. Items were systematically identified for reduction using statistical and theoretical analysis. Correlation testing, means comparison and model fit evaluation confirmed scale equivalence.

RESULTS

Psychometric properties of the 32-item TOPAS were satisfactory and several items were identified for scale reduction. The proposed short-form TOPAS exhibited a high level of internal consistency (α = 0.90) and strong correlation (r = 0.98) to the original scale, while also retaining diversity among items in terms of factor structure and item difficulties.

CONCLUSION

The 32-item and short-form TOPAS' indicated sound validity and reliability to measure resident thriving in the nursing home context.

IMPACT

There is a lack of positive life-world measures for use in nursing homes. The short-form TOPAS indicated sound validity and reliability to measure resident thriving, providing a feasible measure with enhanced functionality for use in aged care research, assessments and care planning for health-promoting purposes in nursing homes.

Monte Carlo Processing on a Chip (MCoaC)-preliminary experiments toward the realization of optimal-hardware for TOPAS/Geant4 to drive discovery

Amongst the scientific frameworks powered by the Monte Carlo (MC) toolkit Geant4 (Agostinelli et al., 2003), the TOPAS (Tool for Particle Simulation) (Perl et al., 2012) is one. TOPAS focuses on providing ease of use, and has significant implementation in the radiation oncology space at present. TOPAS functionality extends across the full capacity of Geant4, is freely available to non-profit users, and is being extended into radiobiology via TOPAS-nBIO (Ramos-Mendez et al., 2018). A current "grand problem" in cancer therapy is to convert the dose of treatment from physical dose to biological dose, optimized ultimately to the individual context of administration of treatment. Biology MC calculations are some of the most complex and require significant computational resources. In order to enhance TOPAS's ability to become a critical tool to explore the definition and application of biological dose in radiation therapy, we chose to explore the use of Field Programmable Gate Array (FPGA) chips to speedup the Geant4 calculations at the heart of TOPAS, because this approach called "Reconfigurable Computing" (RC), has proven able to produce significant (around 90x) (Sajish et al., 2012) speed increases in scientific computing. Here, we describe initial steps to port Geant4 and TOPAS to be used on FPGA. We provide performance analysis of the current TOPAS/Geant4 code from an RC implementation perspective. Baseline benchmarks are presented. Achievable performance figures of the subsections of the code on optimal hardware are presented; Aspects of practical implementation of "Monte Carlo on a chip" are also discussed.

TOPAS Monte Carlo simulation for double scattering proton therapy and dosimetric evaluation

PURPOSE

To construct and commission a double scattering (DS) proton beam model in TOPAS Monte Carlo (MC) code. Dose comparisons of MC calculations to the measured and treatment planning system (TPS) calculated dose were performed.

METHODS

The TOPAS nozzle model was based on the manufacturer blueprints. Nozzle set-up and beam current modulations were calculated using room-specific calibration data. This model was implemented to reproduce pristine peaks, spread-out Bragg peaks (SOBP) and lateral profiles. A stair-shaped target plan in water phantom was calculated and compared to measured data to verify range compensator (RC) modeling.

RESULTS

TOPAS calculated pristine peaks agreed well with measurements, with accuracies of 0.03 cm for range R₉₀ and 0.05 cm for distal dose fall-off (DDF). The calculated SOBP range, modulation width and DDF differences between MC calculations and measurements were within 0.05 cm, 0.5 cm and 0.03 cm respectively. MC calculated lateral penumbra agreed well with measured data, with difference less than 0.05 cm. For RC calculation, TPS underestimated the additional depth dose tail due to the nuclear halo effect. Lateral doses by TPS were 10% lower than measurement outside the target, while maximum difference of MC calculation was within 2%. At deeper depths inside the target volume, TPS overestimated doses by up to 25% while TOPAS predicted the dose to within 5% of measurements.

CONCLUSION

We have successfully developed and commissioned a MC based DS nozzle model. The performance of dose accuracy by TOPAS was superior to TPS, especially for highly inhomogeneous compensator.

Helium CT: Monte Carlo simulation results for an ideal source and detector with comparison to proton CT

PURPOSE

To evaluate the accuracy of relative stopping power and spatial resolution of images reconstructed with simulated helium CT (HeCT) in comparison to proton CT (pCT).

METHODS

A Monte Carlo (MC) study with the TOPAS tool was performed to compare the accuracy of relative stopping power (RSP) reconstruction and spatial resolution of low-fluence HeCT to pCT, both using 200 MeV/u particles. An ideal setup consisting of a flat beam source and a totally absorbing energy-range detector was implemented to estimate the theoretically best achievable RSP accuracy for the calibration and reconstruction methods currently used for pCT. The phantoms imaged included a cylindrical water phantom with inserts of different materials, sizes, and positions, a Catphan phantom with a module containing high-contrast line pairs (CTP528) and a module with cylindrical inserts of different RSP (CTP404), as well as a voxelized 10-year-old female phantom. Dose to the cylindrical water phantom was also calculated. The RSP accuracy was studied for all phantoms except the CTP528 module. The latter was used for the estimation of the spatial resolution, evaluated as the modulation transfer function (MTF) at 10%.

RESULTS

An overall error under 0.5% was achieved for HeCT for the water phantoms with the different inserts, in all cases better than that for pCT, in some cases by a factor 3. The inserts in the CTP404 module were reconstructed with an average RSP accuracy of 0.3% for HeCT and 0.2% for pCT. Anatomic structures (brain, bones, air cavities, etc.) in the digitized head phantom were well recognizable and no artifacts were visible with both HeCT and pCT. The three main tissue materials (soft tissue, brain, and cranium) were well identifiable in the reconstructed RSP-volume distribution with both imaging modalities. Using 360 projection angles, the spatial resolution was 4 lp/cm for HeCT and 3 lp/cm for pCT. Generally, spatial resolution increased with the number of projection angles and was always higher for HeCT than for pCT for the same number of projections. When HeCT and pCT scan were performed to deliver the same dose in the phantom, the resolution for HeCT was higher than pCT.

CONCLUSION

MC simulations were used to compare HeCT and pCT image reconstruction. HeCT images had similar or better RSP accuracy and higher spatial resolution compared to pCT. Further investigation of the potential of helium ion imaging is warranted.

Time-resolved diode dosimetry calibration through Monte Carlo modeling for in vivo passive scattered proton therapy range verification

Purpose

Our group previously introduced an in vivo proton range verification methodology in which a silicon diode array system is used to correlate the dose rate profile per range modulation wheel cycle of the detector signal to the water‐equivalent path length (WEPL) for passively scattered proton beam delivery. The implementation of this system requires a set of calibration data to establish a beam‐specific response to WEPL fit for the selected ‘scout’ beam (a 1 cm overshoot of the predicted detector depth with a dose of 4 cGy) in water‐equivalent plastic. This necessitates a separate set of measurements for every ‘scout’ beam that may be appropriate to the clinical case. The current study demonstrates the use of Monte Carlo simulations for calibration of the time‐resolved diode dosimetry technique.

Methods

Measurements for three ‘scout’ beams were compared against simulated detector response with Monte Carlo methods using the Tool for Particle Simulation (TOPAS). The ‘scout’ beams were then applied in the simulation environment to simulated water‐equivalent plastic, a CT of water‐equivalent plastic, and a patient CT data set to assess uncertainty.

Results

Simulated detector response in water‐equivalent plastic was validated against measurements for ‘scout’ spread out Bragg peaks of range 10 cm, 15 cm, and 21 cm (168 MeV, 177 MeV, and 210 MeV) to within 3.4 mm for all beams, and to within 1 mm in the region where the detector is expected to lie.

Conclusion

Feasibility has been shown for performing the calibration of the detector response for three ‘scout’ beams through simulation for the time‐resolved diode dosimetry technique in passive scattered proton delivery.

TOPAS Simulation of the Mevion S250 compact proton therapy unit

As proton therapy becomes increasingly popular, so does the need for Monte Carlo simulation studies involving accurate beam line modeling of proton treatment units. In this study, the 24 beam configurations of the Mevion S250 proton therapy system installed recently at our institution were modeled using the TOolkit for PArticle Simulation (TOPAS) code. Pristine Bragg peak, spread out Bragg peak (SOBP), and lateral beam profile dose distributions were simulated and matched to the measurements taken during commissioning of the unit. Differences in the range for all Percent Depth Dose (PDD) curves between measured and simulated data agreed to within 0.1 cm. For SOBP scans, the SOBP widths all agreed to within 0.3 cm. With regards to lateral beam profile comparisons between the measured and simulated data, the penumbras differed by less than 1 mm and the flatness differed by less than 1% in nearly all cases. This study shows that Monte Carlo simulation studies involving the Mevion S250 proton therapy unit can be a viable tool in commissioning and verification of the proton treatment planning system.

Dosimetry in MARS spectral CT: TOPAS Monte Carlo simulations and ion chamber measurements

Spectral computed tomography (CT) is an up and coming imaging modality which shows great promise in revealing unique diagnostic information. Because this imaging modality is based on X-ray CT, it is of utmost importance to study the radiation dose aspects of its use. This study reports on the implementation and evaluation of a Monte Carlo simulation tool using TOPAS for estimating dose in a pre-clinical spectral CT scanner known as the MARS scanner. Simulated estimates were compared with measurements from an ionization chamber. For a typical MARS scan, TOPAS estimated for a 30 mm diameter cylindrical phantom a CT dose index (CTDI) of 29.7 mGy; CTDI was measured by ion chamber to within 3% of TOPAS estimates. Although further development is required, our investigation of TOPAS for estimating MARS scan dosimetry has shown its potential for further study of spectral scanning protocols and dose to scanned objects.

The effect of beam purity and scanner complexity on proton CT accuracy

PURPOSE

To determine the dependence of the accuracy in reconstruction of relative stopping power (RSP) with proton computerized tomography (pCT) scans on the purity of the proton beam and the technological complexity of the pCT scanner using standard phantoms and a digital representation of a pediatric patient.

METHODS

The Monte Carlo method was applied to simulate the pCT scanner, using both a pure proton beam (uniform 200 MeV mono-energetic, parallel beam) and the Northwestern Medicine Chicago Proton Center (NMCPC) clinical beam in uniform scanning mode. The accuracy of the simulation was validated with measurements performed at NMCPC including reconstructed RSP images obtained with a preclinical prototype pCT scanner. The pCT scanner energy detector was then simulated in three configurations of increasing complexity: an ideal totally absorbing detector, a single stage detector and a multi-stage detector. A set of 15 cm diameter water cylinders containing either water alone or inserts of different material, size, and position were simulated at 90 projection angles (4° steps) for the pure and clinical proton beams and the three pCT configurations. A pCT image of the head of a detailed digital pediatric phantom was also reconstructed from the simulated pCT scan with the prototype detector.

RESULTS

The RSP error increased for all configurations for insert sizes under 7.5 mm in radius, with a sharp increase below 5 mm in radius, attributed to a limit in spatial resolution. The highest accuracy achievable using the current pCT calibration step phantom and reconstruction algorithm, calculated for the ideal case of a pure beam with totally absorbing energy detector, was 1.3% error in RSP for inserts of 5 mm radius or more, 0.7 mm in range for the 2.5 mm radius inserts, or better. When the highest complexity of the scanner geometry was introduced, some artifacts arose in the reconstructed images, particularly in the center of the phantom. Replacing the step phantom used for calibration with a wedge phantom led to RSP accuracy close to the ideal case, with no significant dependence of RSP error on insert location or material. The accuracy with the multi-stage detector and NMCPC beam for the cylindrical phantoms was 2.2% in RSP error for inserts of 5 mm radius or more, 0.7 mm in range for the 2.5 mm radius inserts, or better. The pCT scan of the pediatric phantom resulted in mean RSP values within 1.3% of the reference RSP, with a range error under 1 mm, except in exceptional situations of parallel incidence on a boundary between low and high density.

CONCLUSIONS

The pCT imaging technique proved to be a precise and accurate imaging tool, rivaling the current x-rays based techniques, with the advantage of being directly sensitive to proton stopping power rather than photon interaction coefficients. Measured and simulated pCT images were obtained from a wobbled proton beam for the first time. Since the in-silico results are expected to accurately represent the prototype pCT, upcoming measurements using the wedge phantom for calibration are expected to show similar accuracy in the reconstructed RSP.