The dosimetric accuracy of a commercial model-based dose calculation algorithm in modeling a six-groove direction modulated brachytherapy tandem applicator.

OBJECTIVE: With advancements in high-dose rate brachytherapy, the clinical translation of intensity modulated brachytherapy (IMBT) innovations necessitates utilization of model-based dose calculation algorithms (MBDCA) for accurate and rapid dose calculations. This study uniquely benchmarks a commercial MBDCA, BrachyVision ACUROSTM (BVA), against Monte Carlo (MI) simulations, evaluating dose distributions for a novel IMBT applicator, termed as the six-groove Direction Modulated Brachytherapy (DMBT) tandem, expanding beyond previous focus on partially shielded vaginal cylinder applicators, through a novel methodology. Approach: The DMBT tandem applicator, made of a tungsten alloy with six evenly spaced grooves, was simulated using the GEANT4 MC code. Subsequently, two main scenarios were created using the BVA and reproduced by the MC simulations: "Source at the Center of the Water Phantom (SACWP)" and "Source at the Middle of the Applicator (SAMA)" for three cubical virtual water phantoms (20 cm)3, (30 cm)3, and (40 cm)3. A track length estimator was utilized for dose calculation and 2D/3D scoring were performed. The difference in isodose surfaces/lines (i.e., coverage) at each voxel, ∆DIsodose Levels/Lines, was thus calculated for relevant normalization points (rref). Results: The coverage was comparable, based on 2D scoring, between the BVA and MC isodose surfaces/lines for the region of clinical relevance, (i.e., within 5 cm radius from the source) with ∆DIsodose Lines (rref: 1 cm from the source) falling within 2% for the two scenarios for all phantom sizes. For the phantom (20 cm)3, ∆DIsodose Level (3D scoring) recorded the range [-3.0% +6.5%] ([-7.4% +7.3%]) for 95% of the voxels of the same scoring volume for the SACWP (SAMA) scenario. Significance: The results indicated that the BVA could render comparable coverage as compared to the MC simulations in the region of clinical relevance for different phantom sizes. ∆DIsodose Lines may offer an advantageous metric for evaluation of MBDCAs in clinical setting. .

Impact of Dose Calculation Algorithms and Radiobiological Parameters on Prediction of Cardiopulmonary Complications in Left Breast Radiation Therapy.

Background: Breast cancer requires evaluating treatment plans using dosimetric and biological parameters. Considering radiation dose distribution and tissue response, healthcare professionals can optimize treatment plans for better outcomes. Objective: This study aimed to evaluate the effects of the different Dose Calculation Algorithms (DCAs) and Biologically Model-Related Parameters (BMRPs) on the prediction of cardiopulmonary complications due to left breast radiotherapy. Material and Methods: In this practical study, the treatment plans of 21 female patients were simulated in the Monaco Treatment Planning System (TPS) with a prescribed dose of 50 Gy in 25 fractions. Dose distribution was extracted using the three DCAs [Pencil Beam (PB), Collapsed Cone (CC), and Monte Carlo (MC)]. Cardiopulmonary complications were predicted by Normal Tissue Complication Probability (NTCP) calculations using different dosimetric and biological parameters. The Lyman-Kutcher-Burman (LKB) and Relative-Seriality (RS) models were used to calculate NTCP. The endpoint for NTCP calculation was pneumonitis, pericarditis, and late cardiac mortality. The ANOVA test was used for statistical analysis. Results: In calculating Tumor Control Probability (TCP), a statistically significant difference was observed between the results of DCAs in the Poisson model. The PB algorithm estimated NTCP as less than others for all Pneumonia BMRPs. Conclusion: The impact of DCAs and BMRPs differs in the estimation of TCP and NTCP. DCAs have a stronger influence on TCP calculation, providing more effective results. On the other hand, BMRPs are more effective in estimating NTCP. Consequently, parameters for radiobiological indices should be cautiously used s to ensure the appropriate consideration of both DCAs and BMRPs.

Neutron flux evaluation algorithm with a combination of Monte Carlo and removal-diffusion calculation methods for boron neutron capture therapy.

