Influence of dental restorative materials on ESR biodosimetry in tooth enamel

Using an experimental model and PENELOPE Monte Carlo simulations, the effects of resin and amalgam on the absorbed doses in tooth enamel were studied to evaluate the feasibility of using restored teeth in electron spin resonance (ESR) dose reconstruction. The model consisted of a phantom containing a plate of these restorative materials placed between powered enamel layers exposed to X rays and a ??Co beam. The experimental results and simulations agreed, showing that the attenuation produced by amalgam and resin with a thickness of 1, 2, and 4 mm is similar to that produced by the enamel itself in the case of the radiation sources employed. For X rays and ??Co γ radiation the attenuation reached almost 100% and 40%, respectively. These results show that for ESR dose reconstruction, the use of all available enamel of a tooth leads to errors in the estimated dose due to attenuation effects in both healthy and restored teeth. Thus the importance of an enamel selection from different sides of the tooth surface to apply ESR dose reconstruction in the case of a practical situation is shown.     

Validation of Monte carlo Geant4 multithreading code for a 6 MV photon beam of varian linac on the grid computing

AimTo evaluate the computation time efficiency of the multithreaded code (G4Linac-MT) in the dosimetry application, using the high performance of the HPC-Marwan grid to determine with high accuracy the initial parameters of the 6 MV photon beam of Varian CLINAC 2100C.BackgroundThe difficulty of Monte Carlo methods is the long computation time, this is one of the disadvantages of the Monte Carlo methods.Materials and methodsCalculations are performed by the multithreaded code G4Linac-MT and Geant4.10.04.p02 using the HPC-Marwan computing grid to evaluate the computing speed for each code. The multithreaded version is tested in several CPUs to evaluate the computing speed according to the number of CPUs used. The results were compared to the measurements using different types of comparisons, TPR_20.10, penumbra, mean dose error and gamma index.ResultsThe results obtained for this work indicate a much higher computing time saving for the G4Linac-MT version compared to the Geant4.10.04 version, the computing time decreases with the number of CPUs used, can reach about 12 times if 64CPUs are used. After optimization of the initial electron beam parameters, the results of the dose simulations obtained for this work are in very good agreement with the experimental measurements with a mean dose error of up to 0.41% on the PDDs and 1.79% on the lateral dose.ConclusionsThe gain in computation time leads us to perform Monte Carlo simulations with a large number of events which gives a high accuracy of the dosimetry results obtained in this work.     

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.     

Monte Carlo modeling of the Siemens Optifocus multileaf collimator

We have developed a new component module for the BEAMnrc software package, called SMLC, which models the tongue-and-groove structure of the Siemens Optifocus multileaf collimator. The ultimate goal is to perform accurate Monte Carlo simulations of the IMRT treatments carried out with Optifocus. SMLC has been validated by direct geometry checks and by comparing quantitatively the results of simulations performed with it and with the component module VARMLC. Measurements and Monte Carlo simulations of absorbed dose distributions of radiation fields sensitive to the tongue-and-groove effect have been performed to tune the free parameters of SMLC. The measurements cannot be accurately reproduced with VARMLC. Finally, simulations of a typical IMRT field showed that SMLC improves the agreement with experimental measurements with respect to VARMLC in clinically relevant cases.PACS number87.55. K-     

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.     

Integration of the M6 Cyberknife in the Moderato Monte Carlo platform and prediction of beam parameters using machine learning

PurposeThis work describes the integration of the M6 Cyberknife in the Moderato Monte Carlo platform, and introduces a machine learning method to accelerate the modelling of a linac.MethodsThe MLC-equipped M6 Cyberknife was modelled and integrated in Moderato, our in-house platform offering independent verification of radiotherapy dose distributions. The model was validated by comparing TPS dose distributions with Moderato and by film measurements. Using this model, a machine learning algorithm was trained to find electron beam parameters for other M6 devices, by simulating dose curves with varying spot size and energy. The algorithm was optimized using cross-validation and tested with measurements from other institutions equipped with a M6 Cyberknife.ResultsOptimal agreement in the Monte Carlo model was reached for a monoenergetic electron beam of 6.75 MeV with Gaussian spatial distribution of 2.4 mm FWHM. Clinical plan dose distributions from Moderato agreed within 2% with the TPS, and film measurements confirmed the accuracy of the model. Cross-validation of the prediction algorithm produced mean absolute errors of 0.1 MeV and 0.3 mm for beam energy and spot size respectively. Prediction-based simulated dose curves for other centres agreed within 3% with measurements, except for one device where differences up to 6% were detected.ConclusionsThe M6 Cyberknife was integrated in Moderato and validated through dose re-calculations and film measurements. The prediction algorithm was successfully applied to obtain electron beam parameters for other M6 devices. This method would prove useful to speed up modelling of new machines in Monte Carlo systems.     

