Simplified method for determining dose to a non-water phantom through the use of ND,w and IAEA TRS 483 for the GammaPod

PurposeGammaPod, a breast stereotactic radiosurgery device, utilizes 25 rotating Co-60 sources to deliver highly conformal dose distributions. The GammaPod system requires that reference dosimetry be performed in a specific vendor-supplied poly-methylmethacrylate (PMMA) phantom. The nonstandard nature of GammaPod dosimetry, in both the phantom material and machine-specific reference (msr), prohibits use of the American Association of Physicists in Medicine Task Group 51 (TG-51) protocol. This study proposes a practical method using TRS 483 to make the reference dosimetry procedure simpler and to reduce overall uncertainties.MethodsThe dose to PMMA (D_PMMA) is determined under msr conditions using TRS 483 with an Exradin A1SL chamber placed in a PMMA phantom. The conversion factor, which converts from the dose-to-water (D_w) in broad-beam Co-60 reference geometry to D_PMMA in the msr small field Co-60 (Q_msr) geometry, is derived using the Monte Carlo simulations and procedure described in TRS 483.ResultsThe new conversion factor value for an Exradin A1SL chamber is 0.974. When combined with N_D,w, D_PMMA differs by 0.5% from the TG-21/N_x method and 0.2% from the IROC values. Uncertainty decreased from 2.2% to 1.6%.ConclusionWe successfully implemented TRS 483 reference dosimetry protocols utilizing N_D,w for the GammaPod in the PMMA phantom. These results show not only agreement between measurements performed with the previously published method and independent thermoluminescent dosimetry measurements but also reductions in uncertainty. This also provides readers with a pathway to develop their own IAEA TRS 483 factor for any new small field machine that may be developed.     

A new cone-beam computed tomography dosimetry method providing optimal measurement positions: A Monte Carlo study

PurposeThe conventional weighted computed tomography dose index (CTDI_w) may not be suitable for cone-beam computed tomography (CBCT) dosimetry because a cross-sectional dose distribution is angularly inhomogeneous owing to partial angle irradiations. This study was conducted to develop a new dose metric (f(0)^CB_w) for CBCT dosimetry to determine a more accurate average dose in the central cross-sectional plane of a cylindrical phantom using Monte Carlo simulations.MethodsFirst, cross-sectional dose distributions of cylindrical polymethyl methacrylate phantoms over a wide range of phantom diameters (8–40 cm) were calculated for various CBCT scan protocols. Then, by obtaining linear least-squares fits of the full datasets of the cross-sectional dose distributions, the optimal radial positions, which represented measurement positions for the average phantom dose, were determined. Finally, the f(0)^CB_w method was developed by averaging point doses at the optimal radial positions of the phantoms. To demonstrate its validity, the relative differences between the average doses and each dose index value were estimated for the devised f(0)^CB_w, conventional CTDI_w, and Haba’s CTDI_w methods, respectively.ResultsThe relative differences between the average doses and each dose index value were within 4.1%, 16.7%, and 11.9% for the devised, conventional CTDI_w, and Haba’s CTDI_w methods, respectively.ConclusionsThe devised f(0)^CB_w value was calculated by averaging four “point doses” at 90° intervals and the optimal radial positions of the cylindrical phantom. The devised method can estimate the average dose more accurately than the previously developed CTDI_w methods for CBCT dosimetry.     

Determination of backscatter factors based on the quality index for diagnostic kilovoltage x-ray beams