BACKGROUND: In Japan, the clinical treatment of boron neutron capture therapy (BNCT) has been applied to unresectable, locally advanced, and recurrent head and neck carcinomas using an accelerator-based neutron source since June of 2020. Considering the increase in the number of patients receiving BNCT, efficiency of the treatment planning procedure is becoming increasingly important. Therefore, novel and rapid dose calculation algorithms must be developed. We developed a novel algorithm for calculating neutron flux, which comprises of a combination of a Monte Carlo (MC) method and a method based on the removal-diffusion (RD) theory (RD calculation method) for the purpose of dose calculation of BNCT. PURPOSE: We present the details of our novel algorithm and the verification results of the calculation accuracy based on the MC calculation result. METHODS: In this study, the "MC-RD" calculation method was developed, wherein the RD calculation method was used to calculate the thermalization process of neutrons and the MC method was used to calculate the moderation process. The RD parameters were determined by MC calculations in advance. The MC-RD calculation accuracy was verified by comparing the results of the MC-RD and MC calculations with respect to the neutron flux distributions in each of the cubic and head phantoms filled with water. RESULTS: Comparing the MC-RD calculation results with those of MC calculations, it was found that the MC-RD calculation accurately reproduced the thermal neutron flux distribution inside the phantom, with the exception of the region near the surface of the phantom. CONCLUSIONS: The MC-RD calculation method is useful for the evaluation of the neutron flux distribution for the purpose of BNCT dose calculation, except for the region near the surface.

Dosimetric verification of four dose calculation algorithms for spine stereotactic body radiotherapy.

The applications of Type B [anisotropic analytical algorithm (AAA) and collapsed cone (CC)] and Type C [Acuros XB (AXB) and photon Monte Carlo (PMC)] dose calculation algorithms in spine stereotactic body radiotherapy (SBRT) were evaluated. Water- and bone-equivalent phantoms were combined to evaluate the percentage depth dose and dose profile. Subsequently, 48 consecutive patients with clinical spine SBRT plans were evaluated. All treatment plans were created using AXB in Eclipse. The prescription dose was 24 Gy in two fractions at a 10 MV FFF on TrueBeam. The doses were then recalculated with AAA, CC and PMC while maintaining the AXB-calculated monitor units and beam arrangement. The dose index values obtained using the four dose calculation algorithms were then compared. The AXB and PMC dose distributions agreed with the bone-equivalent phantom measurements (within ±2.0%); the AAA and CC values were higher than those in the bone-equivalent phantom region. For the spine SBRT plans, PMC, AAA and CC were overestimated compared with AXB in terms of the near minimum and maximum doses of the target and organ at risk, respectively; the mean dose difference was within 4.2%, which is equivalent with within 1 Gy. The phantom study showed that the results from AXB and PMC agreed with the measurements within ±2.0%. However, the mean dose difference ranged from 0.5 to 1 Gy in the spine SBRT planning study when the dose calculation algorithms changed. Users should incorporate a clinical introduction that includes an awareness of these differences.

Difference in target dose distributions between Acuros XB and collapsed cone convolution/superposition and the impact of the tumor locations in clinical cases of stereotactic ablative body radiotherapy for lung cancer.

Objectives: The objective of the study is to analyze the difference in target dose distributions between Acuros XB (AXB) and collapsed cone convolution (CCC)/superposition and the impact of the tumor locations in clinical cases of stereotactic ablative body radiotherapy (SABR) for lung cancer. Materials and Methods: Ninety-six patients underwent SABR for lung cancers Kyushu University Hospital from 2014 to 2017. We recalculated clinical plans originally calculated by AXB using CCC with the identical monitor units (MUs) and beam arrangements. We calculated the following dosimetric parameters: maximum dose (Dmax), minimum dose (Dmin), homogeneity index (HI), conformity index (CI), and D95 of the planning target volume (PTV). We investigated the difference between the results of two calculations and examined the impact of tumor location. Moreover, we determined the target central dose using a thorax phantom and assessed the calculation accuracy of the two algorithms for each fraction. Results: CCC significantly overestimated the dose to PTV, compared to AXB (P < 0.05). The mean differences of Dmax, Dmin, and D95 were 1.17, 1.95, and 1.85 Gy, respectively. The mean differences of HI and CI were 0.02 and - 0.06. Dmin, HI, and D95 had significant correlations with the tumor location, and the difference was greater when the PTV was included the chest wall (P < 0.05). The discrepancy between the calculated and irradiated dose was 2.48% for CCC, whereas it was 0.14% for AXB. Conclusions: We demonstrated that CCC significantly overestimated the dose to PTV relative to AXB in clinical cases of lung SABR.

