Monte Carlo calculation of charge collection efficiencies in ionization chambers.

Collection efficiency is an important quantity in dosimetry with ionization chambers. It can be calculated by solving a hyperbolic system of partial differential equations. This system can be solved only in few, simple, idealized geometries, but for more realistic designs an analytical resolution is no longer possible. In the present work a Monte Carlo scheme that could permit to calculate the collection efficiency for any ionization chamber geometry is proposed. This scheme has been tested against Boag's approach for three chambers with plane-parallel, cylindrical and spherical geometries, operated in the recombination regime. The results obtained in the full Monte Carlo simulation closely agree with the Boag's ones for the three ideal geometries considered. The largest relative difference, ∼0.3%, has been found for the plane-parallel chamber in case of 50 V, the lowest potential difference investigated in this study. Results appear to be stable against changes in the chamber volume, the ion mobility and the recombination constant. The method proposed could be a useful tool to calculate collection efficiencies of ionization chambers, provided the electric field inside them is known.

An on-lattice agent-based Monte Carlo model simulating the growth kinetics of multicellular tumor spheroids.

PURPOSE: To develop an on-lattice agent-based model describing the growth of multicellular tumor spheroids using simple Monte Carlo tools. METHODS: Cells are situated on the vertices of a cubic grid. Different cell states (proliferative, hypoxic or dead) and cell evolution rules, driven by 10 parameters, and the effects of the culture medium are included. About twenty spheroids of MCF-7 human breast cancer were cultivated and the experimental data were used for tuning the model parameters. RESULTS: Simulated spheroids showed adequate sizes of the necrotic nuclei and of the hypoxic and proliferative cell phases as a function of the growth time, mimicking the overall characteristics of the experimental spheroids. The relation between the radii of the necrotic nucleus and the whole spheroid obtained in the simulations was similar to the experimental one and the number of cells, as a function of the spheroid volume, was well reproduced. The statistical variability of the Monte Carlo model described the whole volume range observed for the experimental spheroids. Assuming that the model parameters vary within Gaussian distributions it was obtained a sample of spheroids that reproduced much better the experimental findings. CONCLUSIONS: The model developed allows describing the growth of in vitro multicellular spheroids and the experimental variability can be well reproduced. Its flexibility permits to vary both the agents involved and the rules that govern the spheroid growth. More general situations, such as, e. g., tumor vascularization, radiotherapy effects on solid tumors, or the validity of the tumor growth mathematical models can be studied.

Simple variance reduction in Monte Carlo calculations of specific absorbed fractions: Russian roulette and splitting at the source organ.

PURPOSE: To investigate the capabilities of several variance reduction techniques in the calculation of specific absorbed fractions in cases where the source and the target organs are far away and/or the target organs have a small volume. METHODS: The specific absorbed fractions have been calculated by using the Monte Carlo code PENELOPE and by assuming the thyroid gland as the source organ and the testicles, the urinary bladder, the uterus, and the ovaries as the target ones. A mathematical anthropomorphic phantom, similar to the MIRD-type phantoms, has been considered. Photons with initial energies of 50, 100 and 500 keV were emitted isotropically from the volume of the source organ. Simulations have been carried out by implementing the variance reduction techniques of splitting and Russian roulette at the source organ only and the interaction forcing at the target organs. The influence of the implementation details of those techniques have been investigated and optimal parameters have been determined. All simulations were run with a CPU time of 1.5 · 105 s. RESULTS: Specific absorbed fractions with relative uncertainties well below 10% have been obtained in most cases, agreeing with those used as reference. The best value for the factor defining the application of the Russian roulette technique was r = 0.3. The best value for the splitting number was between s = 3 and s = 10, depending on the specific energies and target organs. CONCLUSIONS: The proposed strategy provides an effective method for computing specific absorbed fractions for the most unfavorable situations, with a computing effort that is considerably reduced with respect to other methodologies.

Computation of the electron beam quality [Formula: see text] factors for the NE2571, NE2571A and NE2581A thimble ionization chambers using PENELOPE.

