The use of dose surface maps as a tool to investigate spatial dose delivery accuracy for the rectum during prostate radiotherapy.

PURPOSE: This study aims to address the lack of spatial dose comparisons of planned and delivered rectal doses during prostate radiotherapy by using dose-surface maps (DSMs) to analyze dose delivery accuracy and comparing these results to those derived using DVHs. METHODS: Two independent cohorts were used in this study: twenty patients treated with 36.25 Gy in five fractions (SBRT) and 20 treated with 60 Gy in 20 fractions (IMRT). Daily delivered rectum doses for each patient were retrospectively calculated using daily CBCT images. For each cohort, planned and average-delivered DVHs were generated and compared, as were planned and accumulated DSMs. Permutation testing was used to identify DVH metrics and DSM regions where significant dose differences occurred. Changes in rectal volume and position between planning and delivery were also evaluated to determine possible correlation to dosimetric changes. RESULTS: For both cohorts, DVHs and DSMs reported conflicting findings on how planned and delivered rectum doses differed from each other. DVH analysis determined average-delivered DVHs were on average 7.1% ± 7.6% (p ≤ 0.002) and 5.0 ± 7.4% (p ≤ 0.021) higher than planned for the IMRT and SBRT cohorts, respectively. Meanwhile, DSM analysis found average delivered posterior rectal wall dose was 3.8 ± 0.6 Gy (p = 0.014) lower than planned in the IMRT cohort and no significant dose differences in the SBRT cohort. Observed dose differences were moderately correlated with anterior-posterior rectal wall motion, as well as PTV superior-inferior motion in the IMRT cohort. Evidence of both these relationships were discernable in DSMs. CONCLUSION: DSMs enabled spatial investigations of planned and delivered doses can uncover associations with interfraction motion that are otherwise masked in DVHs. Investigations of dose delivery accuracy in radiotherapy may benefit from using DSMs over DVHs for certain organs such as the rectum.

Mathematical complexities in radionuclide metabolic modelling: a review of ordinary differential equation kinetics solvers in biokinetic modelling.

Biokinetic models have been employed in internal dosimetry (ID) to model the human body's time-dependent retention and excretion of radionuclides. Consequently, biokinetic models have become instrumental in modelling the body burden from biological processes from internalized radionuclides for prospective and retrospective dose assessment. Solutions to biokinetic equations have been modelled as a system of coupled ordinary differential equations (ODEs) representing the time-dependent distribution of materials deposited within the body. In parallel, several mathematical algorithms were developed for solving general kinetic problems, upon which biokinetic solution tools were constructed. This paper provides a comprehensive review of mathematical solving methods adopted by some known internal dose computer codes for modelling the distribution and dosimetry for internal emitters, highlighting the mathematical frameworks, capabilities, and limitations. Further discussion details the mathematical underpinnings of biokinetic solutions in a unique approach paralleling advancements in ID. The capabilities of available mathematical solvers in computational systems were also emphasized. A survey of ODE forms, methods, and solvers was conducted to highlight capabilities for advancing the utilization of modern toolkits in ID. This review is the first of its kind in framing the development of biokinetic solving methods as the juxtaposition of mathematical solving schemes and computational capabilities, highlighting the evolution in biokinetic solving for radiation dose assessment.

Assessment of intrafractional prostate motion and its dosimetric impact in MRI-guided online adaptive radiotherapy with gating.

PURPOSE: This study aimed to evaluate the intrafractional prostate motion captured during gated magnetic resonance imaging (MRI)-guided online adaptive radiotherapy for prostate cancer and analyze its impact on the delivered dose as well as the effect of gating. METHODS: Sagittal 2D cine-MRI scans were acquired at 4 Hz during treatment at a ViewRay MRIdian (ViewRay Inc., Oakwood Village, OH, USA) MR linac. Prostate shifts in anterior-posterior (AP) and superior-inferior (SI) directions were extracted separately. Using the static dose cloud approximation, the planned fractional dose was shifted according to the 2D gated motion (residual motion in gating window) to estimate the delivered dose by superimposing and averaging the shifted dose volumes. The dose of a hypothetical non-gated delivery was reconstructed similarly using the non-gated motion. For the clinical target volume (CTV), rectum, and bladder, dose-volume histogram parameters of the planned and reconstructed doses were compared. RESULTS: In total, 174 fractions (15.7 h of cine-MRI) from 10 patients were evaluated. The average (±1 σ) non-gated prostate motion was 0.6 ± 1.0 mm in the AP and 0.0 ± 0.6 mm in the SI direction with respect to the centroid position of the gating boundary. 95% of the shifts were within [-3.5, 2.7] mm in the AP and [-2.9, 3.2] mm in the SI direction. For the gated treatment and averaged over all fractions, CTV D98% decreased by less than 2% for all patients. The rectum and the bladder D2% increased by less than 3% and 0.5%, respectively. Doses reconstructed for gated and non-gated delivery were similar for most fractions. CONCLUSION: A pipeline for extraction of prostate motion during gated MRI-guided radiotherapy based on 2D cine-MRI was implemented. The 2D motion data enabled an approximate estimation of the delivered dose. For the majority of fractions, the benefit of gating was negligible, and clinical dosimetric constraints were met, indicating safety of the currently adopted gated MRI-guided treatment workflow.