PurposeThis study aims to investigate the relationship between backscatter factors and Al-half-value-layers (Al-HVL) by making the quality index (QI) a parameter for diagnostic kilovoltage x-ray beams.MethodsBackscatter factors, B_w, for x-ray fluence spectra were calculated from the weighted average of B_w for monoenergetic photons of between 8 and 140 keV with field sizes of 10 cm × 10 cm to 40 cm × 40 cm. The value of B_w for monoenergetic photons was calculated from the ratio of the water kerma at the surface of a water phantom and that at the same point free-in-air using the EGSnrc/cavity code. The weighted averaged backscatter factors were validated by comparing them with those of direct Monte Carlo calculations for the x-ray fluence spectra. The B_w for the x-ray fluence spectra were classified by a QI of 0.35, 0.4, 0.5, 0.6, and 0.7 specified by the ratio of the effective energy and maximum energy. The relationship between B_w and Al-HVL was evaluated for the given QI values. The x-ray fluence spectra were generated for tube voltages of 40–140 kVp with Al-HVLs of 0.5–13.2 mm using the SpekCalc program.ResultsThe weighted averaged backscatter factors for x-ray fluence spectra agreed within 0.7% with those of the direct Monte Carlo calculations. The backscatter factors were represented by the fitting curves of R² > 0.99 with Al-HVL for the given QI values.ConclusionsIt is possible to obtain B_w more accurately by using QI specified by the measured Al-HVL.     

Measuring the dose in bone for spine stereotactic body radiotherapy

PurposeCurrent quality assurance of radiotherapy involving bony regions generally utilises homogeneous phantoms and dose calculations, ignoring the challenges of heterogeneities with dosimetry problems likely occurring around bone. Anthropomorphic phantoms with synthetic bony materials enable realistic end-to-end testing in clinical scenarios. This work reports on measurements and calculated corrections required to directly report dose in bony materials in the context of comprehensive end-to-end dosimetry audit measurements (63 plans, 6 planning systems).Materials and methodsRadiochromic film and microDiamond measurements were performed in an anthropomorphic spine phantom containing bone equivalent materials. Medium dependent correction factors, k_med, were established using 6 MV and 10 MV Linear Accelerator Monte Carlo simulations to account for the detectors being calibrated in water, but measuring in regions of bony material. Both cortical and trabecular bony material were investigated for verification of dose calculations in dose-to-medium (D_m,m) and dose-to-water (D_w,w) scenarios.ResultsFor D_m,m calculations, modelled correction factors for cortical and trabecular bone in film measurements, and for trabecular bone in microDiamond measurements were 0.875(±0.1%), 0.953(±0.3%) and 0.962(±0.4%), respectively. For D_w,w calculations, the corrections were 0.920(±0.1%), 0.982(±0.3%) and 0.993(±0.4%), respectively. In the audit, application of the correction factors improves the mean agreement between treatment plans and measured microDiamond dose from −2.4%(±3.9%) to 0.4%(±3.7%).ConclusionMonte Carlo simulations provide a method for correcting the dose measured in bony materials allowing more accurate comparison with treatment planning system doses. In verification measurements, algorithm specific correction factors should be applied to account for variations in bony material for calculations based on D_m,m and D_w,w.     

Intercomparison of dose enhancement ratio and secondary electron spectra for gold nanoparticles irradiated by X-rays calculated using multiple Monte Carlo simulation codes

PurposeTargeted radiation therapy has seen an increased interest in the past decade. In vitro and in vivo experiments showed enhanced radiation doses due to gold nanoparticles (GNPs) to tumors in mice and demonstrated a high potential for clinical application. However, finding a functionalized molecular formulation for actively targeting GNPs in tumor cells is challenging. Furthermore, the enhanced energy deposition by secondary electrons around GNPs, particularly by short-ranged Auger electrons is difficult to measure. Computational models, such as Monte Carlo (MC) radiation transport codes, have been used to estimate the physical quantities and effects of GNPs. However, as these codes differ from one to another, the reliability of physical and dosimetric quantities needs to be established at cellular and molecular levels, so that the subsequent biological effects can be assessed quantitatively.MethodsIn this work, irradiation of single GNPs of 50 nm and 100 nm diameter by X-ray spectra generated by 50 and 100 peak kilovoltages was simulated for a defined geometry setup, by applying multiple MC codes in the EURADOS framework.ResultsThe mean dose enhancement ratio of the first 10 nm-thick water shell around a 100 nm GNP ranges from 400 for 100 kVp X-rays to 600 for 50 kVp X-rays with large uncertainty factors up to 2.3.ConclusionsIt is concluded that the absolute dose enhancement effects have large uncertainties and need an inter-code intercomparison for a high quality assurance; relative properties may be a better measure until more experimental data is available to constrain the models.     

