Set-up error validation with EPID images: Measurements vs Egs_cbct simulation

AimIn this study, the egs_cbct code’s ability to replicate an electronic portal imaging device (EPID) is explored.BackgroundWe have investigated head and neck (H&N) setup verification on an Elekta Precise linear accelerator. It is equipped with an electronic portal imaging device (EPID) that can capture a set of projection images over different gantry angles.Methods and materialsCone-beam computed tomography (CBCT) images were reconstructed from projection images of two different setup scenarios. Projections of an Anthropomorphic Rando head phantom were also simulated by using the egs_cbct Monte Carlo code for comparison with the measured projections.Afterwards, CBCT images were reconstructed from this data. Image quality was evaluated against a metric defined as the image acquisition interval (IAI). It determines the number of projection images to be used for CBCT image reconstruction.ResultsFrom this results it was established that phantom shifts could be determined within 2 mm and rotations within one degree accuracy using only 20 projection images (IAI = 10 degrees). Similar results were obtained with the simulated data.ConclusionIn this study it is demonstrated that a head and neck setup can be verified using substantially fewer projection images. Bony landmarks and air cavities could still be observed in the reconstructed Rando head phantom. The egs_cbct code can be used as a tool to investigate setup errors without tedious measurements with an EPID system.     

Simulation and measurement of microbeam dose distribution in lung tissue

Microbeam radiation therapy (MRT), a so far preclinical method in radiation oncology, modulates treatment doses on a micrometre scale. MRT uses treatment fields with a few ten micrometre wide high dose regions (peaks) separated by a few hundred micrometre wide low dose regions (valleys) and was shown to spare tissue much more effectively than conventional radiation therapy at similar tumour control rates. While preclinical research focused primarily on tumours of the central nervous system, recently also lung tumours have been suggested as a potential target for MRT.This study investigates the effect of the lung microstructure, comprising air cavities of a few hundred micrometre diameter, on the microbeam dose distribution in lung. In Monte Carlo simulations different models of heterogeneous lung tissue are compared with pure water and homogeneous air–water mixtures. Experimentally, microbeam dose distributions in porous foam material with cavity sizes similar to the size of lung alveoli were measured with film dosimetry at the European Synchrotron Radiation Facility (ESRF) in Grenoble, France.Simulations and experiments show that the microstructure of the lung has a huge impact on the local doses in the microbeam fields. Locally, material inhomogeneities may change the dose by a factor of 1.7, and also average peak and valley doses substantially differ from those in homogeneous material.Our results imply that accurate dose prediction for MRT in lung requires adequate models of the lung microstructure. Even if only average peak and valley doses are of interest, the assumption of a simple homogeneous air–water mixture is not sufficient. Since anatomic information on a micrometre scale are unavailable for clinical treatment planning, alternative methods and models have to be developed.     

Monte Carlo evaluation of occupational exposure during uterine artery embolization

PurposeUterine fibroids affect women mainly of childbearing age, an alternative for the treatment of these fibroids is uterine artery embolization (UAE), a minimally invasive procedure which uses fluoroscopy, providing radiation doses often high, due to the fact that professionals remain in the room throughout the procedure. In this work, equivalent and effective doses were evaluated for the main physician, for the assistant and for the patient during the UAE procedure.MethodsDoses were calculated using computer simulation with the Monte Carlo Method, and virtual anthropomorphic phantoms, in a typical scenario of interventional radiology with field sizes of 20 × 20, 25 × 25 and 32 × 32 cm², tube voltages of 70, 80, 90 and 100 kV, and projections of LAO45, RAO45 and PA.ResultsThe results showed that the highest doses received by the professionals were for the LAO45 projection with 32 × 32 cm² field size and 100 kV tube voltage, which is in accordance with the existing literature. The highest equivalent doses, without the protective equipment, were in the eyes, skin, breast and stomach for the main physician, and for the assistant they were in the eyes, breast, thyroid and skin. When she used the protective equipment, the highest equivalent doses for the main physician were on the skin, brain, bone marrow and bone surface, and for the assistant they were on the skin, brain, red bone marrow and bone surface.ConclusionsEffective doses increased up to 3186% for the main physician, and 2462% for the assistant, without protective equipment, thus showing their importance.     