The quality correction factor [Formula: see text] for electron beams was calculated for three thimble ionization chambers, namely, NE2571, NE2571A and NE2581A. The Monte Carlo code PENELOPE was used to estimate the overall correction factor fc,Q of these chambers for electron beams with nominal energies ranging between 6 and 22MeV, corresponding to a Varian Clinac 2100 C/D. A 60Co beam was used as reference quality Q0. Also eight monoenergetic electron beams reproducing the quality index R50 of the Clinac beams were considered. The [Formula: see text] factors were calculated as the ratio between fc,Q and [Formula: see text] . Those obtained for the NE2571 ionization chamber show a nice agreement with those calculated by Muir and Rogers with EGSnrc. As it occurred to other ionization chambers analyzed in previous works, the [Formula: see text] factors found for the monoenergetic beams are larger (smaller) than those corresponding to the Clinac beams at low (high) R50 values, the differences being slightly above 0.5%. Finally, the [Formula: see text] factors obtained in the case of the NE2571A chamber are systematically ∼0.5% below those of its predecessor chamber, the NE2571.

Air density dependence of the response of the PTW SourceCheck 4pi ionization chamber for 125I brachytherapy seeds.

PURPOSE: To analyze the air density dependence of the response of the new SourceCheck 4pi ionization chamber, manufactured by PTW. METHODS: The air density dependence of three different SourceCheck 4pi chambers was studied by measuring 125I sources. Measurements were taken by varying the pressure from 746.6 to 986.6hPa in a pressure chamber. Three different HDR 1000 Plus ionization chambers were also analyzed under similar conditions. A linear and a potential-like function of the air density were fitted to experimental data and their achievement in describing them was analyzed. RESULTS: SourceCheck 4pi chamber response showed a residual dependence on the air density once the standard pressure and temperature factor was applied. The chamber response was overestimated when the air density was below that under normal atmospheric conditions. A similar dependence was found for the HDR 1000 Plus chambers analyzed. A linear function of the air density permitted a very good description of this residual dependence, better than with a potential function. No significant variability between the different specimens of the same chamber model studied was found. CONCLUSION: The effect of overestimation observed in the chamber responses once they are corrected for the standard pressure and temperature may represent a non-negligible ∼4% overestimation in high altitude cities as ours (700m AMSL). This overestimation behaves linearly with the air density in all cases analyzed.

Monte Carlo calculation of specific absorbed fractions: variance reduction techniques.

The purpose of the present work is to calculate specific absorbed fractions using variance reduction techniques and assess the effectiveness of these techniques in improving the efficiency (i.e. reducing the statistical uncertainties) of simulation results in cases where the distance between the source and the target organs is large and/or the target organ is small. The variance reduction techniques of interaction forcing and an ant colony algorithm, which drives the application of splitting and Russian roulette, were applied in Monte Carlo calculations performed with the code penelope for photons with energies from 30 keV to 2 MeV. In the simulations we used a mathematical phantom derived from the well-known MIRD-type adult phantom. The thyroid gland was assumed to be the source organ and urinary bladder, testicles, uterus and ovaries were considered as target organs. Simulations were performed, for each target organ and for photons with different energies, using these variance reduction techniques, all run on the same processor and during a CPU time of 1.5 · 10(5) s. For energies above 100 keV both interaction forcing and the ant colony method allowed reaching relative uncertainties of the average absorbed dose in the target organs below 4% in all studied cases. When these two techniques were used together, the uncertainty was further reduced, by a factor of 0.5 or less. For photons with energies below 100 keV, an adapted initialization of the ant colony algorithm was required. By using interaction forcing and the ant colony algorithm, realistic values of the specific absorbed fractions can be obtained with relative uncertainties small enough to permit discriminating among simulations performed with different Monte Carlo codes and phantoms. The methodology described in the present work can be employed to calculate specific absorbed fractions for arbitrary arrangements, i.e. energy spectrum of primary radiation, phantom model and source and target organs.

Electron beam quality k(Q,Q0) factors for various ionization chambers: a Monte Carlo investigation with PENELOPE.

