Cancer risk estimation in Digital Breast Tomosynthesis using GEANT4 Monte Carlo simulations and voxel phantoms

The aim of this work was to estimate the risk of radiation induced cancer following the Portuguese breast screening recommendations for Digital Mammography (DM) when applied to Digital Breast Tomosynthesis (DBT) and to evaluate how the risk to induce cancer could influence the energy used in breast diagnostic exams. The organ doses were calculated by Monte Carlo simulations using a female voxel phantom and considering the acquisition of 25 projection images. Single organ cancer incidence risks were calculated in order to assess the total effective radiation induced cancer risk. The screening strategy techniques considered were: DBT in Cranio-Caudal (CC) view and two-view DM (CC and Mediolateral Oblique (MLO)).The risk of cancer incidence following the Portuguese screening guidelines (screening every two years in the age range of 50–80 years) was calculated by assuming a single CC DBT acquisition view as standalone screening strategy and compared with two-view DM. The difference in the total effective risk between DBT and DM is quite low. Nevertheless in DBT an increase of risk for the lung is observed with respect to DM. The lung is also the organ that is mainly affected when non-optimal beam energy (in terms of image quality and absorbed dose) is used instead of an optimal one. The use of non-optimal energies could increase the risk of lung cancer incidence by a factor of about 2.     

Homogeneous vs. patient specific breast models for Monte Carlo evaluation of mean glandular dose in mammography

PurposeTo compare, via Monte Carlo simulations, homogeneous and non-homogenous breast models adopted for mean glandular dose (MGD) estimates in mammography vs. patient specific digital breast phantoms.MethodsWe developed a GEANT4 Monte Carlo code simulating four homogenous cylindrical breast models featured as follows: (1) semi-cylindrical section enveloped in a 5-mm adipose layer; (2) semi-elliptical section with a 4-mm thick skin; (3) semi-cylindrical section with a 1.45-mm skin layer; (4) semi-cylindrical section in a 1.45-mm skin layer and 2-mm subcutaneous adipose layer. Twenty patient specific digital breast phantoms produced from a dedicated CT scanner were assumed as reference in the comparison. We simulated two spectra produced from two anode/filter combinations. An additional digital breast phantom was produced via BreastSimulator software.ResultsWith reference to the results for patient-specific breast phantoms and for W/Al spectra, models #1 and #3 showed higher MGD values by about 1% (ranges [–33%; +28%] and [−31%; +30%], respectively), while for model #4 it was 2% lower (range [−34%; +26%]) and for model #2 –11% (range [−39%; +14%]), on average. On the other hand, for W/Rh spectra, models #1 and #4 showed lower MGD values by 2% and 1%, while for model #2 and #3 it was 14% and 8% lower, respectively (ranges [−43%; +13%] and [−41%; +21%]). The simulation with the digital breast phantom produced with BreastSimulator showed a MGD overestimation of +33%.ConclusionsThe homogeneous breast models led to maximum MGD underestimation and overestimation of 43% and 28%, respectively, when compared to patient specific breast phantoms derived from clinical CT scans.     

Validation of 4D Monte Carlo dose calculations using a programmable deformable lung phantom

PurposeTo validate the accuracy of 4D Monte Carlo (4DMC) simulations to calculate dose deliveries to a deforming anatomy in the presence of realistic respiratory motion traces. A previously developed deformable lung phantom comprising an elastic tumor was modified to enable programming of arbitrary motion profiles. 4D simulations of the dose delivered to the phantom were compared with the measurements.MethodsThe deformable lung phantom moving with irregular breathing patterns was irradiated using static and VMAT beam deliveries. Using the RADPOS 4D dosimetry system, point doses were measured inside and outside the tumor. Dose profiles were acquired using films along the motion path of the tumor (S-I). In addition to dose measurements, RADPOS was used to record the motion of the tumor during dose deliveries. Dose measurements were then compared against 4DMC simulations with EGSnrc/4DdefDOSXYZnrc using the recorded tumor motion.ResultsThe agreements between dose profiles from measurements and simulations were determined to be within 2%/2 mm. Point dose agreements were within 2σ of experimental and/or positional/dose reading uncertainties. 4DMC simulations were shown to accurately predict the sensitivity of delivered dose to the starting phase of breathing motions. We have demonstrated that our 4DMC method, combined with RADPOS, can accurately simulate realistic dose deliveries to a deforming anatomy moving with realistic breathing traces. This 4DMC tool has the potential to be used as a quality assurance tool to verify treatments involving respiratory motion. Adaptive treatment delivery is another area that may benefit from the potential of this 4DMC tool.     