Virtual patient 3D dose reconstruction using in air EPID measurements and a back-projection algorithm for IMRT and VMAT treatments

PurposeAt our institute, a transit back-projection algorithm is used clinically to reconstruct in vivo patient and in phantom 3D dose distributions using EPID measurements behind a patient or a polystyrene slab phantom, respectively. In this study, an extension to this algorithm is presented whereby in air EPID measurements are used in combination with CT data to reconstruct ‘virtual’ 3D dose distributions. By combining virtual and in vivo patient verification data for the same treatment, patient-related errors can be separated from machine, planning and model errors.Methods and materialsThe virtual back-projection algorithm is described and verified against the transit algorithm with measurements made behind a slab phantom, against dose measurements made with an ionization chamber and with the OCTAVIUS 4D system, as well as against TPS patient data. Virtual and in vivo patient dose verification results are also compared.ResultsVirtual dose reconstructions agree within 1% with ionization chamber measurements. The average γ-pass rate values (3% global dose/3 mm) in the 3D dose comparison with the OCTAVIUS 4D system and the TPS patient data are 98.5 ± 1.9%(1SD) and 97.1 ± 2.9%(1SD), respectively. For virtual patient dose reconstructions, the differences with the TPS in median dose to the PTV remain within 4%.ConclusionsVirtual patient dose reconstruction makes pre-treatment verification based on deviations of DVH parameters feasible and eliminates the need for phantom positioning and re-planning. Virtual patient dose reconstructions have additional value in the inspection of in vivo deviations, particularly in situations where CBCT data is not available (or not conclusive).     

EGSnrc application for IMRT planning

The aim of this study was to describe a detailed instruction of intensity modulated radiotherapy (IMRT) planning simulation using BEAMnrc-DOSXYZnrc code system (EGSnrc package) and present a new graphical user interface based on MATLAB code (The MathWorks) to combine more than one. 3ddose file which were obtained from the IMRT plan.This study was performed in four phases: the commissioning of Varian Clinac iX6 MV, the simulation of IMRT planning in EGSnrc, the creation of in-house VDOSE GUI, and the analysis of the isodose contour and dose volume histogram (DVH) curve from several beam angles. The plan paramaters in sequence and control point files were extracted from the planning data in Tan Tock Seng Hospital Singapore (multileaf collimator (MLC) leaf positions – bank A and bank B, gantry angles, coordinate of isocenters, and MU indexes).VDOSE GUI which was created in this study can display the distribution dose curve in each slice and beam angle. Dose distributions from various MLC settings and beam angles yield different dose distributions even though they used the same number of simulated particles. This was due to the differences in the MLC leaf openings in every field. The value of the relative dose error between the two dose ditributions for “body” was 51.23 %. The Monte Carlo (MC) data was normalized with the maximum dose but the analytical anisotropic algorithm (AAA) data was normalized by the dose in the isocenter.In this study, we have presented a Monte Carlo simulation framework for IMRT dose calculation using DOSXYZnrc source 21. Further studies are needed in conducting IMRT simulations using EGSnrc to minimize the different dose error and dose volume histogram deviation.     

In silico investigation of factors affecting the MV imaging performance of a novel water-equivalent EPID

