Depth In Phantom Research Articles

Phase-modulated multi-foci microscopy for rapid 3D imaging

3D imaging technology is pivotal in monitoring the functional dynamics and morphological alterations in living cells and tissues. However, conventional volumetric imaging associated with mechanical z-scanning encounters challenges in limited 3D imaging speed with inertial artifact. Here, we present a unique phase-modulated multi-foci microscopy (PM3) technique to achieve snapshot 3D imaging with the advantages of extended imaging depths and adjustable imaging intervals between each focus in a rapid fashion. To accomplish the tasks, we utilize a spatial light modulator (SLM) to encode the phases of the scattered or fluorescence light emanating from a volumetric sample and then project the multiple-depth images of the sample onto a single charge-coupled device camera for rapid 3D imaging. We demonstrate that the PM3 technique provides ∼55-fold improvement in imaging depth in polystyrene beads phantom compared to the depth of field of the objective lens used. PM3 also enables the real-time monitoring of Brownian motion of fluorescent beads in water at a 15 Hz volume rate. By precisely manipulating the phase of scattered light on the SLM, PM3 can pinpoint the specific depth information in living zebrafish and rapidly observe the 3D dynamic processes of blood flow in the zebrafish trunk. This work shows that the PM3 technique developed is robust and versatile for fast 3D dynamic imaging in biological and biomedical systems.

Experimental comparison of cylindrical and plane parallel ionization chambers for reference dosimetry in continuous and pulsed scanned proton beams

Objective. In this experimental work we compared the determination of absorbed dose to water using four ionization chambers (ICs), a PTW-34045 Advanced Markus, a PTW-34001 Roos, an IBA-PPC05 and a PTW-30012 Farmer, irradiated under the same conditions in one continuous- and in two pulsed-scanned proton beams. Approach. The ICs were positioned at 2 cm depth in a water phantom in four square-field single-energy scanned-proton beams with nominal energies between 80 and 220 MeV and in the middle of 10 × 10 × 10 cm3 dose cubes centered at 10 cm or 12.5 cm depth in water. The water-equivalent thickness (WET) of the entrance window and the effective point of measurement was considered when positioning the plane parallel (PP) ICs and the cylindrical ICs, respectively. To reduce uncertainties, all ICs were calibrated at the same primary standards laboratory. We used the beam quality (k Q ) correction factors for the ICs under investigation from IAEA TRS-398, the newly calculated Monte Carlo (MC) values and the anticipated IAEA TRS-398 updated recommendations. Main results. Dose differences among the four ICs ranged between 1.5% and 3.7% using both the TRS-398 and the newly recommended k Q values. The spread among the chambers is reduced with the newly k Q values. The largest differences were observed between the rest of the ICs and the IBA-PPC05 IC, obtaining lower dose with the IBA-PPC05. Significance. We provide experimental data comparing different types of chambers in different proton beam qualities. The observed dose differences between the ICs appear to be related to inconsistencies in the determination of the k Q values. For PP ICs, MC studies account for the physical thickness of the entrance window rather than the WET. The additional energy loss that the wall material invokes is not negligible for the IBA-PPC05 and might partially explain the low k Q values determined for this IC. To resolve this inconsistency and to benchmark MC values, k Q values measured using calorimetry are needed.

Minimum phantom size for megavoltage photon beam reference dosimetry.

Water phantoms are required to perform reference dosimetry and beam quality measurements but there are no published studies about the size requirements for suchphantoms. To investigate, using Monte Carlo techniques, the size requirements for water phantoms used in reference dosimetry and/or to measure the beam quality specifiers and . The EGSnrc application DOSXYZnrc is used to calculate , the dose per incident fluence at 10cm depth in a water phantom irradiated by incident beams of or 6MV photons. The water phantom dimensions are varied from to and occasionally smaller. The and values are also calculated with care being taken to distinguish results when using Method A (changing depth of water in phantom) and Method B (moving entire phantom). Typical statistical uncertainties are 0.03%. Phantom dimensions have only minor effects for phantoms larger than . A table of corrections to the dose at 10cm depth in beams of or 6MV photons are provided and range from no correction to 0.75% for a beam incident on a phantom. There can be distinct differences in the values measured using Method A or Method B, especially for smaller phantoms. It is explicitly demonstrated that, within 0.15%, values for a phantom measured using Method A or B are independent of source detector distance between 40 and200cm. The phantom sizes recommended in the TG-51 and IAEA TRS-398 reference dosimetry protocols are adequate for accurate reference dosimetry and in some cases are even conservative. Correction factors are necessary for accurate measurement of the dose at 10cm depth in smaller phantoms and these factors are provided. Very accurate beam quality specifiers are not required for reference dosimetry itself, but for specifying beam stability and characteristics it is important to specify phantom sizes and also the method used for measurements.