Effective dose in percutaneous transhepatic biliary drainage examination using PCXMC2.0 and MCNP5 Monte Carlo codes

ObjectivesTo estimate the organ equivalent doses and the effective doses (E) in patient undergoing percutaneous transhepatic biliary drainage (PTBD) examinations, using the MCNP5 and PCXMC2 Monte Carlo-based codes.MethodsThe purpose of this study is to estimate the organ doses to patients undergoing PTBD examinations by clinical measurements and Monte Carlo simulation. Dose area products (DAP) values were assessed during examination of 43 patients undergoing PTBD examination separated into groups based on the gender and the dimensions and location of the beam.ResultsMonte Carlo simulation of photon transport in male and female mathematical phantoms was applied using the MCNP5 and PCXMC2 codes in order to estimate equivalent organ doses. Regarding the PTBD examination the organ receiving the maximum radiation dose was the lumbar spine. The mean calculated H_T for the lumbar spine using the MCNP5 and PCXMC2 methods respectively, was 117.25 mSv and 131.7 mSv, in males. The corresponding doses were 139.45 mSv and 157.1 mSv respectively in females. The H_T values for organs receiving considerable amounts of radiation during PTBD examinations were varied between 0.16% and 73.2% for the male group and between 1.10% and 77.6% for the female group. E in females and males using MCNP5 and PCXMC2.0 was 5.88 mSv and 6.77 mSv, and 4.93 mSv and 5.60 mSv.ConclusionThe doses remain high compared to other invasive operations in interventional radiology. There is a reasonable good coincidence between the MCNP5 and PCXMC2.0 calculation for most of the organs.     

Implementation of new physics models for low energy electrons in liquid water in Geant4-DNA

A new alternative set of elastic and inelastic cross sections has been added to the very low energy extension of the Geant4 Monte Carlo simulation toolkit, Geant4-DNA, for the simulation of electron interactions in liquid water. These cross sections have been obtained from the CPA100 Monte Carlo track structure code, which has been a reference in the microdosimetry community for many years. They are compared to the default Geant4-DNA cross sections and show better agreement with published data.In order to verify the correct implementation of the CPA100 cross section models in Geant4-DNA, simulations of the number of interactions and ranges were performed using Geant4-DNA with this new set of models, and the results were compared with corresponding results from the original CPA100 code. Good agreement is observed between the implementations, with relative differences lower than 1% regardless of the incident electron energy.Useful quantities related to the deposited energy at the scale of the cell or the organ of interest for internal dosimetry, like dose point kernels, are also calculated using these new physics models. They are compared with results obtained using the well-known Penelope Monte Carlo code.     

Low-energy electron penetration range in liquid water

Monte Carlo simulations of electron tracks in liquid water are performed to calculate the energy dependence of the electron penetration range at initial electron energies between 0.2 eV and 150 keV, including the subexcitation electron region (<7.3 eV). Our calculated electron penetration distances are compared with available experimental data and earlier calculations as well as with the results of simulations using newly reported amorphous ice electron scattering cross sections in the range approximately 1-100 eV.     

Fast Monte Carlo codes for occupational dosimetry in interventional radiology