PurposeA Geant4 model of a novel, water-equivalent electronic portal imaging device (EPID) prototype for radiotherapy imaging and dosimetry utilising an array of plastic scintillating fibres (PSFs) has been developed. Monte Carlo (MC) simulations were performed to quantify the PSF-EPID imaging performance and to investigate design aspects affecting performance for optimisation.MethodsUsing the Geant4 model, the PSF-EPID’s imaging performance for 6 MV photon beams was quantified in terms of its modulation transfer function (MTF), noise power spectrum (NPS) and detective quantum efficiency (DQE). Model parameters, including fibre dimensions, optical cladding reflectivity and scintillation yield, were varied to investigate impact on imaging performance.ResultsThe MC-calculated DQE(0) for the reference PSF-EPID geometry employing 30 mm fibres was approximately nine times greater than values reported for commercial EPIDs. When using 10 mm long fibres, the PSF-EPID DQE(0) was still approximately three times greater than that of a commercial EPID. Increased fibre length, cladding reflectivity and scintillation yield produced the greatest decreases in NPS and increases in DQE.ConclusionsThe potential to develop an optimised next-generation water-equivalent EPID with MV imaging performance at least comparable to commercial EPIDs has been demonstrated. Factors most important for optimising prototype design include fibre length, cladding reflectivity and scintillation yield.     

An in-vivo dosimetry procedure for Elekta step and shoot IMRT

PurposeThe aim of this work was to extend an in-vivo dosimetry (IVD) method, previously developed by the authors for 3D-conformal radiotherapy, to step and shoot IMRT treatments for pelvic tumors delivered by Elekta linacs.Materials and methodsThe algorithm is based on correlation functions to convert EPID transit signals into in-vivo dose values at the isocenter point, D_iso. The EPID images were obtained by the so-called “IMRT Dosimetric Weighting” mode as a superposition of many segment fields. This way each integral dosimetric image could be acquired in about 10 s after the end of beam delivery and could be processed while delivering the successive IMRT beams. A specific algorithm for D_iso reconstruction especially featured for step and shoot IMRT was implemented using a fluence inhomogeneity index, FI, introduced to describe the degree of beam modulation with respect to open beams. A γ-analysis of 2D-EPID images obtained day to day, resulted rapid enough to verify the plan delivery reproducibility.ResultsFifty clinical IMRT beams, planned for patients undergoing radiotherapy of pelvic tumors, were used to irradiate a homogeneous phantom. For each beam the agreement between the reconstructed dose, D_iso, and the TPS computed dose, D_iso,TPS, was well within 5%, while the mean ratio R = D_iso/D_iso,TPS resulted for 250 tests equal to 1.006 ± 0.036. The same beams were checked in vivo, i.e. during patient treatment delivery, obtaining 500 tests whose average R ratio resulted equal to 1.011 ± 0.042. The γ-analysis of the EPID images with 5% 3 mm criteria supplied 85% of the tests with pass rates γ_mean ≤ 0.5 and P_γ<1 ≥ 90%.     

An external dosimetry audit programme to credential static and rotational IMRT delivery for clinical trials quality assurance

PurposeExternal dosimetry audits give confidence in the safe and accurate delivery of radiotherapy. The RTTQA group have performed an on-site audit programme for trial recruiting centres, who have recently implemented static or rotational IMRT, and those with major changes to planning or delivery systems.MethodsMeasurements of reference beam output were performed by the host centre, and by the auditor using independent equipment. Verification of clinical plans was performed using the ArcCheck helical diode array.ResultsA total of 54 measurement sessions were performed between May 2014 and June 2016 at 28 UK institutions, reflecting the different combinations of planning and delivery systems used at each institution. Average ratio of measured output between auditor and host was 1.002 ± 0.006. Average point dose agreement for clinical plans was −0.3 ± 1.8%. Average (and 95% lower confidence intervals) of gamma pass rates at 2%/2 mm, 3%/2 mm and 3%/3 mm respectively were: 92% (80%), 96% (90%) and 98% (94%). Moderately significant differences were seen between fixed gantry angle and rotational IMRT, and between combination of planning systems and linac manufacturer, but not between anatomical treatment site or beam energy.ConclusionAn external audit programme has been implemented for universal and efficient credentialing of IMRT treatments in clinical trials. Good agreement was found between measured and expected doses, with few outliers, leading to a simple table of optimal and mandatory tolerances for approval of dosimetry audit results. Feedback was given to some centres leading to improved clinical practice.     

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.     

Development of a four-dimensional Monte Carlo dose calculation system for real-time tumor-tracking irradiation with a gimbaled X-ray head