Evaluation of 133Xe for x-ray fluorescence analysis of cadmium in vivo: a Monte Carlo study

A Monte Carlo model has been developed to support the design of a 180° geometry x-ray fluorescence system for the measurement of cadmium concentration in deep body organs such as the kidney. ¹³³Xe was investigated as the excitation photon source. A total number of 15×10⁶ simulated incident photons were used. Monte Carlo simulations were performed using the EGS4 Monte Carlo code system. The results showed that for distances between the skin and the kidney surface of 30–60 mm, respectively, cadmium concentrations of 15–60 μg/g kidney tissue could easily be detected. The mean skin and kidney doses during such measurements were estimated to be between 8 and 0.9 mGy, respectively. Received: 1 June 1999 / Accepted: 10 February 2000     

Dosimetry for quantitative analysis of the effects of low-dose ionizing radiation in radiation therapy patients

We have developed and validated a practical approach to identifying the location on the skin surface that will receive a prespecified biopsy dose (ranging down to 1 cGy) in support of in vivo biological dosimetry in humans. This represents a significant technical challenge since the sites lie on the patient's surface outside the radiation fields. The PEREGRINE Monte Carlo simulation system was used to model radiation dose delivery, and TLDs were used for validation on phantoms and for confirmation during patient treatment. In the developmental studies, the Monte Carlo simulations consistently underestimated the dose at the biopsy site by approximately 15% (of the local dose) for a realistic treatment configuration, most likely due to lack of detail in the simulation of the linear accelerator outside the main beam line. Using a single, thickness-independent correction factor for the clinical calculations, the average of 36 measurements for the predicted 1-cGy point was 0.985 cGy (standard deviation: 0.110 cGy) despite patient breathing motion and other real-world challenges. Since the 10-cGy point is situated in the region of high-dose gradient at the edge of the field, patient motion had a greater effect, and the six measured points averaged 5.90 cGy (standard deviation: 1.01 cGy), a difference that is equivalent to approximately a 6-mm shift on the patient's surface.     

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.     

A Study on the Dose Distributions in Various Materials from an Ir-192 HDR Brachytherapy Source

Dose distributions of (192)Ir HDR brachytherapy in phantoms simulating water, bone, lung tissue, water-lung and bone-lung interfaces using the Monte Carlo codes EGS4, FLUKA and MCNP4C are reported. Experiments were designed to gather point dose measurements to verify the Monte Carlo results using Gafchromic film, radiophotoluminescent glass dosimeter, solid water, bone, and lung phantom. The results for radial dose functions and anisotropy functions in solid water phantom were consistent with previously reported data (Williamson and Li). The radial dose functions in bone were affected more by depth than those in water. Dose differences between homogeneous solid water phantoms and solid water-lung interfaces ranged from 0.6% to 14.4%. The range between homogeneous bone phantoms and bone-lung interfaces was 4.1% to 15.7%. These results support the understanding in dose distribution differences in water, bone, lung, and their interfaces. Our conclusion is that clinical parameters did not provide dose calculation accuracy for different materials, thus suggesting that dose calculation of HDR treatment planning systems should take into account material density to improve overall treatment quality.     

Towards MR-guided electron therapy: Measurement and simulation of clinical electron beams in magnetic fields

PurposeIn the current era of MRI-linac radiotherapy, dose optimization with arbitrary dose distributions is a reality. For the first time, we present new and targeted experiments and modeling to aid in evaluating the potential dose improvements offered with an electron beam mode during MRI-linac radiotherapy.MethodsSmall collimated (1 cm diameter and 1.5 × 1.5 cm² square) electron beams (6, 12 and 20 MeV) from a clinical linear accelerator (Varian Clinac 2100C) are incident perpendicular and parallel to the strong and localized magnetic fields (0–0.7 T) generated by a permanent magnet device. Gafchromic EBT3 film is placed inside a slab phantom to measure two-dimensional dose distributions. A benchmarked and comprehensive Monte Carlo model (Geant4) is established to directly compare with experiments.ResultsWith perpendicular fields a 5% narrowing of the beam FWHM and a 10 mm reduction in the 15% isodose penetration is seen for the 20 MeV beam. In the inline setup the penumbral width is reduced by up to 20%, and a local central dose enhancement of 100% is observed. Monte Carlo simulations are in agreement with the measured dose distributions (2% or 2 mm).ConclusionA new range of experiments have been performed to offer insight into how an electron beam mode could offer additional choices in MRI-linac radiotherapy. The work extends on historic studies to bring a successful unified experimental and Monte Carlo modeling approach for studying small field electron beam dosimetry inside magnetic fields. The results suggest further work, particularly on the inline magnetic field scenario.     