PurposeInterventional radiology techniques cause radiation exposure both to patient and personnel. The radiation dose to the operator is usually measured with dosimeters located at specific points above or below the lead aprons. The aim of this study is to develop and validate two fast Monte Carlo (MC) codes for radiation transport in order to improve the assessment of individual doses in interventional radiology. The proposed methodology reduces the number of required dosemeters and provides immediate dose results.MethodsTwo fast MC simulation codes, PENELOPE/penEasyIR and MCGPU-IR, have been developed. Both codes have been validated by comparing fast MC calculations with the multipurpose PENELOPE MC code and with measurements during a realistic interventional procedure.ResultsThe new codes were tested with a computation time of about 120 s to estimate operator doses while a standard simulation needs several days to obtain similar uncertainties. When compared with the standard calculation in simple set-ups, MCGPU-IR tends to underestimate doses (up to 5%), while PENELOPE/penEasyIR overestimates them (up to 18%). When comparing both fast MC codes with experimental values in realistic set-ups, differences are within 25%. These differences are within accepted uncertainties in individual monitoring.ConclusionThe study highlights the fact that computational dosimetry based on the use of fast MC codes can provide good estimates of the personal dose equivalent and overcome some of the limitations of occupational monitoring in interventional radiology. Notably, MCGPU-IR calculates both organ doses and effective dose, providing a better estimate of radiation risk.     

Application of Gafchromic EBT2 film for intraoperative radiation therapy quality assurance

PurposeIntraoperative radiation therapy (IORT) using electron beam is commonly done by mobile dedicated linacs that have a variable range of electron energies. This paper focuses on the evaluation of the EBT2 film response in the green and red colour channels for IORT quality assurance (QA).MethodsThe calibration of the EBT2 films was done in two ranges; 0–8 Gy for machine QA by red channel and 8–24 Gy for patient-specific QA by green channel analysis. Irradiation of calibration films and relative dosimetries were performed in a water phantom. To evaluate the accuracy of the film response in relative dosimetry, gamma analysis was used to compare the results of the Monte Carlo simulation and ionometric dosimetry. Ten patients with early stage breast cancer were selected for in-vivo dosimetry using the green channel of the EBT2 film.ResultsThe calibration curves were obtained by linear fitting of the green channel and a third-order polynomial function in the red channel (R2 = 0.99). The total dose uncertainty was up to 4.2% and 4.7% for the red and green channels, respectively. There was a good agreement between the relative dosimetries of films by the red channel, Monte Carlo simulations and ionometric values. The mean dose difference of the in-vivo dosimetry by green channel of this film and the expected values was about 1.98% ± 0.75.ConclusionThe results of this study showed that EBT2 film can be considered as an appropriate tool for machine and patient-specific QA in IORT.     

Using matrix summation method for three dimensional dose calculation in brachytherapy

AimThe purpose of this study is to calculate radiation dose around a brachytherapy source in a water phantom for different seed locations or rotation the sources by the matrix summation method.BackgroundMonte Carlo based codes like MCNP are widely used for performing radiation transport calculations and dose evaluation in brachytherapy. But for complicated situations, like using more than one source, moving or rotating the source, the routine Monte Carlo method for dose calculation needs a long time running.Materials and methodsThe MCNPX code has been used to calculate radiation dose around a ¹⁹²Ir brachytherapy source and saved in a 3D matrix. Then, we used this matrix to evaluate the absorbed dose in any point due to some sources or a source which shifted or rotated in some places by the matrix summation method.ResultsThree dimensional (3D) dose results and isodose curves were presented for ¹⁹²Ir source in a water cube phantom shifted for 10 steps and rotated for 45 and 90° based on the matrix summation method. Also, we applied this method for some arrays of sources.ConclusionThe matrix summation method can be used for 3D dose calculations for any brachytherapy source which has moved or rotated. This simple method is very fast compared to routine Monte Carlo based methods. In addition, it can be applied for dose optimization study.     

Comparative dose levels between CT-scanner and slot-scanning device (EOS system) in pregnant women pelvimetry