PurposeTo develop a four-dimensional (4D) dose calculation system for real-time tumor tracking (RTTT) irradiation by the Vero4DRT.MethodsFirst, a 6-MV photon beam delivered by the Vero4DRT was simulated using EGSnrc. A moving phantom position was directly measured by a laser displacement gauge. The pan and tilt angles, monitor units, and the indexing time indicating the phantom position were also extracted from a log file. Next, phase space data at any angle were created from both the log file and particle data under the dynamic multileaf collimator. Irradiation both with and without RTTT, with the phantom moving, were simulated using several treatment field sizes. Each was compared with the corresponding measurement using films. Finally, dose calculation for each computed tomography dataset of 10 respiratory phases with the X-ray head rotated was performed to simulate the RTTT irradiation (4D plan) for lung, liver, and pancreatic cancer patients. Dose-volume histograms of the 4D plan were compared with those calculated on the single reference respiratory phase without the gimbal rotation [three-dimensional (3D) plan].ResultsDifferences between the simulated and measured doses were less than 3% for RTTT irradiation in most areas, except the high-dose gradient. For clinical cases, the target coverage in 4D plans was almost identical to that of the 3D plans. However, the doses to organs at risk in the 4D plans varied at intermediate- and low-dose levels.ConclusionsOur proposed system has acceptable accuracy for RTTT irradiation in the Vero4DRT and is capable of simulating clinical RTTT plans.     

Dosimetric and Monte Carlo verification of jaws-only IMRT plans calculated by the Collapsed Cone Convolution algorithm for head and neck cancers

AimThe aim of this study is to verify the Prowess Panther jaws-only intensity modulated radiation therapy (JO-IMRT) treatment planning (TP) by comparing the TP dose distributions for head-and-neck (H&N) cancer with the ones simulated by Monte Carlo (MC).BackgroundTo date, dose distributions planned using JO-IMRT for H&N patients were found superior to the corresponding three-dimensional conformal radiotherapy (3D-CRT) plans. Dosimetry of the JO-IMRT plans were also experimentally verified using an ionization chamber, MapCHECK 2, and Octavius 4D and good agreements were shown.Materials and methodsDose distributions of 15 JO-IMRT plans of nasopharyngeal patients were recalculated using the EGSnrc Monte Carlo code. The clinical photon beams were simulated using the BEAMnrc. The absorbed dose to patients treated by fixed-field IMRT was computed using the DOSXYZnrc. The simulated dose distributions were then compared with the ones calculated by the Collapsed Cone Convolution (CCC) algorithm on the TPS, using the relative dose error comparison and the gamma index using global methods implemented in PTW-VeriSoft with 3%/3 mm, 2%/2 mm, 1%/1 mm criteria.ResultsThere is a good agreement between the MC and TPS dose. The average gamma passing rates were 93.3 ± 3.1%, 92.8 ± 3.2%, 92.4 ± 3.4% based on the 3%/3 mm, 2%/2 mm, 1%/1 mm criteria, respectively.ConclusionsAccording to the results, it is concluded that the CCC algorithm was adequate for most of the IMRT H&N cases where the target was not immediately adjacent to the critical structures.     

Monte Carlo simulation for Neptun 10 PC medical linear accelerator and calculations of output factor for electron beam

AimExact knowledge of dosimetric parameters is an essential pre-requisite of an effective treatment in radiotherapy. In order to fulfill this consideration, different techniques have been used, one of which is Monte Carlo simulation.Materials and methodsThis study used the MCNP-4C^b to simulate electron beams from Neptun 10 PC medical linear accelerator. Output factors for 6, 8 and 10 MeV electrons applied to eleven different conventional fields were both measured and calculated.ResultsThe measurements were carried out by a Wellhofler-Scanditronix dose scanning system. Our findings revealed that output factors acquired by MCNP-4C simulation and the corresponding values obtained by direct measurements are in a very good agreement.ConclusionIn general, very good consistency of simulated and measured results is a good proof that the goal of this work has been accomplished.     

