Biomedical subjects

N Reynaert

Publications and source records attributed to N Reynaert.

14 recordsLinked to original sources

PubMed2005

Experimental verification of lung dose with radiochromic film: comparison with Monte Carlo simulations and commercially available treatment planning systems.

The purpose of this study was to assess the absorbed dose in and around lung tissue by performing radiochromic film measurements, Monte Carlo simulations and calculations with superposition convolution algorithms. We considered a layered polystyrene phantom of 12 x 12 x 12 cm3 containing a central cavity of 6 x 6 x 6 cm3 filled with Gammex RMI lung-equivalent material. Two field configurations were investigated, a small 1 x 10 cm2 field and a larger 10 x 10 cm2 field. First, we performed Monte Carlo simulations to investigate the influence of radiochromic film itself on the measured dose distribution when the film intersects a lung-equivalent region and is oriented parallel to the central beam axis. To that end, the film and the lung-equivalent materials were modelled in detail, taking into account their specific composition. Next, measurements were performed with the film oriented both parallel and perpendicular to the central beam axis to verify the results of our Monte Carlo simulations. Finally, we digitized the phantom in two commercially available treatment planning systems, Helax-TMS version 6.1A and Pinnacle version 6.2b, and calculated the absorbed dose in the phantom with their incorporated superposition convolution algorithms to compare with the Monte Carlo simulations. Comparing Monte Carlo simulations with measurements reveals that radiochromic film is a reliable dosimeter in and around lung-equivalent regions when the film is positioned perpendicular to the central beam axis. Radiochromic film is also able to predict the absorbed dose accurately when the film is positioned parallel to the central beam axis through the lung-equivalent region. However, attention must be paid when the film is not positioned along the central beam axis, in which case the film gradually attenuates the beam and decreases the dose measured behind the cavity. This underdosage disappears by offsetting the film a few centimetres. We find deviations of about 3.6% between Monte Carlo and the superposition convolution algorithm of Pinnacle behind the lung region, for both field configurations. Pinnacle is quite accurate in the lung region. Deviations up to 5.6% for the small field are found in the lung region between Monte Carlo and the superposition convolution algorithm of Helax-TMS. Behind the lung region, Helax-TMS is in better agreement with Monte Carlo. Radiochromic film measurements or Monte Carlo simulations are reliable methods to establish the dose in and around lung tissue.

Algorithms↗

PubMed2005

The importance of accurate linear accelerator head modelling for IMRT Monte Carlo calculations.

Two Monte Carlo dose engines for radiotherapy treatment planning, namely a beta release of Peregrine and MCDE (Monte Carlo dose engine), were compared with Helax-TMS (collapsed cone superposition convolution) for a head and neck patient for the Elekta SLi plus linear accelerator. Deviations between the beta release of Peregrine and MCDE up to 10% were obtained in the dose volume histogram of the optical chiasm. It was illustrated that the differences are not caused by the particle transport in the patient, but by the modelling of the Elekta SLi plus accelerator head and more specifically the multileaf collimator (MLC). In MCDE two MLC modules (MLCQ and MLCE) were introduced to study the influence of the tongue-and-groove geometry, leaf bank tilt and leakage on the actual dose volume histograms. Differences in integral dose in the optical chiasm up to 3% between the two modules have been obtained. For single small offset beams though the FWHM of lateral profiles obtained with MLCE can differ by more than 1.5 mm from profiles obtained with MLCQ. Therefore, and because the recent version of MLCE is as fast as MLCQ, we advise to use MLCE for modelling the Elekta MLC. Nevertheless there still remains a large difference (up to 10%) between Peregrine and MCDE. By studying small offset beams we have shown that the profiles obtained with Peregrine are shifted, too wide and too flat compared with MCDE and phantom measurements. The overestimated integral doses for small beam segments explain the deviations observed in the dose volume histograms. The Helax-TMS results are in better agreement with MCDE, although deviations exceeding 5% have been observed in the optical chiasm. Monte Carlo dose deviations of more than 10% as found with Peregrine are unacceptable as an influence on the clinical outcome is possible and as the purpose of Monte Carlo treatment planning is to obtain an accuracy of 2%. We would like to emphasize that only the Elekta MLC has been tested in this work, so it is certainly possible that alpha releases of Peregrine provide more accurate results for other accelerators.