Development of a log file analysis tool for proton patient QA, system performance tracking, and delivered dose reconstruction.

PURPOSE/OBJECTIVE(S): To describe a log file-based patient-specific quality assurance (QA) method and develop an in-house tool for system performance tracking and dose reconstruction in pencil-beam scanning proton therapy that can be used for pre-treatment plan review. MATERIALS/METHODS: The software extracts beam-specific information from the treatment delivery log file and automatically compares the monitor units (MU), lateral position, and size of each spot against the intended values in the treatment plan to identify any discrepancies in the beam delivery. The software has been used to analyze 992 patients, 2004 plans, 4865 fields, and more than 32 million proton spots from 2016 to 2021. The composite doses of 10 craniospinal irradiation (CSI) plans were reconstructed based on the delivered spots and compared with the original plans as an offline plan review method. RESULTS: Over the course of 6 years, the proton delivery system has proved stable in delivering patient QA fields with proton energies of 69.4-221.3 MeV and an MU range of 0.003-1.473 MU per spot. The planned mean and standard deviation (SD) of the energy and spot MU were 114.4 ± 26.4 MeV and 0.010 ± 0.009 MU, respectively. The mean and SD of the differences in MU and position between the delivered and planned spots were 9.56 × 10-8 ± 2.0 × 10-4 MU and 0.029/-0.007 ± 0.049/0.044 mm on the X/Y-axis for random differences and 0.005/0.125 ± 0.189/0.175 mm on the X/Y-axis for systematic differences. The mean and SD of the difference between the commissioning and delivered spot sizes were 0.086/0.089 ± 0.131/0.166 mm on the X/Y-axis. CONCLUSION: A tool has been developed to extract crucial information about the performance of the proton delivery and monitor system and provide a dose reconstruction based on delivered spots for quality improvement. Each patient's plan was verified before treatment to ensure accurate and safe delivery within the delivery tolerance of the machine.

Reconstructed lung doses for the million person study cohort of 26,650 Tennessee Eastman corporation workers employed between 1942 and 1947.

Tennessee Eastman Corporation workers were exposed to uranium dust resulting in high-linear energy transfer (LET) irradiation to lung tissue. In this work, radiation lung doses were reconstructed for 26 650 men and women working at the plant between 1942 and 1947. Site air monitoring data of uranium concentrations and payroll records were used to determine the daily inhaled activities and annualized lung doses. Variations in the activity median aerodynamic diameter of the uranium dust, the solubility of particulate matter in the lungs and the sex-specific breathing rate were investigated as part of a sensitivity analysis. Male and female mean lung doses of 18.9 and 32.7 mGy, respectively, from high-LET alpha irradiation, and there was general agreement with evaluations from previously published epidemiological studies. Annual lung dose estimates and sensitivity analysis for the 26 650 workers in the TEC cohort have been archived on the United States Department of Energy Comprehensive Epidemiologic Data Resource.

[Progress in Development of Dose Verification System Software KylinRay-Dose4D].

Advanced radiotherapy technology enables the dose to more accurately conform to the tumor target area of the patient, providing accurate treatment for the patient, but the gradient of the patient's radiation dose at the tumor edge is getting larger, which putting forward higher requirements for radiotherapy dose verification. The dose verification system software KylinRay-Dose4D can verify the patient's pre-treatment plan and the in vivo/on-line dose during the patient's treatment, providing important reference for the physicist to modify the radiotherapy plan and ensuring that the patient receives accurate treatment. This study introduces the overall design and key technologies of KylinRay-Dose4D, and tests the pre-treatment plan dose checking calculation and 2D/3D dose verification through clinical cases. The test results showed that the 2D/3D gamma pass rate (3 mm/3%) of KylinRay-Dose4D reconstructed dose compared with TPS plan dose and measured dose is larger than 95%, which indicating that the reconstructed dose of KylinRay-Dose4D meets the requirement of clinical application.