Dual-energy contrast-enhanced digital mammography: Glandular dose estimation using a Monte Carlo code and voxel phantom

PurposeTo estimate the mean glandular dose of contrast enhanced digital mammography, using the EGSnrc Monte Carlo code and female adult voxel phantom.MethodsAutomatic exposure control of full field digital mammography system was used for the selection of the X-ray spectrum and the exposure settings for dual energy imaging. Measurements of the air-kerma and of the half value layers were performed and a Monte Carlo simulation of the digital mammography system was used to compute the mean glandular dose, for breast phantoms of various thicknesses, glandularities and for different X-ray spectra (low and high energy).ResultsFor breast phantoms of 2.0–8.0 cm thick and 0.1–100% glandular fraction, CC view acquisition, from AEC settings, can result in a mean glandular dose of 0.450 ± 0.022 mGy −2.575 ± 0.033 mGy for low energy images and 0.061 ± 0.021 mGy – 0.232 ± 0.033 mGy for high energy images. In MLO view acquisition mean glandular dose values ranged between 0.488 ± 0.007 mGy – 2.080 ± 0.021 mGy for low energy images and 0.065 ± 0.012 mGy – 0.215 ± 0.010 mGy for high energy images.ConclusionThe low kV part of contrast enhanced digital mammography is the main contributor to total mean glandular breast dose. The results of this study can be used to provide an estimated mean glandular dose for individual cases.     

Skin models and their impact on mean glandular dose in mammography

PurposeTo quantify the influence of different skin models on mammographic breast dosimetry, based on dosimetric protocols and recent breast skin thickness findings.MethodsBy using an adapted PENELOPE (v. 2014) + PenEasy (v. 2015) Monte Carlo (MC) code, simulations were performed in order to obtain the mean glandular dose (MGD), the normalized MGD by incident air Kerma (DgN), and the glandular depth dose (GDD(z)). The geometry was based on a cranio-caudal mammographic examination. Monoenergetic and polyenergetic beams were implemented, for a breast thickness from 2 cm to 9 cm, with different compositions. Seven skin models were used: a 5 mm adipose layer; a skin layer ranging from 5 mm to 1.45 mm, a 1.45 mm skin thickness with a subcutaneous adipose layer of 2 mm and 3.55 mm.ResultsThe differences, for monoenergetic beams, are higher (up to 200%) for lower energies (8 keV), thicker and low glandular content breasts, decreasing to less than 5% at 40 keV. Without a skin layer, the differences reach a maximum of 1240%. The relative difference in DgN values for 1.45 mm skin and 5 mm adipose layers and polyenergetic beams varies from −14% to 12%.ConclusionsThe implemented MC code is suitable for mammography dosimetry calculations. The skin models have major impacts on MGD values, and the results complement previous literature findings. The current protocols should be updated to include a more realistic skin model, which provides a reliable breast dose estimation.     

A computationally efficient Monte-Carlo model for biomedical Raman spectroscopy

Monte Carlo (MC) modeling is a valuable tool to gain fundamental understanding of light-tissue interactions, provide guidance and assessment to optical instrument designs, and help analyze experimental data. It has been a major challenge to efficiently extend MC towards modeling of bulk-tissue Raman spectroscopy (RS) due to the wide spectral range, relatively sharp spectral features, and presence of background autofluorescence. Here, we report a computationally efficient MC approach for RS by adapting the massively-parallel Monte Carlo eXtreme (MCX) simulator. Simulation efficiency is achieved through “isoweight,” a novel approach that combines the statistical generation of Raman scattered and Fluorescence emission with a lookup-table-based technique well-suited for parallelization. The MC model uses a graphics processor to produce dense Raman and fluorescence spectra over a range of 800 − 2000 cm⁻¹ with an approximately 100× increase in speed over prior RS Monte Carlo methods. The simulated RS signals are compared against experimentally collected spectra from gelatin phantoms, showing a strong correlation.     

Three-dimensional quality assurance of IMRT prostate plans using gel dosimetry