Dose-Response Relationship, Radiation↗

PubMed2004

DOSSCORE: an accelerated DOSXYZnrc code with an efficient stepping algorithm and scoring grid.

DOSSCORE is an accelerated version of DOSXYZnrc that allows photons to cross voxel boundaries of the same medium and utilizes a separate scoring grid superimposed on the geometrical grid. Two different stepping algorithms, the hownear method and the scaling method are implemented in DOSSCORE. The hownear method allows particles to travel larger distances in homogeneous regions where there is no interest in the dose deposition of these particles, whilst the scaling method utilizes a stepping algorithm in which particles are only slowed down by the boundaries of the geometrical voxels and not by the boundaries of the scoring voxels. For CT-based phantoms, only photon ray tracing is applied, which results in a rather modest speed gain of factor 1.2 compared to DOSXYZnrc. The hownear method and scaling method do not increase the speed for CT-based phantoms, but only for homogeneous phantoms and phantoms with a limited number of small heterogeneities. In cases where a small number of scoring voxels are needed, the hownear method performs better than the scaling method, whilst the opposite is true for cases when many scoring voxels are needed. The photon transport is accelerated by almost a factor of 2 for all phantoms (homogeneous, heterogeneous with much homogeneity and CT-based phantoms) compared to DOSXYZnrc. For a small number of scoring voxels, the hownear method is up to a factor of 2.6 and 1.9 faster than DOSXYZnrc for homogeneous and heterogeneous phantoms in the case of photon beams. For an electron beam, a speed gain of factor 2.4 is obtained. For a full scoring grid like the one used in DOSXYZnrc, the scaling method is up to a factor of 2.2 and 1.7 faster than DOSXYZnrc for homogeneous and heterogeneous phantoms in the case of photon beams. For an electron beam, a speed gain of factor 2 is obtained. A speed increase without biasing the results is very relevant. The use of two separate grids, the more efficient stepping algorithms and the accelerated photon transport can be applied to every EGS-based or other Monte Carlo code.

Algorithms↗

PubMed2004

MCDE: a new Monte Carlo dose engine for IMRT.

A new accurate Monte Carlo code for IMRT dose computations, MCDE (Monte Carlo dose engine), is introduced. MCDE is based on BEAMnrc/DOSXYZnrc and consequently the accurate EGSnrc electron transport. DOSXYZnrc is reprogrammed as a component module for BEAMnrc. In this way both codes are interconnected elegantly, while maintaining the BEAM structure and only minimal changes to BEAMnrc.mortran are necessary. The treatment head of the Elekta SLiplus linear accelerator is modelled in detail. CT grids consisting of up to 200 slices of 512 x 512 voxels can be introduced and up to 100 beams can be handled simultaneously. The beams and CT data are imported from the treatment planning system GRATIS via a DICOM interface. To enable the handling of up to 50 x 10(6) voxels the system was programmed in Fortran95 to enable dynamic memory management. All region-dependent arrays (dose, statistics, transport arrays) were redefined. A scoring grid was introduced and superimposed on the geometry grid, to be able to limit the number of scoring voxels. The whole system uses approximately 200 MB of RAM and runs on a PC cluster consisting of 38 1.0 GHz processors. A set of in-house made scripts handle the parallellization and the centralization of the Monte Carlo calculations on a server. As an illustration of MCDE, a clinical example is discussed and compared with collapsed cone convolution calculations. At present, the system is still rather slow and is intended to be a tool for reliable verification of IMRT treatment planning in the case of the presence of tissue inhomogeneities such as air cavities.

Adenocarcinoma↗

PubMed2003