Estimation of the relative biological effectiveness for double strand break induction of clinical kilovoltage beams using Monte Carlo simulations.

The Relative Biological Effectiveness (RBE) of kilovoltage photon beams has been previously investigated in vitro and in silico using analytical methods. The estimated values range from 1.03 to 1.82 depending on the methodology and beam energiesexamined. The focus of this work was to independently estimate RBE values for a range of clinically used kilovoltage beams (70-200 kVp) while investigating the suitability of using TOPAS-nBio for thistask. Previously validated spectra of clinical beams were used to generate secondary electron spectra at several depths in a water tank phantom via TOPAS Monte Carlo (MC) simulations. Cell geometry was irradiated with the secondary electrons in TOPAS-nBio MC simulations. The deposited dose and the calculated number of DNA strand breaks were used to estimate RBEvalues. Monoenergetic secondary electron simulations revealed the highest direct and indirect double strand break yield at approximately 20keV. The average RBE value for the kilovoltage beams was calculated to be 1.14. TOPAS-nBio was successfully used to estimate the RBE values for a range of clinical radiotherapy beams. The calculated value was in agreement with previous estimates, providing confidence in its clinical use in thefuture.

Performance evaluation of commercial and non-commercial shear wave elastography implementations for vascular applications

BackgroundShear wave elastography (SWE) is mainly used for stiffness estimation of large, homogeneous tissues, such as the liver and breasts. However, little is known about its accuracy and applicability in thin (∼0.5–2 mm) vessel walls. To identify possible performance differences among vendors, we quantified differences in measured wave velocities obtained by commercial SWE implementations of various vendors over different imaging depths in a vessel-mimicking phantom. For reference, we measured SWE values in the cylindrical inclusions and homogeneous background of a commercial SWE phantom. Additionally, we compared the accuracy between a research implementation and the commercially available clinical SWE on an Aixplorer ultrasound system in phantoms and in vivo in patients. MethodsSWE measurements were performed over varying depths (0–35 mm) using three ultrasound machines with four ultrasound probes in the homogeneous 20 kPa background and cylindrical targets of 10, 40, and 60 kPa of a multi-purpose phantom (CIRS-040GSE) and in the anterior and posterior wall of a homogeneous polyvinyl alcohol vessel-mimicking phantom. These phantom data, along with in vivo SWE data of carotid arteries in 23 patients with a (prior) head and neck neoplasm, were also acquired in the research and clinical mode of the Aixplorer ultrasound machine. Machine-specific estimated phantom stiffness values (CIRS phantom) or wave velocities (vessel phantom) over all depths were visualized, and the relative error to the reference values and inter-frame variability (interquartile range/median) were calculated. Correlations between SWE values and target/vessel wall depth were explored in phantoms and in vivo using Spearman’s correlations. Differences in wave velocities between the anterior and posterior arterial wall were assessed with Wilcoxon signed-rank tests. Intra-class correlation coefficients were calculated for a sample of ten patients as a measure of intra- and interobserver reproducibility of SWE analyses in research and clinical mode. ResultsThere was a high variability in obtained SWE values among ultrasound machines, probes, and, in some cases, with depth. Compared to the homogeneous CIRS-background, this variation was more pronounced for the inclusions and the vessel-mimicking phantom. Furthermore, higher stiffnesses were generally underestimated. In the vessel-mimicking phantom, anterior wave velocities were (incorrectly) higher than posterior wave velocities (3.4–5.6 m/s versus 2.9–5.9 m/s, p ≤ 0.005 for 3/4 probes) and remarkably correlated with measurement depth for most machines (Spearman’s ρ = -0.873–0.969, p < 0.001 for 3/4 probes). In the Aixplorer’s research mode, this difference was smaller (3.3–3.9 m/s versus 3.2–3.6 m/s, p = 0.005) and values did not correlate with measurement depth (Spearman’s ρ = 0.039–0.659, p ≥ 0.002). In vivo, wave velocities were higher in the posterior than the anterior vessel wall in research (left p = 0.001, right p < 0.001) but not in clinical mode (left: p = 0.114, right: p = 0.483). Yet, wave velocities correlated with vessel wall depth in clinical (Spearman’s ρ = 0.574–0.698, p < 0.001) but not in research mode (Spearman’s ρ = -0.080–0.466, p ≥ 0.003). ConclusionsWe observed more variation in SWE values among ultrasound machines and probes in tissue with high stiffness and thin-walled geometry than in low stiffness, homogeneous tissue. Together with a depth-correlation in some machines, where carotid arteries have a fixed location, this calls for caution in interpreting SWE results in clinical practice for vascular applications.