In this work we calculate the beam quality correction factor k(Q,Q0) for various plane-parallel ionization chambers. A set of Monte Carlo calculations using the code PENELOPE/PENEASY have been carried out to calculate the overall correction factor f(c,Q) for eight electron beams corresponding to a Varian Clinac 2100 C/D, with nominal energies ranging between 6 MeV and 22 MeV, for a (60)Co beam, that has been used as the reference quality Q0 and also for eight monoenergetic electron beams reproducing the quality index R50 of the Clinac beams. Two field sizes, 10 × 10 cm(2) and 20 × 20 cm(2) have been considered. The k(Q,Q0) factors have been calculated as the ratio between f(c,Q) and f(c,Q0). Values for the Exradin A10, A11, A11TW, P11, P11TW, T11 and T11TW ionization chambers, manufactured by Standard Imaging, as well as for the NACP-02 have been obtained. The results found with the Clinac beams for the two field sizes analyzed show differences below 0.6%, even in the case of the higher energy electron beams. The k(Q,Q0) values obtained with the Clinac beams are 1% larger than those found with the monoenergetic beams for the higher energies, above 12 MeV. This difference can be ascribed to secondary photons produced in the linac head and the air path towards the phantom. Contrary to what was quoted in a previous work (Sempau et al 2004 Phys. Med. Biol. 49 4427-44), the beam quality correction factors obtained with the complete Clinac geometries and with the monoenergetic beams differ significantly for energies above 12 MeV. Material differences existing between chambers that have the same geometry produce non-negligible modifications in the value of these correction factors.

A method to relate StarTrack(®) measurements to R50 variations in clinical linacs.

The relation between the data recorded with any device for the daily checking of the behavior of a clinical linac and the reference magnitudes to be monitored may be unknown. An experimental method relating the energy stability of the electron beam measured with StarTrack(®) to the R50 beam quality index is proposed. The bending magnet current is varied producing a change in the exit energy window and, therefore, a modification of the R50 value. For different values of this current, the output data of StarTrack(®) and the R50, obtained from depth doses measured in a water phantom are determined. A linear fit between both sets of data allows the identification of the StarTrack(®) output that provides the best way to obtain the quality index R50, for each beam nominal energy. Using these fits, an historical datum series is used to analyze the method proposed in the daily quality control. The ouput data of the StarTrack(®) and the R50 values show a good linear correlation. It is possible to establish a methodology that allows the monitoring of R50 by direct use of the daily quality control data measured with StarTrack(®). A method to monitor R50 in the daily quality control using the StarTrack(®) device has been developed. The method may be applied to similar devices in which the statistical control variable does not show a linear behavior with R50.

Monte Carlo study for designing a dedicated "D"-shaped collimator used in the external beam radiotherapy of retinoblastoma patients.

PURPOSE: Retinoblastoma is the most common intraocular malignancy in the early childhood. Patients treated with external beam radiotherapy respond very well to the treatment. However, owing to the genotype of children suffering hereditary retinoblastoma, the risk of secondary radio-induced malignancies is high. The University Hospital of Essen has successfully treated these patients on a daily basis during nearly 30 years using a dedicated "D"-shaped collimator. The use of this collimator that delivers a highly conformed small radiation field, gives very good results in the control of the primary tumor as well as in preserving visual function, while it avoids the devastating side effects of deformation of midface bones. The purpose of the present paper is to propose a modified version of the "D"-shaped collimator that reduces even further the irradiation field with the scope to reduce as well the risk of radio-induced secondary malignancies. Concurrently, the new dedicated "D"-shaped collimator must be easier to build and at the same time produces dose distributions that only differ on the field size with respect to the dose distributions obtained by the current collimator in use. The scope of the former requirement is to facilitate the employment of the authors' irradiation technique both at the authors' and at other hospitals. The fulfillment of the latter allows the authors to continue using the clinical experience gained in more than 30 years. METHODS: The Monte Carlo code PENELOPE was used to study the effect that the different structural elements of the dedicated "D"-shaped collimator have on the absorbed dose distribution. To perform this study, the radiation transport through a Varian Clinac 2100 C/D operating at 6 MV was simulated in order to tally phase-space files which were then used as radiation sources to simulate the considered collimators and the subsequent dose distributions. With the knowledge gained in that study, a new, simpler, "D"-shaped collimator is proposed. RESULTS: The proposed collimator delivers a dose distribution which is 2.4 cm wide along the inferior-superior direction of the eyeball. This width is 0.3 cm narrower than that of the dose distribution obtained with the collimator currently in clinical use. The other relevant characteristics of the dose distribution obtained with the new collimator, namely, depth doses at clinically relevant positions, penumbrae width, and shape of the lateral profiles, are statistically compatible with the results obtained for the collimator currently in use. CONCLUSIONS: The smaller field size delivered by the proposed collimator still fully covers the planning target volume with at least 95% of the maximum dose at a depth of 2 cm and provides a safety margin of 0.2 cm, so ensuring an adequate treatment while reducing the irradiated volume.