Skin models and their impact on mean glandular dose in mammography

PurposeTo quantify the influence of different skin models on mammographic breast dosimetry, based on dosimetric protocols and recent breast skin thickness findings.MethodsBy using an adapted PENELOPE (v. 2014) + PenEasy (v. 2015) Monte Carlo (MC) code, simulations were performed in order to obtain the mean glandular dose (MGD), the normalized MGD by incident air Kerma (DgN), and the glandular depth dose (GDD(z)). The geometry was based on a cranio-caudal mammographic examination. Monoenergetic and polyenergetic beams were implemented, for a breast thickness from 2 cm to 9 cm, with different compositions. Seven skin models were used: a 5 mm adipose layer; a skin layer ranging from 5 mm to 1.45 mm, a 1.45 mm skin thickness with a subcutaneous adipose layer of 2 mm and 3.55 mm.ResultsThe differences, for monoenergetic beams, are higher (up to 200%) for lower energies (8 keV), thicker and low glandular content breasts, decreasing to less than 5% at 40 keV. Without a skin layer, the differences reach a maximum of 1240%. The relative difference in DgN values for 1.45 mm skin and 5 mm adipose layers and polyenergetic beams varies from −14% to 12%.ConclusionsThe implemented MC code is suitable for mammography dosimetry calculations. The skin models have major impacts on MGD values, and the results complement previous literature findings. The current protocols should be updated to include a more realistic skin model, which provides a reliable breast dose estimation.     

Mobile shielding evaluation on the fetal dose during a breast radiotherapy using Monte Carlo simulation

PurposeEvaluation of the out-of-field dose is an important aspect in radiotherapy. Due to the fetus radiosensitivity, this evaluation becomes even more conclusive when the patient is pregnant. In this work, a linear accelerator Varian Clinac 2100c operating at 6 MV, a pregnant anthropomorphic phantom (Maria), and different shields added above the abdominal region of the phantom were used for the analysis based on MCNPX. Methods: The simulations were performed for the medial and lateral projections, using either an open field collimation (10×16 cm²) or a multileaf collimator. The added shields (M1 and M2) were designed based on models proposed by Stovall et al. [1], intending to reduce the deposited dose on the fetus and related structures. Results: The presence of the shields showed to be effective in reducing the doses on the fetus, amniotic sac, and placenta, for example. A reduction of about 43% was found in the dose on the fetus when M2 was added, using the open field collimation, in comparison with the situation with no shield, being the lateral projection the main responsible for the dose. The use of MLC significatively reduced the doses in different structures, including on the fetus and amniotic sac, for example, in comparison to the open field situation. A slight increment on the dose in organs such as the eyes, thyroid and brain was found in both collimation systems, due to the presence of the shields. The contribution of the leakage radiation from the tube head of the linear accelerator was found to be in the order of µGy, being reduced by the presence of the M2 shield. Conclusion: Using the shields showed to be an essential feature in order to reduce the dose not only on the fetus, but also in important structures responsible to its development.     

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.     

Monte Carlo dose in a prosthesis phantom based on exact geometry vs streak artefact contaminated CT data as benchmarked against Gafchromic film measurements