PurposeTo estimate fetal absorbed doses for pregnant women pelvimetry, a comparative study between EOS imaging system and low-dose spiral CT-scanner was carried out. For this purpose three different studies were investigated: in vivo, in vitro and Monte Carlo calculations.MethodsIn vivo dosimetry was performed, using OSL NanoDot dosimeters, to determine the dose to the skin of twenty pregnant women. In vitro studies were established by using a cubic phantom of water, in order to estimate the out of field doses. In the latter study, OSLDs were placed at depths corresponding to the lowest, average and highest position of the uterus. Monte Carlo calculations of effective doses to high radio-sensitive organs were established, using PCXMC and CTExpo software suites for EOS imaging system and CT-scanner, respectively.ResultsThe EOS imaging system reduces radiation exposure 4 to 8 times compared to the CT-scanner. The entrance skin doses were 74% (p-values <0.01) higher with the CT-scanner than with the EOS system. In the out of field region, the measured doses of the EOS system were reduced by 80% (p-values <0.02).Monte Carlo calculations confirmed that effective doses to organs are less accentuated for EOS than for CT pelvimetry.ConclusionsThe EOS system is less irradiating than the CT exam. The out-of-field dose which is significant, is lower in the EOS than in the CT-scanner and could be reduced even further by optimizing the time used for image acquisition.     

Structure of Mixed Solvent Electrolyte Solutions via Gibbs Ensemble Monte Carlo Simulation

Abstract

Previously reported Gibbs ensemble Monte Carlo simulations of vapor-liquid equilibrium in methanol-water and methanol-water-NaCl mixtures are extended to permit study of the microscopic structure of the liquid phases of these systems. The salt effect in a prototypical mixed solvent electrolyte solution (water-methanol-NaCl) is microscopically interpreted in terms of the structural changes undergone by the solvation shells of the ions in the liquid phase of water-methanol-NaCl systems in vapor-liquid equilibrium at constant pressure.     

Comparisons of calculations with PARTRAC and NOREC: transport of electrons in liquid water

Monte Carlo computer models that simulate the detailed, event-by-event transport of electrons in liquid water are valuable for the interpretation and understanding of findings in radiation chemistry and radiation biology. Because of the paucity of experimental data, such efforts must rely on theoretical principles and considerable judgment in their development. Experimental verification of numerical input is possible to only a limited extent. Indirect support for model validity can be gained from a comparison of details between two independently developed computer codes as well as the observable results calculated with them. In this study, we compare the transport properties of electrons in liquid water using two such models, PARTRAC and NOREC. Both use interaction cross sections based on plane-wave Born approximations and a numerical parameterization of the complex dielectric response function for the liquid. The models are described and compared, and their similarities and differences are highlighted. Recent developments in the field are discussed and taken into account. The calculated stopping powers, W values, and slab penetration characteristics are in good agreement with one another and with other independent sources.     

Monte Carlo simulations and measurement of DNA damage from x-ray-triggered Auger cascades in iododeoxyuridine (IUdR)

We investigated the DNA damage from Auger electrons emitted from incorporated stable iodine (¹²⁷I), following photoelectric absorption of external x-rays. The effectiveness of the Auger electrons in producing DNA double-strand breaks (DSB) was determined theoretically, using Monte Carlo simulations of the radiation physics and chemistry, and was shown to be in reasonable agreement with DNA damage measured using the comet assay. The DSB yields were measured in CHO cells for ⁶⁰Co (as a non-Auger-promoting radiation) and for tungsten-filtered 100 kVp x-rays capable of producing Auger electron emission. The theoretical study showed that on average, 2.5 Auger electrons were emitted for N-shell orbital vacancies and up to 10 Auger electrons were emitted from L₁-shell vacancies. These Auger bursts produced approximately 0.03 DSB per N-shell vacancy and 0.3 DSB per K-shell or L-shell vacancy. The calculated yield of DSB from Auger cascades per unit dose (1 Gy) in water was approximately 1.7 for tungsten-filtered 100 kVp x-rays, assuming 20% IUdR substitution of thymidine. The comet assay yielded an experimental value of 3.6±1.6 per 1 Gy for the same conditions. The Monte Carlo simulations also demonstrated a high complexity of DSB produced by Auger cascades with virtually all DSB from inner shell orbitals (i.e. K, L shells) accompanied by compounded strand breakage and base damage, indicating a difficult lesion to repair. This finding agrees well with comet assay results of DNA repair, where an increase in the DSB yield in IUdR-sensitized cells was shown to persist after a time of 24 h. We conclude that Auger cascades in iodine produce a modest increase in the number of initial strand breaks of the order of 10% but the complex nature of these DSB makes them very difficult to repair or potentially prone to misrepair. The accentuated DNA damage may have major consequences for cell survival and may be exploitable in kilovoltage photon activation therapy (PAT) of tumors sensitized with iodine. Received: 23 October 2000 / Accepted: 26 March 2001     