Development and evaluation of dose calculation algorithm with a combination of Monte Carlo and point-kernel methods for boron neutron capture therapy.

We developed a 'hybrid algorithm' that combines the Monte Carlo (MC) and point-kernel methods for fast dose calculation in boron neutron capture therapy. The objectives of this study were to experimentally verify the hybrid algorithm and to verify the calculation accuracy and time of a 'complementary approach' adopting both the hybrid algorithm and the full-energy MC method. In the latter verification, the results were compared with those obtained using the full-energy MC method alone. In the hybrid algorithm, the moderation process of neutrons is simulated using only the MC method, and the thermalization process is modeled as a kernel. The thermal neutron fluxes calculated using only this algorithm were compared with those measured in a cubic phantom. In addition, a complementary approach was used for dose calculation in a geometry simulating the head region, and its computation time and accuracy were verified. The experimental verification indicated that the thermal neutron fluxes calculated using only the hybrid algorithm reproduced the measured values at depths exceeding a few centimeters, whereas they overestimated those at shallower depths. Compared with the calculation using only the full-energy MC method, the complementary approach reduced the computation time by approximately half, maintaining nearly same accuracy. When focusing on the calculation only using the hybrid algorithm only for the boron dose attributed to the reaction of thermal neutrons, the computation time was expected to reduce by 95% compared with the calculation using only the full-energy MC method. In conclusion, modeling the thermalization process as a kernel was effective for reducing the computation time.

Planning issues on linac-based stereotactic radiotherapy.

This work aims to summarize and evaluate the current planning progress based on the linear accelerator in stereotactic radiotherapy (SRT). The specific techniques include 3-dimensional conformal radiotherapy, dynamic conformal arc therapy, intensity-modulated radiotherapy, and volumetric-modulated arc therapy (VMAT). They are all designed to deliver higher doses to the target volume while reducing damage to normal tissues; among them, VMAT shows better prospects for application. This paper reviews and summarizes several issues on the planning of SRT to provide a reference for clinical application.

Development and validation of a scatter-corrected CBCT image-guided method for cervical cancer brachytherapy.

Background and purpose: Multiple patient transfers have a nonnegligible impact on the accuracy of dose delivery for cervical cancer brachytherapy. We consider using on-site cone-beam CT (CBCT) to resolve this problem. However, CBCT clinical applications are limited due to inadequate image quality. This paper implements a scatter correction method using planning CT (pCT) prior to obtaining high-quality CBCT images and evaluates the dose calculation accuracy of CBCT-guided brachytherapy for cervical cancer. Materials and methods: The CBCT of a self-developed female pelvis phantom and five patients was first corrected using empirical uniform scatter correction in the projection domain and further corrected in the image domain. In both phantom and patient studies, the CBCT image quality before and after scatter correction was evaluated with registered pCT (rCT). Model-based dose calculation was performed using the commercial package Acuros®BV. The dose distributions of rCT-based plans and corrected CBCT-based plans in the phantom and patients were compared using 3D local gamma analysis. A statistical analysis of the differences in dosimetric parameters of five patients was also performed. Results: In both phantom and patient studies, the HU error of selected ROIs was reduced to less than 15 HU. Using the dose distribution of the rCT-based plan as the baseline, the γ pass rate (2%, 2 mm) of the corrected CBCT-based plan in phantom and patients all exceeded 98% and 93%, respectively, with the threshold dose set to 3, 6, 9, and 12 Gy. The average percentage deviation (APD) of D90 of HRCTV and D2cc of OARs was less than 1% between rCT-based and corrected CBCT-based plans. Conclusion: Scatter correction using a pCT prior can effectively improve the CBCT image quality and CBCT-based cervical brachytherapy dose calculation accuracy, indicating promising prospects in both simplified brachytherapy processes and accurate brachytherapy dose delivery.

Quantitative evaluation of dosimetric uncertainties associated with small electron fields.