External dose reconstruction at the shore of the Metlinsky Pond in the former village of Metlino (Techa River, Russia) based on environmental surveys, luminescence measurements and radiation transport modelling.

The cohorts of people formerly living at the Techa River shoreline in the Southern Urals, Russia, are widely studied cohorts for the investigation of low-dose radiation effects to human health. The nuclear facilities of the Mayak Production Association (PA) discharged their radioactive effluents into the nearby Techa River, especially in the first years of operation. Health status of cohort member data is constantly being improved and updated. Consequently, there is a need to also improve and verify the underlying dosimetry, which gives information about the dose of cohort members. For the Techa River population, the dosimetry is handled in the Techa River Dosimetry System (TRDS). The present work shows results of a feasibility study to validate the TRDS at the location of the village of Metlino, a village just 7 km downstream from the Mayak PA. For this settlement there were two sources of external exposure, the contaminated banks of the Techa River and the contaminated shoreline of the nearby Metlinsky Pond. In the present study the north-western wall of a granary was used as a dose archive to validate dose estimates. Measurements of doses in brick accumulated over many decades and measurements of the current dose rate in bricks were combined with dose rate measurements in air above ground in front of the granary, historical contamination data and Monte-Carlo simulations. Air kerma estimates for 1949-1956 significantly different from zero could not be reconstructed for the Metlinsky Pond shoreline near the granary, but an upper dose limit could be estimated. Implications for TRDS-2016 are discussed.

Estimation of the historical radiation dose of Korean radiation workers, 1961-1983.

Accurate dosimetry plays a key role in evaluating the radiation-induced health risks of radiation workers. The National Dose Registry in Korea contains the dose records of radiation workers in nuclear-related occupations since 1984. Thus, radiation doses for workers before 1984 are often sparse or missing. This study aimed to estimate the historical radiation dose before 1984 for radiation workers in Korea based on dose reconstruction models. The dose reconstruction models were derived from the nationwide self-administered questionnaire survey and the personal badge doses for workers in the cohort of the Korean Radiation Worker Study. The mean estimated annual doses between 1984 and 2016 from the dose reconstruction model were 4.67-0.6 mSv, which closely matched the reported doses of 4.51-0.43 mSv. The mean estimated annual doses between 1961 and 1983 based on the exposure scenarios developed by factors associated with radiation doses ranged from 11.08 to 4.82 mSv. The mean estimated annual doses of individuals in the cohort from 1961 to 1983 ranged from 11.15 to 4.88 mSv. Although caution needs to be exercised in the interpretation of these estimations due to uncertainty owed to the nature of extrapolation beyond the range of observed data, this study offers a sense of the radiation doses for workers during Korea's early period of radiation-related activities, which can be a useful piece of information for radiation-induced health risk assessments.

One-Decade-Spanning transgenerational effects of historic radiation dose in wild populations of bank voles exposed to radioactive contamination following the chernobyl nuclear disaster.

The concept of historic radiation doses associated with accidental radioactive releases and their role in leading to radiation-induced non-targeted effects on affected wild animals are currently being evaluated. Previous research studying Fukushima butterfly, Chernobyl bird and fruit fly populations shows that the effects are transgenerational, underlined by the principles of genomic instability, and varied from one species to another. To further expand on the responses of and their sensitivity in different taxonomically distinct groups, the present study sought to reconstruct historic radiation doses and delineate their effects on bank voles (Clethrionomys glareolus) found within a 400-km radius of the Chernobyl Nuclear Power Plant meltdown site. Historic dose reconstruction from the whole-body dose rates for the bank vole samples for their parental generation at the time of radioactive release was performed. Relationships between the historic doses and cytogenetic aberrations and embryonic lethality were examined via graphical presentations. Results suggest that genomic instability develops at the historic dose range of 20-51 mGy while a radioadaptive response develops at the historic dose range of 51-356 mGy. The Linear No-Threshold (LNT) relationship was absent at historic doses of lower than 356 mGy at all generations. However, LNT was apparent when the very high historic dose of 10.28 Gy in one sampling year was factored into the dose response curve for the bank vole generation 21-22. It is worth being reminded that natural mutation accumulation and other environmental stressors outside the realm of dose effects could contribute to the observed effects in a multiple-stressor environment. Nevertheless, the consistent development of genomic instability and radio-adaptive response across generations and sampling sites unearths the utmost fundamental radiobiological principle of transgenerational non-targeted effects. As a result, it calls for better attention and regulation from global governing bodies of environmental health protection.

Otoliths as object of EPR dosimetric research.