Visualization of photothermal therapy by semiconducting polymer dots mediated photoacoustic detection in NIR II

Visualization of photothermal therapy mediated by photothermal transduction agents (PTAs) is important to promote individual treatment of patients with low side effects. Photoacoustic detection has emerged as a promising noninvasive method for the visualization of PTAs distribution but still has limitations in temperature measurement, including poor measurement accuracy and low tissue penetration depth. In this study, we developed biocompatible semiconducting polymer dots (SPD) for in situ coupling of photothermal and photoacoustic detection in the near-infrared II window. SPD has dual photostability under pulsed laser and continuous-wave laser irradiation with a photothermal conversion efficiency of 42.77%. Meanwhile, a strong correlation between the photoacoustic signal and the actual temperature of SPD can be observed. The standard deviation of SPD-mediated photoacoustic thermometry can reach 0.13 °C when the penetration depth of gelatin phantom is 9.49 mm. Preliminary experimental results in vivo show that SPD-mediated photoacoustic signal has a high signal-to-noise ratio, as well as good performance in temperature response and tumor enrichment. Such a study not only offers a new nanomaterial for the visualization of photothermal therapy but will also promote the theranostic platform for clinical applications.

Application of the ESHO-QA guidelines for determining the performance of the LCA superficial hyperthermia heating system

Purpose This study aimed to assess the quality of the lucite cone applicator (LCA), the standard applicator for superficial hyperthermia at the Erasmus MC Cancer Institute, using the most recent quality assurance guidelines, thus verifying their feasibility. Materials and methods The assessment was conducted on each of the six LCAs available for clinical treatments. The temperature distribution was evaluated using an infrared camera across different layers of a fat-muscle mimicking phantom. The maximum temperature increase, thermal effective penetration depth (TEPD), and thermal effective field size (TEFS) were used as quality metrics. The experimental results were validated through comparison with simulated results, using a canonical phantom model and a realistic phantom model segmented from CT imaging. Results A maximum temperature increase above 6 °C at 2 cm depth in the fat-muscle phantom for all the experiments was found. A mean negative difference between simulated and experimental data was of 1.3 °C when using the canonical phantom model. This value decreased to a mean negative difference of 0.4 °C when using the realistic model. Simulated and measured TEPD showed good agreement for both in silico scenarios, while discrepancies were present for TEFS. Conclusions The LCAs passed all QA guidelines requirements for superficial hyperthermia delivery when used singularly or in an array configuration. A further characterization of parameters such as antenna efficiency and heat transfer coefficients would be beneficial for translating experimental results to simulated values. Implementing the QA guidelines was time-consuming and demanding, requiring careful preparation and correct setup of antenna elements.

Monte Carlo simulations of cell survival in proton SOBP.

Objective. The objective of this study is to develop a multi-scale modeling approach that accurately predicts radiation-induced DNA damage and survival fraction in specific cell lines.Approach. A Monte Carlo based simulation framework was employed to make the predictions. The FLUKA Monte Carlo code was utilized to estimate absorbed doses and fluence energy spectra, which were then used in the Monte Carlo Damage Simulation code to compute DNA damage yields in Chinese hamster V79 cell lines. The outputs were converted into cell survival fractions using a previously published theoretical model. To reduce the uncertainties of the predictions, new values for the parameters of the theoretical model were computed, expanding the database of experimental points considered in the previous estimation. Simulated results were validated against experimental data, confirming the applicability of the framework for proton beams up to 230 MeV. Additionally, the impact of secondary particles on cell survival was estimated.Main results. The simulated survival fraction versus depth in a glycerol phantom is reported for eighteen different configurations. Two proton spread out Bragg peaks at several doses were simulated and compared with experimental data. In all cases, the simulations follow the experimental trends, demonstrating the accuracy of the predictions up to 230 MeV.Significance. This study holds significant importance as it contributes to the advancement of models for predicting biological responses to radiation, ultimately contributing to more effective cancer treatment in proton therapy.