INTRODUCTION: Although many studies have investigated small electron fields, there are several dosimetric issues that are not well understood. This includes lack of charged particle equilibrium, lateral scatter, source occlusion and volume averaging of the detectors used in the measurement of the commissioning data. High energy electron beams are also associated with bremsstrahlung production that contributes to dose deposition, which is not well investigated, particularly for small electron fields. The goal of this work has been to investigate dosimetric uncertainties associated with small electron fields using dose measurements with different detectors as well as calculations with eMC dose calculation algorithm. METHODS: Different dosimetric parameters including output factors, depth dose curves and dose profiles from small electron field cutouts were investigated quantitatively. These dosimetric parameters were measured using different detectors that included small ion chambers and diodes for small electron cutouts with diameters ranging from 15-50mm mounted on a 6 × 6cm2 cone with beam energies from 6-20MeV. RESULTS: Large deviations existed between the output factors calculated with the eMC algorithm and measured with small detectors for small electron fields up to 55% for 6MeV. The discrepancy between the calculated and measured doses increased 10%-55% with decreasing electron beam energy from 20 MeV to 6 MeV for 15mm circular field. For electron fields with cutouts 20mm and larger, the measured and calculated doses agreed within 5% for all electron energies from 6-20MeV. For small electron fields, the maximal depth dose shifted upstream and becomes more superficial as the electron beam energy increases from 6-20MeV as measured with small detectors. DISCUSSION: Large dose discrepancies were found between the measured and calculated doses for small electron fields where the eMC underestimated output factors by 55% for small circular electron fields with a diameter of 15 mm, particularly for low energy electron beams. The measured entrance doses and dmax of the depth dose curves did not agree with the corresponding values calculated by eMC. Furthermore, the measured dose profiles showed enhanced dose deposition in the umbra region and outside the small fields, which mostly resulted from dose deposition from the bremsstrahlung produced by high energy electrons that was not accounted for by the eMC algorithm due to inaccurate modeling of the lateral dose deposition from bremsstrahlung. CONCLUSION: Electron small field dosimetry require more consideration of variations in beam quality, lack of charged particle equilibrium, lateral scatter loss and dose deposition from bremsstrahlung produced by energetic electron beams in a comprehensive approach similar to photon small field dosimetry. Furthermore, most of the commercially available electron dose calculation algorithms are commissioned with large electron fields; therefore, vendors should provide tools for the modeling of electron dose calculation algorithms for small electron fields.

A polar-coordinate-based pencil beam algorithm for VMAT dose computation with high-resolution gantry angle sampling.

PURPOSE: Most commercially available treatment planning systems (TPSs) approximate the continuous delivery of volumetric modulated arc therapy (VMAT) plans with a series of discretized static beams for treatment planning, which can make VMAT dose computation extremely inefficient. In this study, we developed a polar-coordinate-based pencil beam (PB) algorithm for efficient VMAT dose computation with high-resolution gantry angle sampling that can improve the computational efficiency and reduce the dose discrepancy due to the angular under-sampling effect. METHODS AND MATERIALS: 6 MV 1 × 1 m m 2 $1 \times 1{\rm{\;m}}{{\rm{m}}^2}$ pencil beams were simulated on a uniform cylindrical phantom under an EGSnrc Monte Carlo (MC) environment. The MC-generated PB kernels were collected in the polar coordinate system for each bixel on a 40 × 40 c m 2 $40 \times 40{\rm{\;c}}{{\rm{m}}^2}$ fluence map and subsequently fitted via a series of Gaussians. The fluence was calculated using a detectors' eye view with off-axis and MLC transmission factors corrected. Doses of VMAT arc on the phantom were computed by summing the convolution results between the corresponding PB kernels and fluence for each bixel in the polar coordinate system. The convolution was performed using fast Fourier transform to expedite the computing speed. The calculated doses were converted to the Cartesian coordinate system and compared with the reference dose computed by a collapsed cone convolution (CCC) algorithm of the TPS. A heterogeneous phantom was created to study the heterogeneity corrections using the proposed algorithm. Ten VMAT arcs were included to evaluate the algorithm performance. Gamma analysis and computation complexity theory were used to measure the dosimetric accuracy and computational efficiency, respectively. RESULTS: The dosimetric comparisons on the homogeneous phantom between the proposed PB algorithm and the CCC algorithm for 10 VMAT arcs demonstrate that the proposed algorithm can achieve a dosimetric accuracy comparable to that of the CCC algorithm with average gamma passing rates of 96% (2%/2mm) and 98% (3%/3mm). In addition, the proposed algorithm can provide better computational efficiency for VMAT dose computation using a PC equipped with a 4-core processor, compared to the CCC algorithm utilizing a dual 10-core server. Moreover, the computation complexity theory reveals that the proposed algorithm has a great advantage with regard to computational efficiency for VMAT dose computation on homogeneous medium, especially when a fine angular sampling rate is applied. This can support a reduction in dose errors from the angular under-sampling effect by using a finer angular sampling rate, while still preserving a practical computing speed. For dose calculation on the heterogeneous phantom, the proposed algorithm with heterogeneity corrections can still offer a reasonable dosimetric accuracy with comparable computational efficiency to that of the CCC algorithm. CONCLUSIONS: We proposed a novel polar-coordinate-based pencil beam algorithm for VMAT dose computation that enables a better computational efficiency while maintaining clinically acceptable dosimetric accuracy and reducing dose error caused by the angular under-sampling effect. It also provides a flexible VMAT dose computation structure that allows adjustable sampling rates and direct dose computation in regions of interest, which makes the algorithm potentially useful for clinical applications such as independent dose verification for VMAT patient-specific QA.