Otoliths are the organs which fish use for hearing and keeping balance. Otoliths are the most calcified tissues in the fish body. In contrast to bones, otoliths are not affected by remodeling and, therefore, they are expected to accumulate any dose from ionizing radiation during lifetime. Therefore, EPR dosimetry with fish otoliths could be an important tool for dose reconstruction in radiobiology and radioecology. It could also provide useful information remediation actions to de-contaminate waterbodies. Consequently, in the present study, otoliths of three contaminated fish species (roach (Rutilus rutilus), pike (Esox lucius) and perch (Perca Fluviatilis)) were examined with Electron Paramagnetic Resonance (EPR) spectroscopy. The fish were caught at storage reservoirs of liquid radioactive waste from Mayak PA and from the upper reach of the Techa River, which have been contaminated with different levels of radionuclide activity concentrations. It is shown that the radiation-induced EPR signal of otolith is stable and characterized by a linear dose response. However, the slope of the calibration curve (corresponding to the radiation sensitivity of the material) is not the same for different species; this may be caused by differences in mineralization. The reconstructed doses were found to be in the range from undetectable (in fish from the upper stream of the Techa River) up to 265 Gy (in roach from the most contaminated waterbody). In parallel, otoliths were measured with ß-counter to detect 90Sr/90Y. Samples were also tested on the presence of alpha-emitters, but no alpha activity above background could be detected. However, a significant activity concentration of 90Sr was detected (from 1 × 101 to 2 × 104 Bq/g). The EPR doses measured correlated with the 90Sr activity concentration measured in the otolith samples.

Comparison of neutron organ and effective dose coefficients for PIMAL stylized phantom in bent postures in standard irradiation geometries.

Neutron dose coefficients for standard irradiation geometries have been reported in International Commission on Radiological Protection (ICRP) Publication 116 for the ICRP Publication 110 adult reference phantoms. In the present work, organ and effective dose coefficients have been calculated for a receptor in both upright and articulated (bent) postures representing more realistic working postures exposed to a mono-energetic neutron radiation field. This work builds upon prior work by Dewji and co-workers comparing upright and bent postures for exposure to mono-energetic photon fields. Simulations were conducted using the Oak Ridge National Laboratory's articulated stylized adult phantom, "Phantom wIth Moving Arms and Legs" (PIMAL) software package, and the Monte Carlo N-Particle (MCNP) version 6.1.1 radiation transport code. Organ doses were compared for the upright and bent (45° and 90°) phantom postures for neutron energies ranging from 1 × 10- 9 to 20 MeV for the ICRP Publication 116 external exposure geometries-antero-posterior (AP), postero-anterior (PA), and left and right lateral (LLAT, RLAT). Using both male and female phantoms, effective dose coefficients were computed using ICRP Publication 103 methodology. The resulting coefficients for articulated phantoms were compared to those of the upright phantom. Computed organ and effective dose coefficients are discussed as a function of neutron energy, phantom posture, and source irradiation geometry. For example, it is shown here that for the AP and PA irradiation geometries, the differences in the organ coefficients between the upright and bent posture become more pronounced with increasing bending angle. In the AP geometry, the brain dose coefficients are expectedly higher in the bent postures than in the upright posture, while all other organs have lower dose coefficients, with the thyroid showing the greatest difference. Overall, the effective dose estimated for the upright phantom is more conservative than that for the articulated phantom, which may have ramifications in the estimation or reconstruction of radiation doses.

[Impact of Respiratory Motion on Dose Variation in Radiation Therapy for Glottic Cancer].

In radiation therapy for glottic cancer, respiratory motion of larynx may change the dose variation in the target. The purpose of this study is to measure the respiratory motion of the larynx, and quantify the impact of the motion on the dose variation. This study included 10 patients treated by opposing portal irradiation for glottic cancer. We acquired fluoroscopy and respiratory waveform of the patients simultaneously and formulated the relationship between the displacement of larynx and the respiratory phase. We divided one field into 39 sub-fields on the basis of control points. Dose distributions accounting for the displacement were calculated by shifting isocenter calculated using the formula in every sub-fields. Dose variations of clinical target volume (CTV) were evaluated by subtracting dose distributions with displacement consideration and dose distributions without it. Average amplitude and the maximum amplitude of respiratory motion were 2.5 and 8.7 mm, respectively. Average of mean dose variation in CTV was 0.1% of the prescribed dose, and maximum of local dose variation was 2.0% of the prescribed dose. Hence, it is realized that dose variation in CTV by respiratory motion was slight.

Comparison of photon organ and effective dose coefficients for PIMAL stylized phantom in bent positions in standard irradiation geometries.