Intensity modulated radiotherapy (IMRT) is one of the most modern radiation therapy treatment techniques. Although IMRT can deliver high and complex conformational doses to the tumor volume, its implementation requires rigorous quality assurance (QA) procedures that include a dosimetric pre-treatment verification of individual patient planning. This verification usually involves measuring a small volume of absolute dose with an ionization chamber and checking bi-dimensional fluency with an array of detectors. The planning technique has tri-dimensional characteristics, but no tridimensional dosimetry has been established in the clinical routine. One strategy to perform three-dimensional dosimetry is to use polymeric gels associated with magnetic resonance imaging to evaluate dose distribution. Here, we have compared the results of conventional QA procedures involving one- and two-dimensional dosimetry to the results of three-dimensional dosimetry conducted with MAGIC-f gel in 10 cases of prostate cancer IMRT planning. More specifically, we used the gamma index (3%/3 mm) to compare the results of three-dimensional dosimetry to the expected dose distributions obtained with the treatment planning system. Except for one IMRT treatment plan, the gel dosimetry results agreed with the conventional quality control and provided an overview of dose distribution in the target volume.     

Monte Carlo calculations of radiotherapy dose in “homogeneous” anatomy

Given the substantial literature on the use of Monte Carlo (MC) simulations to verify treatment planning system (TPS) calculations of radiotherapy dose in heterogeneous regions, such as head and neck and lung, this study investigated the potential value of running MC simulations of radiotherapy treatments of nominally homogeneous pelvic anatomy. A pre-existing in-house MC job submission and analysis system, built around BEAMnrc and DOSXYZnrc, was used to evaluate the dosimetric accuracy of a sample of 12 pelvic volumetric arc therapy (VMAT) treatments, planned using the Varian Eclipse TPS, where dose was calculated with both the Analytical Anisotropic Algorithm (AAA) and the Acuros (AXB) algorithm. In-house TADA (Treatment And Dose Assessor) software was used to evaluate treatment plan complexity, in terms of the small aperture score (SAS), modulation index (MI) and a novel exposed leaf score (ELS/ELA). Results showed that the TPS generally achieved closer agreement with the MC dose distribution when treatments were planned for smaller (single-organ) targets rather than larger targets that included nodes or metastases. Analysis of these MC results with reference to the complexity metrics indicated that while AXB was useful for reducing dosimetric uncertainties associated with density heterogeneity, the residual TPS dose calculation uncertainties resulted from treatment plan complexity and TPS model simplicity. The results of this study demonstrate the value of using MC methods to recalculate and check the dose calculations provided by commercial radiotherapy TPSs, even when the treated anatomy is assumed to be comparatively homogeneous, such as in the pelvic region.     

Dosimetric and bremsstrahlung performance of a single convergent beam for teletherapy device

The present work investigates preliminary feasibility and characteristics of a new type of radiation therapy modality based on a single convergent beam of photons. The proposal consists of the design of a device capable of generating convergent X-ray beams useful for radiotherapy. The main goal is to achieve high concentrated dose delivery. The first step is an analytical approach in order to characterize the dosimetric performance of the hypothetical convergent photon beam. Then, the validated FLUKA Monte Carlo main code is used to perform complete radiation transport to account also for scattering effects. The proposed method for producing convergent X-rays is mainly based on the bremsstrahlung effect. Hence the operating principle of the proposed device is described in terms of bremsstrahlung production. The work is mainly devoted characterizing the effect on the bremsstrahlung yield due to accessories present in the device, like anode material and geometry, filtration and collimation systems among others.The results obtained for in-depth dose distributions, by means of analytical and stochastic approaches, confirm the presence of a high dose concentration around the irradiated target, as expected. Moreover, it is shown how this spot of high dose concentration depends upon the relevant physical properties of the produced convergent photon beam.In summary, the proposed design for producing single convergent X-rays attained satisfactory performance for achieving high dose concentration around small targets depending on beam spot size that may be used for some applications in radiotherapy, like radiosurgery.     

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-