Dosiomic feature comparison between dose-calculation algorithms used for lung stereotactic body radiation therapy.

To evaluate the reproducibility of dose-based radiomic (dosiomic) features between dose-calculation algorithms for lung stereotactic body radiation therapy (SBRT). We analyzed 105 patients with early-stage non-small cell lung cancer who underwent lung SBRT between March 2011 and December 2017. Radiation doses of 48, 60, and 70 Gy were prescribed to the isocenter in 4-8 fractions. Dose calculations were performed using X-ray voxel Monte Carlo (XVMC) on the iPlan radiation treatment planning system (RTPS). Thereafter, the radiation doses were recalculated using the Acuros XB (AXB) and analytical anisotropic algorithm (AAA) on the Eclipse RTPS while maintaining the XVMC-calculated monitor units and beam arrangements. A total of 6808 dosiomic features were extracted without preprocessing (112 shape, 144 first-order, and 600 texture features) or with wavelet filters to eight decompositions (1152 first-order and 4800 texture features). Features with absolute pairwise concordance correlation coefficients-|CCcon|-values exceeding or equaling 0.85 were considered highly reproducible. Subgroup analyses were performed considering the wavelet filters and prescribed doses. The numbers of highly reproducible first-order and texture features were 34.8%, 26.9%, and 31.0% for the XVMC-AXB, XVMC-AAA, and AXB-AAA pairs, respectively. The maximum difference between the mean |CCcon| values was 0.70 and 0.11 for the subgroup analyses of wavelet filters and prescribed dose, respectively. The application of wavelet filter-based dosiomic analyses may be limited when using different types of dose-calculation algorithms for lung SBRT.

Dosimetric study of Varian Universal Multi-Channel Cylinder System for High-Dose-Rate 192Ir Brachytherapy.

PURPOSE: Recently, the Varian multichannel vaginal cylinder (MCVC) set for high-dose-rate 192Ir brachytherapy was commercially released. This MCVC was distinct from our existing MCVC in its peripheral channel layout and tip design. This investigation sought to assess the dosimetric impact of these changes. METHODS AND MATERIALS: The dimensions of the virtual model for each applicator were compared against both physical and radiographic measurements. Volumetric dose distributions were generated in silico using a model-based dose calculation algorithm (MBDCA). To characterize the effects of the new peripheral channel layout on dose to adjacent areas ("dose-spill"), point doses were compared using two sets of applicator-based reference points: at surface or 5 mm radially from surface. To evaluate the dose-shaping capabilities, a dose distribution was generated for the new applicator and assessed against a representative dose distribution for a patient previously treated with existing equipment. RESULTS: Based on both physical and radiographic measurements, virtual models were representative of each applicator within ±1 mm. Commissioning of the MBDCA was benchmarked based on AAPM Working Group on Dose Calculation Algorithms in Brachytherapy. The layout of the new applicator reduced dose-spill to other reference points significantly, as much as a factor of 16.3, compared with the existing equipment. The rounded tip shape and curve of the peripheral channels in the new applicator produced more conformity to its HR-CTV than existing equipment. CONCLUSIONS: Compared with our existing equipment, the design changes in the new Varian MCVC set offered improved control of dose spill and better conformality to HR-CTV.

Quantitative evaluation of dosimetric uncertainties in electron therapy by measurement and calculation using the electron Monte Carlo dose algorithm in the Eclipse treatment planning system.