PurposeIn radiotherapy, accurate calculation of patient radiation dose is very important for good clinical outcome. In the presence of metallic implants, the dose calculation accuracy could be compromised by metal artefacts generated in computed tomography (CT) images of patients. This study investigates the influence of metal-induced CT artefacts on MC dose calculations in a pelvic prosthesis phantom.MethodsA pelvic phantom containing unilateral Ti prosthesis was CT-scanned and accurate Hounsfield unit (HU) values were assigned to known materials of the phantom as opposed to HU values produced through the artefact CT images of the phantom. Using the DOSXYZnrc MC code, dose calculations were computed in the phantom model constructed from the original CT images containing the artefacts and artefact-free images made from the exact geometry of the phantom with known materials. The dose calculations were benchmarked against Gafchromic EBT3 film measurements using 15 MeV electron and 10 MV photon beams.ResultsThe average deviations between film and MC dose data decreased from 3 ± 2% to 1 ± 1% and from about 6 ± 2% to 3 ± 1% for the artefact and artefact-free phantom models against film data for the electron and photon fields, respectively.ConclusionsFor the Ti prosthesis phantom, the presence of metal-induced CT artefacts could cause dose inaccuracies of about 3%. Construction of an artefact-free phantom model made from the exact geometry of the phantom with known materials to overcome the effect of artefacts is advantageous compared to using CT data directly of which the exact tissue composition is not well-known.     

Two-step validation of a Monte Carlo dosimetry framework for general radiology

The Monte Carlo technique is considered gold standard when it comes to patient-specific dosimetry. Any newly developed Monte Carlo simulation framework, however, has to be carefully calibrated and validated prior to its use. For many researchers this is a tedious work. We propose a two-step validation procedure for our newly built Monte Carlo framework and provide all input data to make it feasible for future related application by the wider community. The validation was at first performed by benchmarking against simulation data available in literature. The American Association of Physicists in Medicine (AAPM) report of task group 195 (case 2) was considered most appropriate for our application. Secondly, the framework was calibrated and validated against experimental measurements for trunk X-ray imaging protocols using a water phantom. The dose results obtained from all simulations and measurements were compared. Our Monte Carlo framework proved to agree with literature data, by showing a maximal difference below 4% to the AAPM report. The mean difference with the water phantom measurements was around 7%. The statistical uncertainty for clinical applications of the dosimetry model is expected to be within 10%. This makes it reliable for clinical dose calculations in general radiology. Input data and the described procedure allow for the validation of other Monte Carlo frameworks.     

Development of a deformable phantom for experimental verification of 4D Monte Carlo simulations in a deforming anatomy

PurposeTo verify the accuracy of 4D Monte Carlo (MC) simulations, using the 4DdefDOSXYZnrc user code, in a deforming anatomy. We developed a tissue-equivalent and reproducible deformable lung phantom and evaluated 4D simulations of delivered dose to the phantom by comparing calculations against measurements.MethodsA novel deformable phantom consisting of flexible foam, emulating lung tissue, inside a Lucite external body was constructed. A removable plug, containing an elastic tumor that can hold film and other dosimeters, was inserted in the phantom. Point dose and position measurements were performed inside and outside the tumor using RADPOS 4D dosimetry system. The phantom was irradiated on an Elekta Infinity linac in both stationary and moving states. The dose delivery was simulated using delivery log files and the phantom motion recorded with RADPOS.ResultsReproducibility of the phantom motion was determined to be within 1 mm. The phantom motion presented realistic features like hysteresis. MC calculations and measurements agreed within 2% at the center of tumor. Outside the tumor agreements were better than 5% which were within the positional/dose reading uncertainties at the measurement points. More than 94% of dose points from MC simulations agreed within 2%/2 mm compared to film measurements.ConclusionThe deformable lung phantom presented realistic and reproducible motion characteristics and its use for verification of 4D dose calculations was demonstrated. Our 4DMC method is capable of accurate calculations of the realistic dose delivered to a moving and deforming anatomy during static and dynamic beam delivery techniques.     