Computational phantoms with articulated arms and legs have been constructed to enable the estimation of radiation dose in different postures. Through a graphical user interface, the Phantom wIth Moving Arms and Legs (PIMAL) version 4.1.0 software can be employed to articulate the posture of a phantom and generate a corresponding input deck for the Monte Carlo N-Particle (MCNP) radiation transport code. In this work, photon fluence-to-dose coefficients were computed using PIMAL to compare organ and effective doses for a stylized phantom in the standard upright position with those for phantoms in realistic work postures. The articulated phantoms represent working positions including fully and half bent torsos with extended arms for both the male and female reference adults. Dose coefficients are compared for both the upright and bent positions across monoenergetic photon energies: 0.05, 0.1, 0.5, 1.0, and 5.0 MeV. Additionally, the organ doses are compared across the International Commission on Radiological Protection's standard external radiation exposure geometries: antero-posterior, postero-anterior, left and right lateral, and isotropic (AP, PA, LLAT, RLAT, and ISO). For the AP and PA irradiation geometries, differences in organ doses compared to the upright phantom become more profound with increasing bending angles and have doses largely overestimated for all organs except the brain in AP and bladder in PA. In LLAT and RLAT irradiation geometries, energy deposition for organs is more likely to be underestimated compared to the upright phantom, with no overall change despite increased bending angle. The ISO source geometry did not cause a significant difference in absorbed organ dose between the different phantoms, regardless of position. Organ and effective fluence-to-dose coefficients are tabulated. In the AP geometry, the effective dose at the 45° bent position is overestimated compared to the upright phantom below 1 MeV by as much as 27% and 82% in the 90° position. The effective dose in the 45° bent position was comparable to that in the 90° bent position for the LLAT and RLAT irradiation geometries. However, the upright phantom underestimates the effective dose to PIMAL in the LLAT and RLAT geometries by as much as 30% at 50 keV.

External dose reconstruction for the former village of Metlino (Techa River, Russia) based on environmental surveys, luminescence measurements, and radiation transport modelling.

In the first years of its operation, the Mayak Production Association, a facility part of the Soviet nuclear weapons program in the Southern Urals, Russia, discharged large amounts of radioactively contaminated effluent into the nearby Techa River, thus exposing the people living at this river to external and internal radiations. The Techa River Cohort is a cohort intensely studied in epidemiology to investigate the correlation between low-dose radiation and health effects on humans. For the individuals in the cohort, the Techa River Dosimetry System describes the accumulated dose in human organs and tissues. In particular, organ doses from external exposure are derived from estimates of dose rate in air on the Techa River banks which were estimated from measurements and Monte Carlo modelling. Individual doses are calculated in accordance with historical records of individuals' residence histories, observational data of typical lifestyles for different age groups, and age-dependent conversion factors from air kerma to organ dose. The work here describes an experimentally independent assessment of the key input parameter of the dosimetry system, the integral air kerma, for the former village of Metlino, upper Techa River region. The aim of this work was thus to validate the Techa River Dosimetry System for the location of Metlino in an independent approach. Dose reconstruction based on dose measurements in bricks from a church tower and Monte Carlo calculations was used to model the historic air kerma accumulated in the time from 1949 to 1956 at the shoreline of the Techa River in Metlino. Main issues are caused by a change in the landscape after the evacuation of the village in 1956. Based on measurements and published information and data, two separate models for the historic pre-evacuation geometry and for the current geometry of Metlino were created. Using both models, a value for the air kerma was reconstructed, which agrees with that obtained in the Techa River Dosimetry System within a factor of two.

Comparison of DVH-based plan verification methods for VMAT: ArcCHECK-3DVH system and dynalog-based dose reconstruction.