In the electron beam radiation therapy, customized blocks are mostly used to shape treatment fields to generate conformal doses. The goal of this study is to investigate quantitatively dosimetric uncertainties associated with heterogeneities, detectors used in the measurement of the beam data commissioning, and modeling of the interactions of high energy electrons with tissue. These uncertainties were investigated both by measurements with different detectors and calculations using electron Monte Carlo algorithm (eMC) in the Eclipse treatment planning system. Dose distributions for different field sizes were calculated using eMC and measured with a multiple-diode-array detector (MapCheck2) for cone sizes ranging from 6 to 25 cm. The dose distributions were calculated using the CT images of the MapCheck2 and water-equivalent phantoms. In the umbra region (<20% isodose line), the eMC underestimated dose by a factor of 3 for high energy electron beams due to lack of consideration of bremsstrahlung emitted laterally that was not accounted by eMC in the low dose region outside the field. In the penumbra (20%-80% isodose line), the eMC overestimated dose (40%) for high energy 20 MeV electrons compared to the measured dose with small diodes in the high gradient dose region. This was mainly due to lack of consideration of volume averaging of the ion chamber used in beam data commissioning which was input to the eMC dose calculation algorithm. Large uncertainties in the CT numbers (25%) resulted from the image artifacts in the CT images of the MapCheck2 phantom due to metal artifacts. The eMC algorithm used the electron and material densities extracted from the CT numbers which resulted large dosimetric uncertainties (10%) in the material densities and corresponding stopping power ratios. The dose calculations with eMC are associated with large uncertainties particularly in penumbra and umbra regions and around heterogeneities which affect the low dose level that cover nearby normal tissue or critical structures.

Comparison of calculation algorithms to predict the IQM detector response for various modulation degrees of VMAT treatment plans on linear accelerator equipped with the HD120 MLC.

PURPOSE: The calculation model for the integral quality monitor (IQM) system does not take into account the characteristics of the HD120 multileaf collimator (MLC), which some Varian accelerators are equipped with. Some treatment plans prepared with this collimator are characterized by a high level of modulation. The aim of the work was to prepare a model for that collimator and to determine the influence of modulation on results of the verification carried out with the use of IQM system. METHODS: The short and long stabilities of the IQM detector response were verified by measuring the signal for a 6 MV flattening filter-free (FFF) beam with the static field of 10 × 10 cm2 size. The obtained results were compared with the measurements performed with the PTW Farmer chamber. Next, the signals for 35 static square fields 4 × 4 cm2 , covering the whole field 38 × 20 cm2 , were measured with the IQM. Based on the results of these measurements, the original calculation model has been changed in order to achieve the smallest differences between calculations and measurements. While tuning the model, the characteristics of the HD120 MLC were included. Measurements were performed for 30 clinical plans (86 arcs) prepared with 6 MV FFF beams. Among those 30 plans, there were were multitarget plans with single isocenter. For each plan, the modulation complexity score (MCS) was calculated. The measurement results were compared with the calculation results performed with the original and authors' calculation model. RESULTS: Very good stability of the short and long stabilities of the IQM detector response was obtained. Measurements performed for 35 static fields revealed that for the manufacturer's and for the authors' models, the deviation exceeded 3% for 12 and five of the 35 static fields, respectively. The differences for the manufacturer's and authors' algorithms were in the range of ±2% for the 15 and 26 of the fields, respectively. For original and the authors' models, the differences between measured and calculated signals (starting with the segment number 40) were within the range of ±3.5% for 87.6% and 96.7% of all arcs for the respective models. For both models, the dependence of the compliance of measurements and calculations on the MCS was observed. For most of the very modulated arcs, the measured signal was at least 3% lower than the calculated one. The largest differences between measurements and calculations were obtained for single-isocenter multitarget plans. CONCLUSIONS: The signal predicted by an algorithm taking into account the real geometry of the collimating system of the Edge accelerator (equipped with the HD120 MLC) made it possible to obtain greater consistency between the measurements and calculations. We characterized the dependence between the MCS of each arc and the compliance of the measurements and calculations. Much worse results were obtained for single-isocenter multitarget plans.

Dose calculation validation of a convolution algorithm in a solid water phantom.