Clinical verification of treatment planning dose calculation in lung SBRT with GATE Monte Carlo simulation code

PurposeThis study aims to use GATE/Geant4 simulation code to evaluate the performance of dose calculations with Anisotropic Analytical Algorithm (AAA) in the context of lung SBRT for complex treatments considering images of patients.MethodsFour cases of non-small cell lung cancer treated with SBRT were selected for this study. Irradiation plans were created with AAA and recalculated end to end using Monte Carlo (MC) method maintaining field configurations identical to the original plans. Each treatment plan was evaluated in terms of PTV and organs at risk (OARs) using dose-volume histograms (DVH). Dosimetric parameters obtained from DVHs were used to compare AAA and MC.ResultsThe comparison between the AAA and MC DVH using gamma analysis with the passing criteria of 3%/3% showed an average passing rate of more than 90% for the PTV structure and 97% for the OARs. Tightening the criteria to 2%/2% showed a reduction in the average passing rate of the PTV to 86%. The agreement between the AAA and MC dose calculations for PTV dosimetric parameters (V₁₀₀; V₉₀; Homogeneity index; maximum, minimum and mean dose; CI_Paddick and D_2cm) was within 18.4%. For OARs, the biggest differences were observed in the spinal cord and the great vessels.ConclusionsIn general, we did not find significant differences between AAA and MC. The results indicate that AAA could be used in complex SBRT cases that involve a larger number of small treatment fields in the presence of tissue heterogeneities.     

A depth dose study between AAA and AXB algorithm against Monte Carlo simulation using AIP CT of a 4D dataset from a moving phantom

Aim

To identifying depth dose differences between the two versions of the algorithms using AIP CT of a 4D dataset.

Assessment of skin dose in breast cancer radiotherapy: on-phantom measurement and Monte Carlo simulation

AimThe main purpose of the present study is assessment of skin dose in breast cancer radiotherapy.BackgroundAccurate assessment of skin dose in radiotherapy can provide useful information for clinical considerations.Materials and methodsA RANDO phantom was irradiated using a 6 MV Siemens Primus linac with medial and tangential radiotherapy fields for simulating breast cancer treatment. Dosimetry was also performed on various positions across the fields using an EBT3 radiochromic film. Similar conditions of measurement on the RANDO phantom including field size, irradiation angle, number of fields, etc. were subsequently simulated via the Monte Carlo N-Particle Transport code (MCNP). Ultimately, dose values for corresponding points from both methods were compared.ResultsConsidering dosimetry using radiochromic films on the RANDO phantom, there were points having underdose and overdose based on the prescribed dose and skin tolerance levels. In this respect, 81.25% and 18.75% of the points had underdose and overdose, respectively. In some cases, several differences were observed between the measurement and the MCNP simulation results associated with skin dose.ConclusionBased on the results of the points which had underdose, it was suggested that a bolus should be used for the given points. With regard to overdose points, it was advocated to consider skin tolerance dose in treatment planning. Differences between the measurement and the MCNP simulation results might be due to voxel size of tally cells in simulations, effect of beam’s angle of incidence, validation time of linac’s head, lack of electronic equilibrium in the build-up region, as well as MCNP tally type.     

Assessment of out-of-field doses in radiotherapy treatments of paediatric patients using Monte Carlo methods and measurements