Comparison and assessment of electron cross sections for Monte Carlo track structure codes.

The purpose of this study was to make an intercomparison and assessment of cross sections for electrons in water used in electron track structure codes. This study is intended to shed light on the extent to which the differences between the input data and physical and chemical assumptions influence the outcome in biophysical modeling of radiation effects. Ionization cross sections and spectra of secondary electrons were calculated by various theories. The analyses were carried out for water vapor cross sections, as these are more abundant and readily available. All suitable published experimental total ionization cross sections were fitted by an appropriate function and used for generation of electron tracks. Three sets of compiled data were used for comparison of total excitation cross sections and mean excitation energy. The tracks generated by a Monte Carlo track code, using various combinations of cross sections, were compared in terms of radial distributions of interactions and point kernels. The spectrum of secondary electrons emitted by the ionization process was found to be the factor that has the most influence on these quantities. A different set of cross sections for excitation and elastic scattering did not affect the electron track structure as much as did ionization cross sections. It is concluded that all codes, using different cross sections and in different phase, currently used for biophysical modeling exhibit close similarities for energy deposition in larger size targets while appreciable differences are observed in B-DNA-size targets. We recommend fitted functions to all available suitable experimental data for the total ionization and elastic cross sections. We conclude that most codes produce tracks in reasonable agreement with the macroscopic quantities such as total stopping power and total yield of strand breaks. However, we predict differences in frequencies of clustering in tracks from the different models.     

Comparison of dosimetric characteristics of physical wedge and enhanced dynamic wedge in inhomogeneous medium using Monte Carlo simulations

BackgroundWidely used physical wedges in clinical radiotherapy lead to beam intensity attenuation as well as the beam hardening effect, which must be considered. Dynamic wedges devised to overcome the physical wedges (PWs) problems result in dosimetry complications due to jaw movement while the beam is on. This study was aimed to investigate the usability of physical wedge data instead of enhanced dynamic wedge due to the enhanced dynamic wedge (EDW) dosimetry measurement hardships of Varian 2100CD in inhomogeneous phantom by Monte Carlo code as a reliable method in radiation dosimetry.Materials and methodsA PW and EDW-equipped-linac head was simulated using BEAMnrc code. DOSXYZnrc was used for three-dimensional dosimetry calculation in the CIRS phantom.ResultsBased on the isodose curves, EDW generated a less scattered as well as lower penumbra width compared to the PW. The depth dose variations of PWs and EDWs were more in soft tissue than the lung tissue. Beam profiles of PW and EDW indicated good coincidence in all points, except for the heel area.ConclusionResults demonstrated that it is possible to apply PW data instead of EDW due to the dosimetry and commissioning hardships caused by EDW in inhomogeneous media.     

Monte Carlo Dosimetry of the 60Co BEBIG High Dose Rate for Brachytherapy

Introduction

The use of high-dose-rate brachytherapy is currently a widespread practice worldwide. The most common isotope source is ¹⁹²Ir, but ⁶⁰Co is also becoming available for HDR. One of main advantages of ⁶⁰Co compared to ¹⁹²Ir is the economic and practical benefit because of its longer half-live, which is 5.27 years. Recently, Eckert & Ziegler BEBIG, Germany, introduced a new afterloading brachytherapy machine (MultiSource^®); it has the option to use either the ⁶⁰Co or ¹⁹²Ir HDR source. The source for the Monte Carlo calculations is the new ⁶⁰Co source (model Co0.A86), which is referred to as the new BEBIG ⁶⁰Co HDR source and is a modified version of the ⁶⁰Co source (model GK60M21), which is also from BEBIG.