Quality indexes based on water measurements for low and medium energy x-ray beams: a theoretical study with PENELOPE.

PURPOSE: To study the use of quality indexes based on ratios of absorbed doses in water at two different depths to characterize x-ray beams of low and medium energies. METHODS: A total of 55 x-ray beam spectra were generated with the codes XCOMP5R and SPEKCALC and used as input of a series of Monte Carlo simulations performed with PENELOPE, in which the percentage depth doses in water and the kQ,Q0 factors, defined in the TRS-398 protocol, were determined for each beam. Some of these calculations were performed by simulating the ionization chamber PTW 30010. RESULTS: The authors found that the relation between kQ,Q0 and the ratios of absorbed doses at two depths is almost linear. A set of ratios statistically compatible with that showing the best fit has been determined. CONCLUSIONS: The results of this study point out which of these ratios of absorbed doses in water could be used to better characterize x-ray beams of low and medium energies.

Retinoblastoma external beam photon irradiation with a special 'D'-shaped collimator: a comparison between measurements, Monte Carlo simulation and a treatment planning system calculation.

Retinoblastoma is the most common eye tumour in childhood. According to the available long-term data, the best outcome regarding tumour control and visual function has been reached by external beam radiotherapy. The benefits of the treatment are, however, jeopardized by a high incidence of radiation-induced secondary malignancies and the fact that irradiated bones grow asymmetrically. In order to better exploit the advantages of external beam radiotherapy, it is necessary to improve current techniques by reducing the irradiated volume and minimizing the dose to the facial bones. To this end, dose measurements and simulated data in a water phantom are essential. A Varian Clinac 2100 C/D operating at 6 MV is used in conjunction with a dedicated collimator for the retinoblastoma treatment. This collimator conforms a 'D'-shaped off-axis field whose irradiated area can be either 5.2 or 3.1 cm(2). Depth dose distributions and lateral profiles were experimentally measured. Experimental results were compared with Monte Carlo simulations' run with the penelope code and with calculations performed with the analytical anisotropic algorithm implemented in the Eclipse treatment planning system using the gamma test. penelope simulations agree reasonably well with the experimental data with discrepancies in the dose profiles less than 3 mm of distance to agreement and 3% of dose. Discrepancies between the results found with the analytical anisotropic algorithm and the experimental data reach 3 mm and 6%. Although the discrepancies between the results obtained with the analytical anisotropic algorithm and the experimental data are notable, it is possible to consider this algorithm for routine treatment planning of retinoblastoma patients, provided the limitations of the algorithm are known and taken into account by the medical physicist and the clinician. Monte Carlo simulation is essential for knowing these limitations. Monte Carlo simulation is required for optimizing the treatment technique and the dedicated collimator.

Low-dose radiation hyper-radiosensitivity in multicellular tumour spheroids.