Uncertainty-aware spot rejection rate as quality metric for proton therapy using a digital tracking calorimeter

Objective. Proton therapy is highly sensitive to range uncertainties due to the nature of the dose deposition of charged particles. To ensure treatment quality, range verification methods can be used to verify that the individual spots in a pencil beam scanning treatment fraction match the treatment plan. This study introduces a novel metric for proton therapy quality control based on uncertainties in range verification of individual spots. Approach. We employ uncertainty-aware deep neural networks to predict the Bragg peak depth in an anthropomorphic phantom based on secondary charged particle detection in a silicon pixel telescope designed for proton computed tomography. The subsequently predicted Bragg peak positions, along with their uncertainties, are compared to the treatment plan, rejecting spots which are predicted to be outside the 95% confidence interval. The such-produced spot rejection rate presents a metric for the quality of the treatment fraction. Main results. The introduced spot rejection rate metric is shown to be well-defined for range predictors with well-calibrated uncertainties. Using this method, treatment errors in the form of lateral shifts can be detected down to 1 mm after around 1400 treated spots with spot intensities of 1 × 107 protons. The range verification model used in this metric predicts the Bragg peak depth to a mean absolute error of 1.107 ± 0.015 mm. Significance. Uncertainty-aware machine learning has potential applications in proton therapy quality control. This work presents the foundation for future developments in this area.

Improvement of neutron sensitivity for lithium formate EPR dosemeters: a Monte Carlo analysis.

This work presents the computational analysis of the sensitivity improvements that could be achieved in lithium formate monohydrate (LFM) electron paramagnetic resonance (EPR) dosemeters exposed to neutron beams. Monte Carlo (MC) simulations were performed on LFM pellets exposed to neutron beams with different energy spectra at various depths inside a water phantom. Various computations were carried out by considering different enrichments of 6Li inside the LFM matrix as well as addition of different amounts of gadolinium oxide inside the pellet blend. The energy released per unit mass was calculated with the aim of predicting the increase in dose achievable by the addition of sensitizers inside the pellets. As expected, a larger amount of 6Li induces an increase of energy released because of the charged secondary particles (i.e. 3H ions and α-particles) produced after neutron capture. For small depths in water phantom and low-energy neutron spectra the dose increase due to 6Li enrichment is high (more than three orders of magnitude with respect to the case of with 7Li). In case of epithermal neutron beams the energy released in 6Li-enriched LFM compound is smaller but larger than in the case of fast neutron beams. On the other hand, the computational analysis evidenced that gadolinium is less effective than 6Li in improving neutron sensitivity of the LFM pellets. Discussion based on the features of MC transport code is provided. This result suggests that 6Li enrichment of LFM dosemeters would be more effective for neutron sensitivity improvement and these EPR dosemeters could be tested for dosimetric applications in Neutron Capture Therapy.

Diamond based integrated detection system for dosimetric and microdosimetric characterization of radiotherapy ion beams.

Ion beam therapy allows for a substantial sparing of normal tissues and higher biological efficacy. Synthetic single crystal diamond is a very good material to produce high-spatial-resolution and highly radiation hard detectors for both dosimetry and microdosimetry in ion beam therapy. The aim of this work is the design, fabrication and test of an integrated waterproof detector based on synthetic single crystal diamond able to simultaneously perform dosimetric and microdosimetric characterization of clinical ion beams. The active elements of the integrated diamond device, that is, dosimeter and microdosimeter, were both realized in a Schottky diode configuration featured by different area, thickness, and shape by means of photolithography technologies for the selective growth of intrinsic and boron-doped CVD diamond. The cross-section of the sensitive volume of the dosimetric element is 4mm2 and 1μm-thick, while the microdosimetric one has an active cross-sectional area of 100×100μm2 and a thickness of about 6.2μm. The dosimetric and microdosimetric performance of the developed device was assessed at different depths in a water phantom at the MedAustron ion beam therapy facility using a monoenergetic uniformly scanned carbon ion beam of 284.7MeV/u and proton beam of 148.7MeV. The particle flux in the region of the microdosimeter was 6·107 cm2 /s for both irradiation fields. At each depth, dose and dose distributions in lineal energy were measured simultaneously and the dose mean lineal energy values were then calculated. Monte Carlo simulations were also carried out by using the GATE-Geant4 code to evaluate the relative dose, dose averaged linear energy transfer (LETd ), and microdosimetric spectra at various depths in water for the radiation fields used, by considering the contribution from the secondary particles generated in the ion interaction processes as well. Dosimetric and microdosimetric quantities were measured by the developed prototype with relatively low noise (∼2keV/μm). A good agreement between the measured and simulated dose profiles was found, with discrepancies in the peak to plateau ratio of about 3% and 4% for proton and carbon ion beams respectively, showing a negligible LET dependence of the dosimetric element of the device. The microdosimetric spectra were validated with Monte Carlo simulations and a good agreement between the spectra shapes and positions was found. Dose mean lineal energy values were found to be in close agreement with those reported in the literature for clinical ion beams, showing a sharp increase along the Bragg curve, being also consistent with the calculated LETd for all depths within the experimental error of 10%. The experimental indicate that the proposed device can allow enhanced dosimetry in particle therapy centers, where the absorbed dose measurement is implemented by the microdosimetric characterization of the radiation field, thus providing complementary results. In addition, the proposed device allows for the reduction of the experimental uncertainties associated with detector positioning and could facilitate the partial overcoming of some drawbacks related to the low sensitivity of diamond microdosimeters to low LET radiation.