The purpose of this study was comparing dose-volume histogram (DVH)-based plan verification methods for volumetric modulated arc therapy (VMAT) pretreatment QA. We evaluated two 3D dose reconstruction systems: ArcCHECK-3DVH system (Sun Nuclear corp.) and Varian dynalog-based dose reconstruction (DBDR) system, developed in-house. Fifteen prostate cancer patients (67.6 Gy/26 Fr), four head and neck cancer patient (66 Gy/33 Fr), and four esophagus cancer patients (60 Gy/30 Fr) treated with VMAT were studied. First, ArcCHECK measurement was performed on all plans; simultaneously, the Varian dynalog data sets that contained the actual delivered parameters (leaf positions, gantry angles, and cumulative MUs) were acquired from the Linac control system. Thereafter, the delivered 3D patient dose was reconstructed by 3DVH software (two different calculating modes were used: High Sensitivity (3DVH-HS) and Normal Sensitivity (3DVH-NS)) and in-house DBDR system. We evaluated the differences between the TPS-calculated dose and the reconstructed dose using 3D gamma passing rates and DVH dose index analysis. The average 3D gamma passing rates (3%/3 mm) between the TPS-calculated dose and the reconstructed dose were 99.1 ± 0.6%, 99.7 ± 0.3%, and 100.0 ± 0.1% for 3DVH-HS, 3DVH-NS, and DBDR, respectively. For the prostate cases, the average differences between the TPS-calculated dose and reconstructed dose in the PTV mean dose were 1.52 ± 0.50%, -0.14 ± 0.55%, and -0.03 ± 0.07% for 3DVH-HS, 3DVH-NS, and DBDR, respectively. For the head and neck and esophagus cases, the dose difference to the TPS-calculated dose caused by an effect of heterogeneity was more apparent under the 3DVH dose reconstruction than the DBDR. Although with some residual dose reconstruction errors, these dose reconstruction methods can be clinically used as effective tools for DVH-based QA for VMAT delivery.

Analysis of EPR and FISH studies of radiation doses in persons who lived in the upper reaches of the Techa River.

Waterborne radioactive releases into the Techa River from the Mayak Production Association in Russia during 1949-1956 resulted in significant doses to about 30,000 persons who lived in downstream settlements. The residents were exposed to internal and external radiation. Two methods for reconstruction of the external dose are considered in this paper, electron paramagnetic resonance (EPR) measurements of teeth, and fluorescence in situ hybridization (FISH) measurements of chromosome translocations in circulating lymphocytes. The main issue in the application of the EPR and FISH methods for reconstruction of the external dose for the Techa Riverside residents was strontium radioisotopes incorporated in teeth and bones that act as a source of confounding local exposures. In order to estimate and subtract doses from incorporated (89,90)Sr, the EPR and FISH assays were supported by measurements of (90)Sr-body burdens and estimates of (90)Sr concentrations in dental tissues by the luminescence method. The resulting dose estimates derived from EPR to FISH measurements for residents of the upper Techa River were found to be consistent: The mean values vary from 510 to 550 mGy for the villages located close to the site of radioactive release to 130-160 mGy for the more distant villages. The upper bound of individual estimates for both methods is equal to 2.2-2.3 Gy. The EPR- and FISH-based dose estimates were compared with the doses calculated for the donors using the most recent Techa River Dosimetry System (TRDS). The TRDS external dose assessments are based on the data on contamination of the Techa River floodplain, simulation of air kerma above the contaminated soil, age-dependent lifestyles and individual residence histories. For correct comparison, TRDS-based doses were calculated from two sources: external exposure from the contaminated environment and internal exposure from (137)Cs incorporated in donors' soft tissues. It is shown here that the TRDS-based absorbed doses in tooth enamel and muscle are in agreement with EPR- and FISH-based estimates within uncertainty bounds. Basically, this agreement between the estimates has confirmed the validity of external doses calculated with the TRDS.

Spectrum learning for super-resolution tomographic reconstruction.

Objective. Computed Tomography (CT) has been widely used in industrial high-resolution non-destructive testing. However, it is difficult to obtain high-resolution images for large-scale objects due to their physical limitations. The objective is to develop an improved super-resolution technique that preserves small structures and details while efficiently capturing high-frequency information.Approach. The study proposes a new deep learning based method called spectrum learning (SPEAR) network for CT images super-resolution. This approach leverages both global information in the image domain and high-frequency information in the frequency domain. The SPEAR network is designed to reconstruct high-resolution images from low-resolution inputs by considering not only the main body of the images but also the small structures and other details. The symmetric property of the spectrum is exploited to reduce weight parameters in the frequency domain. Additionally, a spectrum loss is introduced to enforce the preservation of both high-frequency components and global information.Main results. The network is trained using pairs of low-resolution and high-resolution CT images, and it is fine-tuned using additional low-dose and normal-dose CT image pairs. The experimental results demonstrate that the proposed SPEAR network outperforms state-of-the-art networks in terms of image reconstruction quality. The approach successfully preserves high-frequency information and small structures, leading to better results compared to existing methods. The network's ability to generate high-resolution images from low-resolution inputs, even in cases of low-dose CT images, showcases its effectiveness in maintaining image quality.Significance. The proposed SPEAR network's ability to simultaneously capture global information and high-frequency details addresses the limitations of existing methods, resulting in more accurate and informative image reconstructions. This advancement can have substantial implications for various industries and medical diagnoses relying on accurate imaging.

Cardiorespiratory motion characteristics and their dosimetric impact on cardiac stereotactic body radiotherapy.