PurposeTo assess out-of-field doses in radiotherapy treatments of paediatric patients, using Monte Carlo methods to implement a new model of the linear accelerator validated against measurements and developing a voxelized anthropomorphic paediatric phantom.MethodsCT images of a physical anthropomorphic paediatric phantom were acquired and a dosimetric planning using a TPS was obtained. The CT images were used to perform the voxelization of the physical phantom using the ImageJ software and later implemented in MCNP. In order to validate the Monte Carlo model, dose measurements of the 6 MV beam and Linac with 120 MLC were made in a clinical setting, using ionization chambers and a water phantom. Afterwards TLD measurements in the physical anthropomorphic phantom were performed in order to assess the out-of-field doses in the eyes, thyroid, c-spine, heart and lungs.ResultsThe Monte Carlo model was validated for in-field and out-of-field doses with average relative differences below 3%. The average relative differences between TLD measurements and Monte Carlo is 14,3% whilst the average relative differences between TLD and TPS is 55,8%. Moreover, organs up to 22.5 cm from PTV center show TLD and MCNP6 relative differences and TLD and TPS relative differences up to 21.2% and 92.0%, respectively.ConclusionsOur study provides a novel model that could be used in clinical research, namely in dose evaluation outside the treatment fields. This is particularly relevant, especially in pediatric patients, for studying new radiotherapy treatment techniques, since it can be used to estimate the development of secondary tumours.     

An EGS4 based Monte Carlo code for the calculation of organ equivalent dose to a modified Yale voxel phantom.

Organ or tissue equivalent dose, the most important quantity in radiation protection, cannot be measured directly. Therefore it became common practice to calculate the quantity of interest with Monte Carlo methods applied to so-called human phantoms, which are virtual representations of the human body. The Monte Carlo computer code determines conversion coefficients, which are ratios between organ or tissue equivalent dose and measurable quantities. Conversion coefficients have been published by the ICRP (Report No. 74) for various types of radiation, energies and fields, which have been calculated, among others, with the mathematical phantoms ADAM and EVA. Since then progress of image processing, and of clock speed and memory capacity of computers made it possible to create so-called voxel phantoms, which are a far more realistic representation of the human body. Voxel (Volume pixel) phantoms are built from segmented CT and/or MRI images of real persons. A complete set of such images can be joined to a 3-dimensional representation of the human body, which can be linked to a Monte Carlo code allowing for particle transport calculations. A modified version of the VOX_TISS8 human voxel phantom (Yale University) has been connected to the EGS4 Monte Carlo code. The paper explains the modifications, which have been made, the method of coupling the voxel phantom with the code, and presents results as conversion coefficients between organ equivalent dose and kerma in air for external photon radiation. A comparison of the results with published data shows good agreement.     

Monte Carlo study of the influence of energy spectra,mesh size,high Z element on dose and PVDR based on 1-D and 3-D heterogeneous mouse head phantom for Microbeam Radiation Therapy

PurposeTo evaluate the influence of energy spectra, mesh sizes, high Z element on dose and PVDR in Microbeam Radiation Therapy (MRT) based on 1-D analogy-mouse-head-model (1-D MHM) and 3-D voxel-mouse-head-phantom (3-D VMHP) by Monte Carlo simulation.MethodsA Microbeam-Array-Source-Model was implemented into EGSnrc/DOSXYZnrc. The microbeam size is assumed to be 25 μm, 50 μm or 75 μm in thickness and fixed 1 mm in height with 200 μm c-t-c. The influence of the energy spectra of ID17@ESRF and BMIT@CLS were investigated. The mesh size was optimized. PVDR in 1-D MHM and 3-D VMHP was compared with the homogeneous water phantom. The arc influence of 3-D VMHP filled with water (3-D VMHWP) was compared with the rectangle phantom.ResultsPVDR of the lower BMIT@CLS spectrum is 2.4 times that of ID17@ESRF for lower valley dose. The optimized mesh is 5 µm for 25 µm, and 10 µm for 50 µm and 75 µm microbeams with 200 µm c-t-c. A 500 μm skull layer could make PVDR difference up to 62.5% for 1-D MHM. However this influence is limited (<5%) for the farther homogeneous media (e.g. 600 µm). The peak dose uniformity of 3-D VMHP at the same depth could be up to 8% for 1.85 mm × 1 mm irradiation field, whereas that of 3-D VMHWP is <1%. The high Z element makes the dose uniformity enhance in target. The surface arc could affect the superficial PVDR (from 44% to 21% in 0.2 mm depth), whereas this influence is limited for the more depth (<1%).ConclusionAn accurate MRT dose calculation algorithm should include the influence of 3-D heterogeneous media.     