The value of radiochromic film dosimetry around air cavities: experimental results and Monte Carlo simulations.

In this study we investigate radiochromic film dosimetry around air cavities with particular focus on the perturbation of the dose distribution by the film when the film is parallel to the beam axis. We considered a layered polystyrene phantom containing an air cavity as a model for the air-soft tissue geometry that may occur after surgical resection of a paranasal sinus tumour. A radiochromic film type MD-55 was positioned within the phantom so that it intersected the cavity. Two phantom set-ups were examined. In the first case, the air cavity is at the centre of the phantom, thus the film is lying along the central beam axis. In the second case, the cavity and film are located 2 cm offset from the phantom centre and the central beam axis. In order to examine the influence of the film on the dose distribution and to interpret the film-measured results, Monte Carlo simulations were performed. The film was modelled rigorously to incorporate the composition and structure of the film. Two field configurations, a 1 x 10 cm2 field and a 10 x 10 cm2 field, were examined. The dose behind the air cavity is reduced by 6 to 7% for both field configurations when a film that intersects the cavity contains the central beam axis. This is due to the attenuation exerted by the film when photons cross the cavity. Offsetting the beam to the cavity and the film by 2 cm removes the dose reduction behind the air cavity completely. Another result was that the rebuild-up behind the cavity for the 10 x 10 cm2 field, albeit less significant than for the 1 x 10 cm2 field, could only be measured by the film that was placed offset with respect to the central beam axis. Although radiochromic film is approximately soft-tissue equivalent and energy independent as compared to radiographic films, care should be taken in the case of inhomogeneous phantoms when the film intersects air cavities and contains the beam central axis. Errors in dose measurement can be expected distal to the air cavity due to attenuation in the film itself. This attenuation would not occur in the absence of the film. Both experiments and Monte Carlo computations support this conclusion.

Air↗

PubMed2003

Monte Carlo model of the Elekta SLiplus accelerator: validation of a new MLC component module in BEAM for a 6 MV beam.

A new component module (CM), called MLCE, has been implemented in the BEAM program. The CM takes into account the particular 'tongue-and-groove' design of the Elekta multi-leaf collimator (MLC) and the air gap between the leaves. The model was validated by two series of measurements and simulations. The first benchmarking series focuses on the interleaf leakage and the intraleaf transmission. The measurement showed a total transmission through the MLC of 1.42% of the open field dose. Two Monte Carlo (MC) simulations were made, the first with the new CM MLCE (inclusive of air gap) and the second with the CM MLCQ (exclusive of air gap), which is available in the BEAM distribution. When the air gap between the leaves was determined by varying the parameters of the leaf geometry within tolerance limits on the technical drawing, the total measured transmission of 1.42% was well reproduced by the CM MLCE. In contrast, MC simulations with MLCQ showed that the transmission through the MLC calculated without the interleaf leakage is only 44% of the total transmitted radiation. The relevance of the detailed MLC modelling was demonstrated also by studying the 'adjacent' tongue-and-groove effect, where two adjacent (not opposing) leaves are complementary, opened or closed. The two complementary leaf settings were simulated both with the CM MLCE and MLCQ. A comparison with measurements was made. In regions covered by two or more leaves, the dose increased by 14% for two leaves and by 40% for more than two leaves when the interleaf leakage was included in the transmission. The tongue-and-groove effect was perfectly reproduced by the MLCE module.

Algorithms↗

PubMed2002

Underdosage of the upper-airway mucosa for small fields as used in intensity-modulated radiation therapy: a comparison between radiochromic film measurements, Monte Carlo simulations, and collapsed cone convolution calculations.