Influence of the Presence of Tissue Expanders on Energy Deposition for Post-Mastectomy Radiotherapy

An increasing number of studies have shown that post-mastectomy radiotherapy presents benefits associated with the patients survival and a significant fraction of the treated patients makes use of tissue expanders for breast reconstruction. Some models of tissue expanders have a magnetic disk on their surface that constitutes heterogeneity in the radiation field, which can affect the dose distribution during the radiotherapy treatment. In this study, the influence of a metallic heterogeneity positioned in a breast tissue expander was evaluated by means of Monte Carlo simulations using the MCNPX code and using Eclipse treatment planning system. Deposited energy values were calculated in structures which have clinical importance for the treatment. Additionally, the effect in the absorbed energy due to backscattering and attenuation of the incident beam caused by the heterogeneity, as well as due to the expansion of the prosthesis, was evaluated in target structures for a 6 MV photon beam by simulations. The dose distributions for a breast treatment were calculated using a convolution/superposition algorithm from the Eclipse treatment planning system. When compared with the smallest breast expander volume, underdosage of 7% was found for the largest volume of breast implant, in the case of frontal irradiation of the chest wall, by Monte Carlo simulations. No significant changes were found in dose distributions for the presence of the heterogeneity during the treatment planning of irradiation with an opposed pair of beams. Even considering the limitation of the treatment planning system, the results obtained with its use confirm those ones found by Monte Carlo simulations for a tangent beam irradiation. The presence of a heterogeneity didńt alters the dose distributions on treatment structures. The underdosage of 7% observed with Monte Carlo simulations were found for irradiation at 0°, not used frequently in a clinical routine.     

Assessment of the neutron dose field around a biomedical cyclotron: FLUKA simulation and experimental measurements

In the planning of a new cyclotron facility, an accurate knowledge of the radiation field around the accelerator is fundamental for the design of shielding, the protection of workers, the general public and the environment.Monte Carlo simulations can be very useful in this process, and their use is constantly increasing. However, few data have been published so far as regards the proper validation of Monte Carlo simulation against experimental measurements, particularly in the energy range of biomedical cyclotrons.In this work a detailed model of an existing installation of a GE PETtrace 16.5 MeV cyclotron was developed using FLUKA. An extensive measurement campaign of the neutron ambient dose equivalent H¹(10) in marked positions around the cyclotron was conducted using a neutron rem-counter probe and CR39 neutron detectors. Data from a previous measurement campaign performed by our group using TLDs were also re-evaluated.The FLUKA model was then validated by comparing the results of high-statistics simulations with experimental data. In 10 out of 12 measurement locations, FLUKA simulations were in agreement within uncertainties with all the three different sets of experimental data; in the remaining 2 positions, the agreement was with 2/3 of the measurements.Our work allows to quantitatively validate our FLUKA simulation setup and confirms that Monte Carlo technique can produce accurate results in the energy range of biomedical cyclotrons.     

Monte Carlo studies of a model for lipid-gramicidin A bilayers.

This paper presents results of Monte Carlo simulations of a full bilayer of 200 lipid chains and one gramicidin A dimer. Simulations are described for systems with lipid chains of 14, 16, and 18 carbons, respectively. Using accepted potential functions to calculate interactions between all non-hydrogen atoms a Monte Carlo configuration sampling is generated from which order parameter profiles are calculated and specific configurations are displayed. Results are compared with experimental data for lipid-gramicidin bilayers.     

Towards breast cancer rotational radiotherapy with synchrotron radiation