Monte Carlo simulation for Neptun 10 PC medical linear accelerator and calculations of output factor for electron beam

AimExact knowledge of dosimetric parameters is an essential pre-requisite of an effective treatment in radiotherapy. In order to fulfill this consideration, different techniques have been used, one of which is Monte Carlo simulation.Materials and methodsThis study used the MCNP-4C^b to simulate electron beams from Neptun 10 PC medical linear accelerator. Output factors for 6, 8 and 10 MeV electrons applied to eleven different conventional fields were both measured and calculated.ResultsThe measurements were carried out by a Wellhofler-Scanditronix dose scanning system. Our findings revealed that output factors acquired by MCNP-4C simulation and the corresponding values obtained by direct measurements are in a very good agreement.ConclusionIn general, very good consistency of simulated and measured results is a good proof that the goal of this work has been accomplished.     

Normalized glandular dose coefficients in mammography,digital breast tomosynthesis and dedicated breast CT

PurposeTo provide mean glandular dose (MGD) estimates via Monte Carlo (MC) simulations as a function of the breast models and scan parameters in mammography, digital breast tomosynthesis (DBT) and dedicated breast CT (BCT).MethodsThe MC code was based on GEANT4 toolkit. The simulated compressed breast was either a cylinder with a semi-circular section or ad hoc shaped for oblique view (MLO). In DBT we studied the influence of breast models and exam parameters on the T-factors (i.e. the conversion factor for the calculation of the MGD in DBT from that for a 0-degree projection), and in BCT we investigated the influence on the MGD estimates of the ion chamber volume used for the air kerma measurements.ResultsIn mammography, a model representative of a breast undergoing an MLO view exam did not produce substantial differences (0.4%) in MGD estimates, when compared to a conventional cranio-caudal (CC) view breast model. The beam half value layer did not present a significant influence on T-factors in DBT (<0.8%), while the skin model presented significant influence on MGD estimates (up to 3.3% at 30 degrees scan angle), increasing for larger scan angles. We derived a correction factor for taking into account the different ion chamber volume used in MGD estimates in BCT.ConclusionsA series of MC code modules for MGD estimates in 2D and 3D breast imaging have been developed in order to take into account the most recent advances in breast models.     

Radiation dose around a PET scanner installation: Comparison of Monte Carlo simulations,analytical calculations and experimental results

PurposeMonte Carlo study of radiation transmission around areas surrounding a PET room.MethodsAn extended population of patients administered with ¹⁸F-FDG for PET-CT investigations was studied, collecting air kerma rate and gamma ray spectra measurements at a reference distance. An MC model of the diagnostic room was developed, including the scanner and walls with variable material and thickness. MC simulations were carried out with the widely used code GEANT4.ResultsThe model was validated by comparing simulated radiation dose values and gamma ray spectra produced by a volumetric source with experimental measurements; ambient doses in the surrounding areas were assessed for different combinations of wall materials and shielding and compared with analytical calculations, based on the AAPM Report 108.In the range 1.5–3.0 times of the product between the linear attenuation coefficient and thickness of an absorber (μ x), it was observed that the effectiveness of different combinations of shielding is roughly equivalent. An extensive tabulation of results is given in the text.ConclusionsThe validation tests performed showed a satisfactory agreement between the simulated and expected results. The simulated dose rates incident on, and transmitted by the walls in our model of PET scanner room, are generally in good agreement with analytical estimates performed using the AAPM Publication No. 108 method. This provides an independent confirmation of AAPM's approach. Even in this specific field of application, GEANT4 proved to be a relevant and accurate tool for dosimetry estimates, shielding evaluation and for general radiation protection use.     

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.