An important requirement in radiation therapy is a fast and accurate treatment planning system. This system, using computed tomography (CT) data, direction, and characteristics of the beam, calculates the dose at all points of the patient's volume. The two main factors in treatment planning system are accuracy and speed. According to these factors, various generations of treatment planning systems are developed. This article is a review of the Fast Monte Carlo treatment planning algorithms, which are accurate and fast at the same time. The Monte Carlo techniques are based on the transport of each individual particle (e.g., photon or electron) in the tissue. The transport of the particle is done using the physics of the interaction of the particles with matter. Other techniques transport the particles as a group. For a typical dose calculation in radiation therapy the code has to transport several millions particles, which take a few hours, therefore, the Monte Carlo techniques are accurate, but slow for clinical use. In recent years, with the development of the ‘fast’ Monte Carlo systems, one is able to perform dose calculation in a reasonable time for clinical use. The acceptable time for dose calculation is in the range of one minute. There is currently a growing interest in the fast Monte Carlo treatment planning systems and there are many commercial treatment planning systems that perform dose calculation in radiation therapy based on the Monte Carlo technique.

The purpose of this work is to revisit the impediments and characteristics of fast Monte Carlo techniques for applications in radiation therapy treatment planning using new methods of utilizing pregenerated electron tracks. The limitations of various techniques for the improvement of speed and accuracy of electron transport have been evaluated. A method is proposed that takes advantage of large available memory in current computer hardware for extensive generation of precalculated data. Primary tracks of electrons are generated in the middle of homogeneous materials (water, air, bone, lung) and with energies between 0.2 and 18 MeV using the EGSnrc code. Secondary electrons are not transported, but their position, energy, charge, and direction are saved and used as a primary particle. Based on medium type and incident electron energy, a track is selected from the precalculated set. The performance of the method is tested in various homogeneous and heterogeneous configurations and the results were generally within 2% compared to EGSnrc but with a 40-60 times speed improvement. In a second stage the authors studied the obstacles for further increased speed-ups in voxel geometries by including ray-tracing and particle fluence information in the pregenerated track information. The latter method leads to speed increases of about a factor of 500 over EGSnrc for voxel-based geometries. In both approaches, no physical calculation is carried out during the runtime phase after the pregenerated data has been stored even in the presence of heterogeneities. The precalculated data are generated for each particular material and this improves the performance of the precalculated Monte Carlo code both in terms of accuracy and speed. Precalculated Monte Carlo codes are accurate, fast, and physics independent and therefore applicable to different radiation types including heavy-charged particles.

The Monte Carlo method is the most accurate method for simulation of radiation therapy equipment. The linear accelerators (linac) are currently the most widely used machines in radiation therapy centers. In this work, a Monte Carlo modeling of the Siemens ONCOR linear accelerator in 6 MV and 18 MV beams was performed. The results of simulation were validated by measurements in water by ionization chamber and extended dose range (EDR2) film in solid water. The linac's X-ray particular are so sensitive to the properties of primary electron beam. Square field size of 10 cm × 10 cm produced by the jaws was compared with ionization chamber and film measurements. Head simulation was performed with BEAMnrc and dose calculation with DOSXYZnrc for film measurements and 3ddose file produced by DOSXYZnrc analyzed used homemade MATLAB program. At 6 MV, the agreement between dose calculated by Monte Carlo modeling and direct measurement was obtained to the least restrictive of 1%, even in the build-up region. At 18 MV, the agreement was obtained 1%, except for in the build-up region. In the build-up region, the difference was 1% at 6 MV and 2% at 18 MV. The mean difference between measurements and Monte Carlo simulation is very small in both of ONCOR X-ray energy. The results are highly accurate and can be used for many applications such as patient dose calculation in treatment planning and in studies that model this linac with small field size like intensity-modulated radiation therapy technique.

An important requirement for proton therapy is a software for dose calculation. Monte Carlo is the most accurate method for dose calculation, but it is very slow. In this work, a method is developed to improve the speed of dose calculation. The method is based on pre-generated tracks for particle transport. The MCNPX code has been used for generation of tracks. A set of data including the track of the particle was produced in each particular material (water, air, lung tissue, bone, and soft tissue). This code can transport protons in wide range of energies (up to 200 MeV for proton). The validity of the fast Monte Carlo (MC) code is evaluated with data MCNPX as a reference code. While analytical pencil beam algorithm transport shows great errors (up to 10%) near small high density heterogeneities, there was less than 2% deviation of MCNPX results in our dose calculation and isodose distribution. In terms of speed, the code runs 200 times faster than MCNPX. In the Fast MC code which is developed in this work, it takes the system less than 2 minutes to calculate dose for 106 particles in an Intel Core 2 Duo 2.66 GHZ desktop computer.

Dose escalation with high-energy X rays of medical linear accelerators (linacs) in radiotherapy offers several distinct advantages over the lower energy photons. However, owing to photoneutron reactions, interaction of high-energy photons (>8 MV) with various high-Z nuclei of the materials in the linac head components produces unavoidable neutrons. The aim of this study was to evaluate the photoneutron dose equivalent per unit therapeutic X-ray dose of 18 MV, GE Saturne 20 linac in the treatment room using Monte Carlo (MC) MCNP linac head full simulation as well as thermoluminescence dosemeter measurements. This machine is one of the old linac models manufactured by General Electric Company; however, it is widely used in the developing countries because of low cost and simple maintenance for radiotherapy applications. The results showed a significant photoneutron dose from Saturne 20 linac head components especially at distances near the linac head (<150 cm). Results of this work could be used in several applications, especially designing bunker and entrance door shielding against neutrons produced by photoneutron reactions in GE Saturne 20. However, a detailed cost optimisation for a specific room would require a dedicated calculation.