New operational quantities (RBE × D p lens/Φ) for eye lens neutron dosimetry as a function of energy and angle of incidence in the ICRU 95 formalism

This paper is a continuation of a study published recently by the authors. It presents and discusses computed personal absorbed dose in the lens of the eye (D p lens/Φ), and a relative biological effectiveness (RBE)-weighted absorbed dose (in terms of an newly proposed operational quantity RBE × D p lens/Φ), conversion coefficients for the lens of the eye for neutron exposure at incident energies from thermal to ∼20 MeV and at angles of incidence from 0° to 90° in 15° increments, at 180° and for rotational incidence irradiation geometry (from 0° to 360° in 5° increments). These conversion coefficients were obtained from a simulation model developed for this study that contains the stylised eye model, embedded in the adult UF-ORNL mathematical phantom, whereby the previously stated RBE-weighted absorbed dose was obtained using the proposed RBE versus neutron energy distribution compiled in a previous paper by the authors. The simulations carried out for this study using the Monte Carlo N-Particle transport code version 6.2, were conducted in a realistic human eye model, for the left and right sensitive and whole volume of the lens of the eye, considering the recent proposed redefinition of the operational quantities for external radiation exposure in International Commission on Radiation Units and Measurements (ICRU) report 95. A comprehensive set of tabulated data for neutron fluence-to-dose conversion coefficients (D p lens/Φ in pGy cm2) and RBE-weighted absorbed dose (RBE × D p lens/Φ in pGy cm2) conversion coefficients is included in this paper as a function of incident neutron energy and angle of incidence. Data for D p lens/Φ (pGy cm2) are compared to similar data from the literature for validation of our model. Data for RBE × D p lens/Φ (in pGy cm2), were also compared to the equivalent operational quantity H p(3,α)/Φ (in pSv cm2) conversion coefficients calculated at 3 mm depth in a cylindrical phantom for different incident neutron energies and angles of incidence from 0° to 75° in 15° increments to demonstrate the relevance of this newly proposed operational quantity for doses resulting in tissue reactions (deterministic effects) which should be quoted in Gray (RBE-weighted absorbed dose, RBE × D (Gy)), rather than Sievert (Sv) which is reserved for stochastic effects. The current neutron weighted absorbed dose (RBE × D p lens) is proposed for the tissue reactions in the eye-lens for neutron radiation as per National Council on Radiation Protection and Measurements report 180 and in line with the recent proposal for the review and revision of the System of Radiological Protection to Keeping the International Commission on Radiological Protection (ICRP) recommendations fit for purpose. This method would bring better alignment between the dose limits in ICRP 118 and the new operational quantity consistent with the units of the new eye-lens dose limits without being overly conservative. The utilization of the proposed new operational quantities, as outlined in ICRU 95, has the potential to address the ongoing challenge in enforcing regulatory limits for neutron eye dose, specifically the use of Gy instead of Sv. It should be noted that the applicability of this will vary from country to country as in many countries the legislation is likely to mandate the use of H p(3) until the regulation is amended. This approach can serve as an interim solution while awaiting the issuance of the new ICRP general recommendations, which is expected to take several years. Implementing the new operational quantities can contribute to enhancing the accuracy and effectiveness of neutron eye dose limit enforcement.