Objective and Methods

The purpose of this work is to obtain the dosimetry parameters in accordance with the AAPM TG-43U1 formalism with Monte Carlo calculations regarding the BEBIG ⁶⁰Co high-dose-rate brachytherapy to investigate the required treatment-planning parameters. The geometric design and material details of the source was provided by the manufacturer and was used to define the Monte Carlo geometry. To validate the source geometry, a few dosimetry parameters had to be calculated according to the AAPM TG-43U1 formalism. The dosimetry studies included the calculation of the air kerma strength S _k, collision kerma in water along the transverse axis with an unbounded phantom, dose rate constant and radial dose function. The Monte Carlo code system that was used was EGSnrc with a new cavity code, which is a part of EGS++ that allows calculating the radial dose function around the source. The spectrum to simulate ⁶⁰Co was composed of two photon energies, 1.17 and 1.33 MeV. Only the gamma part of the spectrum was used; the contribution of the electrons to the dose is negligible because of the full absorption by the stainless-steel wall around the metallic ⁶⁰Co. The XCOM photon cross-section library was used in subsequent simulations, and the photoelectric effect, pair production, Rayleigh scattering and bound Compton scattering were included in the simulation. Variance reduction techniques were used to speed up the calculation and to considerably reduce the computer time. The cut-off energy was 10 keV for electrons and photons. To obtain the dose rate distributions of the source in an unbounded liquid water phantom, the source was immersed at the center of a cube phantom of 100 cm³. The liquid water density was 0.998 g/cm³, and photon histories of up to 10¹⁰ were used to obtain the results with a standard deviation of less than 0.5% (k = 1). The obtained dose rate constant for the BEBIG ⁶⁰Co source was 1.108±0.001 cGyh^-1U^-1, which is consistent with the values in the literature. The radial dose functions were compared with the values of the consensus data set in the literature, and they are consistent with the published data for this energy range.     

Validation of the MC-GPU Monte Carlo code against the PENELOPE/penEasy code system and benchmarking against experimental conditions for typical radiation qualities and setups in interventional radiology and cardiology

IntroductionInterventional procedures are associated with potentially high radiation doses to the skin. The 2013/59/EURATOM Directive establishes that the equipment used for interventional radiology must have a device or a feature informing the practitioner of relevant parameters for assessing patient dose at the end of the procedure. Monte Carlo codes of radiation transport are considered to be one of the most reliable tools available to assess doses. However, they are usually too time consuming for use in clinical practice. This work presents the validation of the fast Monte Carlo code MC-GPU for application in interventional radiology.MethodologiesMC-GPU calculations were compared against the well-validated Monte Carlo simulation code PENELOPE/penEasy by simulating the organ dose distribution in a voxelized anthropomorphic phantom. In a second phase, the code was compared against thermoluminescent measurements performed on slab phantoms, both in a calibration laboratory and at a hospital.ResultsThe results obtained from the two simulation codes show very good agreement, differences in the output were within 1%, whereas the calculation time on the MC-GPU was 2500 times shorter. Comparison with measurements is of the order of 10%, within the associated uncertainty.ConclusionsIt has been verified that MC-GPU provides good estimates of the dose when compared to PENELOPE program. It is also shown that it presents very good performance when assessing organ doses in very short times, less than one minute, in real clinical set-ups. Future steps would be to simulate complex procedures with several projections.     

The Detailed Balance Energy-scaled Displacement Monte Carlo Algorithm

Abstract

The Detailed Balance Energy-scaled Displacement Monte Carlo method that stems from the previously published Energy Scaled Displacement Monte Carlo method is presented. The results of tests performed on a dense Lennard-Jones liquid and on two particles in one dimension are reported.