OBJECTIVE: We propose and study a new model aimed at describing the low-dose hyper-radiosensitivity phenomenon appearing in the survival curves of different cell lines. METHODS: The model uses the induced repair assumption, considering that the critical dose at which this mechanism begins to act varies from cell to cell in a given population. The model proposed is compared with the linear-quadratic model and the modified linear-quadratic model, which is commonly used in literature and in which the induced repair is taken into account in a heuristic way. The survival curve for the MCF-7 line of human breast cancer is measured at low absorbed doses and the uncertainties in these doses are estimated using thermoluminiscent dosemeters. RESULTS: It is shown that these multicellular spheroids present low-dose hyper-radiosensitivity. The new model permits an accurate description of the data of two human cell lines (previously published) and of the multicellular spheroids of the MCF-7 line here measured. CONCLUSION: The model shows enough flexibility to account for data with very different characteristics and considers in a faithful way the hypothesis of the repair induction.

Neutron dosimetry in organs of an adult human phantom using linacs with multileaf collimator in radiotherapy treatments.

PURPOSE: To calculate absorbed doses due to neutrons in 87 organs/tissues for anthropomorphic phantoms, irradiated in position supine (head first into the gantry) with orientations anteroposterior (AP) and right-left (RLAT) with a 18 MV accelerator. Conversion factors from monitor units to µGy per neutron in organs, equivalent doses in organs/tissues, and effective doses, which permit to quantify stochastic risks, are estimated. METHODS: MAX06 and FAX06 phantoms were modeled with MCNPX and irradiated with a 18 MV Varian Clinac 2100C/D accelerator whose geometry included a multileaf collimator. Two actual fields of a pelvic treatment were simulated using electron-photon-neutron coupled transport. Absorbed doses due to neutrons were estimated from kerma. Equivalent doses were estimated using the radiation weighting factor corresponding to an average incident neutron energy 0.47 MeV. Statistical uncertainties associated to absorbed doses, as calculated by MCNPX, were also obtained. RESULTS: Largest doses were absorbed in shallowest (with respect to the neutron pathway) organs. In µGyMU(-1), values of 2.66 (for penis) and 2.33 (for testes) were found in MAX06, and 1.68 (for breasts), 1.05 (for lenses of eyes), and 0.94 (for sublingual salivary glands) in FAX06, in AP orientation. In RLAT, the largest doses were found for bone tissues (leg) just at the entrance of the beam in the body (right side in our case). Values, in µGyMU(-1), of 1.09 in upper leg bone right spongiosa, for MAX06, and 0.63 in mandible spongiosa, for FAX06, were found. Except for gonads, liver, and stomach wall, equivalent doses found for FAX06 were, in both orientations, higher than for MAX06. Equivalent doses in AP are higher than in RLAT for all organs/tissues other than brain and liver. Effective doses of 12.6 and 4.1 µSvMU(-1) were found for AP and RLAT, respectively. The organs/tissues with larger relative contributions to the effective dose were testes and breasts, in AP, and breasts and red marrow, in RLAT. Equivalent and effective doses obtained for MAX06/FAX06 were smaller (between 2 and 20 times) than those quoted for the mathematical phantoms ADAM/EVA in ICRP-74. CONCLUSIONS: The new calculations of conversion coefficients for neutron irradiation in AP and RLAT irradiation geometries show a reduction in the values of effective dose by factors 7 (AP) and 6 (RLAT) with respect to the old data obtained with mathematical phantoms. The existence of tissues or anatomical regions with maximum absorbed doses, such as penis, lens of eyes, fascia (part of connective tissue), etc., organs/tissues that classic mathematical phantoms did not include because they were not considered for the study of stochastic effects, has been revealed. Absorbed doses due to photons, obtained following the same simulation methodology, are larger than those due to neutrons, reaching values 100 times larger as the primary beam is approached. However, for organs far from the treated volume, absorbed photon doses can be up to three times smaller than neutron ones. Calculations using voxel phantoms permitted to know the organ dose conversion coefficients per MU due to secondary neutrons in the complete anatomy of a patient.

Neutron and photon spectra in LINACs.