Head-and-neck tumors are often situated at an air-tissue interface what may result in an underdosage of part of the tumor in radiotherapy treatments using megavoltage photons, especially for small fields. In addition to effects of transient electronic disequilibrium, for these small fields, an increased lateral electron range in air will result in an important extra reduction of the central axis dose beyond the cavity. Therefore dose calculation algorithms need to model electron transport accurately. We simulated the trachea by a 2 cm diameter cylindrical air cavity with the rim situated 2 cm beneath the phantom surface. A 6 MV photon beam from an Elekta SLiplus linear accelerator, equipped with the standard multileaf collimator (MLC), was assessed. A 10 x 2 cm2 and a 10 x 1 cm2 field, both widthwise collimated by the MLC, were applied with their long side parallel to the cylinder axis. Central axis dose rebuild-up was studied. Radiochromic film measurements were performed in an in-house manufactured polystyrene phantom with the films oriented either along or perpendicular to the beam axis. Monte Carlo simulations were performed with BEAM and EGSnrc. Calculations were also performed using the pencil beam (PB) algorithm and the collapsed cone convolution (CCC) algorithm of Helax-TMS (MDS Nordion, Kanata, Cahada) version 6.0.2 and using the CCC algorithm of Pinnacle (ADAC Laboratories, Milpitas, CA, USA) version 4.2. A very good agreement between the film measurements and the Monte Carlo simulations was found. The CCC algorithms were not able to predict the interface dose accurately when lateral electronic disequilibrium occurs, but were shown to be a considerable improvement compared to the PB algorithm. The CCC algorithms overestimate the dose in the rebuild-up region. The interface dose was overestimated by a maximum of 31% or 54%, depending on the implementation of the CCC algorithm. At a depth of 1 mm, the maximum dose overestimation was 14% or 24%.

Air↗

PubMed2002

Parameter dependence of the MCNP electron transport in determining dose distributions.

In this paper, a detailed study of the electron transport in MCNP is performed, separating the effects of the energy binning technique on the energy loss rate, the scattering angles, and the sub-step length as a function of energy. As this problem is already well known, in this paper we focus on the explanation as to why the default mode of MCNP can lead to large deviations. The resolution dependence was investigated as well. An error in the MCNP code in the energy binning technique in the default mode (DBCN 18 card = 0) was revealed, more specific in the updating of cross sections when a sub-step is performed corresponding to a high-energy loss. This updating error is not present in the ITS mode (DBCN 18 card = 1) and leads to a systematically lower dose deposition rate in the default mode. The effect is present for all energies studied (0.5-10 MeV) and depends on the geometrical resolution of the scoring regions and the energy grid resolution. The effect of the energy binning technique is of the same order of that of the updating error for energies below 2 MeV, and becomes less important for higher energies. For a 1 MeV point source surrounded by homogeneous water, the deviation of the default MCNP results at short distances attains 9% and remains approximately the same for all energies. This effect could be corrected by removing the completion of an energy step each time an electron changes from an energy bin during a sub-step. Another solution consists of performing all calculations in the ITS mode. Another problem is the resolution dependence, even in the ITS mode. The higher the resolution is chosen (the smaller the scoring regions) the faster the energy is deposited along the electron track. It is proven that this is caused by starting a new energy step when crossing a surface. The resolution effect should be investigated for every specific case when calculating dose distributions around beta sources. The resolution should not be higher than 0.85*(1-EFAC)*CSDA, where EFAC is the energy loss per energy step and CSDA a continuous slowing down approximation range. This effect could as well be removed by determining the cross sections for energy loss and multiple scattering at the average energy of an energy step and by sampling the cross sections for each sub-step. Overall, we conclude that MCNP cannot be used without a caution due to possible errors in the electron transport. When care is taken, it is possible to obtain correct results that are in agreement with other Monte Carlo codes.

Algorithms↗

PubMed2001

On the accuracy of monomer/polymer gel dosimetry in the proximity of a high-dose-rate 192Ir source.