PURPOSE: The dose calculated using a convolution algorithm should be validated in a simple homogeneous water-equivalent phantom before clinical use. The dose calculation accuracy within a solid water phantom was investigated. METHODS: The specific Gamma knife design requires a dose rate calibration within a spherical solid water phantom. The TMR10 algorithm, which approximates the phantom material as liquid water, correctly computes the absolute dose in water. The convolution algorithm, which considers electron density miscalculates the dose in water as the phantom Hounsfield units were converted into higher electron density when the original CT calibration curve was used. To address this issue, the electron density of liquid water was affected by modifying the CT calibration curve. The absolute dose calculated using the convolution algorithm was compared with that computed by the TMR10. The measured depth dose profiles were also compared to those computed by the convolution and TMR10 algorithms. A patient treatment was recalculated in the solid-water phantom and the delivery quality assurance was checked. RESULTS: The convolution algorithm and the TMR10 calculate an absolute dose within 1% when using the modified CT calibration curve. The dose depth profile calculated using the convolution algorithms was superimposed on the TMR10 and measured dose profiles when the modified CT calibration curve was applied. The Gamma index was better than 93%. CONCLUSIONS: Dose calculation algorithms, which consider electron density, require a CT calibration curve adapted to the phantom material to correctly compute the dose in water.

The influence of basic plan parameters on calculated small field output factors - A multicenter study.

PURPOSE: The influence of basic plan parameters such as slice thickness, grid resolution, algorithm type and field size on calculated small field output factors (OFs) was evaluated in a multicentric study. METHODS AND MATERIALS: Three computational homogeneous water phantoms with slice thicknesses (ST) 1, 2 and 3 mm were shared among twenty-one centers to calculate OFs for 1x1, 2x2 and 3x3 cm2 field sizes (FSs) (normalized to 10x10 cm2 FS), with their own treatment planning system (TPS) and the energy clinically used for stereotactic body radiation therapy delivery. OFs were calculated for each combination of grid resolution (GR) (1, 2 and 3 mm) and ST and finally compared with the OFs measured for the TPS commissioning. A multivariate analysis was performed to test the effect of basic plan parameters on calculated OFs. RESULTS: A total of 509 data points were collected. Calculated OFs are slightly higher than measured ones. The multivariate analysis showed that Center, GR, algorithm type, and FS are predictive variables of the difference between calculated and measured OFs (p < 0.001). As FS decreases, the spread in the difference between calculated and measured OFs became larger when increasing the GR. Monte Carlo and Analytical Anisotropic Algorithms, presented a dependence on GR (p < 0.01), while Collapsed Cone Convolution and Acuros did not. The effect of the ST was found to be negligible. CONCLUSIONS: Modern TPSs slightly overestimate the calculated small field OFs compared with measured ones. Grid resolution, algorithm, center number and field size influence the calculation of small field OFs.

Validation of the collapsed cone algorithm for HDR liver brachytherapy against Monte Carlo simulations.

PURPOSE: To validate the collapsed cone (CC) algorithm against Monte Carlo (MC) simulations for model-based dose calculations in high-dose-rate (HDR) liver brachytherapy. METHODS AND MATERIALS: Doses for liver brachytherapy treatment plans of 10 cases were retrospectively recalculated with a model-based approach using Monte Carlo n-Particle Code (MCNP) 6 (Dm,m-MC) and Oncentra Brachy ACE (Dm,m-ACE). Tissue segmentation consisted of assigning uniform compositions and mass densities to predefined Hounsfield Unit (HU) thresholds. Resulting doses were compared according to dose volume histogram parameters typical for clinical routine. These included the percentage liver volume receiving 5 Gy (V5Gy) or 10 Gy (V10Gy), the maximum dose to one cubic centimeter (D1cc) of organs at risk, the clinical target volume (CTV) fractions receiving 150% (V150), 100% (V100), 95% (V95) and 90% (V90) of the prescribed dose and the absolute doses to 95% (D95) and 90% (D90) of the CTV volumes. RESULTS: Doses from Oncentra Brachy ACE agreed well with MC simulations. Differences were seen far from the source, in low-density regions and bone structures. Median percentage deviations were 1.1% for the liver V5Gy and 0.4% for the liver V10Gy, with deviations of largest magnitude amounting to 2.2% and 1.0%, respectively. Organs at risk had median deviations ranging from 0.3% to 1.5% for D1cc, with outliers ranging up to 4.6%. CTV volume parameter deviations ranged between -1.5% and 0.5%, dose parameter deviations ranged mostly between -2% and 1%, with two outliers at -4.0% and -3.4% for a small CTV.

Clinically-implementable template plans for multidwell treatments using Leipzig-style applicators in 192Ir surface brachytherapy.

