Search PubMed⌕ Search

Biomedical subjects

M Fippel

Publications and source records attributed to M Fippel.

7 recordsLinked to original sources
PubMed2000

Investigation of variance reduction techniques for Monte Carlo photon dose calculation using XVMC.

Several variance reduction techniques, such as photon splitting, electron history repetition, Russian roulette and the use of quasi-random numbers are investigated and shown to significantly improve the efficiency of the recently developed XVMC Monte Carlo code for photon beams in radiation therapy. It is demonstrated that it is possible to further improve the efficiency by optimizing transpon parameters such as electron energy cut-off, maximum electron energy step size, photon energy cut-off and a cut-off for kerma approximation, without loss of calculation accuracy. These methods increase the efficiency by a factor of up to 10 compared with the initial XVMC ray-tracing technique or a factor of 50 to 80 compared with EGS4/PRESTA. Therefore, a common treatment plan (6 MV photons, 10 x 10 cm2 field size, 5 mm voxel resolution, 1% statistical uncertainty) can be calculated within 7 min using a single CPU 500 MHz personal computer. If the requirement on the statistical uncertainty is relaxed to 2%, the calculation time will be less than 2 min. In addition, a technique is presented which allows for the quantitative comparison of Monte Carlo calculated dose distributions and the separation of systematic and statistical errors. Employing this technique it is shown that XVMC calculations agree with EGSnrc on a sub-per cent level for simulations in the energy and material range of interest for radiation therapy.

Algorithms↗
PubMed1999

Experimental investigation of a fast Monte Carlo photon beam dose calculation algorithm.

An experimental verification of the recently developed XVMC code, a fast Monte Carlo algorithm to calculate dose distributions of photon beams in treatment planning, is presented. The treatment head is modelled by a point source with energy distribution (primary photons) and an additional head scatter contribution. Utility software is presented, allowing the determination of the parameters for this model using a single measured depth dose curve in water. The simple beam model is considered to be a starting point for more complex models being planned for future versions of the code. This paper is mainly focused on the influence of the different techniques on variance reduction and material property determination for dose distributions. It is demonstrated that XVMC and the simple beam model reproduce measured (by a diamond detector) relative dose distributions with an accuracy of better than +/-2% in various homogeneous and inhomogeneous phantoms. Furthermore, relative dose distributions in solid state phantoms have been measured by film. Also for these cases, measured and calculated dose distributions agree within experimental uncertainty. The short calculation time (depending on voxel resolution, statistical accuracy, field size and energy, a span of 1 min to 1 h using a present-day personal computer) and an interface to a commercial planning system will allow the implementation of the code for routine treatment planning of clinical electron and photon beams.

Algorithms↗
PubMed1999

Fast Monte Carlo dose calculation for photon beams based on the VMC electron algorithm.

A new Monte Carlo algorithm for 3D photon dose calculation in radiation therapy is presented, which is based on the previously developed Voxel Monte Carlo (VMC) for electron beams. The main result is that this new version of VMC (now called XVMC) is more efficient than EGS4/PRESTA photon dose calculation by a factor of 15-20. Therefore, a standard treatment plan for photons can be calculated by Monte Carlo in about 20 min. on a "normal" personal computer. The improvement is caused mainly by the fast electron transport algorithm and ray tracing technique, and an initial ray tracing method to calculate the number of electrons created in each voxel by the primary photon beam. The model was tested in comparison to calculations by EGS4 using several fictive phantoms. In most cases a good coincidence has been found between both codes. Only within lung substitute dose differences have been observed.

Algorithms↗
PubMed1997

Electron beam dose calculations with the VMC algorithm and the verification data of the NCI working group.

The accuracy of the Monte Carlo algorithm for fast electron dose calculation, VMC, is demonstrated by comparing calculations with measurements performed by a working group of the National Cancer Institute (NCI) of the USA. For both energies investigated, 9 and 20 MeV, the measurements in water are taken to determine the energy spectra of the Varian Clinac 1800 accelerator. For the majority of the experiments a good agreement is found between measurements and VMC calculations. However, in some cases deviations have been observed, which could be explained by the incompletely known geometry on the one hand and by inconsistent data on the other hand. As a reference, dose distributions calculated by the MDAH pencil-beam algorithm are also presented. It is shown that, especially near low- or high-density inhomogeneities, large dose overestimations and underestimations are calculated by using a pencil-beam approach, whereas VMC is able to reproduce the correct doses for these cases also.

Algorithms↗
PubMed1996

3D electron dose calculation using a Voxel based Monte Carlo algorithm (VMC).

A new model for calculating electron beam dose has been developed. The algorithm is based on a two- or three-dimensional geometry defined by computerized tomography (CT) images. The Monte Carlo technique was used to solve the electron transport equation. However, in contrast to conventional Monte Carlo models (EGS4) several approximations and simplifications in the description of elementary electron processes were introduced reducing in this manner the computational time by a factor of about 35 without significant loss in accuracy. The Monte Carlo computer program does not need any precalculated data. The random access memory required is about 16 Mbytes for a 128(2) X 50 matrix, depending on the resolution of the CT cube. The Voxel Monte Carlo model (VMC) was tested in comparison to calculations by EGS4 and the "Hogstrom algorithm" (MDAH) using several fictive phantoms. In all cases a good coincidence has been found between EGS4 and VMC, especially near tissue inhomogeneities, whereas the MDAH algorithm has produced dose underestimations of up to 40%.

Algorithms↗