Assembly of fluorophore J-aggregates with nanospacer onto mesoporous nanoparticles for enhanced photoacoustic imaging

Many fluorophores, such as indocyanine green (ICG), have poor photostability and low photothermal efficiency hindering their wide application in photoacoustic (PA) tomography. In the present study, a supramolecular assembly approach was used to develop the hybrid nanoparticles (Hy NPs) of ICG and porous silicon (PSi) as a novel contrast agent for PA tomography. ICG was assembled on the PSi NPs to form J-aggregates within 30 min. The Hy NPs presented a red-shifted absorption, improved photothermal stability, and enhanced PA performance. Furthermore, 1-dodecene (DOC) was assembled into the NPs as a ‘nanospacer’, which enhanced non-radiative decay for increased thermal release. Compared to the Hy NPs, adding DOC into the Hy NPs (DOC-Hy) increased the PA signal by 83%. Finally, the DOC-Hy was detectable in PA tomography at 1.5 cm depth in tissue phantom even though its concentration was as low as 6.25 µg/mL, indicating the potential for deep tissue PA imaging.

Evaluation of proton and carbon ion beam models in TReatment Planning for Particles 4D (TRiP4D) referring to a commercial treatment planning system

PurposeTo investigate the accuracy of the treatment planning system (TPS) TRiP4D in reproducing doses computed by the clinically used TPS SyngoRT. MethodsProton and carbon ion beam models in TRiP4D were converted from SyngoRT. Cubic plans with different depths in a water-tank phantom (WP) and previously treated and experimentally verified patient plans from SyngoRT were recalculated in TRiP4D. The target mean dose deviation (ΔDmean,T) and global gamma index (2%–2 mm for the absorbed dose and 3%–3mm for the RBE-weighted dose with 10% threshold) were evaluated. ResultsThe carbon and proton absorbed dose gamma passing rates (γ-PRs) were ≥99.93% and ΔDmean,T smaller than −0.22%. On average, the RBE-weighted dose Dmean,T was −1.26% lower for TRiP4D than SyngoRT for cubic plans. In TRiP4D, the faster analytical ‘low dose approximation’ (Krämer, 2006) was used, while SyngoRT used a stochastic implementation (Krämer, 2000). The average ΔDmean, T could be reduced to −0.59% when applying the same biological effect calculation algorithm. However, the dose recalculation time increased by a factor of 79–477. ΔDmean,T variation up to −2.27% and −2.79% was observed for carbon absorbed and RBE-weighted doses in patient plans. The γ-PRs were ≥93.92% and ≥91.83% for patient plans, except for one proton beam with a range shifter (γ-PR of 64.19%). ConclusionThe absorbed dose between TRiP4D and SyngoRT were identical for both proton and carbon ion plans in the WP. Compared to SyngoRT, TRiP4D underestimated the target RBE-weighted dose; however more efficient in RBE-weighted dose calculation. Large variation for proton beam with range shifter was observed. TRiP4D will be used to evaluate doses delivered to moving targets. Uncertainties inherent to the 4D-dose reconstruction calculation are expected to be significantly larger than the dose errors reported here. For this reason, the residual differences between TRiP4D and SyngoRT observed in this study are considered acceptable.The study was approved by the Institutional Research Board of Shanghai Proton and Heavy Ion Center (approval number SPHIC-MP-2020-04, RS).

Bragg curve, dose distribution, and target fragmentation for thyroid proton therapy

Information on energy deposition in human tissues and organs is of critical significance for proton radiation treatment planning. In the present paper, we studied the Bragg curves, dose distributions, and target fragmentations of proton therapy for a newly developed thyroid phantom. The phantom is considered to consist of 1.2 mm skin, 4.0 mm adipose, 6.0 mm skeletal muscle, 14.0 mm thyroid gland and 2.0 mm skeletal muscle. The Bragg curves were determined and compared by using three different simulations platforms namely Geant4, FLUKA, and PHITS at several proton energies. The obtained results showed that around 30 MeV, no energy deposition is observed at thyroid phantom and then the energy deposition rates of the layers at deeper (shallower) sides increases (decreases) as the proton energy increases. Moreover, we studied the production cross sections, due to the inelastic interactions caused by proton beam in the thyroid phantom, to quantify the contribution of target fragmentation. The fluence of each fragments at various depths in a thyroid phantom is also investigated and accurately scored. As a result, it is concluded that the target fragmentation is more important before the Bragg peak.