PURPOSE: Multiple dwell positions ("multidwell") within a Leipzig-style applicator can be used to increase dose uniformity and treatment area. Model-based dose calculation algorithms (MBDCAs) are necessary for accurate calculations involving these applicators because of their nonwater equivalency and complex geometry. The purpose of this work was to create template plans from MBDCA calculations and present their dwell times and positions for users of these applicators without access to MBDCAs. METHODS AND MATERIALS: The Leipzig-style solid applicator model within our treatment planning system was used to design template plans. Five template plans, normalized to 0.3 cm depth within a water phantom, were optimized using the treatment planning system MBDCA. Each template plan contained unique dwell positions, times, and active lengths (0.5-1.5 cm). A single-dwell distribution was optimized for comparison. The stem of this applicator stops within the shell; therefore, one template plan contained an intrafraction rotation to determine the largest dose distribution achievable. Effects of imperfect applicator rotation were quantified by inserting rotational offsets and comparing the V100%, D95%, and minimum dose coverage for planning target volumes created from 80%/90% isodose lines. RESULTS: The 90% (80%) isodose line dimensions at 0.3 cm depth for single-dwell increased from 0.94 × 0.97 (1.53 × 1.57) cm2 to 2.09 × 1.24 (2.75 × 1.88) cm2 in the largest template plan. Manually inserted angular offsets up to ±10° for the template plan requiring rotation preserved V100%, D95%, and minimum dose within 2.0%, 1.9%, and 8.0%, respectively. CONCLUSION: A set of template plans was created to provide accessibility to the multidwell methodology, even for users without access to MBDCAs. Each template plan should be reviewed before clinical implementation.

Ipsilateral lung normal tissue complication probability parameters for different dose calculation algorithms in radiotherapy of breast cancer.

PURPOSE: Different dose calculation algorithms (DCAs) predict different dose distributions for the same treatment. Awareness of optimal model parameters is vital for estimating normal tissue complication probability (NTCP) for different algorithms. The aim is to determine the NTCP parameter values for different DCAs in left-sided breast radiotherapy, using the Lyman-Kutcher-Burman (LKB) model. MATERIALS AND METHODS: First, the methodology recommended by International Atomic Energy Agency TEC-DOC 1583 was used to establish the accuracy of dose calculations of different DCAs including: Monte Carlo (MC) and collapsed cone algorithms implemented in Monaco, pencil beam convolution (PBC) and analytical anisotropic algorithm (AAA) implemented in Eclipse, and superposition and Clarkson algorithms implemented in PCRT3D treatment planning systems (TPSs). Then, treatment planning of 15 patients with left-sided breast cancer was performed by the mentioned DCAs and NTCP of the left-lung normal tissue were calculated for each patient individually, using the LKB model. For the PB algorithm, the NTCP parameters were taken from previously published values and new model parameters obtained for each DCA, using the iterative least squares methods. RESULTS: For all cases and DCAs, NTCP computation with the same model parameters resulted in >15% deviation in NTCP values. The new NTCP model parameters were classified according to the algorithm type. Thus, the discrepancy of NTCP computations was reduced up to 5% after utilizing adjusted model parameters. CONCLUSIONS: This paper confirms that the NTCP values for a given treatment type are different for the different DCAs. Thus, it is essential to introduce appropriate NTCP parameter values according to DCA adopted in TPS, to obtain a more precise estimation of lung NTCP. Hence, new parameter values, classified according to the DCAs, must be determined before introducing NTCP estimation in clinical practice.

[Analysis of the internationalized highly cited segment of national surgery publications]. / Analiz internatsionalizirovannogo vysokotsitiruemogo segmenta otechestvennykh publikatsii po khirurgii.

The internationalized segment of manuscripts on surgery published in top-rated journals indexed in the Web of Science Core Collection (WoS CC) was analyzed. Presence of national publications in surgical journals of different quartiles was estimated. It was shown that 59.6% of surgical articles of Russian authors, indexed in the WoS CC, were published in high-ranking journals (Q1 and Q2). The number of such articles has been increased by more than 3 times (from 23 in 2010 to 86 in 2018) on the background of no growth in the fourth quartile. We formed a rating of national educational and academic organizations, which have made the greatest contribution to a portfolio of publications on surgery, indexed in Q1-Q2 quartiles of WoS CC, as well as in the segment of highly cited publications.