The aim of this work was to investigate the applicability of MR-based polymer gel dosimetry to measure the absorbed dose distribution at short distance from an iridium-192 brachytherapy point source. In this paper, different methodological problems that may result in significant errors in the measured dose distribution are discussed. First of all the extent to which physicochemical mechanisms alter the dose response is discussed. The most important among these are the oxygen permeability of the catheter material and monomer-diffusion-related effects during irradiation. The effect of oxygen on the dose-R2 curve has been determined quantitatively and an oxygen map is performed using a well-defined external irradiation beam. The effect of diffusion of monomer during irradiation has been simulated. Another contribution of errors is related to magnetic susceptibility differences between the catheter and the gel during scanning the irradiated gel. The magnetic field distortion has been mapped by using both an experimental MRI technique and by simulation. Moreover, in constructing a dose-versus-distance curve by polar averaging, the sensitivity to the exact distance between source and point of measurement and to partial volume effects is illustrated. An optimization program is proposed to determine the location of the source on a sub-pixel scale.

Acrylamide↗

PubMed2001

Dosimetry of 192Ir sources used for endovascular brachytherapy.

An in-phantom calibration technique for 192Ir sources used for endovascular brachytherapy is presented. Three different source lengths were investigated. The calibration was performed in a solid phantom using a Farmer-type ionization chamber at source to detector distances ranging from 1 cm to 5 cm. The dosimetry protocol for medium-energy x-rays extended with a volume-averaging correction factor was used to convert the chamber reading to dose to water. The air kerma strength of the sources was determined as well. EGS4 Monte Carlo calculations were performed to determine the depth dose distribution at distances ranging from 0.6 mm to 10 cm from the source centre. In this way we were able to convert the absolute dose rate at 1 cm distance to the reference point chosen at 2 mm distance. The Monte Carlo results were confirmed by radiochromic film measurements, performed with a double-exposure technique. The dwell times to deliver a dose of 14 Gy at the reference point were determined and compared with results given by the source supplier (CORDIS). They determined the dwell times from a Sievert integration technique based on the source activity. The results from both methods agreed to within 2% for the 12 sources that were evaluated. A Visual Basic routine that superimposes dose distributions, based on the Monte Carlo calculations and the in-phantom calibration, onto intravascular ultrasound images is presented. This routine can be used as an online treatment planning program.

Angioplasty, Balloon, Coronary↗

PubMed2001

Self-absorption correction for 32P, 198Au and 188Re stents: dose point kernel calculations versus Monte Carlo.

Monte Carlo simulations of dose distributions around radioactive stents are very time intensive. Thus, in order to calculate the dose distribution around a 188Re stent, we chose to test a point kernel method, a method which is known to be faster but the accuracy of which has not been established for this application. The dose point kernel method, which takes into account absorption in the strut material (=self-absorption), was based on different beta-emitting point source distributions in water by itself and surrounded by steel spheres of different thicknesses. This information was input into an integration routine that modeled either a Palmaz-Schatz or Multilink stent. The dose distributions around 198Au and 32P stents calculated with the dose point kernel method were compared to those calculated using EGS4 and MCNP 4B Monte Carlo methods. The resulting correction for self-absorption in steel was distance dependent and averaged 1.12 for 32P and 1.25 for 198Au stents. The dose point kernel method gave nearly identical results to these full Monte Carlo simulations and was thus used to calculate the dose distributions around a 188Re stent. Although 188Re has a half-life of only 17 hours, it is posited to be useful for radioactive restenosis prevention, given that a recently developed rapid electrodeposition procedure allows stents to be made radioactive, at predetermined activities, within 15 minutes. The dose point kernel calculations of a 188Re-coated Multilink stent were compared to its radiochromic film measurements. The dose fall-off agreed with the calculations within 5% over 0.4 to 3.5 mm from the stent surface. The dose point kernel method is a valuable tool to determine depth dose distributions around activated stents taking into account the detailed geometry and the self-absorption in the struts. It not only requires much less processing time than Monte Carlo methods, but also allows the use of higher resolutions in modeling the geometry, which leads to more accurate self-absorption correction factors.

Angioplasty, Balloon, Coronary↗

PubMed2000

Verification of absolute dose rates for intravascular brachytherapy beta sources.