A Monte Carlo calculation, using the MCNPX code, was carried out in order to estimate the photon and neutron spectra in two locations of two linacs operating at 15 and 18 MV. Detailed models of both linac heads were used in the calculations. Spectra were estimated below the flattening filter and at the isocenter. Neutron spectra show two components due to evaporation and knock-on neutrons. Lethargy spectra under the filter were compared to the spectra calculated from the function quoted by Tosi et al. that describes reasonably well neutron spectra beyond 1 MeV, though tends to underestimate the energy region between 10(-6) and 1 MeV. Neutron and the Bremsstrahlung spectra show the same features regardless of the linac voltage.

Quality control for system count rate performance with scatter in gamma cameras.

We aimed to analyze the optimal conditions to carry out the periodical statistical control tests of the gamma camera count rate performance. First we focused in reproducing the actual R(20) value of the gamma camera response. Second we studied the variability of this parameter in the statistical control test. We performed a reference measurement, which consisted of the determination of the complete curve relating observed and incident count rates, the counting model describing it and the reference R(20). This reference determined the conditions for the statistical control tests and the way to analyze the results obtained. Results from three different gamma cameras were studied. Each gamma camera showed a different behavior and required specific data analysis. The optimal conditions to perform the statistical control test were determined in each case. Our procedure provides the information necessary to correlate the average value of R(20) obtained in the quality control test with the reference one. The critical requirement to perform any statistical control test, that is to have a reduced variability of the control variable, can be fulfilled in this case only for relatively high activities.

Specific absorbed fractions in thyroid diagnostics and treatment: Monte Carlo calculation with PENELOPE.

In nuclear medicine, diagnostic and therapy procedures in which a certain radiopharmaceutical is administered to a patient are performed. An important point is the determination of the dose absorbed by the important organs of the patient due to these procedures. This dose depends on the particular radionuclide used and the so-called specific absorbed fractions. In this work, by means of Monte Carlo (MC) simulation, the specific absorbed fractions in case the thyroid gland acts as a source organ and for photon energies between 30 keV and 2 MeV have been determined. The computer code PENELOPE has been used as well as the adult male mathematical phantom provided with the distribution of this code. Three different simulation types were carried out. In one of them, only photon transport was considered. In the other two, electron transport was included, doing a detailed and a mixed simulation. In general, the fractions were estimated with uncertainties <9 %, for the mixed and detailed simulations, and <3 %, for the simulation in which only photons are included. For some target organs and, especially for energies <100 keV, the uncertainties found were larger. The results obtained here have been compared with those obtained by other authors using other MC codes. A good agreement has been found in 80 % of the cases.

Ambient neutron dose equivalent outside concrete vault rooms for 15 and 18 MV radiotherapy accelerators.

In this work, the ambient dose equivalent, H*(10), due to neutrons outside three bunkers that house a 15- and a 18-MV Varian Clinac 2100C/D and a 15-MV Elekta Inor clinical linacs, has been calculated. The Monte Carlo code MCNPX (v. 2.5) has been used to simulate the neutron production and transport. The complete geometries including linacs and full installations have been built up according to the specifications of the manufacturers and the planes provided by the corresponding medical physical services of the hospitals where the three linacs operate. Two of these installations, those lodging the Varian linacs, have an entrance door to the bunker while the other one does not, although it has a maze with two bends. Various treatment orientations were simulated in order to establish plausible annual equivalent doses. Specifically anterior-posterior, posterior-anterior, left lateral, right lateral orientations and an additional one with the gantry rotated 30° have been studied. Significant dose rates have been found only behind the walls and the door of the bunker, near the entrance and the console, with a maximum of 12 µSv h(-1). Dose rates per year have been calculated assuming a conservative workload for the three facilities. The higher dose rates in the corresponding control areas were 799 µSv y(-1), in the case of the facility which operates the 15-MV Clinac, 159 µSv y(-1), for that with the 15-MV Elekta, and 21 µSv y(-1) for the facility housing the 18-MV Varian. A comparison with measurements performed in similar installations has been carried out and a reasonable agreement has been found. The results obtained indicate that the neutron contamination does not increase the doses above the legal limits and does not produce a significant enhancement of the dose equivalent calculated. When doses are below the detection limits provided by the measuring devices available today, MCNPX simulation provides an useful method to evaluate neutron dose equivalents based on a detailed description of linac, patient and bunker.