PurposeWe performed the first investigations, via measurements and Monte Carlo simulations on phantoms, of the feasibility of a new technique for synchrotron radiation rotational radiotherapy for breast cancer (SR³T).MethodsA Monte Carlo (MC) code based on Geant4 toolkit was developed in order to simulate the irradiation with the SR³T technique and to evaluate the skin sparing effect in terms of centre-to-periphery dose ratio at different energies in the range 60–175 keV. Preliminary measurements were performed at the Australian Synchrotron facility. Radial dose profiles in a 14-cm diameter polyethylene phantom were measured with a 100-mm pencil ionization chamber for different beam sizes and compared with the results of MC simulations. Finally, the dose painting feasibility was demonstrated with measurements with EBT3 radiochromic films in a phantom and collimating the SR beam at 1.5 cm in the horizontal direction.ResultsMC simulations showed that the SR³T technique assures a tumour-to-skin absorbed dose ratio from about 7:1 (at 60 keV photon energy) to about 10:1 (at 175 keV), sufficient for skin sparing during radiotherapy. The comparison between the results of MC simulations and measurements showed an agreement within 5%. Two off-centre foci were irradiated shifting the rotation centre in the horizontal direction.ConclusionsThe SR³T technique permits to obtain different dose distributions in the target with multiple rotations and can be guided via synchrotron radiation breast computed tomography imaging, in propagation based phase-contrast conditions. Use of contrast agents like iodinated solutions or gold nanoparticles for dose enhancement (DE-SR³T) is foreseen and will be investigated in future work.     

Benchmark of penelope for low and medium energy X-rays

The Monte Carlo code penelope is benchmarked for X-ray beams with energies between 30 and 300 keV. The results of different simulations performed with penelope are compared with those obtained with a semi-empirical computational model and with experimental measurements. Half-value layer indexes obtained from the attenuation curves for Al and Cu and depth dose curves in water have been considered for this comparison. A good agreement is reached on what guarantees the feasibility of the code.     

Ion passage pathways and thermodynamics of the amphotericin B membrane channel

Amphotericin B is a polyene macrolide antibiotic used to treat systemic fungal infections. Amphotericin B's chemotherapeutic action requires the formation of transmembrane channels, which are known to transmit monovalent ions. We have investigated the ion passage pathways through the pore of a realistic model structure of the channel and computed the associated thermodynamic properties. Our calculations combined the free energy computations using the Poisson equation with a continuum solvent model and the molecular simulations in which solvent molecules were present explicitly. It was found that there are no substantial structural barriers to a single sodium or chloride ion passage. Thermodynamic free energy calculations showed that the path along which the ions prefer to move is off center from the channel's central axis. In accordance with experiments, Monte Carlo molecular simulations established that sodium ions can pass through the pore. When it encounters a chloride anion in the channel, the sodium cation prefers to form a solvent-bridged pair configuration with the anion.     

Experimental measurements and Monte Carlo calculations for 103Pd dosimetry of the 12 mm COMS eye plaque

Monte Carlo simulations and TLD dosimetry have been performed to determine the dose distributions along the central axis of the 12 mm COMS eye plaques loaded with IRA1-¹⁰³Pd seeds. Several simulations and measurements have been employed to investigate the effect of Silastic insert and air in front of the eye on dosimetry results along the central axis of the plaque and at some critical ocular structures. Measurements were performed using TLD-GR200A circular chip dosimeters in a PMMA eye phantom. The central axis TLD chips locations were arranged in one central column of eye phantom, in 3 mm intervals. The off-axis TLD chips locations were arranged in three off-axis columns around the central axis column. Version 5 of the MCNP code was also used to evaluate the dose distribution around the plaque. The presence of the Silastic insert results in dose reduction of 14% at 5 mm; also about 7% dose reduction appears at the interface point, due to the air presence and lack of the scattering condition. The overall dosimetric parameters for the COMS eye plaque loaded with new palladium seeds are similar to a commercial widely used seed such as Theragenics200. As the dose calculations under TG-43 assumptions do not consider the effect of the plaque backing and Silastic insert for accurate dosimetry, it's suggested to apply the effect of the eye plaque materials and air on dosimetry results along the central axis of the plaque and at some critical ocular structures.     

Lower confidence bounds for prediction accuracy in high dimensions via AROHIL Monte Carlo

MOTIVATION: Implementation and development of statistical methods for high-dimensional data often require high-dimensional Monte Carlo simulations. Simulations are used to assess performance, evaluate robustness, and in some cases for implementation of algorithms. But simulation in high dimensions is often very complex, cumbersome and slow. As a result, performance evaluations are often limited, robustness minimally investigated and dissemination impeded by implementation challenges. This article presents a method for converting complex, slow high-dimensional Monte Carlo simulations into simpler, faster lower dimensional simulations. RESULTS: We implement the method by converting a previous Monte Carlo algorithm into this novel Monte Carlo, which we call AROHIL Monte Carlo. AROHIL Monte Carlo is shown to exactly or closely match pure Monte Carlo results in a number of examples. It is shown that computing time can be reduced by several orders of magnitude. The confidence bound method implemented using AROHIL outperforms the pure Monte Carlo method. Finally, the utility of the method is shown by application to a number of real microarray datasets.