BACKGROUND: Cardiac stereotactic body radiotherapy (CSBRT) is an emerging and promising noninvasive technique for treating refractory arrhythmias utilizing highly precise, single or limited-fraction high-dose irradiations. This method promises to revolutionize the treatment of cardiac conditions by delivering targeted therapy with minimal exposure to surrounding healthy tissues. However, the dynamic nature of cardiorespiratory motion poses significant challenges to the precise delivery of dose in CSBRT, introducing potential variabilities that can impact treatment efficacy. The complexities of the influence of cardiorespiratory motion on dose distribution are compounded by interplay and blurring effects, introducing additional layers of dose uncertainty. These effects, critical to the understanding and improvement of the accuracy of CSBRT, remain unexplored, presenting a gap in current clinical literature. PURPOSE: To investigate the cardiorespiratory motion characteristics in arrhythmia patients and the dosimetric impact of interplay and blurring effects induced by cardiorespiratory motion on CSBRT plan quality. METHODS: The position and volume variations in the substrate target and cardiac substructures were evaluated in 12 arrhythmia patients using displacement maximum (DMX) and volume metrics. Moreover, a four-dimensional (4D) dose reconstruction approach was employed to examine the dose uncertainty of the cardiorespiratory motion. RESULTS: Cardiac pulsation induced lower DMX than respiratory motion but increased the coefficient of variation and relative range in cardiac substructure volumes. The mean DMX of the substrate target was 0.52 cm (range: 0.26-0.80 cm) for cardiac pulsation and 0.82 cm (range: 0.32-2.05 cm) for respiratory motion. The mean DMX of the cardiac structure ranged from 0.15 to 1.56 cm during cardiac pulsation and from 0.35 to 1.89 cm during respiratory motion. Cardiac pulsation resulted in an average deviation of -0.73% (range: -4.01%-4.47%) in V25 between the 3D and 4D doses. The mean deviations in the homogeneity index (HI) and gradient index (GI) were 1.70% (range: -3.10%-4.36%) and 0.03 (range: -0.14-0.11), respectively. For cardiac substructures, the deviations in D50 due to cardiac pulsation ranged from -1.88% to 1.44%, whereas the deviations in Dmax ranged from -2.96% to 0.88% of the prescription dose. By contrast, the respiratory motion led to a mean deviation of -1.50% (range: -10.73%-4.23%) in V25. The mean deviations in HI and GI due to respiratory motion were 4.43% (range: -3.89%-13.98%) and 0.18 (range: -0.01-0.47) (p < 0.05), respectively. Furthermore, the deviations in D50 and Dmax in cardiac substructures for the respiratory motion ranged from -0.28% to 4.24% and -4.12% to 1.16%, respectively. CONCLUSIONS: Cardiorespiratory motion characteristics vary among patients, with the respiratory motion being more significant. The intricate cardiorespiratory motion characteristics and CSBRT plan complexity can induce substantial dose uncertainty. Therefore, assessing individual motion characteristics and 4D dose reconstruction techniques is critical for implementing CSBRT without compromising efficacy and safety.

Quantifying dose uncertainties resulting from cardiorespiratory motion in intensity-modulated proton therapy for cardiac stereotactic body radiotherapy.

Background: Cardiac stereotactic body radiotherapy (CSBRT) with photons efficaciously and safely treats cardiovascular arrhythmias. Proton therapy, with its unique physical and radiobiological properties, can offer advantages over traditional photon-based therapies in certain clinical scenarios, particularly pediatric tumors and those in anatomically challenging areas. However, dose uncertainties induced by cardiorespiratory motion are unknown. Objective: This study investigated the effect of cardiorespiratory motion on intensity-modulated proton therapy (IMPT) and the effectiveness of motion-encompassing methods. Methods: We retrospectively included 12 patients with refractory arrhythmia who underwent CSBRT with four-dimensional computed tomography (4DCT) and 4D cardiac CT (4DcCT). Proton plans were simulated using an IBA accelerator based on the 4D average CT. The prescription was 25 Gy in a single fraction, with all plans normalized to ensure that 95% of the target volume received the prescribed dose. 4D dose reconstruction was performed to generate 4D accumulated and dynamic doses. Furthermore, dose uncertainties due to the interplay effect of the substrate target and organs at risk (OARs) were assessed. The differences between internal organs at risk volume (IRV) and OARreal (manually contoured on average CT) were compared. In 4D dynamic dose, meeting prescription requirements entails V25 and D95 reaching 95% and 25 Gy, respectively. Results: The 4D dynamic dose significantly differed from the 3D static dose. The mean V25 and D95 were 89.23% and 24.69 Gy, respectively, in 4DCT and 94.35% and 24.99 Gy, respectively, in 4DcCT. Eleven patients in 4DCT and six in 4DcCT failed to meet the prescription requirements. Critical organs showed varying dose increases. All metrics, except for Dmean and D50, significantly changed in 4DCT; in 4DcCT, only D50 remained unchanged with regards to the target dose uncertainties induced by the interplay effect. The interplay effect was only significant for the Dmax values of several OARs. Generally, respiratory motion caused a more pronounced interplay effect than cardiac pulsation. Neither IRV nor OARreal effectively evaluated the dose discrepancies of the OARs. Conclusions: Complex cardiorespiratory motion can introduce dose uncertainties during IMPT. Motion-encompassing techniques may mitigate but cannot entirely compensate for the dose discrepancies. Individualized 4D dose assessments are recommended to verify the effectiveness and safety of CSBRT.