Development of a compact linear accelerator to generate ultrahigh dose rate high-energy X-rays for FLASH radiotherapy applications.

In recent years, the FLASH effect, in which ultrahigh dose rate (UHDR) radiotherapy (RT) can significantly reduce toxicity to normal tissue while maintaining antitumor efficacy, has been verified in many studies and even applied in human clinical cases. This work evaluates whether a room-temperature radio-frequency (RF) linear accelerator (linac) system can produce UHDR high-energy X-rays exceeding a dose rate of 40Gy/s at a clinical source-surface distance (SSD), exploring the possibility of a compact and economical clinical FLASH RT machine suitable for most hospital treatmentrooms. A 1.65m long S-band backward-traveling-wave (BTW) electron linac was developed to generate high-current electron beams, supplied by a commercial klystron-based power source. A tungsten-copper electron-to-photon conversion target for UHDR X-rays was designed and optimized with Monte Carlo (MC) simulations using Geant4 and thermal finite element analysis (FEA) simulations using ANSYS. EBT3 and EBT-XD radiochromic films, which were calibrated with a clinical machine Varian VitalBeam, were used for absolute dose measurements. A PTW ionization chamber detector was used to measure the relative total dose and a plane-parallel ionization chamber detector was used to measure the relative normalized dose of each pulse. The BTW linac generated 300-mA-pulse-current 11 MeV electron beams with 29kW mean beam power, and the conversion target could sustain this high beam power within a maximum irradiation duration of 0.75 s. The mean energy of the produced X-rays was 1.66 MeV in the MC simulation. The measured flat-filter-free (FFF) maximum mean dose rate of the room-temperature linac exceeded 80Gy/s at an SSD of 50cm and 45Gy/s at an SSD of 67.9cm, both at a 2.1cm depth of the water phantom. The FFF radiation fields at 50cm and 67.9cm SSD at a 2.1cm depth of the water phantom showed Gaussian-like distributions with 14.3 and 20cm full-width at half-maximum (FWHM) values, respectively. This work demonstrated the feasibility of UHDR X-rays produced by a room-temperature RF linac, and explored the further optimization of system stability. It shows that a simple and compact UHDR X-ray solution can be facilitated for both FLASH-RT scientific research and clinical applications.

Spectral analysis of Monte Carlo calculated fluence correction and cema conversion factors for high-energy photon beams at different depths

Purpose: The aim of this study is to investigate the depth-dependent detector response of detailed thimble air-filled ionization chambers by calculating spectral charged particle fluence correction factors at different depths in water. Those spectral correction factors will help to understand, how the detector response varies at different depths and what kind of influences disparate effects have on the spectral detector response.Methods: The cema-approach can be used to obtain spectral charged particle fluence-based correction factors for various measurement conditions by substituting the commonly well-known dose conversion factor with a conversion factor based on the dosimetric quantity cema (“converted-energy per unit mass”). The resulting spectral fluence correction factors were calculated with the EGSnrc software toolkit and analyzed for two air-filled cylindrical ionization chambers (PTW type 31021 Semiflex 3D, SNC125c™) at different depths in a water phantom irradiated with a 6 MV linear accelerator x-ray spectrum. The ionization chamber models have been stepwise decomposed to investigate the perturbation caused by internal and external effects on the fluence distribution within the detector.Results: Monte Carlo calculated fluence-based perturbation correction factors revealed that for all investigated detectors, considerable fluence disturbances occur, especially in the build-up region of depth-dose curves. Our results have shown that even slight variations in depth can have major consequences on the differential charged particle fluence within the ionization chamber, mainly due to internal cavity-specific effects. Furthermore, the results showed that in the case of relative dose measurements, the depth-depending detector response can significantly differ from unity in a range of 1.4%–2.8% depending on the ionization chamber design.Conclusion: The complexity of different effects on the fluence disturbance could be broken down with regard to their influence on the spectral fluence distribution in the sensitive volume of the investigated detectors. It could be demonstrated, that the displacement of water is a depth-depending effect, which can not be compensated or corrected ideally for each investigated water depth by the shift of the effective point of measurement. Generally, the spectral analysis of those energy-dependent correction factors serves to a deeper understanding of the detector response under various conditions.

Technical note: In silico benchmarking of the linear energy transfer-based functionalities for carbon ion beams in a commercial treatment planning system.