Thermal drift reduction with multiple bias current for MOSFET dosimeters.

New thermal compensation methods suitable for p-channel MOSFET (pMOS) dosimeters with the usual dose readout procedure based on a constant drain current are presented. Measuring the source-drain voltage shifts for two or three different drain currents and knowing the value of the zero-temperature coefficient drain current, I(ZTC), the thermal drift of source-drain or threshold voltages can be significantly reduced. Analytical expressions for the thermal compensation have been theoretically deduced on the basis of a linear dependence on temperature of the parameters involved. The proposed thermal modelling has been experimentally proven. These methods have been applied to a group of ten commercial pMOS transistors (3N163). The thermal coefficients of the source-drain voltage and the threshold voltage were reduced from -3.0 mV °C(-1), in the worst case, down to -70 µV °C(-1). This means a thermal drift of -2.4 mGy °C(-1) for the dosimeter. When analysing the thermal drifts of all the studied transistors, in the temperature range from 19 to 36 °C, uncertainty was obtained in the threshold voltage due to a thermal drift of ±9 mGy (2 SD), a commonly acceptable value in most radiotherapy treatments. The procedures described herein provide thermal drift reduction comparable to that of other technological or numerical strategies, but can be used in a very simple and low-cost dosimetry sensor.

Monte Carlo simulation of the dynamic micro-multileaf collimator of a LINAC Elekta Precise using PENELOPE.

Micro-multileaf collimators are devices that are added to LINAC heads for stereotactic radiosurgery. In this work, the performance of an Elekta Precise LINAC with a dynamic micro-multileaf collimator manufactured by 3D-line has been studied. Monte Carlo simulations based on PENELOPE code and measurements with three different detectors (PTW Semiflex 31010 chamber, PTW PinPoint 31016 chamber and PTW Diode 60008) have been carried out. Simulations were tuned by reproducing the experimental TPR(20, 10) quality index, providing a nice description of both the PDD curve and the transverse profiles at the two depths measured. The geometry of the micro-multileaf collimator was tested by calculating the transmission through it, and it was needed to significantly reduce the leaf separation indicated by the manufacturer to reproduce the experimental results. An approximate simulation in which the transport of the particles traversing the dynamic micro-multileaf collimator was described in a simplified way was analyzed, providing good agreement with the full simulations. With the MC model fixed, output factors for various field sizes were calculated and compared to the experimental ones, obtaining good agreement. Percentage depth doses (PDDs) and transverse profiles at two depths measured with the diode for small fields were well reproduced by the simulation, while the measurements performed with the PinPoint chamber showed differences in the PDDs, at large depths, and transverse profiles, at the penumbra. Monte Carlo simulations and Semiflex and diode measurements, performed for a 7.0 cm × 7.0 cm field, were in nice agreement, while those obtained with the PinPoint chamber showed differences that increased with the depth in water. At the phantom entrance, all measurements showed non-negligible differences that made Monte Carlo a good option to estimate the absorbed dose in this region.

Neutron dose equivalent and neutron spectra in tissue for clinical linacs operating at 15, 18 and 20 MV.

In this work, the dose equivalent due to photoneutrons and the neutron spectra in tissue was calculated for various linacs (Varian Clinac 2100C, Elekta Inor, Elekta SL25 and Siemens Mevatron KDS) operating at energies between 15 and 20 MV, using the Monte Carlo code MCNPX (v. 2.5). The dose equivalent in an ICRU tissue phantom has been calculated for anteroposterior treatments with a detailed simulation of the geometry of the linac head and the coupled electron-photon-neutron transport. Neutron spectra at the phantom entrance and at 1-cm depth in the phantom, depth distribution of the neutron fluence in the beam axis and dose distributions outside the beam axis at various depths have also been calculated and compared with previously published results. The differences between the neutron production of the various linacs considered has been analysed. Varian linacs show a larger neutron production than the Elekta and Siemens linacs at the same operating energy. The dose equivalent due to neutrons produced by medical linacs operating at energies >15 MeV is relevant and should not be neglected because of the additional doses that patients can receive.