The absolute dose rate around beta sources for intravascular brachytherapy was measured with thin thermoluminescent (TL) dosimeters and radiochromic film. The calibration of these dosimeters was performed at the dose maximum in a 6 MeV electron beam of a linear accelerator. Compared to the values given by the manufacturers, higher dose rates of about 20% for 90Sr/90Y source trains and 20 to 40% for 90Y line sources were found. We also measured the dose rate for a planar 90Sr/90Y source, which was calibrated at NIST (National Institute of Standards and Technology). This source serves as a reference for the determination of the dose rates by the supplier of the 90Y sources. For this source we measured a 13% higher dose rate compared to the value given in the calibration report by NIST. Our observation can have important consequences for the interpretation and comparison of the dose-related clinical results of intravascular brachytherapy with these and other beta sources. Also the comparison of the clinical results obtained with beta sources and gamma sources, i.e. 192Ir, depends largely on the accuracy of the absolute dosimetry of the beta sources.

Blood Vessels↗

PubMed1999

Monte Carlo calculations of dose distributions around 32P and 198Au stents for intravascular brachytherapy.

3D dose distributions are calculated for a 32P impregnated stent and a 198Au stent for intravascular brachytherapy with the EGS4 Monte Carlo simulation code. The stents were modeled as a combination of eight helicoidal struts. This allowed investigation of the effect of the stent geometry and the electron absorption in the strut material on the dose distributions. Absorbed dose to water was calculated at radial distances ranging from 50 microm to 5 mm from the stent surface. The dose distributions around the stents are compared to the dose distribution around an intravascular brachy-therapy 192Ir source, also calculated with the EGS4 Monte Carlo code. The dose profiles near the struts show hot spots. At 50 microm distance a peak to valley ratio of 3 for 32P and 6 for 198Au in the dose distribution is obtained. For both the isotopes the inhomogeneities decrease with distance and at a radial depth of 350 microm the effect becomes negligible. The calculations showed the importance of the effect of the absorption in the stent material as this leads to a dose decrease to 67% for the 198Au stent and to 77% for 32P near the stent at a distance of 2 mm from the stent axis. It is concluded that from the dosimetric point of view, the 198Au stent is inferior to the 32P stent and the 192Ir source. Application of the 198Au stent in clinical practice requires further investigation of the importance of the adventitia in the restenosis process, and the tolerance dose of the intima.

Angioplasty, Balloon, Coronary↗

PubMed1998

In-water calibration of PDR 192Ir brachytherapy sources with an NE2571 ionization chamber.

An ionometric calibration procedure for 192Ir PDR brachytherapy sources in terms of dose rate to water is presented. The calibration of the source is performed directly in a water phantom at short distances (1.0, 2.5 and 5.0 cm) using an NE2571 Farmer type ion chamber. To convert the measured air-kerma rate in water to dose rate to water a conversion factor (CF) was calculated by adapting the medium-energy x-ray dosimetry protocol for a point source geometry (diverging beam). The obtained CF was verified using two different methods. Firstly, the CF was calculated by Monte Carlo simulations, where the source-ionization chamber geometry was modelled accurately. In a second method, a combination of Monte Carlo simulations and measurements of the air-kerma rate in water (at 1.0, 2.5 and 5.0 cm distance) and in air (1 m distance) was used to determine the CF. The obtained CFs were also compared with conversion factors calculated with the adapted dosimetry protocol for high-energy photons introduced by Tölli. All calculations were done for a Gammamed PDR 192Ir source-NE2571 chamber geometry. The conversion factors obtained with the four different methods agree to within 1% at the three distances of interest. We obtained the following values (medium-energy x-ray protocol): CF(1 cm) = 1.458; CF(2.5 cm) = 1.162; CF(5.0 cm) = 1.112 (1 sigma = 0.7% for the three distances of interest). The obtained results were checked with TLD measurements. The values of the specific dose rate constant and the radial dose function calculated in this work are in accordance with the literature data.

Brachytherapy↗