The increasing number of studies dealing with linear energy transfer (LET)-based evaluation and optimization in the field of carbon ion radiotherapy (CIRT) indicates the rising demand for LET implementation in commercial treatment planning systems (TPS). Benchmarking studies could play a key role in detecting (and thus preventing) computation errors prior implementing such functionalities in aTPS. This in silico study was conducted to benchmark the following two LET-related functionalities in a commercial TPS against Monte Carlo simulations: (1) dose averaged LET (LETd ) scoring and (2) physical dose filtration based on LET for future LET-based treatment plan evaluation and optimization studies. The LETd scoring and LET-based dose filtering (in which the deposited dose can be separated into the dose below and above the user specified LET threshold) functionalities for carbon ions in the research version RayStation (RS) 9A-IonPG TPS (RaySearch Laboratories, Sweden) were benchmarked against GATE/Geant4 simulations. Pristine Bragg peaks (BPs) and cuboid targets, positioned at different depths in a homogeneous water phantom and a setup with heterogeneity were used for thisstudy. For all setups (homogeneous and heterogeneous), the mean absolute (and relative) LETd difference was less than 1 keV/ m (3.5%) in the plateau and target and less than 2 keV/ m (8.3%) in the fragmentation tail. The maximum local differences were 4 and 6 keV/ m, respectively. The mean absolute (and relative) physical dose differences for both low- and high-LET doses were less than 1 cGy (1.5%) in the plateau, target and tail with a maximum absolute difference of 2cGy. No computation error was found in the tested functionalities except for LETd in lateral direction outside the target, showing the limitation of the implemented monochrome model in the tested TPSversion.

Variability of Transrectal Shear Wave Elastography in a Phantom Model.

This study aimed to assess the variability of transrectal shear wave elastography (SWE) using a designed phantom. In a phantom, the SWE values were examined by two radiologists using agarose and emulsion silicone of different sizes (1, 2, and 3 cm) and shapes (round, cubic) at three depths (1, 2, and 3 cm), two region of interest (ROI) and locations (central, peripheral) using two ultrasound machines (A, B from different vendors). Variability was evaluated using the coefficient of variation (CV). The CVs decreased with increasing phantom size. Significant changes in SWE values included; agarose phantom at 3 cm depth (p < 0.001; machine A), 1 cm depth (p = 0.01; machine B), emulsion silicone at 2 cm depth (p = 0.047, p = 0.020; both machines). The CVs increased with increasing depth. Significant changes in SWE values included; 1 cm agarose (p = 0.037, p = 0.021; both machines) and 2 cm agarose phantom (p = 0.047; machine A). Significant differences in SWE values were observed between the shapes for emulsion silicone phantom (p = 0.032; machines A) and between ROI locations on machine B (p ≤ 0.001). The SWE values differed significantly between the two machines (p < 0.05). The intra-/inter-operator agreements were excellent (intraclass correlation coefficient > 0.9). The phantom size, depth, and different machines affected the variability of transrectal SWE.

Dosimetric Properties of Brass Mesh Bolus for High-energy Photon Beam

To investigate fundamental dosimetric properties of surface dose, exit dose, and beam profile of the brass mesh bolus for 4, 6, and 10 MV high-energy photon beams in radiation therapy. Surface dose and exit dose in the water-equivalent phantom were measured, and percent depth doses (PDDs) were calculated with no bolus, one layer of brass mesh, two layers of brass mesh bolus, three layers of brass mesh bolus, and 0.5 cm tissue-equivalent (TE) bolus. Exit dose was measured at a phantom thickness of 10 cm. Beam profiles were measured at phantom depths of 0 cm and 10 cm. All dosimetry was performed for 4, 6, and 10 MV photon beams using a linear accelerator. The surface dose at a phantom depth of 0 cm increased to 37.3%, 36.3%, and 31.0% for 4, 6, and 10 MV, respectively, with the brass mesh bolus compared to the case of no bolus. The surface dose decreased with one layer of brass mesh bolus compared to that with the 0.5 cm TE bolus. On the other hand, the exit dose increased to 22.0%, 23.1%, and 22.8% for 4, 6, and 10 MV, respectively, with the brass mesh bolus compared to the case of no bolus. The beam profile at the depth of 0 cm showed oscillations, and the difference between the maximum and minimum doses was up to 13.1% with one layer of brass mesh bolus. It was suggested that the brass mesh bolus not only increases the surface dose but also has different properties from the conventional TE bolus.