Quantifying dose perturbations in high-risk prostate radiotherapy due to translational and rotational motion of prostate and pelvic lymph nodes.

BACKGROUND: Radiotherapy of the prostate and the pelvic lymph nodes (LN) is a part of the standard of care treatment for high-risk prostate cancer. The independent translational and rotational (i.e., six-degrees-of-freedom, [6DoF]) motion of the prostate and LN target during and between fractions can perturb the dose distribution. However, no standard dose reconstruction method accounting for differential 6DoF target motion is available. PURPOSE: We present a framework for monitoring motion-induced dose perturbations for two independently moving target volumes in 6DoF. The framework was used to determine the dose perturbation for the prostate and the LN target caused by differential 6DoF motion for a cohort of high-risk prostate cancer patients. As a potential first step toward real-time dose-guided high-risk prostate radiotherapy, we furthermore investigated if the dose reconstruction was fast enough for real-time application for both targets. METHODS: Twenty high-risk prostate cancer patients were treated with 3-arc volumetric modulated arc therapy (VMAT). Kilovoltage intrafraction monitoring (KIM) with triggered kilovoltage (kV) images acquired every 3 throughout 7-10 fractions per patient was used for retrospective 6DoF intrafraction prostate motion estimation. The 6DoF interfraction LN motion was determined from a pelvic bone match between the planning CT and a post-treatment cone beam CT (CBCT). Using the retrospectively extracted motion, real-time 6DoF motion-including dose reconstruction was simulated using the in-house developed software DoseTracker. A data stream with the 6DoF target positions and linac parameters was broadcasted at a 3-Hz frequency to DoseTracker. In a continuous loop, DoseTracker calculated the target dose increments including the specified motion and, for comparison, without motion. The motion-induced change in D99.5% for the prostate CTV (ΔD99.5%) and in D98% for the LN CTV (ΔD98%) was calculated using the final cumulative dose of each fraction and averaged over all imaged fractions. The real-time reconstructed dose distribution of DoseTracker was benchmarked against a clinical treatment planning system (TPS) and it was investigated whether the calculation speed was fast enough to keep up with the incoming data stream. RESULTS: Translational motion was largest in cranio-caudal (CC) direction (prostate: [-5.9, +8.4] mm; LN: [-9.9; +11.0] mm) and anterior-posterior (AP) direction (prostate:[-5.6; +6.9] mm; LN: [-9.6; +11.0] mm). The pitch was the largest rotation (prostate: [-22.5; +25.2] deg; LN: [-3.9; +5.5] deg). The prostate CTV ΔD99.5% was [-16.2; +2.5]% for single fractions and [-3.0; +1.7]% when averaged over all imaged fractions. The LN CTV ΔD98% was [-19.8; +1.2]% for single fractions and [-3.1; +0.9]% after averaging. Mean (Standard deviation) absolute dose errors in DoseTracker of 107.8% (Std: 1.9%) for the prostate and 105.5% (Std:1.4%) for the LN were corrected during dose reconstruction by automatically calculated normalization factors. It resulted in accurate calculation of the motion-induced dose errors with relative differences between DoseTracker and TPS dose calculations of -0.1% (Std: 0.5%) (prostate CTV ΔD99.5%) and -0.2% (Std: 0.5%) (LN CTV ΔD98%). The DoseTracker calculation was fast enough to keep up with the incoming inputs for all but two out of 107 184 dose calculations. CONCLUSION: Using the developed framework for dose perturbation monitoring, we found that the differential 6DoF target motion caused substantial dose perturbation for individual fractions, which largely averaged out after several fractions. The framework was shown to provide reliable dose calculations and a sufficiently high-dose reconstruction speed to be applicable in real-time.